有限元与有限差分法基础超详细版本.pptVIP

单元平衡方程组装过程再组装单元②总体刚度方程K称为总体刚度矩阵U称为位移向量F称为载荷向量第94页，课件共162页，创作于2023年2月总体刚度矩阵K的特性对称性奇异性稀疏性非零元素带状分布第95页，课件共162页，创作于2023年2月约束处理过程为什么要约束处理?总体平衡方程组是奇异的消除无限制的刚体运动使总体平衡方程组存在唯一一组解第96页，课件共162页，创作于2023年2月约束处理过程—边界条件边界条件分类力(载荷)边界条件位移边界条件集中载荷力表面分布力自重力热交换引起的温度载荷固定位移约束强制位移约束关联位移约束第97页，课件共162页，创作于2023年2月约束处理过程—模型简化xy第98页，课件共162页，创作于2023年2月约束处理过程—模型简化yxxy第99页，课件共162页，创作于2023年2月约束处理过程—约束方程123456789101112yx第100页，课件共162页，创作于2023年2月约束处理过程—约束处理方法位移约束处理方法赋0赋1法乘大数法第101页，课件共162页，创作于2023年2月约束处理过程—赋0赋1法强制位移约束条件处理U4=C第102页，课件共162页，创作于2023年2月约束处理过程—赋0赋1法强制位移约束条件处理U4=C第103页，课件共162页，创作于2023年2月约束处理过程—赋0赋1法有6个方程，5个未知数，如果约束方程可以消除有限元平衡方程组的奇异性，则取任意5个方程联立求解，都会得到方程组的唯一一组解。系数矩阵由原来的对称的变成了非对称的，这对于大规模有限元方程组求解是十分不利的，采用相同的求解方法，在求解时间和矩阵存贮容量方面都增加了一倍。第104页，课件共162页，创作于2023年2月约束处理过程—赋0赋1法为了保证系数矩阵的对称性，去掉方程组第4行第105页，课件共162页，创作于2023年2月约束处理过程—赋0赋1法引入强制位移约束方程U4=C，使方程组求解时直接将自由度U4求出第106页，课件共162页，创作于2023年2月约束处理过程—赋0赋1法固定位移约束条件处理U4=0第107页，课件共162页，创作于2023年2月约束处理过程—赋0赋1法基本原理利用初等变换对求解方程组进行相同的行列变换，既保证方程组解不会改变，又可以保持方程组系数矩阵的对称性。在进行初等变换时，只要保证对方程组系数矩阵做相同的行列变换，就可以保持方程组系数矩阵的对称性。第108页，课件共162页，创作于2023年2月约束处理过程—乘大数法乘大数法基本原理利用矩阵的初等变换不改变方程组解的思想。第109页，课件共162页，创作于2023年2月约束处理过程—乘大数法强制位移边界条件第110页，课件共162页，创作于2023年2月约束处理过程—乘大数法强制约束方程A是一个大数，是系数矩阵中对角线元素K44的1010倍量级以上为什么要乘以大数A？放大位移约束方程的优势第111页，课件共162页，创作于2023年2月约束处理过程—乘大数法强制位移边界条件第112页，课件共162页，创作于2023年2月约束处理过程—乘大数法固定位移边界条件C=0约束后的方程组简化为第113页，课件共162页，创作于2023年2月约束处理过程—两种方法比较赋0赋1法在约束处理过程中是严格精确的，而乘大数法是一种近似约束处理方法，它的精度取决于所乘大数A值两种方法都可以消除有限元平衡方程的奇异性，得到符合实际边界条件的唯一一组解。但两种方法还是有很大的区别第114页，课件共162页，创作于2023年2月单元模型构造—整体坐标系法1.假设插值多项式2.首先，利用节点值求?0、?1和?2二维3节点三角形单元第62页，课件共162页，创作于2023年2月单元模型构造—整体坐标系法A为单元面积第63页，课件共162页，创作于2023年2月单元模型构造—整体坐标系法3.将?0、?1和?2代入插值多项式，按u1、u2、u3合并同类项第64页，课件共162页，创作于2023年2月单元模型构造—整体坐标系法4.同理可得第65页，课件共162页，创作于2023年2月单元模型构造—整体坐标系法5.单元插值多项式为第66页，课件共162页，创作于2023年2月单元模型构造—整体坐标系法6.单元插值多项式写