板壳问题分析有限元分析作业

5.案例分析5.1问题描述与分析货柜的结构如下图所示，各层隔板厚度为0.015m，两块侧立板厚度为0.01m，背板的厚度为0.01m，材料的弹性模量为1.2X109N/M2,泊松比为0.3，密度为800Kg/m3。

假设货柜底板完全固定，中间三块隔板上存放货物，货物产生的均匀分布应力分布压力作用在上面，压力大小等于200pa，计算货柜的变形和等效应力。

试计算轴承座结构的变形及von Mises 应力分布。

采用自顶向下和自底向上相结合的建模方法建立本实例的三维模型，选择shell63单元划分网格，单元长度控制为0.1。

图1 货柜的结构图5.2定义单元类型和材料属性在单元类型库里面选择Structutal shell ,shell63 ,相应单元类型号设为1。

材料类型定义为：Structutal/Linear/Elastic/Isotropic,并将其弹性常数设为1.2X109N/M2，泊松比设为0.3，密度为800Kg/m3。

5.3创建几何模型与网格划分在几何模型创建中采用Main menu-preprocessor-modeling-areas创建薄板，然后通过Main menu-preprocessor-modeling-copy-areas复制货柜相同的薄板，最后经过通过Main menu-preprocessor-modeling-operate-booleans- partition-volume进行粘接各面得模型图，如图1所示。

用网格划分器MeshTool对几何模型进行单元划分。

选择Main menu-preprocessor-MeshTool命令，弹出MeshTool对话框，选中element attributes下拉列表中的areas，将网格单元长度设置为0.1，划分所得的网格如图2所示。

图2 货柜的单元划分结果图5.4施加边界条件与载荷选择Main menu-solution-define loads-apply-structural-displacement-on Areas对货柜底部横隔板面进行约束；选择Main menu-solution-define loads-apply-structural-pressure-on Areas在货柜中间三个横隔板上施加压力载荷为200 Pa，如图3所示。

基于面积坐标性质，上述位移模式还可改写为：
薄板三角形单元
进行一系列的推导后，可得到
形函数的具体计算式为：
薄板三角形单元
建立了位移模式之后，那么剩下的工作就并不复杂：位移模式→ 应变离散→ 应力离散→刚度矩阵→载荷向量→约束处理→求解。
但位移模式是建立在面积坐标上的，相关的计算怎么进行？
导数间的关系为
单元数 （1/4板）
2×2 4×4 6×6 理论解
四边固定
板中心挠度 wD/PL2
边中点弯矩 M/P
0.00614
-0.1178
0.00580
-0.1233
0.00571
-0.1245
0.00560
-0.1257
薄板三角形单元
三角形单元能较好地适应斜边界，实 际中广泛应用。单元的结点位移仍然 为结点处的挠度wi和绕x，y轴的转角
壳体：壳体的变形除了横向弯曲变形外，同时存在中面变形。 因此可以认为壳体是平面应力问题和平板弯曲问题的组合。当 然，对于厚壳结构，仍需要横向剪切变形的影响。
薄板结构有限元
薄板基础理论知识
薄平板，取其中性面为坐标面，z轴垂直于中性面。其中 t 为 板厚。当板受有垂直于板中性面的外力时，板的中性面将发 生弯扭变形，从而变成一个曲面。板变形的同时，在板的横 截面上将存在内力——弯矩和扭矩。
Mxy=Myx
内力列向量为
薄板基础理论知识
内力可以根据应力进行计算得到
使用记号
平面应力问题 中的弹性矩阵
于是
薄板基础理论知识
进行反向回代，可以得到
在板的上、下表面处，z=±0.5t，于是应力为
薄板基础理论知识
如果薄板在z方向承受分布荷载
此时薄板内部产生应力 则可以采用虚功原理

T i
1 1 2 h 1 1 2
h

BiT DB j abd d dz
(6.17)
21 /44
薄板问题的有限元法
代入 D 、 Bi 和 B j 于是有
D 1 1 b2 T kij N i , N j , uN iT, N T, uN iT, N T, j j 1 1 a 2 ab +2(1- )N
2
24 /44
薄板问题的有限元法
k23 15H ab(i j )(i j ) b2 b2 k31 3Ha (2 3 5 2 ) j0 15 2 j 5i0 a a k32 15H ab(i j )(i j )
23 /44
薄板问题的有限元法
其中
b2 a2 a2 b2 k11 3H 0 15( 2 0 2 0 ) (14 4 5 2 5 2 ) 00 b b a a a2 a2 k12 3Hb (2 3 5 2 ) 0i 15 2 i 5 0i b b b2 b2 k13 3Ha (2 3 5 2 )i0 15 2 i 50 j a a a2 a2 k21 3Hb (2 3 5 2 ) 0 j 15 2 j 5 0i b b a2 k22 Hb 2(1 ) 0 (3 50 ) 5 2 (3 0 )(3 0 ) b
1 E D 2 1 0
薄板问题的有限元法
图 6.2 平板内力
10 /44
薄板问题的有限元法
设 M x 、 M y 和 M xy 表示单位宽度上的内力矩，于是有
2w 2 x Mx h h3 2 w h3 M M y h2 z dz D DC D 'C (6.5) 2 12 y 12 2 M xy 2w 2 xy

其中为薄板的弯曲刚度9899薄板的弹性曲面微分方程薄板横截面上的内力称为薄板横截面上的内力称为薄板内力薄板内力是指薄板横截面的单是指薄板横截面的单位宽度上由应力合成的位宽度上由应力合成的主矢量主矢量和和主矩主矩
平板弯曲问题的有限元分析(1) Kirchhoff弹性薄板理论
参考文献： “弹性力学(下册)”第13章。徐芝纶
x
2w
2 (z2
2
2
)dz 4
E 3 12(1 2 )
x
2w
(c)
同样，在y为常量的截面上，每单位宽度内的 y , yx , yz
也分别合成如下的弯矩，扭矩，和横向剪力：
M y
2 2
z
y dz
E
12(1
3
2
)
(
2w y2
2w x2
)
(d)
M yx
2
2
z yxdz
E 3 12(1 2 )
(9-6)
( z )z q
(f)
2
将(9-6)式代入薄板上板面的边界条件：
得：
E
12(1
3
2
)
4
w
q
(9-7)
或 D4w q， (9-8)
其中
D
E
12(1
3
2
)
(9-9)
薄板的弹性曲面微分方程
为薄板的弯曲刚度
§9-３ 薄板横截面上的内力
► 薄板横截面上的内力，称为薄板内力，是指薄板横截面的单 位宽度上，由应力合成的主矢量和主矩。
对z积分，得到： z
2(1 2 )
2
( 4
z
z2 )4w 3
F3 (x,

3、简述结点力和结点载荷的差别。 结点力：单元与单元间通过结点的相互作用力。 结点载荷：作用于结点上的外载荷。
4、列表给出有限元几类基本单元的图形、结点数、结点自由度数和单元总自由 度数(包括杆单元、梁单元、平面三角形单元、平面四边形单元、轴对称问题三 角形单元、四边形壳单元、四面体单元)。
单元 类型 杆单
（1）单元的类型和形式 为了扩大有限元法的应用领域，新的单元类型和形式不断涌现(等参元,梁板 壳,复合材料) （2）有限元法的理论基础和离散格式 将 Hellinger-Reissner、Hu—Washizu(多场变量变分原理)应用于有限元分析， 发展了混合模型、杂交型的有限元表达格式，应研究了各自的收敛条件；将加权 余量法用于建立有限元的表达格式；进一步研究发展有限元解的后验误差估计和 应力磨平方法。 （3）有限元方程的解法(大型复杂工程结构问题——静态, 特征值, 瞬态等) （4）有限元法的计算机软件(专用软件, 通用软件)
弹性力学中的虚功原理可表达为：在外力作用下处于平衡状态的弹性体，如
果发生了虚位移，那么所有的外力在虚位移上的虚功(外力功)等于整个弹性体内
应力在虚应变上的虚功(内力功)。
根据虚功原理得到 ( εT uT F )d uTTd 0
p
(1 T uT F)d 2
uT
Td
0
其中的 p 即为总势能泛函。由上面变分为零式表明：在所有区域内满足几 何关系，在边界上满足给定位移条件的可能位移中，真实位移使系统的总势能取 驻值(可证明此驻值为最小值)。此即总势能泛函的极值条件。
10, 0
3 2, 0
解：根据拉格朗日插值基函数：
u(x, y) l1(x, y)u1 l2 (x, y)u2 l3(x, y)u3 l4 (x, y)u4

04
此外，随着大数据和人工智能技术的快速发展，有限元分析可以与这 些技术相结合，实现更加智能化、自动化的工程设计和管理。
THANKS
感谢观看
离散化
将连续的物理系统划分为有限个离散的单元， 每个单元具有一定的形状和大小。
集成
将所有单元的数学方程集成为一个整体的有 限元方程组。
单元分析
对每个离散单元进行数学建模，建立单元的 数学方程。
求解
通过求解有限元方程组，得到物理系统的近 似解。
有限元的数学基础
线性代数
01
有限元方法涉及大量的线性代数运算，如矩阵运算、线性方程
定不变的载荷作用下的响应。
它主要关注的是结构的平衡状态 和位移，而不考虑时间因素和动
态效应。
静力分析广泛应用于工程领域， 如建筑、机械、航空航天等，用 于评估结构的强度、刚度和稳定
性。
静力分析的基本步骤
建立数学模型
首先需要建立结构的数学模型，包括对结构的离散化、选 择合适的单元类型和确定边界条件等。
该方法基于离散化的思想，将 复杂的结构分解为简单的、相 互连接的单元，通过求解每个 单元的平衡方程来获得结构的
整体响应。
有限元静力分析在工程领域中 广泛应用于结构强度、刚度、 稳定性等方面的分析，为结构 设计提供了重要的理论依据和 实践指导。
随着计算机技术的发展，有限 元分析软件不断涌现，为工程 师提供了更加高效、精确的数 值分析工具。
施加载荷
根据实际工况，在结构上施加相应的载荷，包括重力、外 部力、压力等。
求解平衡方程
通过有限元方法，将连续的结构离散为有限个单元，并建 立平衡方程组。然后使用数值方法求解这个方程组，得到 各节点的位移和应力等结果。

机电工程学院
第五章 板壳问题的有限元法
章节内容： 5.1 薄板弯曲的基本理论 5.2 薄板单元：矩形单元和三角形单元 5.3 薄壳有限元分析的简介
车辆工程教研室
机电工程学院
5.1 薄板弯曲的基本理论
5.1.1 薄板（thin plate）



工程实际中，存在大量的板壳构件（plate and shell） 几何特点：厚度远远小于其它两个方向的尺寸。 薄板：t/b < 1/15 中面：平分板厚度的平面 坐标系oxyz ：xy轴在中面上，z轴垂直于中面 z 载荷 作用于中面内的载荷：平面应力问题 垂直于中面的载荷：板弯曲
其中
车辆工程教研室
机电工程学院
5.5 薄壳有限元分析

局部坐标系
局部坐标系对整体 坐标系的方向余弦 矩阵（从整体坐标 到局部坐标）

局部坐标系与整体坐标系的关系
车辆工程教研室
机电工程学院
5.5 薄壳有限元分析

坐标变换矩阵
车辆工程教研室
机电工程学院
5.5 薄壳有限元分析

单元刚度矩阵
转换矩阵：
3.
应力
引起的形变很小，在计算变形时可以忽略。
车辆工程教研室
机电工程学院
5.1 薄板弯曲的基本理论
5.1.2 位移

位移分量：薄板中面的挠度 w 根据挠度，可以计算：在x和y轴方向上的位移分量和绕x和y轴方 向的转角。
y
z
b
o
车辆工程教研室
t
x
机电工程学院
5.1 薄板弯曲的基本理论
5.1.3 应变及几何方程
机电工程学院

5.1.5 平衡方程

w S1 w ，
w n

S1
（ 4-1-5）
其中 n 表示边界的法线方向。特例情况下， S1 为固支边，则 w S 0 ，
1
w 0 n S1
146
2 、 在边界 S 2 上，给定位移 w 和力矩 M n ，即
w S w ， M n
2
S2
M n
（ 4-1-6）
特例情况下， S 2 为简支边，则 w S 0 ， M n
2 2
（常曲率和常扭率）项，因为将它们代入式（ 4-1-1）可以得
2 w 2 2 4 x 2 w 2 2 6 y
2
2 w 2 5 x y
（ 4-2-3）
因此，在板弯曲单元的挠度函数中存在常数项、一次项和二次项，就可以满足完备性条件。 （2 ）协调条件： 以单元 1-2 边为例，该边上 y 为常数，挠度 w 是 x 的三次函数
如平板的表面上作用有 z 向的分布荷载 q ， 则从以上各式可以得到经典薄板理论的系统总位能泛函 表达式
1 w T D qw dxdy Q wdS M n dS b n S3 S2 S3 2 n
{M } [Db ]{ }
式中， [ Db ] 为板弯曲弹性矩阵，对于各向同性材料有
（4-1-3）
1 0 Et [D b ] 1 0 2 12(1 ) 1 0 0 2
3
（4-1-4）
为建立平板弯曲问题的能量泛函，还要考虑荷载和边界条件。关于边界条件有三种情况： 1 、 在边界 S1 上，给定位移 w 和截面转角 ，即
2 4 w T e ] z [B i ]{ z[ B ]{ } i} x y i 1

桁架结构梁单元与板壳单元结果对比分析与讨论作者：暂无来源：《智能制造》 2015年第5期撰文/ 大连华锐重工集团股份有限公司设计研究院张守云一、前言在结构有限元分析过程中，为了提高工作效率并方便处理，工程人员习惯将桁架结构简化成梁单元进行处理，以期达到事半功倍的效果。

然而，这种处理方法往往会忽略掉结构的局部细节，尤其是不同杆件连接处的节点板和加强筋板等，应用梁单元很难进行模拟。

而这些局部细节问题处理的适当与否有时也会成为影响结构强度、刚度、屈曲以及疲劳等问题分析计算的关键因素。

因此，为了更加准备地模拟结构的实际情况，有时需要将桁架结构处理成板壳单元进行有限元分析，但如此处理会在建模、求解及后处理上耗费大量的时间和精力，在产品设计和生产周期相对较为紧迫的情况下，一般很难满足快速准确地对结构进行设计、优化或验证的要求。

本文通过对同一桁架结构分别应用梁单元和板壳单元进行有限元模拟分析比较，讨论两种处理方法有限元分析结果的相同与不同之处，以期找到一种分析处理问题相对较为全面而且切实可行的解决方案。

二、桁架结构梁单元及板壳单元有限元模拟分析1. 结构有限元模型本文采用一种大型堆料机的俯仰钢结构作为分析研究对象，对其分别进行梁单元和板壳单元建模。

由于在此只是分析梁单元和板壳单元最终分析结果异同点，所以分析中所用载荷并不完全考虑结构的实际承载大小，而是进行一个假定，即结构悬臂梁部分只承受物料载荷100 吨。

2. 整体刚度分析结果比较由图1 不难看出，基于板壳单元的模型整体刚度要大于基于梁单元模型的整体刚度（前者在相同载荷作用下位移较小）。

其主要原因在于，油缸端部支承处的结构形式经梁单元简化后（图2），无论抗弯还是抗剪刚度都减小不少。

3. 强度分析结果比较（1）悬臂梁处综合应力结果比较。

图3 左侧是悬臂梁梁单元结果，最大综合应力位于悬臂梁中部的上弦杆处，最大值为83MPa。

右侧是悬臂梁板壳单元结果，最大综合应力无论数值、位置和梁单元计算结果都不一致，其最大综合应力出现在斜腹杆的节点板处，最大值为139MPa。

2、分析类型静力分析3、问题描述板壳问题4、ANSYS单位m N Pa5、单元SHELL63：4节点（每个节点6个自由度）6、材料弹性模量和泊松比7、实常数厚度8、建模取1/4模型41进入ANSYS程序→ANSYSED 9.0→ANSYS Product Launcher→change the working directory into yours →input Initial jobname:board→Run2 设置计算类型ANSYS Main Menu: Preferences →select Structural →OK3 定义材料参数ANSYS Main Menu: Preprocessor →Material Props→Material Models →Structural →Linear →Elastic→Isotropic →input EX:2e11, PRXY:0.3 →OK54 选择单元类型ANSYS Main Menu:Preprocessor →Element Type→Add/Edit/Delete→Add →select Shell Elastic 4node 63→OK (back to Element Types window) →Close (the Element Type window)65、定义实常数ANSYS Main Menu: Preprocessor →Real Constants…→Add…→select Type 1→OK→input Shell thickness at node I :0.01→OK →Close (the Real Constants Window)76、创建矩形ANSYS Main Menu:Preprocessor →Modeling →Create →Areas→Rectangle→By Dimensions →依次按下图输入：→OK87、划分单元ANSYS Main Menu: Preprocessor →Meshing →Mesh Tool →OK(Size Controls) Global: Set →input NDIV:20→Mesh Tool →Mesh :select Areas→Shape:Quad→Free→Mesh →Pick All →Close( the Mesh Tool window)98、施加固定边界ANSYS Main Menu: Solution→Define Loads →Apply→Structural →Displacement→On Lines→拾取上边和右边→OK→select Lab2:ALL DOF →OK109、施加对称约束边界ANSYS Main Menu: Solution→Define Loads →Apply→Structural →Displacement →Symmetry B.C →On Lines→拾取下边和左边→OK1111 分析计算ANSYS Main Menu: Solution →Solve →Current LS→OK(to closethe solve Current Load Step window) →OK12、挠度ANSYS Main Menu: General Postproc→Plot Results→Contour Plot→Nodal Solu→select: DOF solution, Z-Component of displacement→OK13按右侧的等侧视图141617板壳单元的压力载荷面1 (I-J-K-L)，（底部，+Z方向）面2 (I-J-K-L)，（顶部，-Z方向）面3 (J-I)，面 4 (K-J)，面 5 (L-K)，面 6 (I-L)19显示壳单元的厚度ANSYS命令菜单栏：PlotCtrls>Style >Size and Style→弹出对话框中（Display of element项设置为On)控制壳单元的结果输出ANSYS Main Menu: General Postproc→Options for OutpTop layer: 顶面Middle layer:中面Bottom layer：底面2021作业三要求1、属于力学的那类问题？2、单位制；单元类型；单元描述；实常数；材料参数3、划分网格的方案，施加载荷和边界条件4、计算结果的体现：MISIS 应力和位移，需要知道最大值以及位置。

ANSYS上机实验报告实验三：板壳的有限元分析班级：姓名：学号：一、实验题目图示正方形平板，承受垂直于板面的均布载荷作用P=20KN/m*m，板厚t=0.1m，平板外缘各边采用固定约束方式，材料选用低碳钢，弹性模量E=210GPa，u=0.33。

二、实验过程1、确定所采用的单位制：N，m，Pa。

2、问题类型：板壳问题。

3、利用ANSYS构造实体模型：1/4模型（正对称）和整体分析。

4、网格划分1）、定义材料属性：Main Menu: Preprocessor →Material Props →Material Models →Structural →Linear →Elastic →Isotropic →input EX: 210e9, PRXY: 0.33 →OK2）、定义单元类型：Main Menu: Preprocessor →Element Type→Add/Edit/Delete →Add →select Shell Elastic 8node 63 →OK (back to Element Types window)3）、定义实常数（厚度）：Main Menu: Preprocessor →Real Constants… →Add… →select Type 1→OK→input TK（I）: 0.1 ，TK（J）: 0.1 ，TK（K）: 0.1 ，TK（L）: 0.1 →OK→Close (the Real Constants Window)4）、划分网格：在size element edge length （单元边长值）处输入0.25、加载及求解。

加载（整体）过程：Main Menu: Solution →Define Loads →Apply →Structural →Displacement →On Nodes →拾取四个边线→OK →select Lab2:ALL OFF →OKMain Menu: Solution →Define Loads →Apply →Structural →pressure →On Areas →拾取面→OK →Value: 20000→OK加载（1/4）过程：Main Menu: Solution →Define Loads →Apply →Structural →Displacement →On Nodes →拾取右边线和上边线→OK →select Lab2:ALL OFF →Apply→拾取左边线→OK→select Lab2:UX→Apply→拾取下边线→OK→select Lab2:UY→OKMain Menu: Solution →Define Loads →Apply →Structural →pressure →On Areas →拾取面→OK →Value: 20000→OK求解：Main Menu: Solution →Solve →Current LS →OK(to close the solve Current Load Step window) →OK6、分析变形、位移和应力状况并抓图。

16
协调性要求 协调单元 满足协调性要求的单元称为 满足协调性要求的单元称为协调单元 收敛的充要条件 w = α1 + α 2 x + α 3 y + α 4 x 2 + α 5 xy + α 6 y 2 + α 7 x 3
+ α 8 x 2 y + α 9 xy 2 + α10 y 3 + α11 x 3 y + α12 xy 3
− 2
h
M xy = ∫ h2 τ xy zdz
− 2
h
{M } = ∫
2 −h 2
h
h {σ }zdz = [D p ]{κ } = [D ]{κ } 12
薄板弯曲的弹性矩阵
11
3
薄板弯曲的应变能 弹性应变能 T 1 1 U = ∫ (σ xε x + σ yε y + τ xyγ xy )dV = ∫ {ε } {σ }dV 2V 2V ⎧ ∂2w ⎫ ⎪ − 2 ⎪ x ⎪ ⎪ ∂2 ⎪ ∂ w ⎪ {σ } = D p {ε } = D p {κ }z {ε } = z ⎨ − 2 ⎬ = z{κ } ∂y ⎪ ⎪ T 1 ∂2w ⎪ U = ∫ {κ } [D p ]{κ }z 2 dV ⎪ 2V ⎪− 2 ∂x∂y ⎪ ⎩ ⎭ T 1 = ∫ {κ } [D ]{κ }dS
∂w 法向导数θ x = ∂y 是x的三次函数，假定
θx = γ1 + γ 2x + γ 3x + γ 4x
2
3
由节点1和节点2处只能提供 θx1,θx2 两个相邻单元在边界上的法向导数的连续性 不能保证。 这种位移函数的矩形单元为非协调单元。

西安工程大学
计算机辅助工程 CAE 讲稿 第 8 章 膜、板、壳结构的有限元法
王益轩编著 2005 年 8 月
99
k y
uk vk wk θkx θkx
o z
ui vi wi
uj vj wj θjx θjy
图 8­4
平面 3 节点 15 个自由度的三角形壳单元
8.4
SHELL63 单元描述
（8­5）
Y
O L
t
X
图 8­1 梁受弯曲力矩作用的变形（a）
西安工程大学
计算机辅助工程 CAE 讲稿 第 8 章 膜、板、壳结构的有限元法
王益轩编著 2005 年 8 月
95
y x
y u
θ
局部放 大
M
M
中性层的挠曲 线 f(x)
图 8­1 梁受弯曲力矩作用的变形 （b）
M
z x y
M
M
（a）
z M y x L M t
8.3 壳结构单元基础理论
壳单元能承受拉伸应力与弯矩力，也就是膜单元与板单元合并后的单元，以三角形单元为例来说， 单元的应变位移关系式（几何关系）与三角形板单元完全相同，只不过每个节点具有5 个节点位移（三 个平动位移，两个转角） ，即 5 个自由度，单元共有 15 个自由度，假设位移模式中应包含 15 个任意常 数。单元图形如图 8­4 所示。推导过程从略，由读者自己推导。
ANSYS 结构分析中的板壳单元 SHELL63
8.4.1
图 8­5 SHELL63 板壳单元
SHELL63 称为弹性壳，因为它只支持线性弹性的材料模式； ANSYS 另有其它 shell 单元可以支持更 广泛的材料模式。SHELL63 有 4 个节点（I, J, K, L） ，每个节点有 6 个自由度：3 个位移（UX, UY, UZ） 及 3 个转角（ROTX, ROTY, ROTZ） ，所以一个单元共有 24 个自由度。若 K、L 两个节点重迭在一起时， 它就退化成一个三角形，如图 8­5 右图所示。I­J­K­L 四个节点假设是共平面，若不共平面则以一最接近 的平面来修正这四个节点。注意，这种修正当然会引进一些误差，所以对那种曲率很大的板壳结构而 言， 必须使用较细的单元。 SHELL63 的单元坐标系统表示在图 8­5 中， 原点是在 I 节点上， X 轴和 I­J 边可以有一角度差 （THETA，

有限元分析论文写作范文（专业推荐6篇）车架作为汽车的承载基体，安装着发动机、传动系、转向系、悬架、驾驶室、货厢等有关部件和总成，承受着传递给它的各种力和力矩。

车架工作状态比较复杂，无法用简单的数学方法对其进行准确的分析计算，而采用有限元方法可以对车架的静动态特性进行较为准确的分析，从而使车架设计从经验设计进入到科学设计阶段。

以下是我们为你准备的6篇有限元分析论文，希望对你有帮助。

有限元分析论文范文第一篇：油罐运输车的有限元分析及优化摘要：为验证油罐运输车的结构强度是否满足使用要求，运用有限元仿真分析方法分别建立其弯曲、扭转、紧急制动3种工况的模型并进行了最大应力分析。

结果显示，罐体结构的应力小于材料的屈服应力，在满足使用要求的基础上，采用尺寸优化分析方法减薄罐体的厚度可实现轻量化。

关键词：油罐运输车；有限元分析；尺寸优化伴随着世界经济持续发展，石油、天然气的需求逐步增加，油罐车作为短途运输交通工具发挥着重要的作用。

存在部分结构不合理和整车质量过重现象及潜在运输的危险性，同时使得运输成本增加。

因此基于CAD/CAE技术对整车进行结构分析与轻量化设计,可以提高产品的科技含量，为企业以后的生产提供设计指导。

1罐车有限元模型的建立1.1单元类型的选择罐体单元主要采用单元类型中的壳单元来划分网格，车架部分由于用梁单元不能分析应力集中问题，所以同样采用壳单元来划分车架网格，这样可以准确地得出分析结果。

罐体的单元选用四边形壳单元（QUAD4），在几何形状复杂的位置可以采用少量的三角形单元（TRIA3）来过渡，以满足总体网格质量的要求，通常要求三角形单元占总单元数的比例不超过5%【2】.罐体以及车架的单元全部为10mm尺寸单元。

1.2罐体与车架连接方式罐体与前后封头、罐体与防波板以及加强板与相应连接部件之间用节点耦合的方式模拟焊接。

大梁与副车架之间的连接采用ACM单元。

ACM单元模拟的是一种特殊的焊接方法（AreaContactMethod），不同于刚性单元结点连接的方法。

板栗真空破壳力学特性有限元分析袁越锦;袁月定;党新安;徐英英;董继先;徐君臣【摘要】Taking the hemispherical chestnut as the study object, a finite element model for the mechanical properties of vacuum shelling chestnut was established by applying the finite element method and transport-process theory. The simulation and experimental results indicated that the finite element model could explain the mechanical properties of vacuum shelling chestnut well. The moisture distribution of chestnut shell was higher in the bottom center area and lower in the shell edge. The value of stress and strain was decreasing from the shell edge to the bottom center and arc peak. The wet stress played an important role in the vacuum shelling process, and was about five times larger than the pressure difference stresses in the initial and middle stage, but the wet stress decreased gradually at the end of the shelling process. The vacuum degree had a great effect on the vacuum shelling process, and the increasing vacuum degree did not only increase the pressure difference stresses but also the wet stresses.%以半球形板栗为研究对象,运用有限元方法和传递过程原理建立了板栗真空破壳的有限元模型,并采用间接顺序耦合法对模型进行了求解.模拟与试验结果表明:运用该模型模拟分析板栗真空破壳力学特性可行；真空破壳过程中应力和应变主要沿栗壳边缘分别向弧顶和壳底中心递减分布；湿应力是真空破壳的主导因素；绝对压力对破壳的影响很大,绝对压力的降低不仅使压差应力增大,同时也使湿应力增大.【期刊名称】《农业机械学报》【年(卷),期】2011(042)005【总页数】6页(P136-141)【关键词】板栗;真空破壳;力学特性;热质传递;有限元【作者】袁越锦;袁月定;党新安;徐英英;董继先;徐君臣【作者单位】陕西科技大学机电工程学院,西安710021;宜春学院数学与计算机科学学院,宜春336000;陕西科技大学机电工程学院,西安710021;陕西科技大学机电工程学院,西安710021;陕西科技大学机电工程学院,西安710021;陕西科技大学机电工程学院,西安710021【正文语种】中文【中图分类】TS255.36引言板栗脱壳是板栗深加工的头道工序,其难点主要是板栗外壳韧硬坚实且没有结合缝。