第五章 三角形单元有限元法.pptVIP

将单元位移???代入式（1-28），其中[S]可由式（1-30）、（1-31）求得，从而可求出单元应力。 （1-28） 习题：算出式（1-28）的结果. 剖分结构为有限个单元，对节点、单元编号 组集单元刚度矩阵和单元等价节点力向量 组集结构刚度矩阵并引入约束 组集结构等价节点荷载向量 从节点平衡方程解未知节点位移 计算结构内力、应力 有限元分析之大腕版 一定要选最变态的题目，什么材料非线性啊，几何非线性啊，接触非线性啊，多物理场耦合啊，都给他弄进去。是个模型就几百万个单元，上千万个节点，画个剖面图就要十几个小时。再整一并行机群，TOP500的，张口就是 High Performance Computation，一口地道的伦敦腔，倍儿有面子。题目一扔进去就跑个把月，你要是一个星期以内出结果，你都不好意思和别人打招呼。你说这样一趟算下来要发多少Paper? 10篇？10篇？就1篇！你还别嫌少，说不定人家还发在会议上。你得琢磨牛人的心理啊，有能耐算这样题目的人，根本就不在乎多发一篇两篇文章。什么叫大牛知道么？大牛就是不求灌水，但求经典。 式(1-48）表明：在ij边上受均布面力的平面问题三角形单元，其等价节点力等于将均布面力合力之半简单地简化到i、j节点上，方向与分布力方向相同。m节点上为零。 (1-48） i j m x y qsx Fs1 Fs3 i j m x y qsy Fs2 Fs4 （3）线性分布面力 i j m 图1-11 x y s 表面力集度在i点为[qsx qsy]T，而在j点为0。设坐标轴s的原点取在j点，沿ji为正向， 。 ij边上任一点的面力集度?qs? s qsi qs i j m 图1-12 x y s l 在ij边上有： 将?qs?和上式代入式(1-46），有 由形函数的性质3： (1-49） 式(1-49）表明：ij边受线性分布面力： i点为[qsx, qsy]T，j点为0 时，其等价节点力可将总载荷的2/3分配给i点，1/3分配给j点，m点为零得出。 x y i j m qsi qs 体积力和表面力向节点的移置符合静力等效原理的前提条件是：线性位移模式。 5.7 系统分析 5.7.1 坐标系 研究各离散单元集合成整体结构，集合整体结构的平衡和变形协调，建立整体结构平衡方程。 单元分析时采用的坐标系成为局部坐标或单元坐标（单元刚度矩阵的通用性）。而结构系统分析时，必须在统一的坐标系内进行（各力学量才能叠加），称为“结构坐标”或“整体坐标”，如图1-13所示。 单元坐标系下，单元位移、单元力、单元刚度矩阵表示为： 整体坐标系下，单元位移、单元力、单元刚度矩阵表示为： X Y X Y P ○ ○ ○ ○ ○ P 图（1-13） (a) 平面桁架(杆件单元) 悬臂深梁 (平面三角形单元) x y x y x y x y 如何从单元坐标转化为结构坐标将在第4章中讨论。 5.7.2 整体刚度矩阵 假设整体结构被划分为ne个单元和n个节点，在整体坐标系下，对于每个单元均有： 将上述这些方程集合起来（整体坐标下叠加），便可得到整个结构的平衡方程。为此，需要将[k]、{δ}、{F}体积膨胀，分别扩大为n1×n1、n1×1和n1×1的矩阵才能相加。膨胀后，原有节点号对应位置的元素不变，而其它元素均为零。 组装方法：建立一个体积为n1×n1的方阵，按单元序号依次把结构坐标单元刚度矩阵的元素放入该方阵中。 放入方法：（1）按单元节点编码对号入座； （2）同位置元素累加。 式中：[K]为整体刚度矩阵，{Δ}为整体节点位移列阵；{P}为整体等价节点荷载列阵。如下： （1-50） i j m i j m 例：平面三角单元 双行双列 5.7.3 结构刚度矩阵特性 1、结构刚度矩阵元素的力学意义 把方程(1-50）写开， =1 =0 =0 =0 =0 =0 (1-51） 2、结构刚度矩阵是对称矩阵 已知单元刚度矩阵是对称矩阵（5.7节），用单元刚度矩阵组集 结构刚度矩阵的 过程又没有破坏 其对称性，结构 刚度矩阵必然也 是对称的。当然， 对称性也可以通 过虚功原理得到证明。 结构刚度矩阵中的任一元素kij是?j为单位位移（ ?j =1），其它位移为零时的Pi。 3、结构刚度矩阵主对角线上的元素恒为正值 由性质(1)可知，任一主对角线上元素kii是使节点位移?i为一单位位移，其它节点位移为零时必须在第i