有限元程序课程设计

重庆大学本科学生课程设计任务书课程设计题目有限元程序设计学院资源及环境科学学院专业工程力学年级2010级参数和设计要求：1.独立完成有限元程序设计。2.独立选择计算算例，并能通过算例判断程序的正确性。3.独立完成程序设计报告，报告内容包括理论公式、程序框图、程序本体、计算算例，算例结果分析、结论等。学生应完成的工作：复习掌握有限单元法的根本原理。掌握弹性力学平面问题3节点三角形单元或4节点等参单元有限元方法的计算流程，以及单元刚度矩阵、等效节点载荷、节点应变、节点应力和高斯积分等的计算公式。用Fortran语言编写弹性力学平面问题3节点三角形单元或4节点等参单元的有限元程序。在VisualFortran程序集成开发环境中完成有限元程序的编辑和调试工作。利用编写的有限元程序，计算算例，分析计算结果。撰写课程设计报告。目前资料收集情况〔含指定参考资料〕：1.王勖成，有限单元法，北京：高等教育出版社，2002。2.O.C.Zienkiewicz,R.L.Taylor,FiniteElementMethod,5thEition,McGraw-HallBookCompanyLimited,2000。3.张汝清，董明，结构计算程序设计，重庆：重庆大学出版社，1988。课程设计的工作方案：第1周星期一上午：教师讲解程序设计方法，程序设计要求和任务安排。第1周星期一至星期二完成程序框图设计。第1周星期三至第2周星期四完成程序设计。第2周星期五完成课程设计报告。任务下达日期2013年6月6日完成日期2013年07月03日指导教师〔签名〕学生〔签名〕一、前言有限单元法是在当今工程分析中获得最广泛应用的数值计算方法，由于其通用性和有效性，受到工程技术界的高度重视。伴随着计算机科学和技术的快速开展，现已成为计算机辅助设计和计算机辅助制造的重要组成局部。由于有限元法是通过计算机实现的，因此有限元程序的编制及相关软件的研发就变得尤为重要。从二十世纪五十年代以来，有限元软件的开展按目的和用途可分为专用软件和大型通用商业软件，而且软件往往集成了网格自动划分，结果分析和显示等前后处理功能，而且随着时间的开展，大型通用商业软件的功能由线性扩展到非线性，由结构扩展到非结构等等，这一系列强大功能的实现与运用都要求我们对有限元法的根底理论知识有较为清楚的认识以及对程序编写的根本能力有较好掌握。1.元分析程序的理论根底作为弹性力学微分方程的等效积分形式，虚位移原理和虚应原理分别是平衡方程与力的边界条件和几何方程与位移边界条件的等效积分形式。将物理方程引入虚位移原理和虚应力原理可以分别导出最小位能原理和最小余能原理，它们本质上和等效积分的伽辽金“弱”形式相一致，这是建立弹性力学有限原方程一般表达格式的理论根底。对于通过弹性力学变分原理建立的弹性力学问题有限元方法，其未知场变量是位移，以节点位移为根本未知量，并以最小位能原理为根底建立的有限单元为位移元。弹性力学平面问题有限元分析表达格式的建立步骤一般为：a．的模型进行单元离散，一般采用三角形单元或四边形单元，不过由于四边形单元具有更高精度，应用更为普遍，一般而言一次或二次单元已足以满足精度要求。离散后对单元进行编号，并且给定单元各节点的整体编码以及局部坐标系下按逆时针局部编码。b．坐标系下对各单元构造形函得到由单元各节点位移表示的单元位移形式，进而得到单元刚度矩阵，利用等参元性质和雅阁比矩阵进行组集，建立整体刚度矩阵。c．等效结点载荷列阵并组集成结构节点等效载荷列阵得到单元各节点位移与结构节点等效载荷列阵的线代方程组，求解得到位移，进而可得应力应变。根据以上一般步骤，编制相应的计算机程序并采用数值积分方法处理有关数学计算就可以得到完整的弹性力学问题有限元分析程序。问题有限元分析教学程序本程序可对二维弹性力学问题行有限元分析计算。计算采用的单元形式为四边形四节点单元或者四边形八节点单元，对于对称和非对称矩阵采用变带宽存储方法，最终输出结果为单元各节点的位移。二、平面4、8节点有限元公式及计算原理2.1根本公式有限元平衡方程单元刚度矩阵〔3〕单元等效结点载荷2.2单元位移插值及坐标变换通过Serendipity四边形单元格式构造插值函数。对于4节点单元，插值函数为：对于8节点单元，插值函数为：其中：〔i=1,2,3,4〕位移插值坐标变换2.3单元刚度的计算〔1〕弹性矩阵平面应力问题：；平面应变问题：；式中E和v分别为材料Young氏模量和Poisson比。应变矩阵形状函数对局部坐标偏导：，形状函数对整体坐标偏导：Jacobi矩阵：Jacobi矩阵的逆：其中：单元刚度矩阵计算单元刚度矩阵Gauss积分：单元刚度矩阵子矩阵：〔i=1，2〕2.4等效结点载荷集中载荷直接施加在结点上体积力产生的等效结点载荷向量体积力产生的等效结点载荷Gauss积分：体积力产生的等效结点载荷子向量：其中：假设体积力作用方向于整体坐标系的y轴正方向角度，那么：2.5单元应变和应力计算单元应变单元中任意点的应变：单元中任意点的应变：其中可作为Gauss积分点或单元结点的整体坐标，其于局部坐标间的关系为：,单元应力三、程序结构〔1〕有限元程序系统的组成及分析过程我们本次设计主要针对二维弹性力学静力学问题。有限元程序系统通常包括处理、有限元程序本体和后处理三局部。有限元软件系统的组成及分析过程如下列图所示：CAD模型CAD模型、材料参数、边界条件、分析类型定义、网格划分等。有限元分析本体程序前处理后处理核心：各种计算方法的实现。以图形、曲线和表格等方式表达、分析计算结果。〔2〕有限元计算程序框图输入离散模型数据输入离散模型数据型数据计算单元刚度阵组集结构刚度矩阵计算单元等效结点载荷组集结构结点载荷列阵引入位移边界条件求解线性方程组其它辅助计算输出结果结束单元循环形成K形成P消除K的奇异求解Ka=P得结点位移a计算应力、应变等〔3〕有限元程序主要编程技巧a.结构模块化主程序：组织调用子程序，过程清晰，不涉及复杂计算。子模块：实现特定功能，可用多个子程序组成。子程序：实现简单功能，便于调用和扩展。b.覆盖技术在主程序中定义一个整数组和一个实型数组，分别存放整型和实型边界数组。设计动态数组表。动态数组覆盖技术。〔4〕有限元程序主程序翻开输入数据，调用PMESH，调用PROFIL计算K的对角元地址和方程号，分配数组大小，调用ASSEM组集整体K，调用PLOAD组集载荷列向量，调用DATRI对刚度矩阵进行三角分解，调用DASOL解方程，调用PRTDIS输出结点位移，调用TRESSS计算GAUSS点和结点应变和应力。〔5〕变量说明NUMNP：节点总数NUMEL：单元总数NUMMAT：材料种类NLOAD：载荷数NDM：每个节点坐标方向数NDF：节点自由度NEN：每个单元节点数。ID〔K，J〕：节点J的第K个自由度约束编号X〔K，J〕：节点J的第K个坐标IX〔K，J〕：J单元中K节点的全局节点编号F〔K，J〕：结点J第K方向上的集中载荷F〔NDF，NLOAD〕：集中载荷X〔NDM，NUMNP〕：结点坐标IX〔NEN，NUMEL〕：单元定义F(K,J)：节点J在K方向上集中载荷大小ALFL=TRUE：非对称矩阵集合ALFL=FALSE：对称矩阵集合EE：杨氏弹性模量XNU：泊松比ITYPE：问题类型Ｓ：单元刚度矩阵Ｐ：载荷和内力向量AD：对角元素AU：上三角元素AL：下三角元素JP：上〔下〕三角元素每行〔列〕最后元素LD：一个单元中各自由度方程数JD：列高ID：边界条件，=0：无约束；=1：固定约束。〔6〕数据输入格式整数：I5格式小数：F10.4格式〔7〕程序本体〔见附录〕程序验证算例算例：平面板受轴向拉伸长3m宽0.4m的平面板左端铰接固定，右端受到水平方向向右、大小为10N的均布载荷，模型如下列图所示：F=10NF=10NF=10N1284519F=10NF=10NF=10N1284519其中结构参数：F=10N,L=3m,h=0.4m材料参数：E=200GPa,ν=0.31、建立坐标系，将板划分成6个如上图所示的四边形四节点单元，进行单元编号，均布载荷处理为右端三节点每个节点受水平向右大小10N的力，竖直方向不受力。1.1程序input文件中输入：example(NEN=4)126132241110.00000.00002001.00000.00003002.00000.00004003.00000.00005110.00000.20006001.00000.20007002.00000.20008003.00000.20009110.00000.400010001.00000.400011002.00000.400012003.00000.400011265223763348745610956711106781211200.0E+90.30001410.00000.0000810.00000.00001210.00000.00001.2调试程序得到结果为：nodaldisplacements10.0000E+000.0000E+0020.3638E-090.2813E-1030.7369E-090.2030E-1040.1126E-080.2437E-1050.0000E+000.0000E+0060.3609E-09-0.2981E-2370.7429E-09-0.1095E-2280.1104E-08-0.2132E-2290.0000E+000.0000E+00100.3638E-09-0.2813E-10110.7369E-09-0.2030E-10120.1126E-08-0.2437E-102、建立坐标系，将板划分成6个如上图所示的四边形八节点单元，进行单元编号，均布载荷处理为右端三节点每个节点受水平向右大小10N的力，竖直方向不受力。其中结构参数：F=10N,L=3m,h=0.4m材料参数：E=200GPa,ν=0.32.1程序input文件中输入：example(NEN=8)296132281110.00000.00002001.00000.00003002.00000.00004003.00000.00005110.00000.20006001.00000.20007002.00000.20008003.00000.20009110.00000.400010001.00000.400011002.00000.400012003.00000.400013000.50000.000014001.50000.000015002.50000.000016110.00000.100017001.00000.100018002.00000.100019003.00000.100020000.50000.200021001.50000.200022002.50000.200023110.00000.300024001.00000.300025002.00000.300026003.00000.300027000.50000.400028001.50000.400029002.50000.400011265131720162237614182117334871519221845610920242723567111021252824678121122262925200.0E+90.30001410.00000.0000810.00000.00001210.00000.00002.2调试程序得到结果为：nodaldisplacements10.0000E+000.0000E+0020.3734E-090.2492E-1030.7549E-090.4768E-1040.1185E-080.6785E-1050.0000E+000.0000E+0060.3770E-09-0.1562E-2270.7371E-09-0.3614E-2280.1086E-08-0.7600E-2290.0000E+000.0000E+00100.3734E-09-0.2492E-10110.7549E-09-0.4768E-10120.1185E-08-0.6785E-10130.1849E-090.2673E-10140.5567E-090.1616E-10150.9277E-090.1568E-11160.0000E+000.0000E+00170.3737E-090.1232E-10180.7517E-090.2430E-10190.1114E-080.3177E-10200.1810E-09-0.5215E-23210.5559E-09-0.2565E-22220.9462E-09-0.5238E-22230.0000E+000.0000E+00240.3737E-09-0.1232E-10250.7517E-09-0.2430E-10260.1114E-08-0.3177E-10270.1849E-09-0.2673E-10280.5567E-09-0.1616E-10290.9277E-09-0.1568E-113.数据分析其理论计算过程如下，对于平面应力问题有：QUOTEQUOTE其中QUOTE：代入计算得任一距原点水平距离为x的点其位移为：理论解与数值解比照图如下：可见接近自由端处两曲线很接近，当X=3m时，u=15x10E-9m，对应输出文件中的节点3、7、11的输出值和3、7、11点有限元解较为接近，误差为很小。同时可以发现靠近固定端处理论解与数值解误差变大。分别对x=2处的理论解和有限元解进行比拟发现误差较小。考虑到实际情况固定端的铰支约束条件对实际位移的影响以及离散单元数目较少等可能造成的误差，可以认为该程序是正确的。另外数值由可以看出同的对称性，符合根本力学理论结果。程序可用于平面问题的求解，可获得较好的位移分布。程序适用于平面应力、平面应变状态结构4、8节点单元有限元分析，可求出结构的位移、应力以及应变。五、算例1、算例：选择超静定梁中点受集中载荷作用的结构，如下列图所示：其中结构参数：F=10kN,L=4m,h=0.4m材料参数：E=200GPa,ν=0.32、单元网格划分及规定对上述结构采用10×1,10×2三种单元网格划分形式，分别采用四结点和八结点四边形单元。规定单元编号先沿y正方向，再沿x正方向。为了方便比拟，所有编号位置固定，例如第12结点，不管是10×1单元四结点，还是10×1单元八结点，其位置不变。其不同之处在于八结点增加的结点编号在四结点的根底上插入，而不是替换原有的结点编号。特殊点如下所示：1161121722116112172212331162817123221233116281712322(1).10X1单元四结点input文件中输入：WSJ_10X1(NEN=4)2210112241110.00000.00002000.40000.00003000.80000.00004001.20000.00005001.60000.00006002.00000.00007002.40000.00008002.80000.00009003.20000.000010003.60000.000011114.00000.000012114.00000.400013003.60000.400014003.20000.400015002.80000.400016002.40000.400017002.00000.400018001.60000.400019001.20000.400020000.80000.400021000.40000.400022110.00000.400011221222232021334192044518195561718667161777815168891415991013141010111213200.0E+90.30001170.0000-10.0E+310X1单元四结点input文件中输出：nodaldisplacements10.0000E+000.0000E+002-0.2029E-06-0.2673E-063-0.3048E-06-0.8380E-064-0.3055E-06-0.1509E-055-0.2059E-06-0.2083E-0560.2994E-21-0.2330E-0570.2059E-06-0.2083E-0580.3055E-06-0.1509E-0590.3048E-06-0.8380E-06100.2029E-06-0.2673E-06110.0000E+000.0000E+00120.0000E+000.0000E+0013-0.2015E-06-0.2671E-0614-0.3018E-06-0.8376E-0615-0.3012E-06-0.1510E-0516-0.1986E-06-0.2077E-0517-0.1853E-21-0.2365E-05180.1986E-06-0.2077E-05190.3012E-06-0.1510E-05200.3018E-06-0.8376E-06210.2015E-06-0.2671E-06220.0000E+000.0000E+00(2).10X1单元八结点input文件中输入：WSJ_10X1(NEN=8)5310112281110.00000.00002000.40000.00003000.80000.00004001.20000.00005001.60000.00006002.00000.00007002.40000.00008002.80000.00009003.20000.000010003.60000.000011114.00000.000012114.00000.400013003.60000.400014003.20000.400015002.80000.400016002.40000.400017002.00000.400018001.60000.400019001.20000.400020000.80000.400021000.40000.400022110.00000.400023000.20000.000024000.60000.000025001.00000.000026001.40000.000027001.80000.000028002.20000.000029002.60000.000030003.00000.000031003.40000.000032003.80000.000033114.00000.200034003.60000.200035003.20000.200036002.80000.200037002.40000.200038002.00000.200039001.60000.200040001.20000.200041000.80000.200042000.40000.200043110.00000.200044000.20000.400045000.60000.400046001.00000.400047001.40000.400048001.80000.400049002.20000.400050002.60000.400051003.00000.400052003.40000.400053003.80000.40001122122234244432232021244145423341920254046414451819263947405561718273848396671617283749387781516293650378891415303551369910131431345235101011121332335334200.0E+90.30001170.0000-10.0E+310X1单元八结点input文件中输出：nodaldisplacements10.0000E+000.0000E+002-0.2896E-06-0.3882E-063-0.4446E-06-0.1198E-054-0.4464E-06-0.2169E-055-0.2954E-06-0.2984E-056-0.1067E-19-0.3336E-0570.2954E-06-0.2984E-0580.4464E-06-0.2169E-0590.4446E-06-0.1198E-05100.2896E-06-0.3882E-06110.0000E+000.0000E+00120.0000E+000.0000E+0013-0.2881E-06-0.3880E-0614-0.4419E-06-0.1198E-0515-0.4415E-06-0.2170E-0516-0.2928E-06-0.2992E-05170.1160E-19-0.3406E-05180.2928E-06-0.2992E-05190.4415E-06-0.2170E-05200.4419E-06-0.1198E-05210.2881E-06-0.3880E-06220.0000E+000.0000E+0023-0.1582E-06-0.1219E-0624-0.3874E-06-0.7546E-0625-0.4638E-06-0.1684E-0526-0.3906E-06-0.2617E-0527-0.1666E-06-0.3254E-05280.1666E-06-0.3254E-05290.3906E-06-0.2617E-05300.4638E-06-0.1684E-05310.3874E-06-0.7546E-06320.1582E-06-0.1219E-06330.0000E+000.0000E+00340.7195E-09-0.3635E-06350.1556E-08-0.1195E-05360.2019E-08-0.2174E-05370.3347E-08-0.3008E-05380.7962E-21-0.3388E-0539-0.3347E-08-0.3008E-0540-0.2019E-08-0.2174E-0541-0.1556E-08-0.1195E-0542-0.7195E-09-0.3635E-06430.0000E+000.0000E+00440.1574E-06-0.1217E-06450.3852E-06-0.7543E-06460.4599E-06-0.1684E-05470.3863E-06-0.2615E-05480.1608E-06-0.3252E-0549-0.1608E-06-0.3252E-0550-0.3863E-06-0.2615E-0551-0.4599E-06-0.1684E-0552-0.3852E-06-0.7543E-0653-0.1574E-06-0.1217E-06(3).10X2单元四结点input文件中输入：WSJ_10X2(NEN=4)3320112241110.00000.00002000.40000.00003000.80000.00004001.20000.00005001.60000.00006002.00000.00007002.40000.00008002.80000.00009003.20000.000010003.60000.000011114.00000.000012114.00000.400013003.60000.400014003.20000.400015002.80000.400016002.40000.400017002.00000.400018001.60000.400019001.20000.400020000.80000.400021000.40000.400022110.00000.400023110.00000.200024000.40000.200025000.80000.200026001.20000.200027001.60000.200028002.00000.200029002.40000.200030002.80000.200031003.20000.200032003.60000.200033114.00000.2000112242322325243342625445272655628276672928778302988931309910323110101133321123242122122425202113252619201426271819152728171816282916171729301516183031141519313213142032331213200.0E+90.30001170.0000-10.0E+310X2单元四结点input文件中输出：nodaldisplacements10.0000E+000.0000E+002-0.2089E-06-0.2816E-063-0.3174E-06-0.8637E-064-0.3188E-06-0.1558E-055-0.2117E-06-0.2146E-056-0.8065E-20-0.2403E-0570.2117E-06-0.2146E-0580.3188E-06-0.1558E-0590.3174E-06-0.8637E-06100.2089E-06-0.2816E-06110.0000E+000.0000E+00120.0000E+000.0000E+0013-0.2074E-06-0.2813E-0614-0.3146E-06-0.8635E-0615-0.3136E-06-0.1558E-0516-0.2072E-06-0.2139E-05170.5719E-20-0.2439E-05180.2072E-06-0.2139E-05190.3136E-06-0.1558E-05200.3146E-06-0.8635E-06210.2074E-06-0.2813E-06220.0000E+000.0000E+00230.0000E+000.0000E+0024-0.7077E-09-0.2647E-0625-0.1592E-08-0.8609E-0626-0.1668E-08-0.1561E-0527-0.5226E-08-0.2158E-0528-0.1099E-20-0.2429E-05290.5226E-08-0.2158E-05300.1668E-08-0.1561E-05310.1592E-08-0.8609E-06320.7077E-09-0.2647E-06330.0000E+000.0000E+00(4).10X2单元八结点input文件中输入：WSJ_10X2(NEN=8)8520112281110.00000.00002000.40000.00003000.80000.00004001.20000.00005001.60000.00006002.00000.00007002.40000.00008002.80000.00009003.20000.000010003.60000.000011114.00000.000012114.00000.400013003.60000.400014003.20000.400015002.80000.400016002.40000.400017002.00000.400018001.60000.400019001.20000.400020000.80000.400021000.40000.400022110.00000.400023110.00000.200024000.40000.200025000.80000.200026001.20000.200027001.60000.200028002.00000.200029002.40000.200030002.80000.200031003.20000.200032003.60000.200033114.00000.200034000.20000.000035000.60000.000036001.00000.000037001.40000.000038001.80000.000039002.20000.000040002.60000.000041003.00000.000042003.40000.000043003.80000.000044114.00000.100045003.60000.100046003.20000.100047002.80000.100048002.40000.100049002.00000.100050001.60000.100051001.20000.100052000.80000.100053000.40000.100054110.00000.100055000.20000.200056000.60000.200057001.00000.200058001.40000.200059001.80000.200060002.20000.200061002.60000.200062003.00000.200063003.40000.200064003.80000.200065114.00000.300066003.60000.300067003.20000.300068002.80000.300069002.40000.300070002.00000.300071001.60000.300072001.20000.300073000.80000.300074000.40000.300075110.00000.300076000.20000.400077000.60000.400078001.00000.400079001.40000.400080001.80000.400081002.20000.400082002.60000.400083003.00000.400084003.40000.400085003.80000.400011224233453555422325243552565333426253651575244527263750585155628273849595066729283948604977830294047614888931304146624799103231424563461010113332434464451123242122557476751224252021567377741325261920577278731426271819587179721527281718597080711628291617606981701729301516616882691830311415626783681931321314636684672032331213646585662.00E+110.30001170.0000-10.0E+310X2单元八结点input文件中输出：nodaldisplacements10.0000E+000.0000E+002-0.3021E-06-0.3948E-063-0.4523E-06-0.1231E-054-0.4533E-06-0.2217E-055-0.3066E-06-0.3052E-056-0.6288E-19-0.3413E-0570.3066E-06-0.3052E-0580.4533E-06-0.2217E-0590.4523E-06-0.1231E-05100.3021E-06-0.3948E-06110.0000E+000.0000E+00120.0000E+000.0000E+0013-0.3006E-06-0.3946E-0614-0.4493E-06-0.1231E-0515-0.4483E-06-0.2217E-0516-0.2990E-06-0.3059E-05170.5240E-19-0.3488E-05180.2990E-06-0.3059E-05190.4483E-06-0.2217E-05200.4493E-06-0.1231E-05210.3006E-06-0.3946E-06220.0000E+000.0000E+00230.0000E+000.0000E+0024-0.7495E-09-0.3869E-0625-0.1469E-08-0.1227E-0526-0.2015E-08-0.2219E-0527-0.7048E-09-0.3060E-0528-0.6477E-20-0.3443E-05290.7048E-09-0.3060E-05300.2015E-08-0.2219E-05310.1469E-08-0.1227E-05320.7495E-09-0.3869E-06330.0000E+000.0000E+0034-0.1663E-06-0.1340E-0635-0.3955E-06-0.7793E-0636-0.4707E-06-0.1723E-0537-0.3970E-06-0.2668E-0538-0.1670E-06-0.3315E-05390.1670E-06-0.3315E-05400.3970E-06-0.2668E-05410.4707E-06-0.1723E-05420.3955E-06-0.7793E-06430.1663E-06-0.1340E-06440.0000E+000.0000E+00450.1448E-06-0.3846E-06460.2216E-06-0.1228E-05470.2224E-06-0.2219E-05480.1473E-06-0.3061E-0549-0.3274E-19-0.3433E-0550-0.1473E-06-0.3061E-0551-0.2224E-06-0.2219E-0552-0.2216E-06-0.1228E-0553-0.1448E-06-0.3846E-06540.0000E+000.0000E+0055-0.3552E-09-0.1062E-0656-0.1107E-08-0.7656E-0657-0.1899E-08-0.1724E-0558-0.2739E-08-0.2682E-0559-0.7310E-08-0.3343E-05600.7310E-08-0.3343E-05610.2739E-08-0.2682E-05620.1899E-08-0.1724E-05630.1107E-08-0.7656E-06640.3552E-09-0.1062E-06650.0000E+000.0000E+0066-0.1433E-06-0.3845E-0667-0.2186E-06-0.1228E-0568-0.2181E-06-0.2219E-0569-0.1418E-06-0.3065E-05700.2054E-19-0.3470E-05710.1418E-06-0.3065E-05720.2181E-06-0.2219E-05730.2186E-06-0.1228E-05740.1433E-06-0.3845E-06750.0000E+000.0000E+00760.1655E-06-0.1338E-06770.3933E-06-0.7791E-06780.4670E-06-0.1723E-05790.3924E-06-0.2668E-05800.1669E-06-0.3311E-0581-0.1669E-06-0.3311E-0582-0.3924E-06-0.2668E-0583-0.4670E-06-0.1723E-0584-0.3933E-06-0.7791E-0685-0.1655E-06-0.1338E-063.变形图用ANSYS软件划分网格绘制变形图如下：A.10X1单元.10X1单元四结点(2).10X1单元八结点B.10X2单元，同理可以画出变形图4.有限元程序Fortran计算结果比拟10X1单元四结点、10X1单元八结点、10X2单元四结点、10X2单元八结点的output比拟如下〔选取局部结点比拟〕：(1).10X1单元四结点和10X1单元八结点比拟〔选取结点6和结点17〕10X1单元四结点：60.2994E-21-0.2330E-0517-0.1853E-21-0.2365E-0510X1单元八结点：6-0.1067E-19-0.3336E-05170.1160E-19-0.3406E-05(2).10X1单元四结点和10X2单元四结点比拟〔选取结点6和结点17〕10X1单元四结点：60.2994E-21-0.2330E-0517-0.1853E-21-0.2365E-0510X2单元四结点：6-0.8065E-20-0.2403E-05170.5719E-20-0.2439E-05(3).10X1单元四结点和10X2单元八结点比拟〔选取结点6和结点17〕10X1单元四结点：60.2994E-21-0.2330E-0517-0.1853E-21-0.2365E-0510X2单元八结点：6-0.6288E-19-0.3413E-05170.5240E-19-0.3488E-05求解理论解此算例的弯矩图如下列图，横截面尺寸如下列图：由材料力学知识可以得到弯矩图，由图乘计算可以求出l/2处位移为：其中：，，，，，代入可求得：分析因结构是两端固定，所以在中间位移变形最大，取梁中间上下两点〔分别是17结点和6结点〕进行结果分析。分别用Fortran、ANSYS、理论解表示。由于17结点和6结点两点处于中间位置，对称，所以X方向位移为零。结点X位移类型10X1-410X1-810X2-410X2-817〔Fortran〕-0.2365E-05-0.3406E-05-0.2439E-05-0.3488E-0517〔ANSYS〕-0.236E-05-0.341E-0517〔理论解〕-0.3125E-05-0.3125E-05-0.3125E-05-0.3125E-05由上表可看出，Fortran和ANSYS计算结果几乎一样，所以只要分析Fortran的计算结果和理论解之间的误差。百分误差如下表：10X1-410X1-810X2-410X2-817〔Fortran〕-0.2365E-05-0.3406E-05-0.2439E-05-0.3488E-0517理论解-0.3125E-05-0.3125E-05-0.3125E-05-0.3125E-0517结点百分误差24.3%8.99%21.9%11.6