- 1、本文档共23页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
结构有限元刚度方程求解基础
结构有限元法刚度阵的特征:
1)一般维数很大
2 )“0”元素很多
3 )非“0”元素多集中在矩阵的对角线附近
4 )一般来说矩阵是对称正定的(特殊的如结构流体耦合问题则刚度阵非对称)
有限元的灵魂:
求解 方程组
已知力 F 求位移 δ 的情况,则 ,实际计算中,我们不直接求 [K]
的逆矩阵,而把 [K] 转化为上三角或下三角矩阵,最后回代求解方程。
线性方程组的解法主要有直接法和迭代法,以下介绍直接法。
1三角方程组
以下矩阵均只考虑方阵。
1) 向前消去
考虑下面方程组
如果 l11l22 ≠0, 则未知数可确定如下 :
这就是所谓的向前消去 , 其一般形式为
算法 1.1 (向前消去:行形式) Ax=b 的解覆盖 b
b(1)=b(1)/L(1,1)
for i=2~n
b(i)=(b(i)-L(i,1~i-1)b(1~i-1))/L(i,i)
end
2) 向后消去
解上三角方程组 的类似算法叫向后消去法 ,
算法 1.2 (向后消去:行形式) Ux=b 的解覆盖 b
b(n)=b(n)/L(n,n)
for i=n-1,1,-1
b(i)=(b(i)-U(i,i+1~n)b(i+1~n))/U(i,i)
end
向后消去法的算法实现 (Fortran):
注 : 数组 du(i) 表示向量 ,whlf(i) 表示向量 , 二维数组 whlk(i,j) 表示矩
阵 ,
变量 neqns 表示方程数 , 临时变量 temp
du(neqns)=whlf(neqns)/whlk(neqns,neqns)
do 100 i=neqns-1, 1, -1
temp=0.
do 200 j=i+1, neqns
temp = temp + whlk(i,j)*du(j)
200 enddo
du(i) = (whlf(i)-temp)/whlk(i,i)
100 enddo
基于列的形式:交换循环顺序可得到以上算法的列形式,考虑向前消去,一旦
x1 解出来,该变量可以从第 2~n 个方程中去掉,我们可只考虑缩小后的方程组
然后我们算出 x2 ,并且从第 3~n 个方程中去掉 x2 ,依次类推,例如:
我们有 x1=3,那么我们处理 2x2 方程组
算法 1.3 (向前消去:列形式) Lx=b 的解覆盖 b
for j=1~n-1
b(j)=b(j)/L(j,j)
b(j+1~n)=b(j+1~n)-b(j)L(j+1~n,j)
end
b(n)=b(n)/L(n,n)
算法 1.4 (向后消去:列形式) Ux=b 的解覆盖 b
for j=n,2,-1
b(j)=b(j)/U(j,j)
b(1~j-1)=b(1~j-1)-b(j)U(1~j-1,j)
end
b(1)=b(1)/U(1,1)
2LU 分解
如前面所见,三角方程比较容易求解,那么我们可以把刚度阵 [K] 经过适当的线
性变换等价成一个三角方程组— Gauss消去法
例:在方程组
中,将第一个方程乘以 2,并且在第二个方程中减去它则:
矩阵形式:
这就是 n=2 的 Gauss 消去法
此算法线性变换过程也可写成矩阵形式, 则最终求出一个下三角矩阵 L 和一个上
三角矩阵 U,使得 A=LU,
然后原始问题 Ax=b 可通过两个三角求解过程求得:
Ly=b,Ux=y = Ax=LUx=Ly=b
在 n=2 的水平上,如果 x1 ≠0且 τ=x2/x1, 那么
更一般的,设 且 ,如果
定义:
则
形如 Mk 的矩阵的前 k 个分量都为 0 时,一般称高斯变换,这样的矩阵是下三角
的,元素 称为乘子,向量 称为高斯向量。
上三角化,设
您可能关注的文档
- 经济学基础曲线大全.pdf
- 经纬度计算和转换工具说明.pdf
- 结构力学第十章习题集.pdf
- 给排水计算软件使用说明.pdf
- 绝缘栅场效应管(IGFET)的基本知识.pdf
- 维修改造工程竣工验收汇报材料施工总结报告.pdf
- 2025届河南省周口市扶沟高级中学高考化学全真模拟密押卷含解析.doc
- 陕西省延安中学2025届高考考前模拟化学试题含解析.doc
- 2025届林芝市重点中学高考化学四模试卷含解析.doc
- 云南省云天化中学2025届高考历史考前最后一卷预测卷含解析.doc
- 四川省德阳市罗江中学2025届高三考前热身化学试卷含解析.doc
- 山东省枣庄现代实验学校2025届高三下学期第五次调研考试化学试题含解析.doc
- 吉林省长春市十一高中等九校教育联盟2025届高三一诊考试生物试卷含解析.doc
- 2025届江苏省盐城市伍佑中学高考仿真模拟化学试卷含解析.doc
- 2025届广西贺州中学高考冲刺押题(最后一卷)生物试卷含解析.doc
- 安徽省池州市贵池区2025届高三第一次模拟考试生物试卷含解析.doc
- 宁夏银川一中2025届高三(最后冲刺)化学试卷含解析.doc
- 广东省广州市增城区四校联考2025届高考压轴卷化学试卷含解析.doc
- 2025届邯郸市第一中学高考生物必刷试卷含解析.doc
- 2025届安徽省安庆市石化第一中学高考仿真卷化学试卷含解析.doc
最近下载
- DBJ50T-136-2012建筑地基基础检测技术规范(高清版).pdf
- 【中职-情景模拟】24.学生、家长不认同劳动教育,对卫生值日、劳动实践等有抵触心理.docx VIP
- 小学三年级《生命安全教育》全册教案(湖北版).pdf VIP
- 《Unit 1 Teenage Life Reading for Writing》教案(附导学案)1.docx
- 桩基础土方开挖施工方案.docx VIP
- 恶性骨肿瘤患者的护理ppt课件.pptx
- 苏教版四年级上册同步奥数培优 第十六讲 等量代换.pdf VIP
- 2024年中国兵器工业集团招聘3774人公开引进高层次人才和急需紧缺人才笔试参考题库(共500题)答案详解版.docx VIP
- 老年冠心病慢病管理指南.pptx VIP
- 人教版普通高中地理选择性必修2区域发展.pdf
文档评论(0)