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

《ANSYS分析实例及其在工程中的应用-精华》由会员分享，可在线阅读，更多相关《ANSYS分析实例及其在工程中的应用-精华（182页珍藏版）》请在装配图网上搜索。

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长柱形天然气罐168编号：时间：2021年x月x日书山有路勤为径，学海无涯苦作舟页码：第182页 共182页1 第一章1.1 !1.1 铰接杆在外力作用下的变形! 本程序来源于邢静忠等编著的ANSYS分析实例与工程应用，机械工业出版社（2004年）FINISH/CLEAR,NOSTART! (1) 设置工程的选项,

8、定义几何参数/FileName, EX1.1/TITLE, EX1.1, DEFLECTION 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) 进入求解模块,开始求解/SOLUS

9、OLVEFINISH! (6) 进入后处理模块,显示节点位移和杆件内力/POST1PLDISP,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

10、,VALUE(1,2),DATA,STRSS,DISP*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.2 人字形

11、屋架的静力分析! 本程序来源于邢静忠等编著的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, LINK1MP

12、, 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, 8, FY,

13、 -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,1SOLVEFINISH

15、! (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,VALUE(1,3),DATA,ABS(REAC_1/900

16、) ,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 平行杆件与刚性梁连接的热应力问题! 本程序来源于邢静忠等编著的ANSYS分析实例与工程应用，机械工业出版社（2004年）FINISH/CLEA

17、R,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,-1N,6,1,-1FILL! (3) 定义单元信息E,1,4E,3,6MAT,2E,2,5CP,1,UY,5,4,

18、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 (STEEL_N)COPPER_E = ENEARN (COPPER_N)ETABLE,STRS_ST,LS,1ETA

19、BLE,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_COLABEL(1,2) = (MPa) , (MPa) *VFILL,VALUE(1,1),DATA,34.59,13.47*VFILL,VALUE(1,2),DATA,STRSS_ST/1e6,STRSS_CO/1e

20、6*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),VALUE(1,2),VALUE(1,3)(1X,A8,A8, ,F10.2, ,F10.2, ,1F5.3)/COM,-/OUTFINISH*LIST,EX1_4,out1.5 端部有间隙的杆

21、的热膨胀! 本程序来源于邢静忠等编著的ANSYS分析实例与工程应用，机械工业出版社（2004年）FINISH/CLEAR,START! (1) 定义分析类型,单元类型,截面面积和材料的力学,热学参数/FILNAME, EX1.5/PREP7/TITLE, EX1.5, THERMAL EXPANSION TO CLOSE A GAP AT A 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

22、,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,UX,2,1OUTPR,BASIC,LASTFINISH! (5) 在求解模块开始求解/SOLUSOLVEFINISH! (6) 在一般后处理模块中,定义水平应力和铅直应力单元表,并提取3号单元的应力结果/POST1ETABLE,STRSX,S,XETAB

23、LE,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, (MLABEL(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

24、.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,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

25、年）FINISH/CLEAR/FILNAME, EX2-1T/UNITS, CGS/TITLE, EX2-1T.txt, 3 NODES SIMULATION OF OVER-DETERMINED BEAM 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, 1NL

26、ISTDSYS, 1NLISTET, 1, BEAM3MP, EX, 1, 207e9MP, NUXY, 1, 0.3R, 1, 1, 0.020833, 0.5E, 1, 2E, 2, 3 $ E, 3, 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,

27、 1, PRES, 60, 120 $ SFBEAM, 3, 1, PRES, 120, 180SFBEAM, 4, 1, PRES, 180, 240 $ SFBEAM, 5, 1, PRES, 240, 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,IM

28、OMENT, SMISC, 6ETABLE,JMOMENT, SMISC, 12ETABLE, ISHEAR, SMISC, 2ETABLE, JSHEAR, SMISC, 8PRETAB/TITLE, SHEAR FORCE DIAGRAMPLLS, ISHEAR, JSHEAR/TITLE, BENDING MOMENT DIAGRAMPLLS, IMOMENT, JMOMENTFINISH2.2 !2.2 工字形截面外伸梁的平面弯曲! 本程序来源于邢静忠等编著的ANSYS分析实例与工程应用，机械工业出版社（2004年）Finish/CLEAR, NOSTART/PREP7/FileNam

29、e, EX2.2/TITLE, EX2.2, BEAM STRESSES AND DEFLECTIONSC* STR. OF MATL., TIMOSHENKO, PART 1, 3RD ED., PAGE 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

30、,UYSFBEAM,1,1,PRES,5000SFBEAM,4,1,PRES,5000FINISH! (5) 进入求解模块，定义求解选项并求解/SOLUOUTPR,BASIC,1SOLVEFINISH! (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,

31、LS,3*GET,STRSS,ELEM,MID_ELM,ETAB,STRS! (8) 将计算结果写入到文本文件*DIM,LABEL,CHAR,2*DIM,VALUE,2,3LABEL(1) = STRS_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 COM

32、PARISON -/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, EX2_2, out2.3 !2.3 矩形截面梁的纵横弯曲分析! 本程序来源于邢静忠等编著的ANSYS分析实例与工程应用，机械工业出版社（2004年）FINISH/CLEAR, NOSTART/PREP7/FileName, EX2.3/TITLE, EX2.3, TIE ROD WITH LA

33、TERAL LOADING, NO STREES STIFFENINGC* STR. OF MATLS., 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,

34、ALL,UY,ROTX,ROTZD,1,UZNSEL,S,5DSYM,SYMM,XNSEL,ALL! (5) 施加荷载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

35、,U,Z*GET,SLOPE_C2,NODE,LFT_END,ROT,YFINISH! (8) 在时间历程后处理器中，叠加处理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) 进入求解模块，制定最小收

36、敛误差后，开始非线性求解/SOLUCNVTOL,F,.0001,1SOLVEFINISH! (11) 进入后处理模块，显示位移计算结果，提取左边支座节点的转角和跨中挠度/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! (1

37、3) 生成计算结果文本文件EX21.out*DIM,LABEL,CHAR,3,2*DIM,VALUE_C1,3,3*DIM,VALUE_C2,3,3LABEL(1,1) = UZ MAX ,SLOPE ,MOMENT MLABEL(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),ABS(SLO

38、PE_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 | ANSY

39、S | RATIO/COM,/COM,RESULTS FOR F0 (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,2),VALUE_C2(1,3)(1X,A8,A8, ,F17.7, ,F17.7, ,

40、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单元计算悬臂梁的组合变形C* USING SOLID5 HEXAHEDRONSANTYPE,STATICET,1,SOLID5,2MP,EX,1,69E3MP,NUXY,

41、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,10ESIZE,25.4V,1,2,4,3,5,6,8,7VMESH,1! (4) 选择左侧截面处的所有节点，并固定这些节点。在右侧的所有节点上施加集中力NSE

42、L,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,FZ,4.4483/4*REPEAT,4,2SOLVEFINISH! (8) 进入一般后处理模块/POST1CSYS,0! (9) 定义宏程序，完成

43、自由端位移计算结果的处理*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,NODE,N1,U,X! (11) 处理第2个荷载步的计算结果*USE,MAC,2,(1/2.4712)SET,2,1LCSEL,S,0,0*GET,UY1

44、,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) = 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),

45、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 CANTILEVER BEAM by SOLID92 ELEMENT.smrt,offMOPT,VMESH,MAINMOPT,AMESH,ALTE! (15) 使用10节点四面体单元SOLID92ANTYPE,STATICET,1,SOLID92MP,EX,1,69E3MP,PRXY,1

46、,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,ALL,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

47、*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.4NSEL,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) 给上面选择

48、出的这些节点施加沿着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! (19) 对这种荷载情况进行计算/SOLUSOLVE! (20) 改变荷载，给上面选择出的这些节点施加沿着Y方向的轴向荷载，重新计算FDELE,ALLF,RGHT,FY, 4.4483/6F,LFT,FY, 4.4483/6F,TOPP,FY, 4.4483/6F,

49、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/6F,MIDD,FZ, 4.4483/3SOLVEFINISH! (22) 在后处理模块中，显示位移计算结果和理论结果的比值/POST1/COM * USE MACRO TO PROCESS ALL 3 LOADCASES *! (23) 利用前面定义的宏，处理前面计算的

50、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 L

51、ABEL(1,2) = in,in,in*VFILL,VALUE(1,1),DATA,7.6E-4,2.4712,10.9646*VFILL,VALUE(1,2),DATA,UX2,UY2,UZ2*VFILL,VALUE(1,3),DATA,ABS(UX2/(7.6E-4),ABS(UY2/2.4712),ABS(UZ2/10.9646)SAVE, TABLE_2FINISH! (25) 改变为10节点四面体耦合场单元SOLID98，重新计算悬臂梁在端部集中力作用下的双向弯曲/CLEAR, NOSTART/PREP7/TITLE, EX2.4(3), STRAIGHT CANTILEVER B

52、EAM by SOLID98 ELEMENT.smrt,offMOPT,VMESH,MAINMOPT,AMESH,ALTE/COM * REPEAT USING SOLID98 TETRAHEDRONS *RESUMEET,1,SOLID98,2! (26) 选择边界上的某些节点，施加约束和荷载/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,NDMX!NSEL,S,LOC,X,152.4NSEL,R,LOC,Y,5.08NSEL

53、,R,LOC,Z,1.27*GET,TOPP,NDMX!NSEL,S,LOC,X,152.4NSEL,R,LOC,Y,0.0NSEL,R,LOC,Z, 1.27*GET,BOTT,NDMX!NSEL,S,LOC,X, 152.4NSEL,R,LOC,Y, 2.54NSEL,R,LOC,Z, 2.54*GET,LFT,NDMX!NSEL,S,LOC,X,152.4NSEL,R,LOC,Y, 2.54NSEL,R,LOC,Z,.00*GET,RGHT,NDMXNSEL,ALL! (27) 给上面选择出的这些节点施加沿着X方向的轴向荷载/COM * APPLY LOADS TO PARAMETRIC

54、 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! (28) 改变荷载，给上面选择出的这些节点施加沿着Y方向的轴向荷载，重新计算/SOLUSOLVEFDELE,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/3! (29) 改变荷载，给上面选择出的这些节点施加沿着Z方向的轴

55、向荷载，重新计算SOLVEFDELE,ALLF,RGHT,FZ, 4.4483/6F,LFT,FZ, 4.4483/6F,TOPP,FZ, 4.4483/6F,BOTT,FZ, 4.4483/6F,MIDD,FZ, 4.4483/3SOLVEFINISH! (30) 在后处理模块中，显示位移计算结果和理论结果的比值/POST1/COM * USE MACRO TO PROCESS ALL 3 LOADCASES *! (31) 利用前面定义的宏，处理前面计算的3中荷载工况*USE,MAC,1,(1/7.6E-4)SET,1,1LCSEL,S,0,0N1=NODE(152.4,0,0)*GET,UX3,NODE,N1,U,X!*USE,MAC,2,(1/2.4712)SET,2,1LCSEL,S,0,0*GET,UY3,NODE,N1,U,Y!*USE,MAC,3,(1/.432)SET,3,1LCSEL,S,0,0*GET,UZ3,NODE,N1,U,Z! (32) 定义字符变量，比较理论结果和有限元结果的比值*DIM,LABEL,CHAR,3,