有限元法基础9杆系结构力学问题.pptxVIP

第九章 杆件结构力学问题 9.1 结构单元概论9.2 等截面直杆-梁单元9.3 平面杆件9.4 空间杆件9.5 小结有限元法基础9. 杆件结构力学问题 本章要点结构单元的力学假设经典梁单元和Timoshenko梁单元平面杆系和空间杆系的单元矩阵和向总体坐标的转换有限元法基础9. 杆件结构力学问题 关键概念 结构单元 C1连续 经典梁理论Timoshenko梁理论 平面杆系 空间杆系 剪切自锁（shear locking)有限元法基础9.1 结构单元概论结构单元 杆－梁单元 板壳单元特点：1）在一个方向或二个方向上的几何尺度远小于其他方向；2）引入一定的假设，相应简化为一维或二维问题。结构单元有限元法基础9.1 结构单元概论使用实体单元计算结构问题，理论上是可行的引起的问题1）如果划分网格适应结构特点，使单元的一个方向或两个方向的尺度比其他方向小很多，将使单元刚度系数差异过大，影响求解精度2）保持单元各个方向的尺度接近，将导致单元总数过分庞大有限元法基础9.1 结构单元概论例：一由三个弹簧单元组成的系统，单元刚度为K1＝K2＝105，K3＝1，施加于节点的轴力为P2＝0，P1＝P3＝1。利用5为有效数字的算法求解位移。有限元法基础9.1 结构单元概论系统平衡方程有限元法基础9.1 结构单元概论代入数据高斯消去法求解100000＋1＝100000求解失败有限元法基础9.1 结构单元概论解决办法 1）提高计算有效位数 有限元法基础9.1 结构单元概论2）相对自由度方法(relative DOF)令代入有限元方程 有限元法基础9.1 结构单元概论3）主从自由度方法（master-slave DOF） 分析系统的特点，由于 ，近似地认为将u3作为主自由度，其他为从自由度。只需对u3 进行建立平衡条件由于假设产生了误差。有限元法基础9.1 结构单元概论基于主从自由度原理的梁单元和板壳单元 结构力学中的梁理论和板壳理论引入变形假设实际上正是应用了主从自由度的原理，将问题归结为中面位移。基于相对自由度的梁单元和板壳单元 这类单元 本质上是二维、三维实体单元，为使中面法线变形后保持直线，沿厚度方向只设置两个节点，其中一个保持原来意义上的自由度，另一个节点的自由度用它与上述节点的相对位移来代替。 有限元法基础9.2 等截面直杆-梁单元一. 杆单元 承受轴向拉压载荷，并只经历轴向拉压变形的细长构件，称为杆（Rod）。 假设：应力在截面上均匀分布，原来垂直于轴线的截面变形后仍然垂直于轴线。 有限元法基础9.2 等截面直杆-梁单元直杆的基本方程 几何关系 本构方程 平衡方程 端部条件 有限元法基础9.2 等截面直杆-梁单元最小势能原理 有限元法基础9.2 等截面直杆-梁单元单元列式位移插值函数 其中N＝[N1,N2,…,Nn], qe＝[u1,u2,…,un]T坐标变换2节点单元3节点单元 有限元法基础9.2 等截面直杆-梁单元单元矩阵 由变分原理变分取驻值，得有限元方程其中 有限元法基础9.2 等截面直杆-梁单元二.扭转单元 几何关系 本构关系 平衡方程 端部条件 一般 m(x)＝0，在端部给定扭矩M。 有限元法基础9.2 等截面直杆-梁单元最小势能原理设转角的插值函数为单元刚度矩阵 有限元法基础9.2 等截面直杆-梁单元扭转单元的限制 1）采用自由扭转理论，截面在变形后依然保持平面，这只适用于圆截面情形； 2）在约束扭转时，截面将发生翘曲，需引入翘曲修正因子； 3）只有在截面有两个对称轴时，如圆、椭圆矩形，截面才是绕形心转动。 有限元法基础9.2 等截面直杆-梁单元三.梁单元 经典梁理论假设：变形前垂直于中心线的截面，变形后仍保持为平面，且垂直于中心线。有限元法基础9.2 等截面直杆-梁单元基本方程 几何方程 本构关系 平衡方程 端部条件 有限元法基础9.2 等截面直杆-梁单元最小势能原理本质边界条件：位移端部条件自然边界条件：力的端部条件 端部条件本质边界条件自然边界条件简支固定自由有限元法基础9.2 等截面直杆-梁单元单元列式 梁弯曲是C1类问题，采用Hermite插值 有限元法基础9.2 等截面直杆-梁单元将插值函数代入泛函，由 得到有限元方程单元刚度矩阵为 有限元法基础9.2 等截面直杆-梁单元4.考虑剪切变形的梁单元 经典梁理论只适用于细长梁，即梁高远小于跨度。 变形假设 变形前垂直于中线的截面，变形后仍然保持平面，但不再垂直于中线。有限元法基础9.2 等截面直杆-梁单元剪切变形曲率与转角的关系最小势能原理 折剪系数有限元法基础9.2 等截面直杆-梁单元关于折剪系数 由于变形假设，剪应力在截面上为常数，但是，实际情况剪应力