非线性有限元 第三篇 章 连续介质力学.pptVIP

非线性有限元第3章 连续介质力学 ;第2讲 连续介质力学 ;1 引言 ;1 引言 ;1 引言 ;2 变形和运动 ;2 变形和运动 ;2 变形和运动 ; 空间变量 x 和时间 t 的任何函数的材料时间导数可以通过链规则得到;变形梯度－是运动函数的Jacobian矩阵 ;2 变形和运动 ;2 变形和运动 ;2 变形和运动 ;运动条件;2 变形和运动 ;例3.1;三角形3节点线性位移单元的构形 ;将未变形构形中的节点坐标代入上式 ; 在单元中的位移是材料坐标的线性函数，变形梯度仅为时间函数，若给定时间，F 为常数。Jacobian行列式给出为;例3.3;沿着由节点1和2以及节点1和4所定义的边界上位移场为零，运动为 ;例3.4;解：;3 应变度量;一个单元绕着原点转动了θ角。计算线性应变 ;3 应变度量; 线性分析的适用性则在于能够容许误差的量级，最终取决于感兴趣的应变的大小。如果感兴趣的应变量级是10-2，那么1％的误差是能够接受的（几乎总是这样）。如果感兴趣的应变更小，可接受的转动更小，对于10-4量级的应变，为满足1％的误差，转动必须是10-3 弧度量级的。 这些指导数据假设平衡解答是稳定的，即不可能发生屈曲。然而，屈曲是可能的，即使是在很小的应变下，所以当可能发生屈曲时，应该使用能适合应付大变形的度量。 ;3 应变度量;3 应变度量;3 应变度量;第二个运动度量D，称为速度应变, 是变形的率度量, 定义速度梯度 ;没有变形，转动张量和角速度张量相等：W＝Ω。 由速度梯度定义，在刚体运动中变形率D＝0，所以L＝W ，积分;证明在刚体运动中变形率D＝0;3 应变度量;3 应变度量;例3.5 拉伸和转动联合作用下的应变度量 ;得到Green应变张量;为了确定变形率的时间历史，计算变形梯度的时间导数和逆 ; 一个单元经历了图示的变形阶段。在这些阶段之间的运动是时间的线性函数。计算每一阶段的变形率张量D，对于回到未变形构形的整个变形循环，获得变形率的时间积分。 ;得到速度梯度和变形率为 ;例3.7 计算变形率的时间积分 ;例3.7 计算变形率的时间积分 ;4 应力度量;4 应力度量;4 应力度量;4 应力度量;在物体中的每个点都构造了一个坐标系。这个坐标系随着材料或单元一起转动。通过将这些张量表达在一个随材料而转动的坐标系中，很容易处理结构单元和各向异性材料。 ;4 应力度量;旋转Cauchy应力和旋转变形率定义为;例3.8 平面问题 ;例：平面问题 ;5 守恒方程;5 守恒方程;5 守恒方程;5 守恒方程;5 守恒方程;质量守恒方程的几种特殊形式 ;5 守恒方程;5 守恒方程;4 能量守恒 ;5 守恒方程;5 守恒方程;5 守恒方程;5 守恒方程;5 守恒方程;6 Lagrangian守恒方程;6 Lagrangian守恒方程;6 Lagrangian守恒方程;7 极分解和框架不变性;极分解定理 ;7 极分解和框架不变性;7 极分解和框架不变性;在面积坐标的形式下，运动描述为 ;7 极分解和框架不变性;7 极分解和框架不变性;7 极分解和框架不变性; 公式的不足在于它不能解释材料的转动。通过应力张量的客观率可以解释材料的转动；称为框架不变率。 考虑三种客观率： Jaumann率， Truesdell率， Green－Naghdi率。 框架不变性的核心是应力的(变化)材料导数不受刚体位移的影响。 所有这些都应用于当前的有限元软件中，如ABAQUS。还有许多其它的客观率将在后面讨论。;黄先生书描述固体本构大变形给出3种定义： 1 SO(Simo-Ortiz)定义来自于Green－Naghdi率本构模型，只不过将后者从参考构型前推到卸载构形(令温度和结构不变，应力全部卸除后的残余变形，也称为中间构形)和当前构型； MOS(Moran-Ortiz-Shih)本构理论来自于Jaumann率，将变形张量分解为对称(平动)和反对称部分(转动)。在中间构形建立本构关系，把中间构形中的Green应变率定义为弹性变形率D，dE/dt＝D既反映了当前构形、也反映了中间构形的变化。 RH(Rice-Hill)与SO的差别是不分别定义Green应变的弹性和塑性部分，而是分解Green应变率为弹性和塑性部分。 Cauchy应力与Jaumann率构成ABAQUS的核心部分。;Cauchy应力的Jaumann率 ;7 极分解和框架不变性;7 极分解和框架不变性;7 极分解和框架不变性;7 极分解和框架不变性;7 极分解和框架不变性; 在刚体运动中初始应力（或预应力）嵌入在固体中，即在刚体转动中由于变形没有发生变化，观察者所看到的随着物体运动的应力也不应该变化。 单元坐标系：在刚体转动中，Ja