等全参数单元作业002.docxVIP

实用文案 实用文案 标准 标准 有限元等参数单元 学 生姓名 学生学号 专业班级 任课老师 2015 年 4月26日 一 ?等参数单元概述 利用等参数单元解决体复杂的曲线边界或需要布置疏密不匀的 网格，克服了网格划分受边界的形状影响，单元大小可以不相等，是 一种精度高而且应用广泛的单元。本次作业结合无限小单元法和有限 单元法的优点，提出了平面等效桁架模型，分析等效后桁架杆件的初 始刚度、截面面积以及等效弹性模量等特征值。 建立了等效桁架模型 的单刚矩阵、应变矩阵和单元的轴力阵。在钢架结构的应用中 ，给出 了平面等效桁架模型的分离式模型，建立了等参单元平面等效桁架模 型的单元刚度矩阵。 从数学角度上讲，有限单元法是将一个结构或连续体的求解域 离散为若干个子域(单元),并通过边界节点连接求解的一种数值计算 方法。有限元解的精确性直接与位移插值函数的选取有关。 模型复杂, 计算烦琐。因此，本文提出了一种新的单元形式即平面等效桁架模型 ， 无限小微元体宏观化后全部等效为平面二力杆件。 此模型的最大特点 是不需要设形函数，无需考虑介质之间的相互作用，可以全程追踪构 件或结构的整个受力过程。特别是对钢筋混凝土结构，用此单元模型 中二力杆件的破坏可以明确判断钢筋和混凝土的破坏情况。 二、平面等效桁架单元模型 平面等效桁架模型的假定和总体描述 对宏观化之后的平面应力单元做以下基本假定(图1)。 方形微元体x、y方向的长度dx = dy= h, 与之对应的 平面等效桁架模型的外围杆长也为 h。z方向厚度为t。 微元体的材料弹性模量为E,剪切模量为G,在弹性 受力阶段有:G = E /[ 2( 1+ v)] 其中v为材料的泊松比 平面等效桁架模型中水平和竖向杆的初期刚度相同 ，为k eq1; 内部两根斜杆的初期刚度为ke q2。 ⑴平瀟应h单元 〔B)平而等效桁朝单.元 图1平面应力单元与平面等效桁架模型等效图 假定在外力作用下，在计算x方向的变形时，两单元y方向上的 变形相同。这样在计算两单元刚度等效时 ，省去y方向的变形这一 项。反之亦然。由上可知，将一个正方形微元体简化为4个质点和6 根杆组成的平面桁架体系，由于构成桁架单元的各个杆件都是二力 杆，只承受轴力，使计算过程简单化，只要这个桁架体系的力学性能 (线性和非线性)与原来的平面单元相同，则整体就可以看作由轴向拉 压杆组成的集合体。这就必须要求在相同的结点荷载作用下 ，两种单 元模型在相对应的方向上变形相同。 平面等效桁架模型的计算方法 由材料力学的基本原理可知单元在外力作用下会产生变形 ，两单元 的变形如下图2( a)、( b)所示。 (a)轴向等效变形 (a)轴向等效变形 (b)剪切等效变形 图2两种单元的等效荷载-变形关系 设图2中微元体轴向变形和剪切变形为 D1、D2,与之对应的 平面等效桁架模型为Del、Dc2,等效计算的条件为等效结点荷 载作用下的结点位移相等，即： D1 = Del; D2 = De2 由力学原理得桁架体系中杆件初始刚度为：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;即： Aeql = ht /2 Aeq2 = 2K e q2Aeq1 /K e q1 = ht 2 v + ( 1 + v ) 2 + 1 ( 2) 其中keql, 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): Fx ="吒触Z2— <|T;.A//2 F认二 - ^*ht/2 i I t r F =曲?月 f/2+ ? h t/2 F = h t /2 + di ? /? f X2 <J 导 77 T 玛 (4) F込二曲；加Z!? rns：ht/1 Fkv 二曲：族左一