虚功原理移置单元载荷.PPT

第五节 有限元法应注意的问题 1.对计算对象进行离散化 1）根据对象的特点选择合适的单元： 平面问题单元：三节点三角形、四节点矩形、六节点三角形、八节点矩形。 空间单元：四节点四面体单元、八节点六面体单元等。 2）对划分的单元进行编号，对各单元的节点进行编号。 单元编号原则：起始编号，增量，用①、②、③…表示 节点编号原则：起始编号，增量，相邻节点号差值， 用1、2、3…表示 1.对计算对象进行离散化 3）单元划分的疏密对计算精度和计算成本的影响。 1.对计算对象进行离散化 4）边界曲折、应力集中、应力变化大的部位，单元划分应细化，否则，可划分的稀疏一些。单元由细到疏应当逐步过渡。 1.对计算对象进行离散化 5）对三角形单元，三条边长应尽量接近，不应出现钝角，以免计算结果出现大的误差；对矩形单元，长度和宽度不宜相差过大。 1.对计算对象进行离散化 7）如果计算对象具有不同的厚度或不同的弹性系数，则厚度或弹性系数突变之处应是单元的边线，边线上相邻单元具有共同的节点。 1.对计算对象进行离散化 8）应力分布有突变之处，或者受有应力集中载荷处布置节点，其附近的单元也应划分得细些。 1.对计算对象进行离散化 9）对称性的利用 应充分利用受力体结构的对称性（几何形状和支承条件对某轴对称，同时截面和材料性质也对称此轴）。 2.确定单元的位移模式 单元位移模式的概念：以节点位移为基本未知量，选择一个简单的函数近似表示位移随坐标变化的规律，这个函数就叫位移模式。 常用单元的位移模式 假设： a1、a2…….an为待定常数； u为单元中某节点的x位移； v为单元中某节点的y位移； w为单元中某节点的z位移； x、y、z为某节点的x、y、z坐标。 2.确定单元的位移模式 三节点三角形单元的位移模式： u=a1+a2x+a3y v=a4+a5x+a6y 四节点矩形单元的位移模式： u=a1+a2x+a3y+a4xy v=a5+a6x+a7y+a8xy 2.确定单元的位移模式 六节点三角形单元的位移模式： u=a1+a2x+a3y+a4x2+a5xy+a6y2 v=a7+a8x+a9y+a10x2+a11xy+a12y2 八节点矩形单元的位移模式： u=a1+a2x+a3y+a4x2+a5xy+a6y2+a7x2y+a8xy2 v=a9+a10x+a11y+a12x2+a13xy+a14y2+a15x2y+a16xy2 常用单元位移模式 3.载荷向节点的转移 在进行有限元分析时，受力体要受到体力、集中载荷、分布载荷等，这些载荷必须向节点转移。载荷转移的依据是：圣维南原理、虚功原理。 1）体力向节点的转移 重心的性质：重心到 顶点的距离与重心到 对边交点的距离之比 为2：1。 3.载荷向节点的转移 2）集中载荷向节点的转移 3.载荷向节点的转移 3）分布载荷向节点的转移 4.弹性力学问题的基本解法 1）前面已经学过了弹性力学的平衡微分方程（体力和应力的关系，3个）、几何方程（应变和位移的关系，6个）和物理方程（应力和应变的关系，6个），共有15个方程。而弹性力学问题单元的待求未知参数也为15个（μ，ν，ω，σx， σy， σz ，τxy， τyz， τzx， εx， εy ， εz ，γxy ， γyz ， γzx ）。于是，15个方程有15个未知函数，加上边界条件用于确定积分常数，可求解出各未知函数。 * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * 所以结构总方程为： 其中 考虑到边界条件： 返回 用对角元乘大数法消除奇异性后的结构总体方程为： 返回 由以上方程解得的各节点的位移为： 返回 然后将相应的节点位移代入公式，可分别求得各单元的应 变和应力。 对于单元1： 返回 对于单元2： 返回 对于单元3： 返回 对于单元4： 返回 第十节 热载荷与应力计算 当物体温度发生变化时,物体将由于膨胀而产生线应变aT，其中a为材料的线膨胀系数，T表示弹性体内任意点的温度改变值。在平面问