《现代设计方法4-3 三角形三节点平面单元》由会员分享，可在线阅读，更多相关《现代设计方法4-3 三角形三节点平面单元（101页珍藏版）》请在金锄头文库上搜索。
1、4.3 4.3 平面平面问题问题问题问题 的有限的有限单单单单元法元法三角形三节点平面单元结构离散化单元分析整体分析有限元分析的基本步骤：1 结构离散化例：图示一悬臂梁，梁的厚度为t，设泊松 比m =1/3，弹性模量为E，试用三节点三角形单元 进行离散。l2 m1 mAP/2P/23P/2P/2412yiijjmmx2 单元分析l单元分析的主要内容：由节点位移求内部任一点的位移，由节点位移求单元应变、应力和节点力。l单元分析的步骤：节点 位移单元内部 各点位移单元 应变单元 应力节点 力(1)(2)(3)(4)单元分析(1) 单元位移模式l有限元法中，通常采用位移法进行计算，即取节点的位移分量
2、为基本未知量，单元中的位移、应变、应力等物理量，都和基本未知量相关联。l 节点i的位移分量可写成di = ui viTl单元节点位移向量d e可写成 d e = di dj dmT = ui vi uj vi um vmT三角形三节点单元xymjiviuiujvjvmumuv2.1 由节点位移求单元内部任一点位移l六个位移分量需六个待定参数a1，a2，a6 l设单元内任一点位移分量是坐标(x,y)的线性函数：l (4.22)l写成矩阵的形式为：l (4.23)l设节点i，j，m坐标分别是xi，yi；xj，yj；xm，ym。把三个节点的坐标及其水平位移代入式(4.22)中得：l (a)l解得：l
3、 (b)l其中， (2) 由单元节点位移d e求位移参数a l对v同理可列出a4、a5、a6的方程 ：vi = a4+a5 xi+a6 yi i, j, ml解得： a4、a5、a6l为了书写的简便与规格化，引入记号lai, bi, ci分别为行列式A中各行各个元素的代数余子式i j m解出a1 a6结果：la =Ld e为三角形单元面积。将a写成矩阵形式，有a =Ld el由单元节点位移d e求单元内部任一点位移d (x,y)ld (x,y)=m(x,y)a =m(x,y)Ld ei j mlNi, Nj, Nm是坐标的连续函数，反映单元内位移的分布状态，称为位移的形状函数，简称形函数(sh
4、ape function)。矩阵N称为形函数矩阵(shape function matrix) 。形函数物理意义：Ni(x,y) ，节点i单位位移，其它节点位移分量为0，单元内部产生位移分布形状l在单元任一点上，三个形函数之和等于1。l形函数Ni在i点的函数值为1，在j点及m点的函数值为零。l三角形单元i, j, m在ij边上的形函数与第三个顶点的坐标无关。形函数的性质：例：求图示单元和单元的形函数矩阵3P/2P/2412yiijjmmx（a）(b)(0, a)(0, 0)(b,0)yijmxab(b, a)(0, a)yijmx(b,0)分别如上图所示：i j m单元如图所示。设a=1m,
5、b=2m.l(或直接由图形可知其面积)求系数ai, aj, am, bi, bj, bm, ci, cj, cm(0, a)(0, 0)(b,0)yijmxabi j m求形函数矩阵l代入相关常数：将a=1,b=2代入得：(b, a)(0, a)yijmx(b,0)求常数l单元如图所示=ab/2求形函数矩阵将a=1, b=2m代入上式得：作业：l已知三角形三节点单元坐标如图示，设单元中一点A的坐标(0.5,0.2)，已知三角形三节点单元i节点位移(2.0,1.0)， j节点位移(2.1,1.1)， m节点位移(2.15,1.05)，1)写出单元的位移函数；2)求A点的位移分量。i(1,0)j(
6、1,1)m(0,0)xyA2.2 由节点位移求单元的应变几何 方程i j mi j m简记为 e =Bd eB可写成分块的形式： B=Bi Bj BmB称为应变 矩阵，它的元素都只与单元的几何性质有关的常量。这种单元称为平面问题的常应变三角形单元。单元应变与单元节点位移关系(i, j, m)l物理方程s =De 而 e =Bde ls = DBd e= S d e (求应力的表达式) l记 S=DB S应力矩阵: S=Si Sj Sm2.3 由单元节点位移求单元的应力1) 平面应力问题 ：l代入D及Bl得: S =Si Sj Sm.l对于平面应力：l (i, j, m )l将上式中以 E/(1
7、-m2) 代 E，以 m/(1-m) 代 m 则子矩阵：l (i, j, m)2)对于平面应变问题 ：l考虑单元平衡，节点力是作用在单元上的外力，与单元应力平衡。有限元法中以虚功方程代替平衡方程。l节点力列阵及单元内应力列阵：单元节点力是指单元和节点相连接的内力；考虑节点平衡，节点力为外力，与节点外载荷平衡；sysxtxy2.4 由节点位移求单元节点力(求单元刚度矩阵)节点虚位移列阵及虚应变：l令实际受力状态在虚位移状态上做虚功，虚功方程：l由 e =Bd e 知 e *=Bd *el则 e *T=(d *e)TBTl由于d *e中的元素为常量，提至前积分号前，故：l(对于三角形三节点单元，B
8、和s 为常量，单元厚度t也是常量；为三角形单元面积，用表示）单元刚度矩阵物理方程s = De 简记为 Fe = ked eke = BTDBt 或= BTSt几何方程e = Bd e对于平面应力问题：(r =i, j, m; s =i, j, m).将上式中的E 换成m换成,对于平面应变问题 ：求例4.2（p84）单元的单元刚度矩阵(0, a)(0, 0)(b,0 )yijmxab(i, j, m)l解：（1）求矩阵B（2）求矩阵S（3）求矩阵kelke=BTDBt =BTStl代入a=1, b=2m 得：l可算出，当a=b时单元刚度矩阵与尺寸a,b 无关在单元节点力列阵Fe、单元应力列阵s
9、e、单元应变列阵e e和单元节点位移列阵d e的四个列阵之间，存在五个转换关系，可得五个转换矩阵。平衡条 件物理关 系几何关 系单元刚度矩阵的特性：(1)单元刚度矩阵的物理意义：单元刚度矩阵ke表示了单元抵抗变形的能力，即表示了节点位移d e与节点力Fe之间的关系。kij表示节点j发生单位位移时，其它节点位移分量均为零时，在节点i上产生的节点力。(2)分块性质:单元刚度矩阵可以分块运算。(4.51)按节点进行分块，则单 元刚度矩阵的分块形式可写为：(4.52)（3）对称性 单元刚度矩阵是一个对称矩阵式(4.51)中： knl=kln (n=16; l=16)分块形式中： krs=ksrT (r
10、 =i, j, m; s =i, j, m)（4）奇异性 单元刚度矩阵是一个奇异矩阵|ke|=0，表明其逆矩阵不存在。即，如果给定了单元节点位移可以得出唯一的节点力: Fe=ked e。反之，如果给出节点力却无法求出确定的节点位移。因为这时的单元未考虑所受约束时，可能存在不引起单元应力和节点力的刚性位移，这部分刚体位移由节点力是无法唯一确定的。三节点三角形单元每行元素之和为零三节点三角形单元每行元素之和为零例：证明图示单元刚度矩阵: kI=kIIIl证明：由于单元刚度矩阵l ke=BTDBtl可知：当两个三角形单元几何尺寸相同时，t值和单元面积值均相同；当两个单元的材料性质相同时，弹性矩阵D也
11、时相同的。故ke是否相同，取决于矩阵B是否相同。l不难验证 ,I、III单元的上述br, cr (i, j, m)值均相 等。l结论：l两个单元刚度矩阵ke相等的条件为：只要两单元的形状、大小，方向和单元弹性常数均相同，并且编号的方式也相同(如按逆时针方向编号为 i, j, m，直角顶点编号为m)，则两个单元的刚度矩阵时相等的。2.5 单元载荷的移置. （离散时每个单元受载作用于节点上）(a)原则：将单元载荷向节点处移置，按照虚功等效的原则进行。对于变形体(包括弹性体)，虚功等效是指原载荷与节点载荷载在任何虚位移上做的虚功相等。当位移模式确定后，载荷移置(或分解)其结果是唯一的。虚功等效包含了
12、刚体体系的静力等效，当虚位移为刚体位移时，虚功等效即为静力等效，静力等效是虚功等效的特例。离散(b)载荷移置公式(1)集中力设单元i, j, m中任一点M(x, y)处受有集中力P=Px PyT，移置到该单元各节点处载荷列阵为Re=Xi Yi Xj Yj Xm YmT (1) 集中力P=Px PyT(1)集中力l假设该单 元发生一微小虚位移，M点相应的虚位移为 f *，该单元各节点处相应虚位移为d * ，由静力等效原理，载荷与节点等效载荷在虚位移上所作虚功相等：l (d *e)TRe= f *TPl 将 f *=N d *e 代入上式：l有 (d *e)TRe = f *TP=(d *e)TNTPl则 Re=NTPlRe=Xi Yi Xj Yj Xm YmT=NiPx NiPy NjPx NjPy NmPx NmPyTl设单元i, j, m的一边受有分布的面力P=Px Py将微元面积tds上的面力合力Ptds当作集中载荷，可得面力的移置公式:(2)面力(2) 面力(3)体力l设单元i, j, m受有分布体力G=Px PyT将微分体积tdxdy上的体力合力Gtdxdy当作集中载荷dG同理可得：(3) 分布体力G=X Y T(4) 三节点三角形单元上同时有体力、面力和集中力 等(a) 集中力P=Px PyT(b) 分布体力G=X Y
