等参数单元作业002(共9页)

1、湖南大学等参数作业 有限元等参数单元 学生姓名 学生学号 专业班级任课老师 2015 年4 月 26 日一.等参数(cnsh)单元概述(i sh) 利用(lyng)等参数(cnsh)单元解决(jiju)体复杂的曲线边界或需要布置疏密不匀的网格，克服了网格划分受边界的形状影响，单元大小可以不相等，是一种精度高而且应用广泛的单元。本次作业结合无限小单元法和有限单元法的优点, 提出了平面等效桁架模型，分析等效后桁架杆件的初始刚度、截面面积以及等效弹性模量等特征值。建立了等效桁架模型的单刚矩阵、应变矩阵和单元的轴力阵。在钢架结构的应用中, 给出了平面等效桁架模型的分离式模型,建立了等参单元平面等效桁架

2、模型的单元刚度矩阵。 从数学角度上讲, 有限单元法是将一个结构或连续体的求解域离散为若干个子域(单元), 并通过边界节点连接求解的一种数值计算方法。有限元解的精确性直接与位移插值函数的选取有关。模型复杂, 计算烦琐。因此, 本文提出了一种新的单元形式即平面等效桁架模型, 无限小微元体宏观化后全部等效为平面二力杆件。此模型的最大特点是不需要设形函数, 无需考虑介质之间的相互作用, 可以全程追踪构件或结构的整个受力过程。特别是对钢筋混凝土结构, 用此单元模型中二力杆件的破坏可以明确判断钢筋和混凝土的破坏情况。二、平面等效桁架单元模型1.平面等效桁架模型的假定和总体描述对宏观化之后的平面应力单元做以

3、下基本假定(图1)。( 1)方形微元体x、y 方向的长度dx = dy= h, 与之对应的平面等效桁架模型的外围杆长也为h。z 方向厚度为t。( 2)微元体的材料弹性模量为E, 剪切模量为G, 在弹性受力阶段有: G = E / 2( 1+ v) 其中v为材料的泊松比( 3)平面等效桁架模型中水平和竖向杆的初期刚度相同, 为k eq1; 内部两根斜杆的初期刚度为ke q2。图1 平面(pngmin)应力单元与平面等效桁架模型等效图( 4)假定在外力作用下, 在计算x 方向的变形时, 两单元y 方向上的变形相同。这样在计算两单元刚度等效时, 省去y 方向的变形这一项。反之亦然。由上可知, 将一个

4、正方形微元体简化为4个质点和6根杆组成(z chn)的平面桁架体系, 由于构成桁架单元的各个杆件都是二力杆, 只承受轴力, 使计算过程简单化, 只要这个桁架体系的力学性能(线性和非线性)与原来的平面单元相同, 则整体就可以看作由轴向拉压杆组成的集合体。这就必须要求在相同的结点荷载作用下, 两种单元模型在相对应的方向上变形相同。2.平面等效桁架(hngji)模型的计算方法由材料力学的基本原理可知单元在外力作用下会产生变形, 两单元的变形如下图2( a)、( b)所示。( a)轴向等效(dn xio)变形( b)剪切等效(dn xio)变形图2 两种单元的等效(dn xio)荷载- 变形关系设图2

5、中微元体轴向变形和剪切变形为D1、D2, 与之对应的平面等效桁架模型为Dc1、Dc2, 等效计算的条件为等效结点荷载作用下的结点位移相等, 即: D1 = Dc1; D2 = Dc2由力学原理得桁架体系中杆件初始刚度为: K eq1 = Eht 1+ v v+ ( 1+ v ) 2 + 1 K eq2 = 2# h tE / 2( 1+ v) ( 1) 假设等效后桁架模型中桁架杆受拉或受压强度与原平面应力单元对应相同, 由于考虑微元体的弹塑性特性, 在桁架体系中杆件达到抗压强度之前, x 向桁架杆已被拉断, 斜杆也退出工作, 而外力主要由y向杆件承受, 因此可得桁架体系中x 向和y 向杆的等效

6、面积A1 为原微元体截面面积的1 /2, 再假设桁架杆的刚度与其面积成正比, 与长度成反比, 即可得到斜杆的等效面积Ae q2; 即: Aeq1 = ht /2 Aeq2 = 2K e q2Aeq1 /K e q1 = ht 2 v + ( 1 + v ) 2 + 1 ( 2) 其中keq1, ke q2 为式( 1) 中的桁架杆件的等效初始刚度, h 为模型的边长, t为厚度, v为泊松比。 由轴向拉压杆单位长度刚度公式K eq = E eqAeq, 可得桁架杆等效弹性模量: E eq1 = 2E 1 + v v + ( 1 + v) 2 + 1 E eq2 = 2E 1 + v v + (

7、 1 + v) 2 + 1 ( 3) 图3 单元节点(ji din)力等效关系图在有限单元法单元分析(fnx)中, 需将单元边界应力转化为单元等效节点力, 等效关系如图3所示。单元等效节点力的计算方法如式( 4) : 式( 4) 中F ix, F iy, Fj x, Fjy, Fkx, Fky, Fm lx, F ly 表示平面等效桁架单元(dnyun)节点力在x、y 方向上的分量。图4 平面等效桁架单元(dnyun)杆件和结点自由度编号经验证, 平面(pngmin)应力单元和平面等效桁架单元两者弹性形变势能、外力势能和总势能分别相等(xingdng)5 , 证实等效模型合理。3.平面等效桁架

8、模型的计算步骤基于有限元的解题步骤如下。( 1) 由相同外力作用下结点等效荷载相等的等效条件, 求出平面等效桁架模型的初始刚度K eq1, K e q2。( 2) 对平面等效桁架单元中的各杆件编号, 并对结点自由度和节点力编号, 计算各杆的刚度矩阵, 由刚度集成原理得单元的单刚K e, 图4给出了单元的编号和自由度编号(假设节点力编号和结点自由度编号相同)。K e 为桁架单元的单元刚度矩阵, 形式如式( 5)。 ( 3) 平面桁架(hngji)单元的结点力和结点位移关系可表示为: F e = K e D e ( 6) 整体(zhngt)分析得: F = K D ( 7) 由式( 6)、式( 7

9、) 结合边界条件可以(ky)求得整体结点位移。( 4) 根据第( 3) 步骤求得的结点位移, 由应变的几何意义可以求出各杆件应变和结点位移的关系为: Nhx, Nhy 分别表示桁架单元x、y 方向各杆件应变的值, Nik, Njl 分别表示内部斜向杆件的应变。式中 D 称为平面桁架单元的几何矩阵。由本构方程可知: R e = D N e = D B D e ( 9)进一步可求出等效杆件的轴力: N e = A eq D B D e ( 10)由桁架杆的轴力可以很容易的计算出桁架单元的等效结点荷载: 上式中称M 为平衡矩阵。综上可知, 对应于等效(dn xio)桁架某一结点的位移可以由上述步骤(

10、bzhu)计算出结构的响应。三、 平面等效桁架单元在钢筋混凝土结构(jigu)中的应用1.平面等效桁架模型单元刚度矩阵的修正由于钢筋混凝土结构是由钢筋和混凝土复合材料组成, 本文选择钢筋混凝土有限元分析中的分离式模型, 对不规则形状的任意四边形单元这里引入了等参单元, 并给出了等参变换后的带钢筋单元的平面等效桁架的单元刚度(图5)。图5 等参平面桁架模型钢筋结点位移和单元结点位移关系 无钢筋时上述四边形等效桁架模型等参单元的刚度矩阵可表示为: 其中| J | 为雅可比行列式。钢筋单元可按照平面杆件线性单元计算, 其刚度矩阵可按照式( 13) 计算: C, S 分别表示杆件与水平(shupng)轴夹角的余弦、正弦。通过图5的几何关系可以(ky)将钢筋结点位移和四边形单元结点位移通过座标转换阵联系起来。 R 为座标转换(zhunhun)矩阵。 ks = R l ks R ( 15) 故计算程序中钢筋混凝土的刚度矩阵为: K CS = K C + K S ( 16) 通过等参数单元建立等效平面桁架代替平面应力单元处理非规则的四边形结构来分析平面问题。可将复杂的模型简单化, 只要划分的单元