下载文档原格式
有限元上机实验报告[——以Abaqus软件进行的有限元分析]汪健强1008320139实验1——平面问题应力集中分析目的要求：掌握平面问题的有限元分析方法和对称性问题建模的方法。
通过简单力学分析，可以知道本实验问题属于平面应力问题，基于结构和载荷的对称性，可以只取模型的1/4进行分析。
用8节点四边形单元分析X=0截面σx的分布规律和最大值，计算圆孔边的应力集中系数，并与理论解对比。
一、实验过程概述：1、启动ABAQUS/CAE2、创建部件3、创建材料和截面属性4、定义装配件5、设置分析步6、定义边界条件和载荷7、划分网格8、提交分析作业9、后处理10、退出ABAQUS/CAE二、实验结果：（1）边界受力图（1）X方向应力分量σx应力云图：（2）左边界直线与圆弧边交点的σx值为： 2.96714 MPa；（2）左右对称面上的σx曲线：三、实验内容分析：a)模型全局σx应力分布：σx应力集中分布于中心圆孔与x、y轴相交的地方，且与x轴相交处应力为负，与y轴相交处应力为正；沿圆周向周围，σx迅速减小；沿y 方向的σx应力大于沿x方向的σx应力。
b)应力集中系数为 2.92975，小于理论值3.0。
误差来源：有限元分析方法是将结构离散化，网格划分得越稀疏，计算出的结果就越偏离理论值。
分的越密集，结果越接近与理论值。
四、实验小结与体会：通过本次实验，对理论课所学有限元基本方法有了一个更加直观、深入的理解。
通过对Abaqus软件三个步骤：前处理、分析计算、后处理的操作，了解了这款软件的基本应用和它对有限元的一些很好的应用。
试验中，遇到诸多问题，仔细思考，加之请教老师，逐一解决，确实很有收获。
更增加了对有限元的认识，和对其功能之强大有了更深的理解。
实验二平面问题有限元解的收敛性一、实验目的和要求：（1）在ABAQUS软件中用有限元法探索整个梁上σx和σy的分布规律。
（2）计算梁底边中点正应力σx的最大值；对单元网格逐步加密，把σx的计算值与理论解对比，考察有限元解的收敛性。
一，实验描述：1、本作业属于哪类问题在本作业中，根据板的结构特点和受力情况，确定该问题属于平面对称应力问题，定义分析类型为静力学分析。
2、本文采用如何的单位制本题中，长度单位为m，故为方便起见，采用力的单位为N，压强的单位为pa，时间的单位为s，质量的单位为kg。
3、单元类型：对单元描述；材料；实常数单元类型为选取shell Elastic 4node63单元材料：弹性模量为2.1e11pa，泊松比为0.33。
由题意确定时常数，即厚度为0.1m4、划分网格。
网格划分设置。
单元数，节点数。
网格化分设置：设置单元边长值为0.1m，指定单元形状为Areas。
5、加载描述（1）对整体模型，首先对四周固定端加载位移约束，设定其位移值为0；然后，对面施加面载荷，设置面载荷为20000N；接着将施加在实物体上的载荷转换到有限元模型上，并显示施加在有限元模型上的载荷。
加载完后，对该有限元模型进行求解。
（2）对四分之一模型，首先对两边固定端加载位移约束，设定其位移值为0；然后，对两坐标轴所在的边施加对称载荷，最后对面施加面载荷，设置面载荷为20000N；接着将施加在实物体上的载荷转换到有限元模型上，并显示施加在有限元模型上的载荷。
加载完后，对该有限元模型进行求解。
6、后处理：最大MISIS应力和最大位移的位置和大小。
绘制结构的应力和变形图。
二，实验步骤（一）绘制整体实体模型1，在ANSYS中构造实体模型，如下图所示2，根据结构特点及所受载荷地情况，选取shell Elastic 4node63单元，设置材料常数：弹性模量E=2.1e11，u=0.33，单元边值为0.1m对其进行网格剖分，网格划分图如下：3，正确施加载荷和边界条件，结果如下：边界条件施加载荷20000N/m求解以后4，绘制平板的应力和变形图，并给出最大应力和变形的位置及大小：应力图应变图从图中可以看出其中实体边界中点位置的应力最大为0.583e+07pa，最大变形在中间圆弧的位置，0.144e-3m（二）绘制四分之一的实体模型1，在ANSYS中构造实体模型，如下图所示2,根据结构特点及所受载荷地情况，选取shell Elastic 4node63单元，设置材料常数：弹性模量E=2.1e11，u=0.33，单元边值为0.1m对其进行网格剖分，网格划分图如下：3，正确施加载荷和边界条件，注意此处有一个对称载荷的加载，结果如下：4，绘制平板的应力和变形图，并给出最大应力和变形的位置及大小：应力图应变图三，实验小结这次ANSYS上机实验课是使用shell中的Elastic 4node63单元，这使我对ANSYS软件中的单元有了更深的认识，同时对平面问题的静力分析的基本思路和操作步骤更加熟悉。
有限元上机实验报告学生专业学生学号学生姓名实验日期南京理工大学机械工程学院一、实验设备机械工程软件工具包Ansys二、实验主要流程和步骤（1）建立有限元模型的几何、输入模型的物理和材料特性、边界条件和载荷的描述、模型检查的整个过程。
具体操作如下： ①定义文件名 ②建模③选用单元类型 ④设定单元的厚度 ⑤设定材料属性 ⑥离散几何模型 ⑦施加位移约束 ⑧施加压强⑨查看最后的有限元模型（2）对建立的有限元模型选择相应的求解器进行求解运算。
（3）对计算结果进行考察和评估，比如绘制应力、变形图，将结果与失效准则进行比较等。
习题11、已知条件简支梁如图3.1.1所示，截面为矩形，高度h=200mm ，长度L=1000mm ，厚度t=10mm 。
上边承受均布载荷，集度q=1N/mm 2，材料的E=206GPa ，μ=0.29。
平面应力模型。
X 方向正应力的弹性力学理论解如下：)534()4(622223-+-=h y h y q y x L h q x σ2、目的和要求（1）在Ansys 软件中用有限元法探索整个梁上x σ，y σ的分布规律。
（2）计算下边中点正应力x σ的最大值；对单元网格逐步加密，把x σ的计算值与理论解对比，考察有限元解的收敛性。
（3）针对上述力学模型，对比三节点三角形平面单元和4节点四边形平面等参元的求解精度。
3、实验步骤（1） 定义文件名， （2）建模，（3）选用单元类型 （4） 设定单元的厚度 （5） 设定材料属性 （6） 离散几何模型 （7）施加位移约束 （8） 施加压强（9） 查看最后的有限元模型 （10） 提交计算 （11） 查看位移（12） 查看模型X 方向应力（13） 查看X 方向上的应力关于X 轴的位移图模型图1MNMXXY Z0.116E-06.232E-06.348E-06.464E-06.580E-06.696E-06.812E-06.927E-06.104E-05APR 13 201309:15:22NODAL SOLUTION STEP=1SUB =1TIME=1USUM (AVG)RSYS=0DMX =.104E-05SMX =.104E-05位移云图1MNMXXY Z-188808-147068-105329-63589-218501989061629103369145108186848APR 13 201309:23:35NODAL SOLUTION STEP=1SUB =1TIME=1SX (AVG)RSYS=0DMX =.104E-05SMN =-188808SMX =186848应力云图1107.036283.180459.324635.468811.612987.7561163.9001340.0441516.1881692.3321868.479(x10**2) 0.1.2.3.4.5.6.7.8.91DISTAPR 13 201309:32:04POST1STEP=1SUB =1TIME=1PATH PLOT NOD1=1NOD2=2X1X 向应力关于X 轴位移图 三角单元三角单元模型1MNMXXY Z0.964E-07.193E-06.289E-06.386E-06.482E-06.578E-06.675E-06.771E-06.867E-06APR 13 201309:42:17NODAL SOLUTION STEP=1SUB =1TIME=1USUM (AVG)RSYS=0DMX =.867E-06SMX =.867E-06三角单元位移图1MNMXXY Z-129669-100854-72038-43223-14408144084322372038100854129669APR 13 201309:43:16NODAL SOLUTION STEP=1SUB =1TIME=1SX (AVG)RSYS=0DMX =.867E-06SMN =-129669SMX =129669三角单元应力云图1104.842224.027343.212462.397581.582700.767819.952939.1371058.3221177.5071296.688(x10**2) 0.1.2.3.4.5.6.7.8.91DISTAPR 13 201309:46:38POST1STEP=1SUB =1TIME=1PATH PLOT NOD1=1NOD2=2X1三角单元X 向应力关于X 轴位移图1MNMXXY Z-158263-123094-87924-52754-17585175855275487924123094158263APR 13 201309:50:47ELEMENT SOLUTION STEP=1SUB =1TIME=1SX (NOAVG)RSYS=0DMX =.867E-06SMN =-158263SMX =158263X 向应力中间最大两边小，有限元解只是一种数值近似与理论解还是有误差的。
有限元分析上机报告上海电机学院机械学院目录Project1 桁架结构静力有限元分析 (1)Project2 梁结构静力有限元分析 (3)Project3 平面结构静力有限元分析 (5)Project4 实体结构静力有限元分析 (6)有限元分析上机报告Project1 桁架结构静力有限元分析上机地点上机日期指导教师班级学号姓名一、上机目的二、主要仪器设备三、上机结果1、根据任务要求提交图片或者数据；2、提交.log文件；3、验证有限元分析结果。
四、体会与建议有限元分析上机报告Project2 梁结构静力有限元分析上机地点上机日期指导教师班级学号姓名一、上机目的二、主要仪器设备三、上机结果1、根据任务要求提交图片或者数据；2、提交.log文件；3、验证有限元分析结果。
四、体会与建议有限元分析上机报告Project3 平面结构静力有限元分析上机地点上机日期指导教师班级学号姓名一、上机目的二、主要仪器设备三、上机结果1、根据任务要求提交图片或者数据；2、提交.log文件。
四、体会与建议有限元分析上机报告Project4 实体结构静力有限元分析上机地点上机日期指导教师班级学号姓名一、上机目的二、主要仪器设备三、上机结果1、根据任务要求提交图片或者数据；2、提交.log文件。
四、体会与建议出师表两汉：诸葛亮先帝创业未半而中道崩殂，今天下三分，益州疲弊，此诚危急存亡之秋也。
然侍卫之臣不懈于内，忠志之士忘身于外者，盖追先帝之殊遇，欲报之于陛下也。
诚宜开张圣听，以光先帝遗德，恢弘志士之气，不宜妄自菲薄，引喻失义，以塞忠谏之路也。
宫中府中，俱为一体；陟罚臧否，不宜异同。
若有作奸犯科及为忠善者，宜付有司论其刑赏，以昭陛下平明之理；不宜偏私，使内外异法也。
侍中、侍郎郭攸之、费祎、董允等，此皆良实，志虑忠纯，是以先帝简拔以遗陛下：愚以为宫中之事，事无大小，悉以咨之，然后施行，必能裨补阙漏，有所广益。
将军向宠，性行淑均，晓畅军事，试用于昔日，先帝称之曰“能”，是以众议举宠为督：愚以为营中之事，悉以咨之，必能使行阵和睦，优劣得所。
轴流式通风机叶轮与机座有限元分析分析与优化报告书第2 页共47 页目录第一部分机座的有限元分析与优化—-———--—--—--—--———--——---——--——--—- 41。
1 机座分析的已知条件--—--—--—--—-----—-———---—-————--—-—-——-—— 41。
2 材料的力学性能--—--——-—-——--———-——-—--——---—--------—-————--- 41。
3 有限元分析模型——-—-—--—-—--—------——----———-————-———------—-- 41.3.1 分析前的假设--——-——-——---—-———-——-—---———-—---—-————— 41。
3.2 建立分析模型—--—-————--———---—————--—--—-————-——---—— 51。
3.3 建立有限元分析模型—-——-——-————---———--———-----—--—-- 71.4 计算结果——----——----—--—--—--—————---------———-—————————-—---— 71.4.1 变形结果———---—-——-—-—--——-------——-------—-——————-—-—- 71.4.2 应力结果-——-—--————-----——-—-——--—-—--—-——-—--————----— 81.4。
3 路径结果—-——-----——-—----——-—---—-—-—-———--——--————---- 111。
4。
4 分析结果评判-———-----———-----——-———-—-----——--—--—--—- 131.5 机座优化-———-—---—————-—-------——--——--——--——-——-—---——--—---- 141.5。
1 优化参数的确定—-—-—--—---—-——------——--——-----————-—— 141.5。
有限元法基础及应用上机报告南京理工大学2015年12月上机实验一1 实验题目设计一个采用减缩积分线性四边形等参元的有限元模型，通过数值试验来研究网格密度、位移约束条件与总刚度矩阵奇异性、沙漏扩展、求解精度的关系，并验证采用减缩积分时保证总刚度矩阵非奇异的必要条件。
总结出你的研究结论，撰写实验报告。
2 实验目的通过实验来研究减缩积分方案中网格密度和位移约束条件对总体刚度矩阵奇异性和求解精度的影响，以此加深对有限元减缩积分的理解，和对减缩积分中保证总体刚度矩阵非奇异性的认识。
3建模概述先保持位移约束条件不变，研究网格密度对总体刚度矩阵奇异性和求解精度的影响，并验证采用减缩积分时保证总刚度矩阵非奇异的必要条件。
如下图1所示，建立一个简支和链杆的约束条件，然后不断增加网格密度，通过ABAQUS 来计算位移和应力的变化规律。
简支（两个约束）链杆（一个约束）积分点（3个独立关系式）节点（两个自由度）4 计算结果分析讨论与结论 1）1*1单元四边形减缩积分实验载荷 布种/单元应力云图2）2*1单元四边形减缩积分实验载荷 单元应力云图3）4*4单元四边形减缩积分实验载荷布种单元应力云图结果分析5 实验体会与小结单元刚度矩阵的特征：(1)对称性(2)奇异性(3)主元恒正K相同(4)平面图形相似、弹性矩阵D、厚度t相同的单元，eK的分块子矩阵按结点号排列，每一子矩阵代表一个结点，占两行两(5)e列，其位置与结点位置对应。
整体刚度矩阵的特征：(1)对称性(2)奇异性(3)主元恒正(4)稀疏性(5)非零元素呈带状分布。
[K]的物理意义是任意给定结构的结点位移所得到的结构结点力总体上满足力和力矩的平衡。
为消除[K]的奇异性，需要引入边界条件，至少需给出能限制刚体位移的约束条件。
对于一个给定形式的单元，如果采用精确积分，则插值函数中所有项次在|J|=常数的条件下能被精确积分，并能保证刚度矩阵的非奇异性。
如果采用减缩积分，因为插值函数中只有完全多项式的项次能被精确积分，因此需要进行刚度矩阵非奇异必要条件的检查。
有限元作业一上机报告20120430231 余凯峰T1243-2一、题目图示折板上端固定，右侧受力F=1000N，该力均匀分布在边缘各节点上；板厚t=2mm，材料选用低碳钢，弹性模量E=210Gpa，μ=0.33.二、有限元分析的目的1、利用ANSYS构造实体模型；2、根据结构的特点及所受载荷的情况，确定所用单元类型；正确剖分网格并施加外界条件；3、绘制结构的应力和变形图，给出最大应力和变形的位置及大小；并确定折板角点Ａ处的应力和位移；4、研究网格密度对Ａ处角点应力的影响；5、若在Ａ处可用过渡圆角，研究Ａ处圆角半径对Ａ处角点应力的影响三、有限元模型的特点1、作业类型：平面应力2、单位制：N mm MPa3、单元类型：在Preferences选Structural，Preprocessor>ElemmentType>Add/Edit/Delete中定义单元类型为Plane-82，K3设置为Plane strs w/thk材料: EX=2.10E5MPa, PRXY=0.33 实常数：THK=2mm4、划分网格：Preprocessor—Meshing—Mesh Tool—Size Controls下的Global 点击Set—在第一行输入你要定义的网格长度—OK—点击Clear旁边的Mesh—Pick All这一步划分网格5、加载描述：Preprocessor—Loads—Define Loads—Apply—Structural—Displacement—On Lines—拾取模型被固定的上边—ok—All DOF—okPreprocessor—Loads—Define Loads—Apply—Structural—Pressure—On Lines—拾取右边线—Load pres valu输入-1000/60—ok6、后处理：网格密度设置为“3”时应力和位移图如下：网格密度设置为“4”时应力和位移图如下网格密度设置为“5”时应力和位移图如下移也减小但变化不大。
有限元上机作业系所：同济大学土木工程学院专业：建筑与土木工程 .学号： ####### .姓名： ## .指导教师： ### .二零一五年一月二十八日目录第1题 (3)1 问题描述 (3)1.1运用弹性力学方法求解 (3)1.2 Abaqus有限元分析 (4)1.2.1 材料只考虑弹性时 (5)1）创建部件 (5)2）划分网格 (5)3）创建材料和截面属性 (6)4）装配部件 (6)5）设置分析步 (6)6）定义边界条件和荷载 (7)7）提交作业 (8)8）后处理 (8)9）误差分析 (8)10）误差验证 (8)1.2.2材料考虑弹塑性时 (9)1.2.3 对比分析 (11)第2题 (12)2.1 问题描述 (13)2.2 Abaqus建立模型 (13)1）创建部件 (13)2）划分网格 (14)3）创建材料和截面属性 (14)4）装配部件 (14)5）设置分析步 (15)6）定义边界条件和荷载 (15)7）提交作业 (16)2.3结果分析 (16)第3题 (19)3.1 问题描述 (19)3.2 Abaqus建立模型 (20)3.2.1对12边同时施加x和y方向的位移约束 (20)1）创建部件 (20)2）划分网格 (20)3）创建材料和截面属性 (21)4）装配部件 (21)5）设置分析步 (21)6）定义荷载和边界条件 (21)7）提交作业 (22)8）结果分析 (22)3.2.2对12边施加x方向的位移约束,对12边的中间一点施加y方向的位移约束 (22)3.2.3结果对比分析 (23)1）Mises应力分析 (23)2）S11应力分析 (24)3）S22应力分析 (24)4）S12应力分析 (25)5）Max principle应力分析 (25)6）Mid principle应力分析 (26)7）Min principle应力分析 (26)8）结论： (27)第1题（原题号5） 通过Abaqus 有限元软件模拟具有小圆孔的平板均匀拉伸产生的应力集中问题，验证弹性力学中的结论以及对比材料只考虑弹性和考虑弹塑性的区别。
有限元分析基础结课报告任课教师：聂志峰学生姓名：XXX学号：XXXXXXXXXXXX班级：XXXXXXXXXXXX4m 5 m2m水深4m习题1：选用Plane82单元分析如图1所描述的水坝受力情况，设坝体材料的平均密度为2g/cm3，考虑自重影响，材料弹性模量为E=700Mpa, 泊松比为0.3。
按水坝设计规范，在坝体底部不能出现拉应力。
分析坝底的受力情况，是否符合要求。
解：（1）思路：建模和分析过程参考上机指南中的Project2。
（2）建模和分析：从已知条件知，此计算类型为Structural力学类型；由于考虑自重的影响，故需设定密度和施加重力载荷；单元类型选择Solid Quad 4node 42；定义材料参数为EX:2.1e11, PRXY:0.3（根据已知条件）；生成几何模型利用点生成面方式；网格划分参照Project2；模型施加约束，坝体的底部施加x和y的约束，其余部位不施加约束，载荷在坝体的右端施加水的压力载荷，施加方式9800*{4-{y}};最后分析计算，查看应力图和变形图结果，保存数据。
图1 水坝截面图（3）ANSYS软件分析过程：1.1进入ANSYS程序→ANSYSED 10.0 →Interactive →change the working directory into yours →input Initial jobname: dam→Run1.2设置计算类型ANSYS Main Menu: Preferences →select Structural →OK1.3选择单元类型ANSYS Main Menu: Preprocessor →Element Type→Add/Edit/Delete →Add →select Solid Quad 4node 42 →OK (back to Element Types window)→Options… →select K3: Plane Strain →OK→Close (the Element Type window)1.4定义材料参数ANSYS Main Menu: Preprocessor →Material Props →Material Models →Structural→Linear→Elastic→Isotropic→input EX:2.1e11, PRXY:0.3→OKANSYS Main Menu→Preprocessor →Material Props →Material Models→Structural →Density →Dens: 20001.5生成几何模型✓生成特征点ANSYS Main Menu: Preprocessor →Modeling →Create →Keypoints →In Active CS→依次输入四个点的坐标：input:1(0,0,0),2(4,0,0),3(2,4,0),4(0.,5,0)→OK✓生成坝体截面ANSYS Main Menu: Preprocessor →Modeling →Create →Areas →Arbitrary →Through KPS→依次连接四个特征点，1(0,0,0),2(4,0,0),3(2,4,0),4(0.,5,0) →OK1.6网格划分ANSYS Main Menu: Preprocessor →Meshing →Mesh Tool→(Size Controls) lines: Set →依次拾取两条横边:OK→input NDIV: 15 →Apply→依次拾取两条纵边:OK →input NDIV: 20 →OK →(back to the mesh tool window)Mesh: Areas, Shape: Quad, Mapped →Mesh →Pick All (in Picking Menu) →Close( the Mesh Tool window)1.7模型施加约束✓分别给下底边施加x和y方向的约束ANSYS Main Menu: Solution→Define Loads →Apply→Structural →Displacement→On lines →pick the lines(拾取坝体的下底边) →OK→select Lab2:UX, UY →OK✓给竖直边施加x方向的分布载荷ANSYS 命令菜单栏: Parameters→Functions →Define/Edit→1) 在下方的下拉列表框内选择x ,作为设置的变量；2) 在Result窗口中出现{X},写入所施加的载荷函数：9800*{4-{y}}；3) File>Save(文件扩展名：func) →返回：Parameters→Functions →Read from file：将需要的.func文件打开，任给一个参数名，它表示随之将施加的载荷→OK →ANSYS Main Menu: Solution→Define Loads →Apply→Structural →Pressure →On Lines →拾取斜边；OK →在下拉列表框中，选择：Existing table →OK →选择需要的载荷参数名→OK施加重力载荷ANSYS Main Menu →Solution →Define Loads →Apply →Structural →Gravity →ACEL Y: 9.8 →OK1.8 分析计算ANSYS Main Menu: Solution →Solve →Current LS→OK(to close the solve Current Load Step window) →OK1.9 结果显示ANSYS Main Menu: General Postproc →Plot Results→Deformed Shape…→select Def + Undeformed→OK (back to Plot Results window)→Contour Plot→Nodal Solu…→select: DOF solution, UX,UY, Def + Undeformed , Stress ,SX,SY,SZ, Def + Undeformed→OK1.10 退出系统ANSYS Utility Menu: File→Exit…→Save Everything→OK（4）结果分析：从应力图可以看出，无论是在X和Y方向的应力值都是负值，在其方向是受压应力的作用(无拉应力的作用)。
（5）变形图和应力图：图2 变形图（无网格）图3 X方向的应力分布图图4 Y方向的应力分布图图5 XY Shear方向的应力分布图习题2、如图2所示的短圆筒，内半径为0.3m，外半径为0.5m，高度为1m。
假定圆筒内、外壁温度均为200℃，上端面温度为300℃,下端面绝热，导热系数为40w/mc°，计算圆筒的温度场分布。
解：（1）思路：建模和分析过程参考上机指南中的Project4。
（2）建模和分析：从已知条件知，此计算类型为Thermal力学类型；单元类型选择Thermal Solid Quad 4node 55；定义材料参数为KXX:7.5；生成几何模型利用点生成面方式；网格划分参照Project4；模型施加约束，分别给两条内外直边施加约束，Value: 200，上端面的约束设置为300;最后分析计算，查看应力图和变形图结果，保存数据。
R1=03R2=05图6 受温度载荷的圆筒示意图（3）ANSYS软件分析过程：2.1进入ANSYS程序→ANSYSED 10.0 →Interactive →change the working directory into yours →input Initial jobname: cylinder →Run2.2设置计算类型ANSYS Main Menu: Preferences… →select Thermal →OK2.3选择单元类型ANSYS Main Menu: Preprocessor →Element Type→Add/Edit/Delete →Add →select Thermal Solid Quad 4node 55 →OK (back to Element Types window)→Options… →select K3: Axisymmetric →OK →Close (the Element Type window)2.4定义材料参数ANSYS Main Menu: Preprocessor →Material Props →Material Models →Thermal→Conductivity→Isotropic→input KXX:7.5→OK2.5生成几何模型✓生成特征点ANSYS Main Menu: Preprocessor →Modeling →Create →Keypoints →In Active CS→依次输入四个点的坐标：input:1(0.3,0),2(0.5,0),3(0.5,1),4(0.3,1)→OK✓生成圆柱体截面ANSYS Main Menu: Preprocessor →Modeling →Create →Areas →Arbitrary →Through KPS→依次连接四个特征点，1(0.3,0),2(0.5,0),3(0.5,1),4(0.3,1) →OK2.6网格划分ANSYS Main Menu: Preprocessor →Meshing →Mesh Tool→(Size Controls) lines: Set →拾取两条水平边:OK→input NDIV: 5 →Apply→拾取两条竖直边:OK →input NDIV: 15 →OK →(back to the mesh tool window)Mesh: Areas, Shape: Quad, Mapped →Mesh →Pick All (in Picking Menu) →Close( the Mesh Tool window)2.7模型施加约束✓分别给两条直边施加约束ANSYS Main Menu: Solution→Define Loads →Apply→Thermal →Temperature →On Lines →拾取左边, Value: 200 →Apply(back to the window of apply temp on lines)→拾取右边，Value：200 →拾取顶边，Value:300→OK2.8 分析计算ANSYS Main Menu: Solution →Solve →Current LS→OK(to close the solve Current Load Step window) →OK2.9 结果显示ANSYS Main Menu: General Postproc →Plot Results→Deformed Shape…→select Def + Undeformed→OK (back to Plot Results window)→Contour Plot→Nodal Solu…→select: DOF solution, Temperature TEMP →OK(4)热应力图：图7 热应力DOF SOLUTION图图8 热应力的变化幅度图（Thermal gradient X方向）图9 热应力的变化幅度图（Thermal gradient Y方向）图10 热应力的变化总幅度图（Thermal gradient vector方向）习题3、矩形截面超静定粱的受力与约束情况如图13（a）所示，截面如图13（b）所示，b=20mm,h=80mm 材料的弹性模量为Mpa.2⨯=，泊松比为0.3。
