结构总刚方程的特性及其求解方法.pdf

《结构总刚方程的特性及其求解方法.pdf》由会员分享，可在线阅读，更多相关《结构总刚方程的特性及其求解方法.pdf（8页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、第九章 结构总刚方程的特性及其求解方法一、一、 结构总刚矩阵的特性结构总刚矩阵的特性对称性对称性由： 稀疏性稀疏性1T1KTTKTK KT22当结构的离散单元较多时，对一个节点p 的变形能提供刚度仅与该节点相连的那些单元有关称这些单元上的节点为节点 p 的相关联节点 。相关联节点总比总节点数要少的多，故总体刚度矩阵中每一行的非零元素只占该行元素总数的一少部分。 这样的矩阵称为稀疏矩阵。可利用总刚阵的稀疏性，设法只存储总刚阵中与非零元素相关的那些矩阵元素，可大大节省电脑的内存。带状性带状性带状性是指稀疏矩阵中的元素比较集中地位于对角线元附近。当我们进行结构离散化网格，只要对节点编号稍加注意， 就

2、总能做到使每一个节点与相关联节点的编号比较接近。这样，就可使总刚阵中的非零元素较集中于对角线的两侧。奇异性奇异性有限元的分析计算步骤中，最初形成单元刚阵和组集总刚阵时，未考虑分片差值曲线面对边界条件的满足。从能量原理讲，此时无法获得最小能量解。由此，必须置入适当的边界条件。边界条件的置入是在离散化单元组集后的能量变分后进行的 因为此时获得了明确的结构刚度矩阵方程 。从物理上分析，在置入边界条件前，结构总刚阵是奇异的。二、结构总刚的处理方案二、结构总刚的处理方案奇异消除的置大数法及其原理以前学过强加边界条件的划行划列及置“1”法。虽然这些方法是准确的，但都要对矩阵结构做很多处理，显然不十分方便。

3、置大数法是强加边界条件的一种近似方法，该方法对刚阵的处理却十分简单。/ 置大数法步骤1设给定的边界条件：i aii为一个或多个节点位移在结构总刚对应i行的对角线上，取：Kii D KiiD 为一大数，1020402在载荷列阵相应的行上，取：P a Dii/ 验证置大数法的近似性对应第 i 行刚阵，写出其相应刚度方程：ki11ki22 Diki,i1i1ki.nn aiDDi aiDkijjj1ni aij1nkijDjkij Di ai/ 置大数法可用来近似求解支反力在给定节点位移的点上，节点力是未知的注意在位移与载荷边界同点上的处理原则，在这些点上载荷不再考虑 。在结构总刚度方程中，即对应给

4、定节点位移的那一行刚度方程的右端是未知数。即：ki11ki22kiiaiki,i1i1ki,nn ri由置大数我们获得了全部位移解，此时的方程应为：ki11 ki22(D kii)i ki,nn ri用 D 遍除：_j1jinkijkii_ _rij( (1) )iDDD_ri上式近似于：iriiDD但应注意：i一定是给定零位移条件下有较好的精度，否则精度十分差。如：i1 ri D显然差之千里。/罚函数约束变分原理置大数法可归入罚函数约束变分原理当中，以下简介之。考虑对一个弹性体的无约束变分问题，增加一组约束方程组：C(C(u ) ) 0_u(研究域)该约束条件的置入，可采用lagrange乘

5、子法。 T TC(C(u ) )d_但拉氏乘子法变分原理是以扩大计算变量为代价的， 当然， 可通过对约束方程的了解和推导，建立起拉氏乘子的物理形式，从而获得广义变分原理。但总的来说，计算量扩大或变得难以处理。罚函数约束变分原理没有上述缺点，也不过多改变刚度方程性质，但解的性质是随罚函数的增大而趋向好的近似性。思路及方法：TC Cu 0 C C C C 0及C C C C 0T罚函数约束变分原理是把直接约束条件转换成积分约束，即新泛函为： C C C C d 求解时再不考虑C Cu 0的条件_T其中，是一个罚数，如果本身是要取极小值的，则应当是一个正数因为后面的积分大于或等于零 ， 使泛函取驻值

6、而得到的解只是近似地满足约束条件，值越大，约束条件满足得越好。直观上看，要使有极小值，假设当很大时，如：C C C C d为有限值，则不可能取极小值，故只有当很大时，C C C C 0这不是证明，只能说明这种解的趋势 。举一个例子： 2a12a1a2a218a16a2(a1a2)引入单个约束时不需积分关于参数a1,a2的极小值，即：_222T_T_ 0 (42)a1(22)a218 0a1_ 0 (22)a1(22)a26 0a2可以看到随的增大，a2 a11a1=-12.00a2=-13.52-12.00-13.006-12.00-12.4310-12.00-12.27100-12.00-1

7、2.03在有限元中，假设引入一强制边界条件：iai 02能量泛函修正为： (iai)； _1TKF2_由 0_Kiii2iai fifi fi 2ai这正是上面所讲的置大数法。v/ 结构总刚的存贮方案之一：变带宽一维存贮法因为结构总刚是稀疏对称的，形如右图。为节省内存，可采用存贮每一行从第一个非“0”元到对角线元的方法。各行顺序接起来，并用一管理数组记录每一行从第一个非“0”元到对角线元的个数即可。 这就是常用的一种有限元结构总刚的存贮方案。从第一个非 “0” 元到对角线元结构的元素个数叫做每行的带宽，记为K。 maxk叫作矩阵的包络。KDKM /1,n1n1n1+1n2n2+1n3. . .

8、三、结构总刚方程的计算方法三、结构总刚方程的计算方法1LDL分解法T由对称满秩矩阵的三角分解定理，K可作如下唯一分解：KLDLT其中L11L21L Ln1L22下三角阵 Lnn 1L11D总刚方程的计算：1L221 LnnKPLDLPT记：1TUDLL12L111L13L11L23L221Ln1L11Ln2L22 1PLU1记：UZ. 2得：LZp.2是下三角阵，可采用向下回代过程计算2由方程即：Zm (PmLj1m1mjZj)/Lmm1中是上三角单位阵，可采用向上回代过程求得由方程m Zm( (jm1Unjmj) )T剩余的问题是如何将K进行LDL分解，采用对等法：KLU可只计算下三角阵部分

9、元素即可由此，k11 L11k21 L21k22 L21L12/L11 L22kn1 Ln1kn2 Ln1m1j1L12 Ln2L11Lim kimLijLmj/Ljji 1,n；m 1,i仔细研究上述数值计算，可以发现：1在每行的第一个非“0”元素前，分解过程不产生非“0”元。即L ( j i)元素仅处ij于原总刚的带宽之间，这保证了一维存贮方式的有效性。2在分解 K 阵过程中，每一个k ( j i)元素仅用到过一次，这保证了分解好的元素可ij以放在总刚的原位置上。这样不必开设新的内存空间，仅用原总刚存贮分解好的总刚。3仔细研究以下的分解计算步骤：mLijLkmLjjkLijLkmi1mk由

10、公式：Lim kimj1m1LijLmjLjj从上图中，只能看到第一个非零元 看求和号内元素乘积规律，可以发现：a. 从 1 到m 1列，元素都为“0” ，又kij( j 1,.m1)为“0” 。故不会在带宽之外产生非“0”元。b.Lim kim(第一个非“0”元)c.Lik(km) KikL Lijj1m1mj/Ljjd. 再细加程序判断，即可以形成算法程序。作业：研究中间变量Z及的生成算法细图。2子结构方法我们以前多次谈到内节点和内自由度问题， 其实是指不与其它节点或自由度相关联的那些，或仅在局部单元内部，而不与外部发生联系的那些节点或自由度。如图a 中的点 i 或图b 梁元两端的转角自由

11、度。这个概念可以延伸到一个复杂的结构上。 把任何数量毗连的元素拼在一起成为一个新的结构单元，称之为“超单元”或“子结构” ，用它外边界上的节点和其它的子结构以及一般的元素共同组集成整个结构。这就是“子结构法”的思想。梁元i子结构的方法的关键点是把 “超单元”的那些内部节点或自由度通过等价变换 称为缩象 ，用外边界上的自由度来表征包括节点或载荷列阵 。这对于求解大型复杂结构 例如整架飞机的优点是可大大降低对机器的存贮量。具体子结构法步骤：设：内部节点的编号为i，外边界上点的编号为b由此可形成该超单元的刚度矩阵，即：KbbKibKbiUbSbUSKiiii应当注意，这里的 Si是已知的, 即作用在

12、内节点上的等效载荷。用矩阵块乘法得：KbbUb KbiUi SbKibUb KiiUi Si1212式U K( (S K U ) )由iiiiibb3111式得：(K代入bbKbiKiiKib)Ub SbKbiKiiSi KUb Sb44与其它子结构刚度方程组集，通过计算线性方程组可获得U解，再通过3式由方程b获得Ui解。具体计算思路：*11证：Kib KiiKib，Si* KiiSi*KiiKibTSi*KibSi11这样可以利用LDL方法对Kii进行分解，从而利用矩阵分解法获得KiiKib及KiiSi，11而不必计算Kii的逆。这是因为Kii的稀疏性较好，而Kii则可能是满秩的，而利用分解

13、法求解方程，可保持矩阵的稀疏性，而不必用求逆方法计算。四、关于应力计算四、关于应力计算在位移模式的有限元方法中，第一步获得的解是结点位移，由结点位移通过单元的几何矩阵求得应变、应力，即BDB对于应力计算在位移有限元方法中应当有几点概念：1应力的连续性程度较位移的连续性程度要低一个阶次，特别在单元的边界上，可能发生应力的不连续的状态。 如最明显的常应力单元导致问题的最主要根源在于形状函数的选择，特别是分片插值函数的选择。2应力计算的精度较位移计算的精度要低一个阶次一个数量级 ，主要误差来源有二， 一是计算过程的数值误差的传递；二是插值函数的导数精确逼近。3在等参元当中，可以证明 Gauss 点f

14、 (x)f1(x)f2(x)上的应力与其他点上的应力计算精度要高一个阶次。因此，在等参元中一般仅计算 Gauss 点上的应力，代表单元的应力水平。4对于应力不连续的情形可采用多种修正方案来处理，使得应力计算趋于合理化。第十章现代有限元数值计算技术的发展1各种能量原理的分区组合应用，以求解关心域的高精度参数结果。2各种杂交单元的设计，获得高精度的计算结果。3有限元程序设计技术的发展，内存的全稀疏管理，结点编码优化技术。4有限元技术与 CAD 技术的结合，形成多个Commercial 软件，加强了前后置处理功能，减轻了数据准备的工作，增加了计算结果的可视化功能。5有限元方法深入到各种复杂问题及各种领域。 有限元并行计算技术。