您好,欢迎访问三七文档
当前位置:首页 > 临时分类 > 有限元理论与方法-第7讲
1讲授内容备注第7讲(第7周)2.应变矩阵确定了单元位移后,可以很方便地利用几何方程和物理方程求得单元的应变和应力。作为平面问题,单元内具有3个应变分量εx、εy、γxy(各符号的意义见附录1),用矩阵表示为xvyuyvxuxyyxε将(2.1.4)式代入上式中,得到emmjjiimjimjibcbcbccccbbbAδε00000021或eδBε(2-1-7)式中B称为应变矩阵,写为分块形式,即B=[BiBjBm](2-1-8)而其子阵为),,(0021mjibccbAiiiiiB(2-1-9)3节点三角形单元的B是常量阵,所以称为常应变单元。在应变梯度较大(也即应力梯度较大)的部位,单元划分应适当密集,否则将不能反映应变的真实变化而导致较大的误差。上述应变中包括与应力有关的应变和与应力无关的应变两部分,无关的应变ε0又称为初应变0000xyyxεε0由温度变化、收缩、晶体生长等因素引起,对工程结构一般只考虑温度应变,无论线性和非线性温度,计算时可近似地采用平均温度33refTTTTTmji式中,Ti、Tj、Tm分别为节点i、j、m的温度,Tref为参考温度。对于平面应力问题,温度T引起的初始应变为200TTε其中,α为线膨胀系数。由于温度变化在各向同性介质中不引起剪切变形,所以γxy0=0。以后所述问题,除非特别说明,都指各向同性介质。对平面应力问题,温度T引起的初始应变为0)1(0TTε当不考虑温度的影响时,当前温度即为参考温度。以后所述问题,除非特别说明,不考虑温度影响。3.单元应力根据物理方程,对平面应力问题,取应变分量xyxyxyyyxxEEEEE)1(2由上式解出xyxyxyyxyyxxEEEE211)1(2)(1)(1222εDσxyyx(2-1-10)式中,D为弹性矩阵,2100010112ED(2-1-11)取决于弹性常数E和μ。将式(2-1-7)代入式(2-1-10)得eδSσ(2-1-12)mjiSSSBDS(2-1-13)3),,(2121)1(22mjibccbcbAEiiiiiiiS式中,S为应力矩阵,反映了单元应力与单元节点位移之间的关系。由于单元应力和应变分量为常量,所以单元边界上有应力阶越,随单元划分变密,突变将减小。对平面应变问题,有四个应力分量:σx、σy、τxy和σz。取应变分量01)1(211yxzzxyxyzxyyzyxxEEEE由应变分量解出σx、σy、τxy,弹性矩阵为)1(22100011011)21)(1()1(ED(2-1-14)根据物理方程可以求解各应力分量。4.单元刚度矩阵单元节点力为Fe,节点虚位移为(δ*)e,节点虚应变为(ε*)e,平面单元的厚度为t。应用虚位移原理yxtyxtAAeeddddT*T*T*DεεσεFδ将e**)(δBε及eδBε代入上式整理得到eAeyxtδDBBFddT可见单元刚度矩阵为yxtAeddTDBBK(2-1-15)对于三节点三角形单元,面积为A,所取为线性位移模式,单元刚度矩阵为tAeTSBK进一步表示为mmmjmijmjjjiimijiieKKKKKKKKKK对平面应力问题有4mjismjirbbcccbbcbccbccbbAEtsrsrsrsrsrsrsrsrrs,,;,,21212121142K(2-1-16)单元刚度矩阵表达单元抵抗变形的能力,其元素值为单位位移所引起的节点力,与普通弹簧的刚度系数具有同样的物理本质。例如子块Kij22211211ijijijijrskkkkK其中:上标1表示x方向自由度,2表示y方向自由度,后一上标代表单位位移的方向,前一上标代表单位位移引起的节点力方向。如11ijk表示j节点产生单位水平位移时在i节点引起的水平节点力分量,21ijk表示j节点产生单位水平位移时在i节点引起的竖直节点力分量,其余类推。单元刚度矩阵为对称矩阵。由于单元可有任意的刚体位移,给定的节点力不能惟一地确定节点位移,可知单元刚度矩阵不可求逆,具有奇异性。5.等效节点载荷有限单元法分析只采用节点载荷,作用于单元上的非节点载荷都必须移置为等效节点载荷。可依照静力等效原则,即原载荷与等效节点载荷在虚位移上所作的虚功相等,求等效节点载荷。(1)集中力的移置。设单元ijm内坐标为(x,y)的任意一点M受有集中载荷f=[fxfy]T,移置为等效节点载荷Pe=[XiYiXjYjXmYm]T。假想单元发生了虚位移,其中,M点虚位移为u*=NT(δ*)e,其中(δ*)e为单元节点虚位移。按照静力等效原则有fNδfuPδTT*T*T*))(()())((eee则fNPTe(2-1-17)(2)体力的移置。设单元承受有分布体力,单位体积的体力记为q=[qxqy]T,其等效节点荷载为yxtAeddTqNP(2-1-18)(3)面力的移置。设在单元的某一个边界上作用有分布的面力,单位面积上的面力为p=[pxpy]T,在此边界上取微面积tds,对整个边界面积分,得到lestdTPNP(2-1-19)例2-1求单元在以下受力情况下的等效节点荷载:y方向的重力为G、图2-2所示ij边受x方向均布力p、图2-3所示jm边受x方向线性分布力。xyxyOpmlljmijipO图2-2均布力图2-3线性分布力求解说明利用上述公式求等效节点载荷,当原载荷是分布体力或面力时,进行积分运算是比较繁琐的。但在线性位移模式下,可以按照静力学中力的分解原理直接求出等效节点载荷,上述三种情况等效节点荷载分别为51010103GeP00021021ptleP031032002ptleP6.整体分析结构的整体分析就是将离散后的所有单元通过节点连接成原结构物进行分析,分析过程是将所有单元的单元刚度方程组集成总体刚度方程,引进边界条件后求解整体节点位移向量。总体刚度方程实际上就是所有节点的平衡方程,由单元刚度方程组集总体刚度方程应满足以下两个原则:(1)各单元在公共节点上协调地彼此连接,即在公共节点处具有相同的位移。由于基本未知量为整体节点位移向量,这一点已经得到满足。(2)结构的各节点离散出来后应满足平衡条件,也就是说,环绕某一节点的所有单元作用于该节点的节点力之和应与该节点的节点载荷平衡。每一节点统一使用整体节点编号(如图2-4所示),第4单元节点编号i、j、m统一依次改为8、7、5。确定各单元的大域变换矩阵,如第4单元为0000000000000000000000004IIIG其中,I为2×2阶单位矩阵。51234678123456789imj图2-4节点与单元编号求出各单元刚度矩阵,利用大域变换法求出结构整体刚度矩阵[K],引入边界条件,得到结构的节点平衡方程PKδ(2-1-19)进而,求解节点位移、单元应力和应变。例2-2如图2-5所示,一悬臂梁,自由端受合力为P的均布力作用,梁厚t=1,μ=1/3,求节点位移。12221x341Py2/P2/P图2-5结构与离散6求解说明结构为平面应力问题,划分为2个三角形单元①、②,有四个节点1、2、3、4,坐标分别为(0,0)、(2,0)、(2,1)、(0,1)。单元①、②节点顺序分别取3、1、2和1、3、4,刚度矩阵完全一样。对单元①:bi=0,bj=-1,bm=1,ci=2,cj=0,cm=-2。对单元②:bi=0,bj=1,bm=-1,ci=-2,cj=0,cm=2。单元的刚度矩阵为134121224723241210022303201220212024200432321EKK利用大域变换法求出整体刚度矩阵为134120012247230024121301220423072440001221341200244723122041213024402307323EK节点荷载向量为P=[000-P/20-P/2000]T位移向量为δ=[u1v1u2v2u3v3u4v4]T由结构平衡方程求得节点位移为T0042.850.199.888.100EPδ7.平面问题高次单元如前所述,三节点三角形单元因其位移模式是线性函数,应变与应力在单元内都是常量,而弹性体实际的应力场是随坐标而变化的。因此,这种单元在各单元间边界上应力有突变,存在一定误差。为了更好地逼近实际的应变与应力状态,提高单元本身的计算精度,可以增加单元节点而采用更高阶次的位移模式,称为平面问题高次单元。如六节点三角形单元、矩形单元等。这里只介绍六节点三角形单元与矩形单元的位移模式,其它单元的位移模式和具体求解步骤与三节点三角形单元类似,且应用较少,不在赘述。图2-6六节点三角形单元图2-7四节点矩形单元六节点三角形单元如图2-6所示,在三角形单元各边中点处增加一个节点,则每个单元有六个节点,共有12个自由度。位移模式的项数应与自由度数相当,阶次应选得对称以保证几何各向同性。其7位移模式应取完全二次多项式2121121098726524321yxyxyxvyxyxyxu(2-1-20)将节点的坐标和位移代入即可求出广义坐标。四节点矩形单元如图2-7所示,共有8个自由度,取位移模式为xyyxvxyyxu87654321(2-1-21)将节点的坐标和位移代入即可求出广义坐标。
本文标题:有限元理论与方法-第7讲
链接地址:https://www.777doc.com/doc-2375719 .html