有限元法及其在岩体力学中的应用.pdfVIP

第19卷第1期 西昌学院学报·自然科学版 V01．19．No．1 Science Edition 2005年3月 XichangCoUege·Natural Mar一2005 有限元法及其在岩体力学中的应用 伍迪，袁前胜，余明东，胡建春 (西昌学院，四川西昌615013) 【摘 要】本文介绍了有限元法的基本思想，并探讨了用有限元分析法分析解决岩体工程中 计算的理论和步骤，目的在于简化复杂曲计算问题。 【关键词】有限元；岩体力学；计算 【中图分类号】0241．82【文献标识码】A 1有限元法的基本思想 f(xn-。)) 下面利用上述近似计算公式计算定积分： 有限元法的基本思想是用较简单的方法来处理 ，2 2 复杂问题。它将求解域看成若干个有限元的小子域 J。eIdx (单元)构成，对于每一个单元假设一个较简单的近 的近似值。 似解，然后推导求解该域应满足的条件，进而得到问 设精确度s=O．0005，为了达到精度要求，这里采 题的近似解；它能把复杂问题简单化，把无限问题有 用变步长梯形法，先把区间[O，2]分为2等份，若达不 限化；这对于解决许多难以求得精确解的实际问题 到精度要求，再分为4等份，依次下去，直到满足精度 的求解非常有效。有限元法的求解不仅计算精度高， 要求为止。现用计算机来处理，其计算程序如下： 而且能适应各种复杂问题，因而已成为工程领域的 #includemath．h 有效分析手段之一。 main() float double 有限元是那些集合在一起的能够表示成实际连 {inti； a=0，b=2，n=2，h，ep； 续域的离散单元。例如：在定积分应用中，定积分： ex0，exl； ，b h=(b—a)／2；ep=0．0005；exl=0．O； f(x)dx J 8 do 是由曲线y=f(x)与直线x=a、x=b和y=0围成的曲 {ex0=exl；exl=(exp(a：lca)+exp(b宰b))／2； 边梯形的面积，在近似计算中把该曲边梯形等分为n for(i-1；in；i++) 个相邻小曲边梯形，即把区间[a，b]n等份，每个小 曲边梯形的面积用相应的小梯形(或小矩形)的面积 exl=exl宰h；h=h／2；n=n术2；} 代替；第i个小梯形的两个底分别为f(xi一。)和f(x：)，小while(fabs(exl一exO)=ep)； 梯形的高h=(b—a)／n，则得到第i个小梯形的面积 printf(‘‘jifenjinshizhi_％f”，ex1)；} h(f(xH)“(xi))／2其中：a=xox1…xi_1)【i …x。=b；i=1，2，…，n) 这里用的近似计算方法，其实就是用了有限元 于是得到定积分的近似值： 方法的基本