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

4.3 平面问题的有限单元法 三角形三节点平面单元 有限元分析的基本步骤： 1 结构离散化 例：图示一悬臂梁，梁的厚度为t，设泊松比m =1/3，弹性模量为E，试用三节点三角形单元进行离散。 2 单元分析 单元分析的主要内容：由节点位移求内部任一点的位移，由节点位移求单元应变、应力和节点力。 单元分析的步骤： 2.1 由节点位移求单元内部任一点位移 (1) 单元位移模式 有限元法中，通常采用位移法进行计算，即取节点的位移分量为基本未知量，单元中的位移、应变、应力等物理量，都和基本未知量相关联。 节点i的位移分量可写成 {di} = [ui vi]T 单元节点位移向量{d }e可写成 {d }e = [di dj dm]T = [ui vi uj vi um vm]T 六个位移分量需六个待定参数 a1，a2，……，a6 设单元内任一点位移分量是坐标(x,y)的线性函数： (4.22) 写成矩阵的形式为： (4.23) (2) 由单元节点位移{d }e求位移参数{a } 设节点i，j，m坐标分别是xi，yi；xj，yj；xm，ym。把三个节点的坐标及其水平位移代入式(4.22)中得： (a) 解得： (b) 其中， 对v同理可列出a4、a5、a6的方程 ： vi = a4+a5 xi+a6 yi i, j, m 解得： a4、a5、a6 为了书写的简便与规格化，引入记号 ai, bi, ci分别为行列式A中各行各个元素的代数余子式 解出a1~ a6结果： {a }=[L]{d }e 将a写成矩阵形式，有{a }=[L]{d }e 由单元节点位移{d }e求单元内部任一点位移{d (x,y)} {d (x,y)}=[m(x,y)]{a }=[m(x,y)][L]{d }e Ni, Nj, Nm是坐标的连续函数，反映单元内位移的分布状态，称为位移的形状函数，简称形函数(shape function)。矩阵[N]称为形函数矩阵(shape function matrix) 。 形函数物理意义：Ni(x,y) ，节点i单位位移，其它节点位移分量为0，单元内部产生位移分布形状 形函数的性质： 在单元任一点上，三个形函数之和等于1。 形函数Ni在i点的函数值为1，在j点及m点的函数值为零。 三角形单元i, j, m在ij边上的形函数与第三个顶点的坐标无关。 例：求图示单元Ⅰ和单元Ⅱ的形函数矩阵 分别如上图所示： ①单元Ⅰ如图所示。设a=1m, b=2m. (或直接由图形可知其面积) ②求系数ai, aj, am, bi, bj, bm, ci, cj, cm ③求形函数矩阵 代入相关常数： ④求常数 单元Ⅱ如图所示Δ=ab/2 作业： 已知三角形三节点单元坐标如图示，设单元中一点A的坐标(0.5,0.2)，已知三角形三节点单元i节点位移(2.0,1.0)， j节点位移(2.1,1.1)， m节点位移(2.15,1.05)， 1)写出单元的位移函数；2)求A点的位移分量。 2.3 由单元节点位移求单元的应力 物理方程 {s }=[D]{e} 而 {e }=[B]{d}e {s } = [D][B]{d }e = [S] {d }e (求应力的表达式) 记 [S]=[D][B] [S]应力矩阵: [S]=[Si Sj Sm] 1) 平面应力问题： 代入[D]及[B] 得: [S ]=[Si Sj Sm]. 对于平面应力：