固体力学线弹性问题有限元分析课件.pptVIP

第五讲固体力学-线弹性问题有限元分析元计算技术部

性力学作固体力学的一个重要分支，研究性物体在外力和其他外界因素作用下生的形和内力，它是材料力学、构力学、塑性力学和某些交叉学科的基。广泛用在建筑、机械、化工、航天等工程域。本将分支，从其物理模型，有限元弱形式推，以及ELAB.1.0有限元分析、ELAB1.0有限元件公式等各个方面行介。基本方程ELAB1.0模型向导实现有限元脚本文件分析

?线弹性问题的基本方程从静力学、几何学和物理学方面考虑得到线弹性稳态问题对应的数学物理方程为：平衡方程：本构方程：几何方程：其中σ、σ、σ表示直角坐标系下三个方向的正应力，ε、ε、ε表示对应的正应变，yyxxyyzzxxzzσ、σ、σ表示直角坐标系下三个剪应力，ε、ε、ε表示对应的剪应变，u、v、w表示xyxzyzxyxzyz直角坐标系下三个方向的位移，E表示杨氏模量，ν表示泊松比。

边界条件：第一类边界条件:第二类边界条件：第三类边界条件：?有限元分析运用迦辽金有限元法求位移，由上面的平衡方程可得：

其中σu、σv、σw表示三个方向的虚位移。对上式进行分部积分化为弱形式可得：将本构方程带入到上面的弱形式，得到求解位移的最终弱形式表达式：

对于弹性体的应力，采用最小二乘法，由线弹性问题的本构方程可以得到如下的弱形式：线弹性问题属于固体力学中基础的学科分支，在ELAB1.0有限元软件中以公式库的形式提供给大家，因此可以采用【公式库-固体力学-线弹性】直接生成的方式生成程序代码，下面通过一个算例用ELAB1.0公式库来实现。

?工程背景三维工字形部件线弹性体，如下图所示，底面为边长为8m的正方体，上下两部分高度为2m，中间部分高度为10m。该部件的弹性模量为1.0e10N/m2，泊松比为0.3，地面边界固定，上表面施加100N的均布力载荷，分析该部件的位移、应力以及变形情况。几何模型

u固体力学--线弹性ELAB1.0软件实现?工程建模3、“坐系”1、点“工程向”入公式2、“固体力学”研究域

5、“型”6、定义工程名和工程路径，完成工程设置4、“元型”

?定义材料参数点工具“参数置”→“材料参数”，如下所示：材料参数框中定相的材料参数，如下所示：b场材料参数a场材料参数

?前理点击工具栏中“前处理”按钮进入GID，建立该工程的几何模型。注：入GID后要行ELAB1.0的数据化data→problemtype→ELAB几何建模：建立几何模型的具体操作详见《有限元分析基础与应用》相关章节。注：模型建立后，Geometry——Edit——Collapse——models，中所建模型，按鼠中束，将所有的体一体。保没有孤立的点、或者面。

?添加材料参数和边界条件：b场材料施加均布力地面固定边界a场材料

?划分网格：设置划分网格的单元类型(要与工程建模中选择的单元类型一致)，以及网格尺寸，划分网格如下图所示：?求解计算点击工具栏中“求解计算”按钮，完成模型的求解计算。

?后处理点击工具栏中的“后处理”按钮进入GID，查看计算结果。x向位移uy向位移v

z向位移w位移矢量x向应力dxxy向应力dyy

z向应力dzz变形图（放大2.1392e6倍）

?有限元语言描述文件为生成该问题有限元计算的所有程序源代码，针对之前的ELAB1.0有限元分析得到的微分方程弱形式，ELAB1.0软件提供简洁的有限元语言描述文件，包括微分方程描述文件、多物理场描述文件以及求解命令流控制文件。针对该问题的有限元描述文件包括delxyz.fde(求解位移微分方程)，delxyz.fbc(求解位移边界条件)，selxyz.fde(求解应力微分方程)，solid.mdi，solid.gcnü微分方程描述文件delxyz.fde(求解位移微分方程)在delxyz.fde给出单元的待求未知量，涉及到的材料参数，单元的形函数表达式，刚度矩阵表达式和载荷表达式，以及为描述刚度矩阵和载荷向量而自定义的函数。以下给出微分方程描述文件中与微分方程弱形式对应的部分(详细的解析见《有限元分析基础和应用》中相关章节)：

微分方程弱形式：ü求解位移的微分方程描述文件未知变量：DISPuvw未知变量对应微分方程弱形式中的变量(几何方程中)uvw

材料参数：MATEpepvfxfyfzroualpha材料参数行对应微分方程弱形式中的变量Evfxfyfz单元刚度矩阵：dist=+[ev_i;ev_j]*sm_i_j*fact+[ep_i;ep_i]*shear*fact单元刚度矩阵对应微分方程弱形式中的左端项

单元载荷向量：load=+[u_i]*f_i*vol单元载荷向量对应微分方程弱形式中的右端项的体积分项ü求解位移