空间问题的有限单元法.pptVIP

空间问题的有限单元法演示文稿;（优选）空间问题的有限单元法;7.1三维应力状态;由弹性力学知，应变与位移间的几何关系是;弹性体受力作用，内部任意一点的应力状态也是三维的，用列向量表示为;7.2空间结构的离散化;其中，最简单的空间单元是四面体单元。采用四面体单元和线性位移函数处理空间问题，可以看作是平面三角形单元的推广。

如图7-2(f)所示，一个平行6面体可由5个四面体组成，其基本单元仍是四面体。它们分别由如下结点组成：

2→3→4→7，1→2→4→5，2→4→7→5，2→7→6→5，

4→7→5→8。;;由式(7-7)求出、、和，再代回式(7-5)

中，整理后得;为使四面体的体积V不为负值，在右手坐标系中，使右手旋转按着由i-j-m的转向转动时，是向法向n方向前进。

用求位移u的同样方法，可求得;7.3.2单元刚度矩阵;同样，用虚功原理建立结点力和结点位移间的关系式，从而得出简单四面体单元的刚度矩阵。;(r=i、j、m、n，S=i、j、m、n)(7-15);常应变四面体单元;体积坐标;一、体积坐标定义;二、体积坐标与

整体坐标的关系;;高次四面体单元;高次四面体单元;高次四面体单元;高次四面体单元;五面体单元;五面体单元;六面长方体单元;六面长方体单元;六面体等参元;;;(a)直角坐标系与实际单元(b)自然坐标系与基本单元;(7-18);;;;式中第i行j列的子矩阵为;;;A点承受Z方向集中载荷Fz=5000N和Y方向集中载荷

Fy=-5000N；B点承受X方向集中载荷Fx=5000N；C点承受Z方向集中载荷Fz=-5000N；D点承受X方向集中载荷

Fx=-5000N。弹性模量为EX=210GP，泊松比μ=0.3。;（1）选择单元类型

运行PreprocessorElementTypeAdd/Edit/Delete，弹出ElementTypes对话框，如图7-5所示。然后单击

Add，弹出LibraryofElementTypes窗口，如图7-7所

示，选择SOLID45单元，单击OK。;（2）设置材料属性

运行PreprocessorMaterialPropsMaterialModels，弹出如图7-7所示对话框。双击Isotropic，弹出Linear

IsotropicPropertiesforMaterialNumber1对话框，如图7-8所示，在EX选项栏中设置数值2.1e11，在PRXY选项

栏中设置数值0.3。设置完毕单击OK按钮。;（3）建立模型

运??PreprocessorModelingCreateAreasRectangle

By2Corners，弹出如图7-9所示对话框，在WPX选项栏中填写0，在WPY选项栏中填写0，在Width选项栏中填写0.05，在Height选项栏中填写0.2，点击OK。生成如图

7-10所示图形。;将长方形旋转成柱体，运行

PreprocessorModelingOperate

ExtrudeAreasAboutAxis，弹出如图7-11所示拾取框。选择图7中长方形后弹出单击OK，再选择长方形左上角和左下角结点后，单击OK.。弹出如图7-12所示对话框。在ARC选项栏中填入旋转角度360度，设置完毕单击OK按钮，生成如图7-13所示

圆柱体。;运行MeshingSizeCntrlsManualSizeGlobalSize弹出如图7-14所示对话框，设置SIZE选项栏中的数据为0.01。运行MeshingMeshVolumesFree自由划分网格后得到如图7-15所示图形。;（7）施加载荷

显示图形的关键点，运行PlotCtrlsNumbering弹出如图7-17所示对话框，激活KPNumbers后面的选框，使它变成on形式。选择菜单SolutionDefineLoadsApply;图7-17圆柱实体示意图;图7-18云图显示对话框;图7-19位移云图图7-20应力云图;圆筒直径0.4m，高度0.6m，壁厚0.005m；材料Q235，弹性模量E=2.1e11Pa，泊松比μ=0.3；约束：圆筒的下部在轴线方向固定，其它方向自由；载荷：顶部环线上承受轴向线压力P＝-200000N/m。

;