现代设计方法4-3三角形三节点平面单元.pptxVIP

4.3 平面问题的有限单元法三角形三节点平面单元第1页，共101页。有限元分析的基本步骤：结构离散化单元分析整体分析第2页，共101页。P/2A1 ml2 mP/2yP/243imjⅡⅠjxmi12P/21 结构离散化 例：图示一悬臂梁，梁的厚度为t，设泊松比m =1/3，弹性模量为E，试用三节点三角形单元进行离散。第3页，共101页。2 单元分析单元分析的主要内容：由节点位移求内部任一点的位移，由节点位移求单元应变、应力和节点力。单元分析的步骤：节点位移单元内部各点位移单元应变单元应力节点力(4)(1)(2)(3)单元分析第4页，共101页。yvjujjvivuvmuiiummx2.1 由节点位移求单元内部任一点位移(1) 单元位移模式有限元法中，通常采用位移法进行计算，即取节点的位移分量为基本未知量，单元中的位移、应变、应力等物理量，都和基本未知量相关联。 节点i的位移分量可写成{di} = [ui vi]T单元节点位移向量{d }e可写成 {d }e = [di dj dm]T = [ui vi uj vi um vm]T三角形三节点单元第5页，共101页。六个位移分量需六个待定参数a1，a2，……，a6 设单元内任一点位移分量是坐标(x,y)的线性函数： (4.22)写成矩阵的形式为： (4.23)第6页，共101页。(2) 由单元节点位移{d }e求位移参数{a }设节点i，j，m坐标分别是xi，yi；xj，yj；xm，ym。把三个节点的坐标及其水平位移代入式(4.22)中得： (a)解得： (b)其中， 第7页，共101页。i j m对v同理可列出a4、a5、a6的方程 ： vi = a4+a5 xi+a6 yi i, j, m解得：a4、a5、a6为了书写的简便与规格化，引入记号ai, bi, ci分别为行列式A中各行各个元素的代数余子式第8页，共101页。解出a1~ a6结果：{a }=[L]{d }e?为三角形单元面积。第9页，共101页。将a写成矩阵形式，有{a }=[L]{d }e由单元节点位移{d }e求单元内部任一点位移{d (x,y)}{d (x,y)}=[m(x,y)]{a }=[m(x,y)][L]{d }e第10页，共101页。i j m第11页，共101页。Ni, Nj, Nm是坐标的连续函数，反映单元内位移的分布状态，称为位移的形状函数，简称形函数(shape function)。矩阵[N]称为形函数矩阵(shape function matrix) 。形函数物理意义：Ni(x,y) ，节点i单位位移，其它节点位移分量为0，单元内部产生位移分布形状第12页，共101页。形函数的性质：在单元任一点上，三个形函数之和等于1。形函数Ni在i点的函数值为1，在j点及m点的函数值为零。三角形单元i, j, m在ij边上的形函数与第三个顶点的坐标无关。第13页，共101页。yP/243imjⅡⅠjmxi12P/2例：求图示单元Ⅰ和单元Ⅱ的形函数矩阵第14页，共101页。yy(0, a)(b, a)(0, a)jimⅡaⅠjxixbi j mm(b,0)(0, 0)(b,0)(b) （a）分别如上图所示：第15页，共101页。①单元Ⅰ如图所示。设a=1m, b=2m.(或直接由图形可知其面积) 第16页，共101页。y(0, a)jaⅠixi j mbm(0, 0)(b,0)②求系数ai, aj, am, bi, bj, bm, ci, cj, cm第17页，共101页。 ③求形函数矩阵代入相关常数：第18页，共101页。将a=1,b=2代入得：第19页，共101页。y(0, a)(b, a)imⅡjx(b,0)④求常数单元Ⅱ如图所示Δ=ab/2第20页，共101页。⑤求形函数矩阵第21页，共101页。将a=1, b=2m代入上式得：第22页，共101页。yj(1,1)xm(0,0)i(1,0)作业：已知三角形三节点单元坐标如图示，设单元中一点A的坐标(0.5,0.2)，已知三角形三节点单元i节点位移(2.0,1.0)， j节点位移(2.1,1.1)， m节点位移(2.15,1.05)，1)写出单元的位移函数；2)求A点的位移分量。·A第23页，共101页。2.2 由节点位移求单元的应变几何方程 第24页，共101页。i j mi j m第25页，共101页。(i, j, m)简记为{e }=[B]{d }e单元应变与单元节点位移关系 [B]称为应变矩阵，它的元素都只与单元的几何性质有关的常量。这种单元称为平面问题的常应变三角形单元。[B]可写成分块的形式： [B]=[Bi Bj Bm]第26页，共101页。2