有限元法基础-9杆系结构力学问题要点详解.ppt

9.2 等截面直杆-梁单元 注意: 是梁高h的3次函数 是梁高h的1次函数 当 时，并且 是非奇异的 有限元方程 零解 解决办法 有限元法基础 “剪切自锁”现象 （1）减缩积分 （2）假设应变 * 9.2 等截面直杆-梁单元 （1）减缩积分 2节点单元对 使用1点积分 使用1点积分 有限元法基础 矩阵秩为1 * 9.2 等截面直杆-梁单元 产生原因 2节点单元的剪应变 当 ,要求 有限元法基础 第一式表示转角与挠度的关系 第二式导致转角相等，不能表示弯曲变形 * 9.2 等截面直杆-梁单元 （1）假设应变 由于直接导出的剪切应变可能不正确，因此，重新假设剪应变 并令 有限元法基础 插值点取在应变计算比较准确的点上，一般在应力最优点上，如2节点单元去在形心点处，3节点单元取2点Gauss积分点上。 * 9.2 等截面直杆-梁单元 例：计算悬臂梁问题 使用一个单元计算 （1）经典梁单元 有限元法基础 * 9.2 等截面直杆-梁单元 （2）Timoshenko梁单元 纯弯曲 有限元法基础 减缩积分 精确积分 * 9.2 等截面直杆-梁单元 （2）Timoshenko梁单元 剪切力 有限元法基础 减缩积分 精确积分 罚因子 误差25％ * 9.3 平面杆系 特点 （1）平面杆系受力状态 a）横向弯曲 b）沿轴向的拉压 （2）在局部坐标下 平面杆单元＝梁单元＋受轴向力杆单元 （3）平面杆系单元通过坐标转换 单元坐标 － 总体坐标 坐标转换 有限元法基础 * 9.3 平面杆系 单元刚度矩阵 例：2节点单元刚度矩阵 有限元法基础 * 9.3 平面杆系 坐标变换 总体坐标系 单元坐标系 单元坐标与总体坐标的 夹角为 ，方向余弦为 有限元法基础 * 9.3 平面杆系 节点DOF的转换关系 有限元法基础 * 9.3 平面杆系 总体坐标下的刚度矩阵 有限元法基础 * 9.3 平面杆系 铰节点的处理 铰接的特点：在铰接处的所有节点位移相同，转 动不同，不承受弯矩 处理办法：在集成总体刚度矩阵前将相关转动 DOF凝聚，不与其他单元的刚度矩阵集成。 有限元法基础 * 9.3 平面杆系 有限元法基础 * 9.3 平面杆系 如图2＃的一个端点与刚架铰接 单元刚度矩阵为 为与4＃节点连接的铰接DOF 由第二组方程得 有限元法基础 * 9.3 平面杆系 凝聚后的方程组 有限元法基础 为了编程方便，凝聚后的刚度矩阵，仍以零元素补充为6X6的矩阵 * 9.4 空间杆系 空间杆系的受力特点 1）两个方向的弯曲 2）沿轴向的拉压 3）绕轴线的扭转 有限元法基础 力和平移 力矩和转动 单元坐标和总体坐标 坐标转换关系 1）3个角度定位 2）9个方向余弦 * 9.4 空间杆系 空间杆单元 每节点有6个DOF，分别对应3个平移和3个转动 单元刚度矩阵中的截面几何常数： 截面面积 A xz面内的惯性矩 Iy xy面内的惯性矩Iz 绕x轴的扭转惯性矩 Ix 翘曲常数 扭转常数 有限元法基础 * 9.4 空间杆系 空间杆单元 每节点有6个DOF，分别对应3个平移和3个转动 单元刚度矩阵中的截面几何常数： 截面面积 A xz面内的惯性矩 Iy xy面内的惯性矩Iz 绕x轴的扭转惯性矩 Ix 有限元法基础 * 9.4 空间杆系 ANSYS中 BEAM188单元 截面性质 有限元法基础 y z * 9.4 空间杆系 在空间中的位置 由单元端部节点坐标可确认 单元在总体坐