平面桁架的有限元法.ppt

* * 平面桁架的有限元法 一、建立有限元模型 平面桁架是现成的有限元模型， 只有节点力 二、单元分析 位移假设，几何方程，物理方程， 单元刚度矩阵 三、总刚度矩阵、总载荷列阵、引入支撑 每个节点建立平衡方程 集成技术（叠加法、对号入座） 四、求解、结果分析，后处理 一、建立有限元模型 平面桁架是现成的有限元模型， 只有节点力 桁架原型=有限元模型 二、单元分析 平面二力杆单元 注意：一个单元的节点数目（ 2）； 一个节点的自由度数目（2）； 有整体坐标系和局部坐标系之分 二、单元分析（先在局部坐标系中进行） 问题：当节点位移一定，相应的节点力为多大？ 节点力与节点位移之间的比例常数为刚度系数 目标：建立单元刚度矩阵。 方法：结构力学的位移法加虚功原理， 皆用节点位移表示。 1、位移假设： 设各点的位移，即位移场为 其中 为待定系数，由节点坐标和节点位移确定。 将位移写成矩阵形式： 在节点上满足： 写成矩阵形式： 解线性代数方程组，得 代入 得 节点位移与单元内位移的关系 其中 为形函数。 用Machematica推导形函数： 在编程中只关心节点位移 的系数矩阵 2、应变（从位移求应变） 只有轴向应变 3、应力（从应变求应力，本构方程） 4、内力（从应力求内力） 二力杆只有轴向力 用Machematica推导应力矩阵： (EA的两个字母连写作参数名) 用内力表示的物理方程 4、单元刚度矩阵（用虚功原理） 给节点一组虚位移，单元内即产生相应的虚应变，实际节点力（轴向力、横向力、弯矩）在虚位移上所做的功等于实际应力（横向剪应力可忽略不计）在虚应变上产生的应变能。 用Machematica推导单刚矩阵[Ke]： 用Machematica推导，所得结果以矩阵形式显示 三、总刚度矩阵、总载荷列阵、引入支撑 1、坐标变换 前面的推导是在局部坐标系中进行的，而总刚度矩阵、总载荷列阵、引入支撑等要在整体坐标中进行，因此，存在坐标变换问题。 三、总刚度矩阵、总载荷列阵、引入支撑 1、坐标变换 三、总刚度矩阵、总载荷列阵、引入支撑 2、在总结构和整体坐标中讨论问题，所需要的信息： （1）节点号； （2）单元长度； （2）单元倾角(单元起点到终点的有向线段与整体横轴之夹角）； （4）节点力所在的节点号和作用方向； （5）支撑所在的节点号和作用方向。 3、总节点位移向量的排列规律 按节点号和节点自由度的顺序（从小到大）排列。 三、总刚度矩阵、总载荷列阵、引入支撑 建立节点号数组（用Mathematica)： jm=Table[{0, 0}, {i, ne}]; “jm数组名，ne 总单元数;” jm={…，{i节点号， j节点号}，…}; 建立单元长度数组 gc=Table[{0}, {i, ne}]; “gc数组名，ne 总单元数;” gc={L1, L2, … Lne}; 建立单元方向角数组 gj=Table[{0}, {i, ne}]; “gj数组名，ne 总单元数;” gj={af1, af2, … afne}; 3、总节点位移向量的排列规律 jm={{1, 2}，{2, 3}，…，{i节点号， j节点号}，…}; 单元节点位移排列规律： rr =2(k-1)+kk k =1, 2 (单元节点); kk =1, 2, (节点自由度) k=1 k=2 kk=1 kk=2 总节点位移向量{d}的排列规律(从小到大排序）： ss=2(jm[[e, k]]-1)+kk; e=单元号 对每个单元进行节点循环和自由度循环，从单元地址到结构地址搬家 是单元到结构的映射(map) 3、总节点位移向量的排列规律 单元节点位移号：rr =2(k-1)+kk；k =1, 2 (单元节点)； kk =1, 2(节点自由度) 总节点位移号：ss=2(jm[[e, k]]-1)+kk； 例如，已解得结构的总节点位移向量{d}求各单元的应力。 用Mathematica： deg=Table[0, {i, 4}]; “整体坐标系中的单元节点位移向量” ； del=Table[0, {i, 4}]; “局部坐标系中的单元节点位移向量” ；del= se=Table[{0, 0}, {i, ne}]; “单元内力数组