等参数单元作业002.docVIP

有限元等参数单元 学生姓名 学生学号 专业班级 任课老师 2015 年 4 月 26 日 一.等参数单元概述 利用等参数单元解决体复杂的曲线边界或需要布置疏密不匀的网格，克服了网格划分受边界的形状影响，单元大小可以不相等，是一种精度高而且应用广泛的单元。本次作业结合无限小单元法和有限单元法的优点, 提出了平面等效桁架模型，分析等效后桁架杆件的初始刚度、截面面积以及等效弹性模量等特征值。建立了等效桁架模型的单刚矩阵、应变矩阵和单元的轴力阵。在钢架结构的应用中, 给出了平面等效桁架模型的分离式模型,建立了等参单元平面等效桁架模型的单元刚度矩阵。 从数学角度上讲, 有限单元法是将一个结构或连续体的求解域离散为若干个子域(单元), 并通过边界节点连接求解的一种数值计算方法。有限元解的精确性直接与位移插值函数的选取有关。模型复杂, 计算烦琐。因此, 本文提出了一种新的单元形式即平面等效桁架模型, 无限小微元体宏观化后全部等效为平面二力杆件。此模型的最大特点是不需要设形函数, 无需考虑介质之间的相互作用, 可以全程追踪构件或结构的整个受力过程。特别是对钢筋混凝土结构, 用此单元模型中二力杆件的破坏可以明确判断钢筋和混凝土的破坏情况。 二、平面等效桁架单元模型 1.平面等效桁架模型的假定和总体描述 对宏观化之后的平面应力单元做以下基本假定(图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 平面应力单元与平面等效桁架模型等效图 ( 4)假定在外力作用下, 在计算x 方向的变形时, 两单元y 方向上的变形相同。这样在计算两单元刚度等效时, 省去y 方向的变形这一项。反之亦然。由上可知, 将一个正方形微元体简化为4个质点和6根杆组成的平面桁架体系, 由于构成桁架单元的各个杆件都是二力杆, 只承受轴力, 使计算过程简单化, 只要这个桁架体系的力学性能(线性和非线性)与原来的平面单元相同, 则整体就可以看作由轴向拉压杆组成的集合体。这就必须要求在相同的结点荷载作用下, 两种单元模型在相对应的方向上变形相同。 2.平面等效桁架模型的计算方法 由材料力学的基本原理可知单元在外力作用下会产生变形, 两单元的变形如下图2( a)、( b)所示。 ( a)轴向等效变形 ( b)剪切等效变形 图2 两种单元的等效荷载- 变形关系 设图2中微元体轴向变形和剪切变形为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 向杆的等效面积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 + ( 1 + v) 2 + 1 ( 3) 图3 单元节点力等效关系图 在有限单元法单元分析中, 需将单元边界应力转化为单元等效节点力, 等效关系如图3所示。单元等效节点力的计算方法如式( 4) : 式(