第九章__等参单元.pptVIP

写成矩阵形式 第九章 等参数单元 当 r=0, 以 代替 求解是否收敛还需看位移函数是否满足三个条件: 位移函数中包含有刚体位移 位移函数中包含有常应变 位移函数内部连续和边界协调 由于刚体位移和常应变状态可写成下列形式: 第九章 等参数单元 §9-4 等参数单元的收敛性 但上二式可化成下式: 根据形函数的性质: ΣNi=1 第九章 等参数单元 所以构成的位移函数包含刚体位移和常应变 多项式在单元内部无奇异点，所以连续。相邻单元的边界上可由线上节点座标唯一确定，而母单元相邻边界都是直线，且形函数相同，所以交界线上的位移由节点位移唯一确定，相邻单元连续。 则上二式都化成 六面体 20节点等参数单元 母单元是边长为2的立方体，通过映射可得一曲面六面体 第九章 等参数单元 §9-5 空间等参数单元 三维等参数单元的位移函数 第九章 等参数单元 第九章 等参数单元 要 点 采用等参数单元的意义 等参数名称的来由 如何根据形函数的特性建立位移函数 数值积分采用高斯积分法 §9-1 如何提高有限元法计算精度 — 采用等参数单元 第九章 等参数单元 9-1-1 采用等参数单元的意义 1. 找到结构物实际变形的模式，然后取与实际变形模式一致的单元离散结构物。 2. 从三角形单元到矩形单元，随着单元位移函数中变量幂次的增加则单元内应力位移的计算精度提高，所以就需要找一些高阶的多项式作为单元的位移模式建立新的单元。常用的多项式是根据杨辉三角形的次序建立的。 3. 由于项数取得多就必须增加单元的自由度，即增加单元节点，于是就出现了等参数单元。 9-1-2 等参数单元 第九章 等参数单元 每一种等参数单元有二种形式：一种形式在局部座标 ξ η 或 ξ η ζ 中，是以边长为2的正方形单元或边长为2的立方体单元用作计算称为母单元；另一种形式在整体座标 x y或x y z中，通过母单元映射到整体座标上的单元用于离散结构物称为子单元。由于母单元上的位移函数与座标变换的函数具有相同的参数。这就是等参数单元名称的来由。 母单元 的位移函数: 座标变换的映射关系: 第九章 等参数单元 §9-2 平面等参数单元 9-2-1 四边形四节点等参数单元 四边形四节点单元位移模式: 其中 母单元 若以母单元上12边为例，通过映射可得在平面内任一直线，12边的方程为 η= -1,代入 第九章 等参数单元 通过局部座标与整体座标的映射关系把母单元变换到整体座标上成为一个任意四边形用于离散结构物，它能适合于任意曲边的形状。 同理可得 第九章 等参数单元 把以参数 ξ 代表的 x方程和 y方程消去 ξ ，则得 x , y 所组成的直线方程 y=kx+b 所以母单元上的四个边都可以通过映射在x, y座标面上得出一个任意四边形，用该四边形离散结构物。 四边形四节点 单元的应变 第九章 等参数单元 i = 1,2,3,4 其中 其中 第九章 等参数单元 由于形函数 N 是 x, y的函数，现对ξ, η求导，这是复合函数求导 把上二式写成矩阵形式 其中 第九章 等参数单元 称为雅可比(Jacobian)矩阵 为了把 [Bi]矩阵中的Ni,x和 Ni,y化成 Ni,ξ和 Ni,η ，则代入下式 其中雅可比矩阵的逆阵由下式给出 也可把形函数Ni对的ξ, η 偏导数写成通式: i= 1,2,3,4 第九章 等参数单元 四边形四节点 单元的应力 对于平面应力情况 第九章 等参数单元 四边形四节点单元的刚度矩阵[K]e: [K]e可划分成四行四列的子矩阵 : i, j = 1,2,3,4 对于平面应力情况: 第九章 等参数单元 等效节点力计算: 1. 集中力的等效节点力:由于计算复杂，把有集中力处设置为节点。 2. 体枳力的等效节点力: 3. 表面力的等效节点力: 第九章 等参数单元 展开后 i= 1,2,3,4 单元上某一点的节点力 由于表面力的运算是曲线积分，所以常把 qx ,qy化成切线, 第九章 等参数单元 展开后 单元上某一点的节点力 i=1,2,3,4 第九章 等参数单元 把用 σ,τ 表示 qx , qy的代入 i=1 ,2 3. 4 若以母单元上的34边为例，则η =1 i= 3,4 母单元是以2为边长的正方形具有八个节点，由于