第2章有限元理论基础

第2章有限元理论基础以结构计算为例。有限元计算(ANSYS软件)的理论基础－基本方程，边界条件。基本方程：描述应力状态的平衡方程描述应变状态的几何方程－－－－－有限元计算的核心思想。描述应力应变关系的本构方程对应的边界条件。2.1应力状态分析图2.1为单元体的应力状态。图2.1单元体的任一点的应力状态描述：333231232221131211ij，剪应力互等jiij，六个独立分量单元体的静力平衡问题。单元体沿三个坐标轴方向的力的平衡条件和对三个轴的力矩平衡条件。三维力的平衡微分方程：01312111Fzyx0jiijFx0iijx02322212Fzyx0,jiijF0,iij03332313Fzyxj=1,2,3j=1,2,3note:1.11在垂直x轴平面的应力，在X轴的分量。2.F为体力，包括：重力、磁力、惯性力，与物体的质量成正比。Fi为I轴的体力分量。3.物体表面单位面积的面力T三个分量为Tx,Ty,Tz,或T１,T２,T３,应力的三个分量x,y,z或1,2,3应力边界条件：332313232221223121111TTT3,2,1iTTijijii表达作用在物体表面单位面积丧的面力T与物体内的应力分量之间的关系。2.2应变状态分析图2.2为单元体的应变状态。图2.2单元体的一点的应变状态的张量描述：333231232221131211zzzyzxyzyyyxxzxyxxij与应力状态相似。由剪应变互等定律jiij九个应变分量中只有六个是完全独立的。描述一点的应变分量与点的位移分量之间的关系的几何方程。zwyvxuzyxxwzuzvywyuxvxzzxzyyzyxxy如果用张量表示：定义：zxxzzxyzzyyzxyyxxyzzzyyyxxx21,21,21,,,)(21)(21)(21211,33,1313,22,3232,11,22112123,3332,2221,11111uuuuuuxuxuuuuxuzxyzxyzzyyxx)(21,,ijjiijuu2.2应力－应变关系分析－－用于定义材料的性质、设定边界条件。2.3.1弹性材料弹性理论研究的对象是理想的弹性体。符合下述四个基本假设：(1)物体是连续的,由连续介质组成,没有空隙，物体中的应力、应变以及位移等量是连续的，可以用坐标的连续函数表示。(2)物体是均匀的，各部分具有相同的性质，弹性常数不随位置坐标而改变。(3)物体各向同性、弹性常数不随方向而改变。(4)物体是完全弹性的。线弹性问题是以理想弹性体的微小位移和应变为前提条件。(略去高次项成线性方程)。反之为非线性弹性理论。处于(弹)塑性状态的应力－应变关系理论为(弹)塑性理论。线弹性理论的广义虎克定律描述弹性体内任一点的应力－应变分量之间的关系。zxzxyzyzxyxyyyxxzzzzxxzzyyyyzzyyxxxxGGGEEE212121)]([1)]([1)]([1ijijijEE1式中：332211iiij单位张量ijij21E杨氏模量；v泊松比；G为剪切模量)1(2EG上述为物理方程，或本构方程。用应变表示应力的形式：zxzxyzyzxyxyzzzzyyyyxxxxeee222222ijijije2基本方程通解＋定解条件完整的问题初始条件――讨论问题的特定时刻（动态问题）边界条件――弹性体边界上的外力和位移的情况。弹性力学中的边界条件：(1)应力边界条件已知作用在物体内的体力Fi和物体表面处的面力Ti,就已知物体边界上的应力。(2)位移边界条件已知作用在物体内的体力Fi和物体表面处的位移ui,已知物体边界上的位移。(3)混合边界条件已知作用在物体内的体力Fi，物体表面处的位移ui和力已知。解的存在定理和解的唯一行定理：弹性力学问题的解是存在的。在小变形条件下，受一组平衡力系作用的物体的应力和应变的解是唯一的。平面问题：一般都是三维问题，求解、分析是相当复杂和困难的。在特殊情况下，如果能够把它作为平面问题，大大简化了计算和分析的过程。WHAT?平面问题：薄板一类的结构。包括平面应力和平面应变问题。平面问题的基本方程和边界条件平衡方程：00yyyxyxyxxxFyxFyx几何方程：)(21yuxvyvxuyxxyyyxx边界条件：力2121nnTnnTyyxyyyxxxx位移：u=u1v=v1平面应力问题：等厚度、板厚度远小于板的其它两个方向的尺寸，所受的外力平行于板面且不受厚度变化。满足上述条件，应力状态是平面的，二维的；应变状态是空间的，三维的。弹性体平面应力问题：应力状态描述：00000yyxyyxxxij应变状态描述：zzyyxyyxxxij0000本构方程：0121)(1)(1)(1zxyzxyxyxyyyxxzzxxyyyyyyxxxxEGEEE平面应变问题：沿z轴方向的尺寸远大于其它另外两个方向的尺寸。Z轴方向的位移或应变相对于本身尺寸很小，可忽略不计。满足上述条件的应变状态是平面的，二维的；其应力状态是空间三维的。弹性体的平面应变问题：应力状态描述：zzyyxyyxxxij0000应变状态描述00000yyxyyxxxij本构方程：xyxyxxyyyyyyxxxxGEE21)(1)(1''''式中：2'1vEEvvv1'轴对称问题：弹性体的几何形状和所受外力均对称于弹性体的中心轴。典型的有：柱、筒、环等。特征：剪应力为零0rr，只有正应力CrBrArCrBrArrrrr2)ln23()(2)ln21()(22轴对称平面应力问题：平行板面的外力和几何尺寸沿z轴对称，且z轴的尺寸远远小于x,y方向的尺寸，如薄板圆环。轴对称平面应变问题：所受的外力和几何尺寸沿z轴对称，沿z轴方向的尺寸远大于其它另外两个方向的尺寸或两端受约束。Z轴方向的位移或应变相对于本身尺寸很小，可忽略不计。两种情况的应力分量表达式一样：02222rrrrCrACrA平面应力问题的应变分量表达式：02])1(2)1(2)3()1([1])1(2)1(2)31()1[(122rrrrrCvnrBvBvrAvECvnrBvBvrAvE平面应变问题的应变表达式：将E换为21'vEE和v换为vvv1'。2.3.２弹塑性材料应力较小时，材料处于弹性变形阶段，应力－应变线性关系，服从虎克定律，当应力超过材料的屈服极限s时，材料进入塑性变形阶段，总应变ij有弹性应变eij和塑性应变pij两部分组成。pijeijij应力－应变曲线如图。在塑性阶段不仅与应力状态有关，而且与加载路径有关。介绍几种常用的弹塑性理想（简单）模型(1)理想的刚塑性模型材料在屈服前处于刚性无变形状态，一旦屈服后进入不强化的理想塑性流动状态。应力－应变关系见图２.1。图2.1其表达式：s(2)理想的弹塑性模型弹性阶段的应力－应变关系是线性的，屈服后是不强化的理想塑性流动状态。应力－应变关系见图2.2。s0s0s图2.2其表达式：sEss以上为理想塑性模型。(3)刚塑性线性强化模型屈服前为刚性无变形状态，屈服后为线性强化。应力－应变关系如图2.3。图2.3其表达式：1Es(4)弹塑性线性强化模型－－－双直线模型弹性阶段是线性的，屈服后线性强化。应力－应变关系如图2.4。图2.4s0arctgEs0arctgE其表达式：)(1ssEEss(5)幂强化模型应力－应变遵从幂指数关系，nA式中：A0n为强化系数。n=0理想刚塑性材料；n=１理想线弹性材料。见图2.5。图2.5(6)多阶段模型为了描述更加接近实际材料的情况。大多数材料是非理想的材料，其普遍的应力应变关系遵从在弹性段是线性的，符合虎克定律。屈服后发生塑性变形，应力应变规律是复杂非线性的，为研究方便，将塑性阶段分段分析。如图2.6。s0ArctgE1ArctgE2ArctgE3图2.6(7)增量理论模型塑性状态下，应力－应变关系的本质－－－增量关系。这是一种更接近实际塑性材料的模型。Reuss假设：'ijpijdd2.3.３塑性材料塑性问题的求解是非常复杂的。不少实际问题可用刚塑性材料的平面应变问题来近似―――采用的滑移线场理论。塑性平面应变问题，则是指物体内个质点的塑性流动均平行与XOY平面，与z无关。0),,(),,(wyxvvyxuu基本方程应力平衡方程：00)(41222yxyxkyyxyyyxxxyyyxx应变几何方程：0,0,,zxyzxyxxyyxxyvxuyvxu滑移线场理论中常遇到的应力边界条件：（1）应力自由表面00xyyy（2）无摩擦接触面0xy一些特殊情况查书。2.3.４粘弹性材料固体同时显示出弹性和粘性的行为。以线弹性和线粘性为例讨论。当开始加载时，初始材料显示弹性性能，符合虎克定律E。在加载时形变随时间而增大的现象，卸载后形变继续保留下来，应力－应变关系符合牛顿定律。dtd线性粘弹性材料的力学模型：Maxwell模型：一个弹簧和一个阻尼器串联。近似描述粘弹性固体的应力松弛行为。E讨论：当0,0，得：tEet0)(应力松弛――应力随时间减小。当0得：tt0)(蠕变――应力随时间增大。Kelvin模型：一个弹簧和一个阻尼器并联。近似描述粘弹性固体的蠕变行为。E讨论：0)1()(0tEeEt蠕变――应力随时间增大0,0E弹性行为――虎克定律标准线性固体模型：Maxwell模型和另一个弹簧并联。同时描述应力松弛和蠕变现象的模型。)(RM讨论：当0,0tERReMMt)()(000应力松弛当0tERReMMt)()(000应变松弛－蠕变粘弹性固体的动态力学行为内耗：tgWWQsin211标准线性固体的动态力学模型：2211MQMMMRuM前面讲的是由外力引起的应力―――结构应力。下面讲由热――温度变化引起的应力――热应力。热应力引起热冲击破坏（温度骤变）和热疲劳破坏。―――高温场合尤为重要。2.5热应力2.5.１热应力的产生构件的热胀冷缩受到外界约束。（棒两端固定，然后加热，产生压应力）相连的两个构件存在温差，且相互约束。（温差引起热传导，构件热胀冷缩，相互约束引起应力。）构件内部存在温差，且相互约束。（构件加热－冷却不均匀热处理，内部