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

(a)直角坐标系与实际单元(b)自然坐标系与基本单元图7-320结点三维等参单元位移函数和几何坐标的变换式应取为相同的参数，其坐标变换关系可表示为(7-17)则单元的位移函数可写成第7章空间问题的有限单元法第30页，共60页，星期日，2025年，2月5日(7-18)在自然坐标系中，各结点的形状函数可写成如下形式对于8个顶角结点（i＝1，2，……，8）式中xi、yi、zi——结点i的坐标；ui、vi、wi——结点i沿x、y、z方向的位移；Ni——对应于i结点的形状函数。对于的边上点（i＝9，11，13，15）对于的边上点（i＝10，12，14，17）第7章空间问题的有限单元法第31页，共60页，星期日，2025年，2月5日对于的边上点（i＝17，18，19，20）(7-19)7.4.2单元刚度矩阵三维变形状态下，一点的应变与位移的几何关系为(7-20)第7章空间问题的有限单元法第32页，共60页，星期日，2025年，2月5日式中为单元结点位移列向量，而单元应变转换矩阵[B]可按结点分块表示为其中每个子矩阵又可分为上下两块，有（i＝1，2，…20）(7-21)单元的刚度矩阵为第7章空间问题的有限单元法第33页，共60页，星期日，2025年，2月5日(7-22)为便于以下计算，弹性矩阵[D]可分块写为(7-23)令,则，为60×60的方阵，可按结点写为子块形式第7章空间问题的有限单元法第34页，共60页，星期日，2025年，2月5日式中第i行j列的子矩阵为(7-24)将(7-20)、(7-22)分块式代入(7-23)，其被积函数可写为(7-25)式中按坐标变换式(7-17)，应有第7章空间问题的有限单元法第35页，共60页，星期日，2025年，2月5日同样可有(7-27)三维7面体的雅可比矩阵为(7-27)同理可采用三维高斯求积公式计算单元刚度矩阵。即第7章空间问题的有限单元法第36页，共60页，星期日，2025年，2月5日式中，L，M，N为沿、、方向的积分点数目，而积分点坐标及权重可由高斯积分表查得。7.5ANSYS桁架结构计算示例7.5.1问题描述如图7-4所示，一个圆柱实体。柱高0.2m，圆柱横截面直径为0.1m。约束方式：底面全约束。承受载荷：第7章空间问题的有限单元法第37页，共60页，星期日，2025年，2月5日A点承受Z方向集中载荷Fz=5000N和Y方向集中载荷Fy=-5000N；B点承受X方向集中载荷Fx=5000N；C点承受Z方向集中载荷Fz=-5000N；D点承受X方向集中载荷Fx=-5000N。弹性模量为EX=210GP，泊松比μ=0.3。图7-4圆柱实体示意图7.5.2ANSYS求解操作过程打开Ansys软件，在Ansys环境下做如下操作:第7章空间问题的有限单元法第38页，共60页，星期日，2025年，2月5日（1）选择单元类型运行PreprocessorElementTypeAdd/Edit/Delete，弹出ElementTypes对话框，如图7-5所示。然后单击Add，弹出LibraryofElementTypes窗口，如图7-7所示，选择SOLID45单元，单击OK。图7-5单元类型对话框图7-7单元类型库对话框第7章空间问题的有限单元法第39页，共60页，星期日，2025年，2月5日（2）设置材料属性运行PreprocessorMaterialPropsMaterialModels，弹出如图7-7所示对话框。双击Isotropic，弹出LinearIsotropicPropertiesforMaterialNumber1对话框，如图7-8所示，在EX