结构有限元分析.docVIP

第七章 结构有限元分析 引 言 求解具体结构工程中的问题是有限元素法的最终目的，而实际工程结构是复杂多样的，要很好的运用有限元素法还得解决好像坐标变换、对称边界条件运用以及复杂结构连接等问题。本章即为解决有限元方法应用于工程结构中实际问题的算法。 一、杆系或梁系的刚度坐标变换 1、向量的坐标变换公式 一维向量的平面分解 ii) 一维向量的三维空间分解 iii) 平面向量的坐标变换： iv)在三维空间中平面向量的坐标变换 u v w 2、杆元局部系下刚阵与整体系下刚阵的变换 i）局部系下的单元平衡方程： 由坐标变换（对节点力） = 由局部坐标系下的平衡方程 由位移（节点）的坐标变换 代入{P}的表达式： 故 杆系举例： 节点编号 单元编号 形成各单元的总体坐标系下刚阵 单元拼装 求解总体刚度方程 3、平面梁元局部系下刚阵到整体系的坐标变换 i). 梁元局部系下的单元刚度平衡方程 ii) 坐标变换 iii) 空间梁元有更复杂的变换关系 iv）其它单元的坐标变换 Homework: 列出平面弹性问题的刚度矩阵向三维空间的变换 二、节点局部系下的刚度矩阵变换（斜支撑条件的应用） 实际问题 ii）问题：一些节点在总体坐标系下，一些节点是在局部坐标系下，这类问题称为混合坐标架问题，即最终的刚度矩阵是一个混合标架下的形式。我们的任务是从总体系下的刚阵推导混合标架下的刚阵计算。 iii) 推导（以杆元或平面单元为例） ，总体坐标系 若单元节点1取为局部坐标系，即： 一般规律： i 三、不同类型元素的结合和各种坐标变换矩阵 引：工程结构特别是飞行器、船舶等不可避免地要同时采用不同类型的元素，从而引起不同类型元素之间的位移协调性问题。为保证这种协调性，必须引入特殊的变换矩阵以实现之，实际情况繁多，难于一一举例，仅通过几个例子说明处理这类问题的方法。 1、桁架与刚架混合结构使用杆元与梁元（不同节点不同自由度） 以杆元a- b为例，在组建总刚阵时，a- b杆的刚阵要适当变换以适应节点a的自由度，即： 改坐标变换可在组装时实现。 方法2：还可以利用子结构方法，将a、d两点的转角作为内部自由度，而刚架的出口刚阵不包括该自由度（即聚缩掉），因该自由度与邻近元素没有协调要求。 2、梁元素可能与平面应力板元素混合使用 梁元素的节点取在轴线a- b上，但是位移协调关系发生在角节点1-2-3- 4上，且前者含有转角自由度，，而后者仅含有平移自由度。 即： 需要建立转换关系：(a点到1－2点；b点到3－4点；) T可依据初等弯曲理论的基本假设来决定。 如图： 则： Note: 上述变换，仅作了节点位移协调，而不能保证元素边界上各点的位移协调性。实际当中，很难保证它们的相容性，不得不放松要求，仅保证节点上位移一致，相应地产生一定程度上的误差。这样做，仍可保证一定的工程精度。 作业：推导该梁元与平面应力元位移协调的刚阵形式。 偏心梁元： 用1、2两点的广义位移取代中性轴上a、b节点的广 义位移，这样，1、2两点的位移可以与其它单元的节 点协调。 由初等弯曲理论： 同理： 四、约束不足和附加约束 这种问题的发生一般在工程结构上，例如如下结构，杆元在y’局部坐标系下计算没有任何问题，此时的总刚平衡方程中，对应的 节点位移仅有y’方向三个节点的位移。但放在总体 系下，或在x’y’系下计算时，总刚平衡方程中，在 一个节点上有两个位移自由度，而结构本身是不能 提供4节点x’方向的位移约束。从计算总平衡方程看， 总刚阵此时是奇异的（因为缺少约束条件），因从平 衡观点来看，当节点4上作用有x’方向的外力时， 结构内不可能产生内力与之平衡；换句话说，结构 内产生的内力值为无穷大，反过来说，若要节点4沿x’方向发生位移，无需施加任何外力（几何可变（瞬时）结构）。 再如空间平面应力元素组成的薄壁结构，节点a与上面所描述的道理一致，同属于约束不定的问题。 这类问题在总刚矩阵中一般有两类表现： 如所述杆系问题，当坐标系取为x’y’时， 单元和在4节点x’方向对应的刚阵元 素（行和列）为零，从而，总刚矩阵奇异； 当坐标系取为xy时，节点4在xy方向 对应的刚阵元素线性相关，同样，总刚矩阵奇异。 处理方法： 在这类特殊节点处，采用节点坐标系，并将不定的那个自由度抛去不计，这样做，在刚阵中导致奇异的那些行（列）就不再出现了。 这种做法的缺点是，使各节点的自由度数不等，坐标轴不统一，因而，程序设计上使工作复杂化。 将包