有限元基础理论教程

有限元基础理论教程第1页，课件共24页，创作于2023年2月先考虑结点n有水平方向位移约束，与n结点水平方向对应的平衡方程为：在[K]矩阵中，第2n-1行的对角线元素改为1，该行中全部非对角线元素改为0；在{P}中，第2n-1个元素改为0。为了保持[K]矩阵的对称性，将第2n-1列的全部非对角元素也改为0。应该换成下面的方程：第2页，课件共24页，创作于2023年2月如果结点n在垂直方向有位移约束，则（2-29）中的第2n个方程修改为，在[K]矩阵中，第2n行的对角线元素改为1，该行中全部非对角线元素改为0；在{P}中，第2n个元素改为0。为了保持[K]矩阵的对称性，将第2n列的全部非对角元素也改为0。第3页，课件共24页，创作于2023年2月例2.9、结构的位移约束条件如图所示，结构平衡的方程组如下。修改整体刚度矩阵。第4页，课件共24页，创作于2023年2月根据结点4的位移约束，修改该整体刚度矩阵的系数。第5页，课件共24页，创作于2023年2月根据结点1和结点6的位移约束条件继续修改整体刚度矩阵可以得到以下的形式，第6页，课件共24页，创作于2023年2月如果结点n处存在一个已知非零的水平方向位移，这时的约束条件为，在[K]矩阵中，第2n-1行的对角线元素乘上一个大数A，将方程修改为，第7页，课件共24页，创作于2023年2月处理约束条件的意义在于，强迫单元内部的分片近似位移场满足整体结构的位移边界条件，得到整个结构的近似解。在定义位移约束时，要消除结构的刚体位移。第8页，课件共24页，创作于2023年2月1）用相关单元的单元刚度矩阵计算结点力，再由力的平衡关系得到约束反力。与位移约束对应的约束反力如何计算？2）给矩阵相应的对角元素加上一个大数，将载荷列阵的对应元素置为零。相当于用一个刚度很大的弹簧代替位移约束。n结点在x方向位移为零，用刚度系数为A的弹簧代替，由于刚度很大，un是一个很小的值。方程修改为，第9页，课件共24页，创作于2023年2月2.8整体刚度矩阵的特点与存储方法

整体刚度矩阵具有以下几个显著的特点：对称性，稀疏性，非零系数带形分布。1）对称性 由单元刚度矩阵的对称性和整体刚度矩阵的集成规则，可知整体刚度矩阵必为对称矩阵。利用对称性，只保存整体矩阵上三角部分的系数即可。2）稀疏性

单元刚度矩阵的多数元素为零，非零元素的个数只占较小的部分。第10页，课件共24页，创作于2023年2月3）非零元素带形分布整体刚度矩阵的非零元素分布在以对角线为中心的带形区域内，这种矩阵称为带形矩阵。在包括对角线元素的半个带形区域内，每行具有的元素个数叫做半带宽。最大半带宽用d表示，例2.9所示结构的最大半带宽为，第11页，课件共24页，创作于2023年2月二维等带宽存储

设整体刚度矩阵[K]为一个n行、n列的矩阵，最大半带宽为d。利用带形矩阵的特点和对称性，只需要保存以d为固定带宽的上半带的元素，称为二维等带宽存储。 整体刚度矩阵[K]每行中的上半带元素取出，保存在另一个矩阵[K*]的对应行中，得到一个n行、d列矩阵[K*]。 把元素在[K]矩阵中的行、列编码记为r、s，在矩阵[K*]中的行、列编码记为r*、s*，对应关系如下：

r*=r s*=s-r+1第12页，课件共24页，创作于2023年2月[K]矩阵中的对角线元素保存在新矩阵中的第1列中，[K]矩阵中的r行元素仍然保存在新矩阵的r行中，[K]矩阵中的s列元素则按照新的列编码保存在新矩阵的不同列中。仍然会保存一些零元素，但是元素寻址很方便。整体刚度矩阵[K]等带宽矩阵[K*]第13页，课件共24页，创作于2023年2月二维等带宽存储所存储的元素数量取决于最大半带宽d的值，d的值则由单元结点的编码方式决定。相同的有限单元网格按照图2.13(a)的结点编码，最大的半带宽为14；按照图2.13(b)的结点编码，最大的半带宽为18。图2.13(a)图2.13(b)第14页，课件共24页，创作于2023年2月一维变带宽存储按行的一维变带宽存储，按照每一行的半带宽把半带宽内的元素保存到一维数组中。使用辅助数组定义对角元素在一维数组中的位置。需要存储的矩阵元素最少。解方程时的地址计算比较复杂，会带来一些额外的计算量。第15页，课件共24页，创作于2023年2月2.9线性方程组解法线性方程组的解法：直接解法

包括高斯消去法、等带宽高斯消去法、三角分解法，以及适用于大型方程组求解的分块算法和波前法等。迭代解法 雅可比迭代法、高斯-赛德尔迭代法、超松弛迭代法和共轭梯度法等。第16页，课件共24页，创作于2023年2月利用矩阵的对称、稀疏、带状分布等特点提高方程求解效率是关键。在方程组的阶数不是特别高时，通常采用直接解法。当方程组的阶数过高时，为避免舍入误差和消元时有效数损失等对计算精度的影响，可以选择迭代方法。第17页，课件共24页，创作于2023年2月利用矩阵的对称、稀疏、带状分布等特点提高方程求解效率是关键。在方程组的阶数不是特别高时，通常采用直接解法。当方程组的阶数过高时，为避免舍入误差和消元时有效数损失等对计算精度的影响，可以选择迭代方法。第18页，课件共24页，创作于2023年2月高斯消去法的一般公式

第m次消元时，以第m-1次消元后的第m行元素为主元行，对第i行元素(i>m)的消元公式为，若原系数矩阵是对称矩阵，则在消元过程中的待消矩阵仍然保持对称。

如果第19页，课件共24页，创作于2023年2月雅可比迭代的原理A为非奇异矩阵，主元非零D为对角矩阵通过迭代逐步逼近方程组的解，第20页，课件共24页，创作于2023年2月波前法的特点 刚度矩阵K和载荷列阵P不按照自然编号进入内存，而按照参加运算的顺序排列；以集成完毕的自由度作为主元对其它列的元素进行消元修正。不形成整体刚度矩阵。共轭梯度法 求解最优化问题的搜索算法。迭代公式为，搜索方向，线性方程组，第21页，课件共24页，创作于2023年2月ANSYS提供了多种求解器供选择，分为直接解法和迭代解法。直接解法包括：波前法（FrontalSolver）稀疏法（SparseDirectSovler）迭代解法包括：雅可比共轭梯度法（JacobiConjugateGradientSolver，JCG），不完全共轭梯度法（IncompleteCholeskyConjugateGradientSolver，ICCG）预处理共轭梯度法（PreconditionedConjugateGradientSolver，PCG）代数多格法（AlgebraicMultigridSolver，AMG）区域分割法（DistributedDomainSolver，