分析实例及其在工程中的应用精华

1、ADAMS应用实例-精华经典1第一章31.1!1.1 铰接杆在外力作用下的变形31.2!1.2 人字形屋架的静力分析41.3!1.3 超静定拉压杆的反力计算51.4!1.4 平行杆件与刚性梁连接的热应力问题71.5端部有间隙的杆的热膨胀82第二章102.1!2.1 等截面简单超静定梁的平面弯曲分析(单元多)102.2!2.2 工字形截面外伸梁的平面弯曲112.3!2.3 矩形截面梁的纵横弯曲分析132.4!2.4 悬臂梁的双向弯曲152.5!2.5 圆形截面悬臂杆的弯扭组合变形242.6!2.6 悬臂等强度梁的弯曲262.7!2.7 弹性地基半无限长梁在端部力和力偶作用下的变形292.8!2.

2、8 偏心受压杆的大变形分析303第三章313.1!3.1 利用梁单元计算压杆稳定性313.2!3.2 利用实体单元计算压杆稳定性333.3!3.3 悬臂压杆的过曲屈分析343.4!3.4 平面刚架的平面外失稳384第四章424.1!4.1 均布荷载作用下深梁的变形和应力(直接建模)424.2!4.1 均布荷载作用下深梁的变形和应力(实体建模)434.3!4.2 一对集中力作用下的圆环444.4!4.2 一对集中力作用下的圆环(实体建模)454.5!4.3 用实体单元分析变截面杆的拉伸464.6!4.4 用二维实体单元分析等截面悬臂梁的平面弯曲484.7!4.5 变截面悬臂梁在端部集中力作用下的

3、平面静力分析504.8!4.6 纯弯曲悬臂曲梁的二维静力分析524.9!4.7 端部集中力作用的悬臂圆环曲梁平面弯曲的三维分析544.10!4.8 均匀拉力作用下含圆孔板的孔边应力集中594.11!4.9 两端固定的厚壁管道在自重作用下的变形和应力655第五章675.1!5.1 含椭圆孔的椭圆薄膜在外部张力作用下的静力分析675.2!5.2 圆形薄膜大变形静力分析705.3!5.3 柱形容器在内压作用下的静力分析715.4!5.4 圆柱形薄壳在均匀内压作用下的静力分析736第六章746.1!6.1 简支和固支圆板的在不同荷载作用下的弯曲746.2!6.2 悬臂长板的大挠度弯曲776.3!6.3

4、 用壳体单元分析受均布荷载作用的固支圆板大挠度弯曲806.4!6.4 利用拉伸操作建立膨胀弯管模型826.5!6.5 两端简支开口柱壳在自重作用下的静力分析836.6!6.6 圆筒在一对横向集中力作用下的变形856.7!6.7 两边简支开口柱壳在集中力作用下的大变形曲屈877第七章907.1!7.1 单自由度弹簧质量系统的频率计算907.2!7.2 悬索自由振动的频率917.3!7.3 用弹簧单元连接的圆盘的扭转振动937.4!7.4 圆杆连接圆盘的扭转振动947.5!7.5 钻杆的扭转自由振动958第八章978.1!8.1 简支梁的自振频率计算978.2!8.1 简支梁的自振频率计算(增加节

5、点后的情况)998.3!8.2 自由-自由梁的纵向自由振动1018.4!8.3 有轴向压力作用的简支梁的自由振动1028.5!8.4 用壳体单元计算悬臂等强度梁的自由振动1048.6!8.5 矩形截面薄壁悬臂梁的自由振动1059第九章1089.1!9.1 圆形张紧薄膜的自由振动1089.2!9.2 薄膜二维非轴对称自由振动分析1099.3!9.3 薄膜三维非轴对称振动分析1129.4!9.4 悬臂长板的自由振动频率1159.5!9.5 悬臂宽板的模态分析1169.6!9.6 固支圆板的自由振动1179.7!9.7 用实体单元分析圆环的振动1199.8!9.8 机翼模型的振动分析12010第十章

6、12210.1!10.1 带三个圆孔的平面支座分析12210.2!10.2 角支座应力分析12310.3!10.3 立体斜支座的实体建模12410.4!10.4 四分之一车轮实体建模12510.5!10.5 轴承支座的实体建模12611第十一章12811.1!11.4 变截面悬臂梁的外形形状优化12811.2!11.5 平面刚架的优化设计13212第十二章13412.1!12.1 四边简支方形层合板在均布外载作用下的变形13412.2!12.2 均布拉力作用下含裂纹板的应力强度因子计算13913第十三章14713.1!第13章 用APDL实现空间网壳结构参数化建模14714第十四章14914.

7、1!第14章 塔式起重机静动力分析14915第十五章16515.1!15.3 利用轴对称壳单元SHELL51计算长柱形天然气罐16515.2!15.4 利用8节点2D实体单元PLANE82单元计算长柱形天然气罐16715.3!15.5 用20节点3D实体单元solid95计算1/4长柱形天然气罐1681 第一章1.1 !1.1 铰接杆在外力作用下的变形! 本程序来源于邢静忠等编著的ANSYS分析实例与工程应用，机械工业出版社（2004年）FINISH/CLEAR,NOSTART! (1) 设置工程的选项, 定义几何参数/FileName, EX1.1/TITLE, EX1.1, DEFLECT

8、ION OF A HINGED SUPPORT*AFUN,DEGTHETA=30A=10000B=A/2*TAN(THETA)! (2) 进入前处理模块, 定义节点/PREP7N,1N,2,A/2,-BN,3,A! (3) 在前处理模块中,定义单元类型,材料参数和各个单元ET,1,LINK1R,1,1000MP,EX,1,210E3E,1,2E,2,3! (4) 在前处理模块中,定义位移约束D,1,ALL,3,2F,2,FY,-1000OUTPR,1FINISH! (5) 进入求解模块,开始求解/SOLUSOLVEFINISH! (6) 进入后处理模块,显示节点位移和杆件内力/POST1PLD

9、ISP,2MID_NODE = NODE (A/2,-B,0 )*GET,DISP,NODE,MID_NODE,U,YLEFT_EL = ENEARN (MID_NODE)ETABLE,STRS,LS,1*GET,STRSS,ELEM,LEFT_EL,ETAB,STRS! (7) 申明数组,提取计算结果,并比较计算误差*DIM,LABEL,CHAR,2*DIM,VALUE,2,3LABEL(1) = 'STRS_MPa','DEF_mm'*VFILL,VALUE(1,1),DATA,1,-0.05498*VFILL,VALUE(1,2),DATA,STRSS,D

10、ISP*VFILL,VALUE(1,3),DATA,ABS(STRSS /1 ) ,ABS( DISP /0.05498 )/OUT,EX1_1,out/COM/COM,- EX1.1 RESULTS COMPARISON -/COM,/COM, | TARGET | ANSYS | RATIO/COM,*VWRITE,LABEL(1),VALUE(1,1),VALUE(1,2),VALUE(1,3)(1X,A8,' ',F10.3,' ',F10.3,' ',1F5.3)/COM,-/OUTFINISH*LIST,EX1_1,out1.2 !1

11、.2 人字形屋架的静力分析! 本程序来源于邢静忠等编著的ANSYS分析实例与工程应用，机械工业出版社（2004年）FINISH/CLEAR, START! (1) 设置工程的选项/FILNAME,EX1.2/UNITS, SI/TITLE, EX1.2, Plane Roof Tuss Model.! (2) 进入前处理模块, 定义节点/PREP7N, 1, 0, 0N, 2, 2, 0N, 3, 4, 0N, 4, 6, 0N, 5, 8, 0N, 6, 2, 1N, 7, 4, 2N, 8, 6, 1! (3) 在前处理模块中, 定义单元类型LINK1, 材料参数和各个单元ET, 1, L

12、INK1MP, EX, 1, 207E9R, 1, 0.01E, 1, 2E, 2, 3E, 3, 4E, 4, 5E, 1, 6E, 6, 7E, 2, 6E, 2, 7E, 3, 7E, 4, 7E, 4, 8E, 7, 8E, 8, 5/PNUM, NODE,1/NUMBER, 2/PNUM, ELEM, 1EPLOTFINISH! (4) 进入求解模块, 设置求解选项, 设置位移边界条件, 加载并开始求解/SOLUANTYPE, STATICOUTPR, BASIC, ALLD, 1, ALL, 0D, 5, UY, 0F, 6, FY, -1000F, 7, FY, -1000F,

13、8, FY, -1000SOLVEFINISH! (5) 进入一般后处理模块, 显示结构变形图/POST1PLDISP,2! (6) 列表显示节点位移和单元的计算结果PRDISPETABLE, MFORX,SMISC,1ETABLE, SAXL, LS, 1ETABLE, EPELAXL, LEPEL, 1PRETAB/NUMBER, 0PLETAB, MFORXFINISH1.3 !1.3 超静定拉压杆的反力计算! 本程序来源于邢静忠等编著的ANSYS分析实例与工程应用，机械工业出版社（2004年）FINISH/CLEAR,START! (1) 设置工程的选项, 进入前处理模块,声明单元类型

14、,实常数和材料参数/FILNAME, EX1.3/PREP7/TITLE, EX1.3, STATICALLY INDETERMINATE REACTION FORCE ANALYSISANTYPE,STATICET,1,LINK1R,1,1MP,EX,1,210E9! (2) 定义节点N,1N,2,0.4N,3,0.7N,4,1.0! (3) 定义单元E,1,2EGEN,3,1,1! (4) 定义位移约束和荷载D,1,ALL,4,3F,2,FY,-500F,3,FY,-1000FINISH! (5) 进入求解模块SOLU/SOLUOUTPR,BASIC,1OUTPR,NLOAD,1SOLVE

15、FINISH! (6) 进入一般后处理模块POST1/POST1NSEL,S,LOC,Y,1.0FSUM*GET,REAC_1,FSUM,ITEM,FYNSEL,S,LOC,Y,0FSUM*GET,REAC_2,FSUM,ITEM,FY! (7) 申明数组,输出计算结果并比较*DIM,LABEL,CHAR,2*DIM,VALUE,2,3LABEL(1) = 'R1, N ','R2, N '*VFILL,VALUE(1,1),DATA,900,600*VFILL,VALUE(1,2),DATA,ABS(REAC_1),ABS(REAC_2)*VFILL,VALU

16、E(1,3),DATA,ABS(REAC_1/900) ,ABS(REAC_2/600)/OUT,EX1_3,out/COM/COM,- EX1.3 RESULTS COMPARISON -/COM,/COM, | TARGET | ANSYS | RATIO/COM,*VWRITE,LABEL(1),VALUE(1,1),VALUE(1,2),VALUE(1,3)(1X,A8,' ',F10.1,' ',F10.1,' ',1F5.3)/COM,-/OUTFINISH*LIST,EX1_3,out1.4 !1.4 平行杆件与刚性梁连接的热应力问

17、题! 本程序来源于邢静忠等编著的ANSYS分析实例与工程应用，机械工业出版社（2004年）FINISH/CLEAR,START! (1) 定义分析类型,单元类型,截面面积和材料的力学,热学参数/FILNAME,EX1.4/PREP7/TITLE, EX1.4, THERMALLY LOADED SUPPORT STRUCTUREANTYPE,STATICET,1,LINK1R,1,65e-6MP,EX,1,100E9MP,ALPX,1,1.6E-5MP,EX,2,210E9MP,ALPX,2,1.3E-5TREF,70! (2) 定义节点位置信息N,1,-1N,3,1FILLN,4,-1,-1

18、N,6,1,-1FILL! (3) 定义单元信息E,1,4E,3,6MAT,2E,2,5CP,1,UY,5,4,6! (4) 定义位移约束和荷载（外力和温度荷载）D,1,ALL,3F,5,FY,-4000BFUNIF,TEMP,80FINISH! (5) 在求解模块,定义求解选项,并开始求解/SOLUOUTPR,BASIC,1OUTPR,NLOAD,1NSUBST,1SOLVEFINISH! (6) 在一般后处理模块中,寻找特定位置的节点和单元,并提取它们的内力/POST1STEEL_N = NODE (,)COPPER_N = NODE (1,0,0)STEEL_E = ENEARN (ST

19、EEL_N)COPPER_E = ENEARN (COPPER_N)ETABLE,STRS_ST,LS,1ETABLE,STRS_CO,LS,1*GET,STRSS_ST,ELEM,STEEL_E,ETAB,STRS_ST*GET,STRSS_CO,ELEM,COPPER_E,ETAB,STRS_CO! (7) 将计算结果及其与理论解的比较情况输出到文件vm3.txt*DIM,LABEL,CHAR,2,2*DIM,VALUE,2,3LABEL(1,1) = 'STRSS_ST','STRSS_CO'LABEL(1,2) = ' (MPa) ',&

20、#39; (MPa) '*VFILL,VALUE(1,1),DATA,34.59,13.47*VFILL,VALUE(1,2),DATA,STRSS_ST/1e6,STRSS_CO/1e6*VFILL,VALUE(1,3),DATA,ABS(STRSS_ST/34.59e6 ) ,ABS( STRSS_CO/13.47e6 )/COM/OUT,EX1_4,out/COM,- EX1.4 RESULTS COMPARISON -/COM,/COM, | TARGET | ANSYS | RATIO/COM,*VWRITE,LABEL(1,1),LABEL(1,2),VALUE(1,1),

21、VALUE(1,2),VALUE(1,3)(1X,A8,A8,' ',F10.2,' ',F10.2,' ',1F5.3)/COM,-/OUTFINISH*LIST,EX1_4,out1.5 端部有间隙的杆的热膨胀! 本程序来源于邢静忠等编著的ANSYS分析实例与工程应用，机械工业出版社（2004年）FINISH/CLEAR,START! (1) 定义分析类型,单元类型,截面面积和材料的力学,热学参数/FILNAME, EX1.5/PREP7/TITLE, EX1.5, THERMAL EXPANSION TO CLOSE A GAP AT A

22、RIGID SURFACEANTYPE,STATICET,1,CONTAC26R,1,2.625E15ET,2,PLANE42,3R,2,1MP,EX,1,72.4E3MP,ALPX,1,2.25E-5MP,NUXY,1,0! (2) 定义节点信息N,1,2000,1000N,2,3000,1000N,3,3000,4000N,4,2000,4000N,11,1000,998N,12,4000,998! (3) 定义单元信息E,1,11,12E,2,11,12TYPE,2REAL,2E,1,2,3,4TREF,21BFUNIF,TEMP,77! (4) 定义位移约束D,3,ALL,4D,1,U

23、X,2,1OUTPR,BASIC,LASTFINISH! (5) 在求解模块开始求解/SOLUSOLVEFINISH! (6) 在一般后处理模块中,定义水平应力和铅直应力单元表,并提取3号单元的应力结果/POST1ETABLE,STRSX,S,XETABLE,STRSY,S,Y*GET,STRSSX,ELEM,3,ETAB,STRSX*GET,STRSSY,ELEM,3,ETAB,STRSY! (7) 将计算结果输出到结果文件*DIM,LABEL,CHAR,2,2*DIM,VALUE,2,3LABEL(1,1) = 'SIGX, (M','SIGY, (M'LA

24、BEL(1,2) = 'Pa) ','Pa) '*VFILL,VALUE(1,1),DATA,-91.22,-42.96*VFILL,VALUE(1,2),DATA,STRSSX,STRSSY*VFILL,VALUE(1,3),DATA,ABS(STRSSX/91.22) ,ABS(STRSSY/42.96)/COM/OUT,EX1_5, out/COM,- EX1.5 RESULTS COMPARISON -/COM,/COM, | TARGET | ANSYS | RATIO/COM,*VWRITE,LABEL(1,1),LABEL(1,2),VALUE(1

25、,1),VALUE(1,2),VALUE(1,3)(1X,A8,A8,' ',F10.3,' ',F10.3,' ',1F5.3)/COM,-/OUTFINISH*LIST,EX1_5,out2 第二章2.1 !2.1 等截面简单超静定梁的平面弯曲分析(单元多)! 本程序来源于邢静忠等编著的ANSYS分析实例与工程应用，机械工业出版社（2004年）FINISH/CLEAR/FILNAME, EX2-1T/UNITS, CGS/TITLE, EX2-1T.txt, 3 NODES SIMULATION OF OVER-DETERMINED BEAM

26、 WITH MORE NODES/PREP7N, 1, 0, 0N, 2, 0.06, 0$ N, 3, 0.12, 0 $ N, 4, 0.18, 0N, 5, 0.24, 0 $ N, 6, 0.30, 0 $ N, 7, 0.36, 0N, 8, 0.42, 0 $ N, 9, 0.48, 0 $ N, 10, 0.54, 0N, 11, 0.60, 0NPLOT NPLOT, 1NLISTDSYS, 1NLISTET, 1, BEAM3MP, EX, 1, 207e9MP, NUXY, 1, 0.3R, 1, 1, 0.020833, 0.5E, 1, 2E, 2, 3 $ E, 3,

27、 4 $ E, 4, 5 $ E, 5, 6E, 6, 7 $ E, 7, 8 $ E, 8, 9 $ E, 9,10E, 10, 11EPLOTELISTFINISH/SOLUANTYPE, STATICOUTPR, BASIC, ALLD, 1, UX, 0, , , , , UY, ROTZD, 11, UY, 0 SFBEAM, 1, 1, PRES, 0, 60SFBEAM, 2, 1, PRES, 60, 120 $ SFBEAM, 3, 1, PRES, 120, 180SFBEAM, 4, 1, PRES, 180, 240 $ SFBEAM, 5, 1, PRES, 240,

28、 300SFBEAM, 6, 1, PRES, 300, 360 $ SFBEAM, 7, 1, PRES, 360, 420SFBEAM, 8, 1, PRES, 420, 480 $ SFBEAM, 9, 1, PRES, 480, 540SFBEAM, 10, 1, PRES, 540, 600F, 6, FY, -200SOLVEFINISH/POST1PRDISPETABLE,IMOMENT, SMISC, 6ETABLE,JMOMENT, SMISC, 12ETABLE, ISHEAR, SMISC, 2ETABLE, JSHEAR, SMISC, 8PRETAB/TITLE, S

29、HEAR FORCE DIAGRAMPLLS, ISHEAR, JSHEAR/TITLE, BENDING MOMENT DIAGRAMPLLS, IMOMENT, JMOMENTFINISH2.2 !2.2 工字形截面外伸梁的平面弯曲! 本程序来源于邢静忠等编著的ANSYS分析实例与工程应用，机械工业出版社（2004年）Finish/CLEAR, NOSTART/PREP7/FileName, EX2.2/TITLE, EX2.2, BEAM STRESSES AND DEFLECTIONSC* STR. OF MATL., TIMOSHENKO, PART 1, 3RD ED., PAGE

30、 98, PROB. 4! (1) 定义求解类型，实常数和材料参数ANTYPE,STATICET,1,BEAM3KEYOPT,1,9,9R,1,35e-4,4.5E-5,0.20MP,EX,1,210E9MP,PRXY,1,0.3! (2) 定义节点位置N,1N,5,4FILL! (3) 定义单元连接关系E,1,2EGEN,4,1,1! (4) 定义位移约束和荷载条件D,2,UX,UYD,4,UYSFBEAM,1,1,PRES,5000SFBEAM,4,1,PRES,5000FINISH! (5) 进入求解模块，定义求解选项并求解/SOLUOUTPR,BASIC,1SOLVEFINISH! (

31、6) 进入后处理模块列表显示节点位移计算结果和图形显示变形图/POST1SET,1,1PRNSOL,U,COMPPRNSOL,ROT,COMPPLDISP,1! (7) 提取指定位置的节点和单元，定义单元表并获取弯曲应力MID_NODE = NODE (2, )*GET,DISP,NODE,MID_NODE,U,YMID_ELM = ENEARN (MID_NODE)ETABLE,STRS,LS,3*GET,STRSS,ELEM,MID_ELM,ETAB,STRS! (8) 将计算结果写入到文本文件*DIM,LABEL,CHAR,2*DIM,VALUE,2,3LABEL(1) = 'S

32、TRS_MPa','DEF_mm'*VFILL,VALUE(1,1),DATA,-5.5556,0.1323*VFILL,VALUE(1,2),DATA,STRSS/1e6,DISP*1000*VFILL,VALUE(1,3),DATA,ABS(STRSS/1e6/5.5556) ,ABS(DISP/0.1323e-3)/OUT, EX2_2, out/COM/COM,- EX2.2 RESULTS COMPARISON -/COM,/COM, | TARGET | ANSYS | RATIO/COM,*VWRITE,LABEL(1),VALUE(1,1),VALUE(

33、1,2),VALUE(1,3)(1X,A8,' ',F10.3,' ',F10.3,' ',1F5.3)/COM,-/OUTFINISH*LIST, EX2_2, out2.3 !2.3 矩形截面梁的纵横弯曲分析! 本程序来源于邢静忠等编著的ANSYS分析实例与工程应用，机械工业出版社（2004年）FINISH/CLEAR, NOSTART/PREP7/FileName, EX2.3/TITLE, EX2.3, TIE ROD WITH LATERAL LOADING, NO STREES STIFFENINGC* STR. OF MATLS.,

34、 TIMOSHENKO, PART 2, 3RD ED., PAGE 42, ART. 6! (1) 定义求解类型，实常数和材料参数ANTYPE,STATICET,1,BEAM4,1R,1,0.403225e-2, 0.135492e-5, 0.135492e-5, 0.0635, 0.0635MP,EX,1,210E9MP,PRXY,0.3! (2) 定义节点位置N,1N,5,2.54FILL! (3) 定义单元连接关系E,1,2EGEN,4,1,1! (4) 定义位移约束和荷载条件D,ALL,UY,ROTX,ROTZD,1,UZNSEL,S,5DSYM,SYMM,XNSEL,ALL! (5

35、) 施加荷载F,1,FX,-97740SFBEAM,ALL,1,PRES,314FINISH! (6) 进入求解模块，开始求解/SOLUSOLVEFINISH! (7) 进入后处理模块，显示位移计算结果，提取左边支座节点的转角和跨中挠度/POST1NSEL,S,1,5,4PRNSOL,U,ZPRNSOL,ROT,YNSEL,ALLPRRSOLRGHT_END = NODE (2.54,0,0)LFT_END = NODE (0,0,0)*GET,UZ_MX_C2,NODE,RGHT_END,U,Z*GET,SLOPE_C2,NODE,LFT_END,ROT,YFINISH! (8) 在时间历程

36、后处理器中，叠加处理5号节点的弯矩/POST26RFORCE,2,RGHT_END,M,YSTORE*GET,M_MX_C2,VARI,2,EXTREM,VMAXFINISH! (9) 重新进入前处理模块，重新计算考虑应力强化效应的压杆/PREP7/TITLE, EX2.3(2), TIE ROD WITH LATERAL LOADING, STRESS STIFFENING PRESENTSSTIF,ONNSUBST,5AUTOTS,ONFINISH! (10) 进入求解模块，制定最小收敛误差后，开始非线性求解/SOLUCNVTOL,F,.0001,1SOLVEFINISH! (11) 进入

37、后处理模块，显示位移计算结果，提取左边支座节点的转角和跨中挠度/POST1NSEL,S,1,5,4PRNSOL,U,ZPRNSOL,ROT,YPRRSOL*GET,UZ_MX_C1,NODE,RGHT_END,U,Z*GET,SLOPE_C1,NODE,LFT_END,ROT,YFINISH! (12) 在时间历程后处理器中，叠加处理跨中（5号节点的）弯矩/POST26RFORCE,2,RGHT_END,M,YSTORE*GET,M_MX_C1,VARI,2,EXTREM,VMAX! (13) 生成计算结果文本文件EX21.out*DIM,LABEL,CHAR,3,2*DIM,VALUE_C1

38、,3,3*DIM,VALUE_C2,3,3LABEL(1,1) = 'UZ MAX ','SLOPE ','MOMENT M'LABEL(1,2) = '(m) ','(rad) ','AX N-m'*VFILL,VALUE_C1(1,1),DATA,-0.0050274,.00321,-521.52*VFILL,VALUE_C1(1,2),DATA,UZ_MX_C1,SLOPE_C1,M_MX_C1*VFILL,VALUE_C1(1,3),DATA,ABS(UZ_MX_C1/0.0050274),

39、ABS(SLOPE_C1/.00321), ABS(M_MX_C1/521.52)*VFILL,VALUE_C2(1,1),DATA,-0.00957,.006028,-1013*VFILL,VALUE_C2(1,2),DATA,UZ_MX_C2,SLOPE_C2,M_MX_C2*VFILL,VALUE_C2(1,3),DATA,ABS(UZ_MX_C2/0.00957),ABS(SLOPE_C2/.006028), ABS(M_MX_C2/1013)/COM/OUT, EX2_3, out/COM,- EX2.3 RESULTS COMPARISON -/COM,/COM, | TARGET

40、 | ANSYS | RATIO/COM,/COM,RESULTS FOR F<>0 (STIFFENED):/COM,*VWRITE,LABEL(1,1),LABEL(1,2),VALUE_C1(1,1),VALUE_C1(1,2),VALUE_C1(1,3)(1X,A8,A8,' ',F17.7,' ',F17.7,' ',1F5.3)/COM,/COM,RESULTS FOR F=0 (UNSTIFFENED):/COM,*VWRITE,LABEL(1,1),LABEL(1,2),VALUE_C2(1,1),VALUE_C2(1

41、,2),VALUE_C2(1,3)(1X,A8,A8,' ',F17.7,' ',F17.7,' ',1F5.3)/COM,-/OUTFINISH*LIST, EX2_3, out2.4 !2.4 悬臂梁的双向弯曲! 本程序来源于邢静忠等编著的ANSYS分析实例与工程应用，机械工业出版社（2004年）Finish/CLEAR, START/FileName,EX2.4/PREP7/TITLE, EX2.4(1), STRAIGHT CANTILEVER BEAM BY SOLID5 ELEMENT.! (1) 利用SOLID5单元计算悬臂梁的组合

42、变形C* USING SOLID5 HEXAHEDRONSANTYPE,STATICET,1,SOLID5,2MP,EX,1,69E3MP,NUXY,1,0.3! (2) 定义悬臂梁上的8个的关键点K,1K,2,152.4! 用关键点生成命令KGEN, ITIME, NP1, NP2, NINC, DX, DY, DZ, KINC 生成其它关键点KGEN,2,1,2,1,5.08KGEN,2,1,4,1,2.54! 上面两行命令完成了节点3到4的定义和节点5到8的位置定义! (3) 定义一条边线，制定这条线划分单元的个数和一般单元的尺寸L,1,2smrt,offLESIZE,ALL,10ESI

43、ZE,25.4V,1,2,4,3,5,6,8,7VMESH,1! (4) 选择左侧截面处的所有节点，并固定这些节点。在右侧的所有节点上施加集中力NSEL,S,LOC,X,0D,ALL,ALLNSEL,ALLFK,2,FX,4.4483/4*REPEAT,4,2NOORDERFINISH! (5) 求解在自由端拉力作用下悬臂梁的变形，主要是拉压变形/SOLUSOLVE! (6) 改变荷载，计算在Y方向的集中力作用下的情形FKDELE,ALL,FXFK,2,FY,4.4483/4*REPEAT,4,2SOLVE! (7) 改变荷载，计算在Z方向的集中力作用下的情形FKDELE,ALL,FYFK,2

44、,FZ,4.4483/4*REPEAT,4,2SOLVEFINISH! (8) 进入一般后处理模块/POST1CSYS,0! (9) 定义宏程序，完成自由端位移计算结果的处理*CREATE,MACSET,ARG1,1LCDEF,ARG1,ARG1NSEL,S,LOC,X,152.4PRNSOL,U,COMPLCFACT,ARG1,ARG2LCASE,ARG1PRNSOL,U,COMP*END! (10) 使用上面定义的宏程序处理所有的3个工况/COM *USE,MAC,1,(1/7.6E-4)SET,1,1LCSEL,S,0,0N1=NODE(152.4 , 0, 0)*GET,UX1,NOD

45、E,N1,U,X! (11) 处理第2个荷载步的计算结果*USE,MAC,2,(1/2.4712)SET,2,1LCSEL,S,0,0*GET,UY1,NODE,N1,U,Y! (12) 处理第3个荷载步的计算结果*USE,MAC,3,(1/10.9646)SET,3,1LCSEL,S,0,0*GET,UZ1,NODE,N1,U,Z! (13) 定义字符串变量和相应的理论计算结果*DIM,LABEL,CHAR,3,2*DIM,VALUE,3,3LABEL(1,1) = 'DEFL X ','DEFL Y ','DEFL Z 'LABEL(1,2)

46、 = 'mm','mm','mm'*VFILL,VALUE(1,1),DATA,7.6E-4,2.471,10.965*VFILL,VALUE(1,2),DATA,UX1,UY1,UZ1*VFILL,VALUE(1,3),DATA,ABS(UX1/(7.6E-4),ABS(UY1/2.471),ABS(UZ1/10.965)SAVE, TABLE_1FINISH! (14) 改变单元类型，重新分析矩形截面悬臂梁在端部集中力作用下的变形/CLEAR, NOSTART/PREP7/TITLE, EX2.4(2), STRAIGHT CANTILEVE

47、R BEAM by SOLID92 ELEMENT.smrt,offMOPT,VMESH,MAINMOPT,AMESH,ALTE! (15) 使用10节点四面体单元SOLID92ANTYPE,STATICET,1,SOLID92MP,EX,1,69E3MP,PRXY,1,0.3K,1K,2,5.08KGEN,2,1,2,1,2.54KGEN,2,1,4,1,152.4/9ESIZE,25.4V,1,2,4,3,5,6,8,7VMESH,1! (16) 利用第1块体的单元模式，生成其余的8部分单元，定义固定端部的约束VGEN,9,1,1,1,152.4/9,4NSEL,S,LOC,X,0D,AL

48、L,ALLSAVE! (17) 获取相应节点的编号，并在这些节点上施加荷载/COM * GET NODE NUMBERS FOR LOAD APPLICATION *NSEL,S,LOC,X,152.4NSEL,R,LOC,Y,2.54NSEL,R,LOC,Z,1.27*GET,MIDD,NDMXNSEL,S,LOC,X,152.4NSEL,R,LOC,Y,5.08NSEL,R,LOC,Z,1.27*GET,TOPP,NDMXNSEL,S,LOC,X,152.4NSEL,R,LOC,Y,0.0NSEL,R,LOC,Z, 1.27*GET,BOTT,NDMXNSEL,S,LOC,X, 152.4

49、NSEL,R,LOC,Y, 2.54NSEL,R,LOC,Z, 2.54*GET,LFT,NDMXNSEL,S,LOC,X,152.4NSEL,R,LOC,Y, 2.54NSEL,R,LOC,Z,.00*GET,RGHT,NDMXNSEL,ALL! (18) 给上面选择出的这些节点施加沿着X方向的轴向荷载/COM * APPLY LOADS TO PARAMETRIC NODE NUMBERS *F,RGHT,FX, 4.4483/6F,LFT,FX, 4.4483/6F,TOPP,FX, 4.4483/6F,BOTT,FX, 4.4483/6F,MIDD,FX, 4.4483/3FINISH

50、! (19) 对这种荷载情况进行计算/SOLUSOLVE! (20) 改变荷载，给上面选择出的这些节点施加沿着Y方向的轴向荷载，重新计算FDELE,ALLF,RGHT,FY, 4.4483/6F,LFT,FY, 4.4483/6F,TOPP,FY, 4.4483/6F,BOTT,FY, 4.4483/6F,MIDD,FY, 4.4483/3SOLVE! (21) 改变荷载，给上面选择出的这些节点施加沿着Z方向的轴向荷载，重新计算FDELE,ALLF,RGHT,FZ, 4.4483/6F,LFT,FZ, 4.4483/6F,TOPP,FZ, 4.4483/6F,BOTT,FZ, 4.4483/6

51、F,MIDD,FZ, 4.4483/3SOLVEFINISH! (22) 在后处理模块中，显示位移计算结果和理论结果的比值/POST1/COM * USE MACRO TO PROCESS ALL 3 LOADCASES *! (23) 利用前面定义的宏，处理前面计算的3中荷载工况*USE,MAC,1,(1/7.6E-4)SET,1,1LCSEL,S,0,0N1=NODE(152.4,0,0)*GET,UX2,NODE,N1,U,X!*USE,MAC,2,(1/2.4712)SET,2,1LCSEL,S,0,0*GET,UY2,NODE,N1,U,Y!*USE,MAC,3,(1/.432)SET,3,1LCSEL,S,0,0*GET,UZ2,NODE,N1,U,Z! (24) 定义字符变量，比较理论结果和有限元结果的比值*DIM,LABEL,CHAR,3,2*DIM,VALUE,3,3LABEL(1,1) = 'DEFL X ','DEFL Y ','DEFL Z 'LABEL(1,2) = 'in','in','in'*VFILL,VAL