命令流ansys经典实例

合集下载

ansysy命令流例子

题二：按图1-1所示，将截面沿路径扫掠成体。

单位（m）路径截面!画面沿线拖拉成体finiclear.all/prep7k,1,0,0,0k,2,0,-6,0k,3,3,-9,0k,4,5,-9,0k,5,3,-6,0LSTR,1,2 !连接1和2点生成直线1LARC,2,3,5,3 !起始点2结束点3 控制点5半径3画弧线LSTR,3,4 !连接3和4点K,6,-0.4,0,0K,7,0.4,0,0K,8,0,0,-0.4K,9,0,0,0.692820323LSTR,9,6LSTR,9,7LARC,6,7,8,0.4AL,4,5,6!由线生成面VDRAG,1,,, , , ,1, 2, 3, , , !将面1沿着线1 2 3拖拉成体题三：给L=100m的梁定义横截面，横截面如图2-1：（单位m）finiclear.all/prep7ET,1,PLANE82k,1,0,0,0k,2,2,0,0k,3,1,-1,0k,4,1,1.732050808,0LARC,1,2,3,1 !起始点1和终点2 和关键点3生成半径为1的弧线LSTR,1,4!连接点1 4生成线LSTR,2,4AL,1,2,3!将线1 2 3成面LSEL,ALLLESIZE, all, , , 6AMESH,ALLSECWRITE,jiemian1,SECT!将截面保存在D中的USER 中的tangyufengstarfini/CLEAR/prep7ET,1,beam188k,1,0,0,0k,2,100,0,0LSTR,1,2SECTYPE,1,BEAM,MESHSECOFFSET,CENTSECREAD,'jiemian1','SECT','',MESH!将截面读入ESIZE,, 300LMESH,1!在将PlotCTRls中的STYLE中的size and sharp中的Disply of element打开题目四：实体建模，并采用映射网格生成单元。

Ansys命令流大全爆款.doc

1、A,P1,P2,P3,P4,P5,P6,P7,P8,P9此命令用已知的一组关键点点（P1~P9）来定义面（Area），最少使用三个点才能围成面，同时产生转围绕些面的线。

点要依次序输入，输入的顺序会决定面的法线方向。

如果超过四个点，则这些点必须在同一个平面上。

Menu Paths:MainMenu>Preprocessor>Create>Arbitrary>Through KPs2、*ABBR，Abbr,String－－定义一个缩略语．Abbr:用来表示字符串＂String＂的缩略语，长度不超过８个字符．String：将由＂Abbr＂表示的字符串，长度不超过６０个字符．3、ABBRES，Lab，Fname，Ext－从一个编码文件中读出缩略语．Lab：指定读操作的标题，NEW：用这些读出的缩略语重新取代当前的缩略语（默认）CHANGE：将读出的缩略语添加到当前缩略语阵列，并替代现存同名的缩略语．Ext：如果＂Fname＂是空的，则缺省的扩展命是＂ABBR＂．4、ABBSAV，Lab，Fname，Ext－将当前的缩略语写入一个文本文件里Lab：指定写操作的标题，若为ALL，表示将所有的缩略语都写入文件（默认）5、add, ir, ia,ib,ic,name,--,--,facta, factb, factc将ia,ib,ic变量相加赋给ir变量ir, ia,ib,ic：变量号name: 变量的名称6、Adele,na1,na2,ninc,kswp ！kswp=0时只删除掉面积本身，＝1时低单元点一并删除。

7、Adrag, nl1,nl2,nl3,nl4,nl5,nl6,nlp1,nlp2,nlp3,nlp4,nlp5,nlp6 ！面积的建立，沿某组线段路径，拉伸而成。

8、Afillt,na1,na2,rad ！建立圆角面积，在两相交平面间产生曲面，rad为半径。

9、*AFUN，Lab在参数表达式中，为角度函数指定单位．Lab：指定将要使用的角度单位．有３个选项．RAD：在角度函数的输入与输出中使用弧度单位（默认）DEG：在角度函数的输入与输出中使用度单位．STAT：显示该命令当前的设置（即是度还是弧度）．10、Agen,itime,na1,na2,ninc,dx,dy,dz,kinc,noelem,imove ！面积复制命令。

ANSYS土木工程经典实例命令流大全

ANSYS土木工程经典实例命令流大全ANSYS是目前最为领先的工程仿真软件之一，广泛应用于土木工程领域。

本文将介绍一些ANSYS土木工程的经典实例以及相关的命令流，帮助工程师更好地应用该软件进行仿真分析。

1. 桥梁结构分析实例实例简介一座桥梁由多个零部件组成，包括桥墩、桥面、桥拱等。

如何分析这些零部件的受力情况，以便于对桥梁结构进行优化和改进呢？ANSYS提供了一系列的分析工具和命令流，可以帮助我们完成这项任务。

命令流详解首先需要创建一个桥梁模型，并进行网格划分。

然后利用ANSYS的各种分析工具进行仿真分析，得到桥梁各个零部件的受力情况。

在此基础上，可以进行结构优化，最终得到一个强度和稳定性都较好的桥梁结构。

以下是桥梁结构分析实例的一些关键命令流：•创建单元网格：ET, SOLID186•定义材料属性：MP, EX, NU, DENS•定义边界条件：*BOUNDARY，MP，SYM，FIX•加载边界条件：DLOAD，TYPE，P，_LOC•计算位移和应力分布：*POST1，DISPL，NF，S2. 地基基础分析实例实例简介地基基础是土木工程中的重要组成部分，承载着整个工程的重量。

如何对地基基础的承载力进行分析和计算呢？ANSYS也提供了相应的分析工具和命令流，帮助土木工程师完成这项任务。

命令流详解首先需要建立地基基础的三维模型，并进行网格划分。

然后利用ANSYS的各种分析工具进行仿真分析，计算地基基础承载力、变形等相关指标。

在此基础上，可以进行结构优化，最终得到一个承载能力和稳定性都较好的地基基础。

以下是地基基础分析实例的一些关键命令流：•创建单元网格：ET, SOLID186•定义材料属性：MP, EX, NU, DENS•定义边界条件：*BOUNDARY，MP，SYM，FIX•加载边界条件：DLOAD，TYPE，P，_LOC•计算应力分布和变形：*POST1，S，EPTO，ETA3. 挖土工程分析实例实例简介挖土工程是土木工程中的重要环节，需要对不同参数下的挖土工程进行分析和优化。

ansys中典型案例命令流汇总

lfillt,649,650,0.008075 lfillt,650,651,0.008075
FLST,2,6,4 FITEM,2,652 FITEM,2,651 FITEM,2,654 FITEM,2,650 FITEM,2,653 FITEM,2,649 AL,P51X
vext,271,271,0,0,-0.0594,0 vovlap,all
lstr,17,18 lstr,18,19 lstr,19,20
lfillt,25,26,0.0022 lfillt,26,27,0.0022
FLST,2,5,4,ORDE,2 FITEM,2,25
FITEM,2,-29 LCOMB,P51X, ,0
K,21,0,0.00225,-0.012075 k,22,0.020075,0.00225,-0.012075 k,23,0.020075,0.00225,0.012075 K,24,0,0.00225,0.012075
!---------------------------------------------------------
K,,0,0.00225,-0.015675 !生成次外层13匝线圈 k,,0.023675,0.00225,-0.015675 k,,0.023675,0.00225,0.015675 K,,0,0.00225,0.015675
lstr,129,130 lstr,130,131 lstr,131,132
lfillt,186,187,0.0022 lfillt,187,188,0.0022
lstr,125,129 lstr,132,128
FLST,2,12,4 FITEM,2,191 FITEM,2,181 FITEM,2,182 FITEM,2,184 FITEM,2,185 FITEM,2,183 FITEM,2,192 FIT

ANSYS单元特性之命令流算例

EX1.1 （LINK1‎‎）（1）进入后处理‎模块,显示节点位‎移和杆件内‎力MID_N‎O DE = NODE (A/2,-B,0 )! 寻找距离位‎置（A/2,-B,0）最近的点，存入MID‎_NODE‎*GET,DISP,NODE,MID_N‎O DE,U,Y！提取节点M‎I D_NO‎D E上的位‎移UY，若果已知要‎求的节点，直接提取即‎可。

LEFT_‎E L = ENEAR‎N (MID_N‎O DE)! 需找距离节‎点MID_‎N ODE最‎近的单元，存入LEF‎T_ELETABL‎E,STRS,L S,1! 用轴向应力‎S AXL的‎编号“LS,1”定义单元表‎S TRS*GET,STRSS‎,ELEM,LEFT_‎E L,ETAB,STRS! 从单元表S‎‎TRS中提‎取LEFT‎‎_EL单元‎的应力结果‎，存入变量S‎‎TRSS。

注意：提取的轴向‎应力结果具‎体到指定的‎单元。

（2）申明数组,提取计算结‎果,并比较计算‎误差*DIM,LABEL‎,CHAR,2！定义2个元‎素的字符型‎数组LAB‎‎E L*DIM,VALUE‎,,2,3！定义2*3的数值型‎数组V AL‎‎U ELABEL‎(1) = 'STRS_‎M Pa','DEF_m‎m' ! 给字符型数‎组的第1个‎元素赋值*VFILL‎,VALUE‎(1,1),DATA,1,-0.05498‎! 给其他数值‎型数组中的‎元素赋值*VFILL‎,VALUE‎(1,2),DATA,STRSS‎,DISP*VFILL‎,VALUE‎(1,3),DATA,ABS(STRSS‎/1 ) ,ABS( DISP /0.05498‎)/OUT,EX1_1‎,out ！将输出内容‎重定向到文‎件EX1_‎‎1.out/COM ! 以注释形式‎输出内容/COM,------------------- EX1.1 RESUL‎T S COMPA‎R ISON‎---------------------/COM,/COM, | TARGE‎T | ANSYS‎ | RATIO‎/COM,*VWRIT‎E,LABEL‎(1),VALUE‎(1,1),VALUE‎(1,2),VALUE‎(1,3)(1X,A8,' ',F10.3,' ',F10.3,' ',1F5.3)/COM,----------------------------------------------------------------/OUT ! 结束数据重‎定向，关闭输出文‎件FINIS‎H*LIST,EX1_1‎,out ! 列表显示文‎件EX1_‎‎1.out的内‎容EX1.2 （LINK1‎‎）/PNUM, NODE,1！打开节点编‎号显示/NUMBE‎R, 2！只显示编号‎，不使用色彩‎列表显示节‎点位移和单‎元的计算结‎果PRDIS‎P! 列表显示节‎点位移值计‎算结果ETABL‎E, MFORX‎,SMISC‎,1！以杆单元的‎轴力为内容‎，建立单元表‎M F ORX‎‎ETABL‎E, SAXL, LS, 1 ！以杆的轴向‎应力为内容‎，建立单元表‎S A XLETABL‎E, EPELA‎X L, LEPEL‎, 1! 以杆单元的‎轴向应变为‎内容，建立单元表‎E PELA‎‎X LPRETA‎B! 显示单元表‎中的计算结‎果/NUMBE‎R, 0！显示编号，并使用彩色‎PLETA‎B, MFORX‎！用色度图显‎示单元表M‎‎F ORX中‎杆件轴力图‎EX1.3 （LINK1‎‎）NSEL,S,LOC,Y,1.0 ！选择所有位‎于Y=1.0位置上的‎节点FSUM！累计叠加选‎择集中所有‎节点上的反‎力*GET,REAC_‎1,FSUM,,ITEM,FY ！将累加结果‎中的FY（Y方向的力‎）保存到变量‎R E AC_‎‎1中EX1.4 （LINK1‎‎）R,1,65e-6！定义第1类‎实常数，杆件截面面‎积为65m‎‎m^2，在转化为国‎际单位制时‎操作TREF,70 ! 设定参考温‎度为70度‎BFUNI‎F,TEMP,80 ! 温度从原来‎的70度均‎匀上升到8‎‎0度（TREF+10）EX1.5 （PLANE‎‎42 AND CONT A‎‎C26）ETABL‎E,STRSX‎,S,X！定义X方向‎的应力为单‎元表S TR‎‎S X*GET,STRSS‎X,ELEM,3,ETAB,STRSX‎！从单元表S‎‎TRS X中‎提取3号单‎元的X向应‎力，存入STR‎‎S S X。

ansys旋转经典命令流

ansys旋转经典命令流1 旋转摩擦 (1)2. 电磁三d命令流实例（论坛看到） (11)3. 帮助感应加热例⼦induction heating of a solid cylinder billet (15)4. 感应加热温度场的数值模拟（论⽂）inducheat30命令流 (19)5. 如何施加恒定的⾓速度？Simwe仿真论坛 (24)6. 旋转⼀个已经⽣成好的物体 (27)7. 产⽣这样的磁⼒线 (28)8. 旋转摩擦⽣热简单例⼦（⼆维旋转） (32)8.1. 原版 (32)8.2. 部分gui操作 (35)9. VM229 Input Listing (39)10 轴承---耦合+接触分析 (47)11. 板的冲压仿真 (52)1 旋转摩擦FINISH/FILNAME,Exercise24 !定义隐式热分析⽂件名/PREP7 !进⼊前处理器ET,1,SOLID5 !选择单元类型MPTEMP,,,,,,,,MPTEMP,1,0MPDATA,DENS,1,,7800 !定义材料1的密度MPTEMP,,,,,,,,MPTEMP,1,0MPDATA,C,1,,460 !定义材料1的⽐热MPTEMP,,,,,,,,MPTEMP,1,0MPDATA,KXX,1,,66.6 !定义材料1的热传导系数MPTEMP,,,,,,,,MPTEMP,1,0UIMP,1,REFT,,,30 !定义材料1的热膨胀系数的参考温度MPDATA,ALPX,1,,1.06e-5 !定义材料1的热膨胀系数MPTEMP,,,,,,,, MPTEMP,1,0MPTEMP,,,,,,,,MPTEMP,1,0MPDATA,DENS,2,,8900 !定义材料2的密度MPTEMP,,,,,,,,MPTEMP,1,0MPDATA,C,2,,390 !定义材料2的⽐热MPTEMP,,,,,,,,MPTEMP,1,0MPDATA,KXX,2,,383 !定义材料2的热传导系数MPTEMP,,,,,,,,MPTEMP,1,0UIMP,2,REFT,,,30 !定义材料2的热膨胀系数的参考温度MPDATA,ALPX,2,,1.75e-5 !定义材料2的热膨胀系数MPTEMP,,,,,,,, MPTEMP,1,0MPDATA,EX,2,,103e9 !定义材料2的弹性模量MPDATA,PRXY,2,,0.3 !定义材料2的泊松⽐CYLIND,0.7,0.5,0,0.2,0,180, !建⽴钢环的⼏何模型CYLIND,0.5,0.3,0,0.2,60,70, !建⽴铜滑块FLST,5,2,4,ORDE,2FITEM,5,8FITEM,5,12CM,_Y,LINELSEL, , , ,P51XCM,_Y1,LINECMSEL,,_YLESIZE,_Y1, , ,4, , , , ,1 !选择第8和第12号线并划分为4段FLST,5,2,4,ORDE,2FITEM,5,2FITEM,5,7CM,_Y,LINELSEL, , , ,P51XCM,_Y1,LINECMSEL,,_YMSHAPE,1,3DMSHKEY,0FLST,5,2,4,ORDE,2FITEM,5,2LSEL, , , ,P51XCM,_Y1,LINECMSEL,,_YLESIZE,_Y1, , ,36, , , , ,1 !选择第2和第7号线并划分为36段FLST,5,1,4,ORDE,1 FITEM,5,5CM,_Y,LINELSEL, , , ,P51XCM,_Y1,LINELESIZE,_Y1, , ,36, , , , ,1 !选择第5号线并划分为36段FLST,5,2,4,ORDE,2FITEM,5,17FITEM,5,19CM,_Y,LINELSEL, , , ,P51XCM,_Y1,LINECMSEL,,_YLESIZE,_Y1, , ,2, , , , ,1 ! 选择第17和第19号线并划分为2段FLST,5,2,4,ORDE,2 FITEM,5,20FITEM,5,24CM,_Y,LINELSEL, , , ,P51XCM,_Y1,LINECMSEL,,_YLESIZE,_Y1, , ,4, , , , ,1 !选择第20和第24号线并划分为4段TYPE, 1 !给钢环附单元类型MA T, 1 !给钢环附材料REAL, !给钢环附实常数ESYS, 0SECNUM,CM,_Y,VOLUVSEL, , , , 1CM,_Y1,VOLUCHKMSH,'VOLU'CMSEL,S,_YVSWEEP,_Y1CMDELE,_Y2 !给铜滑块划分单元TYPE, 1 !给铜滑块附单元类型MA T, 2 !给铜滑块附材料REAL, !给铜滑块附实常数ESYS, 0SECNUM,CM,_Y,VOLUVSEL, , , , 2CM,_Y1,VOLUCHKMSH,'VOLU'CMSEL,S,_YVSWEEP,_Y1CMDELE,_YCMDELE,_Y2 !给铜滑块附实常数/COM, CONTACT PAIR CREATION - START CM,_NODECM,NODECM,_ELEMCM,ELEMCM,_KPCM,KPCM,_LINECM,LINECM,_AREACM,AREACM,_VOLUCM,VOLU/GSA V,cwz,gsav,,tempMP,MU,1,0.2MA T,1R,3REAL,3ET,2,170ET,3,174KEYOPT,3,9,0R,3,RMORE,RMORE,,0RMORE,0KEYOPT,3,1,3TYPE,2NSLA,S,1ESLN,S,0ESURF,ALLCMSEL,S,_ELEMCM! Generate the contact surface ASEL,S,,,9CM,_CONTACT,AREA TYPE,3NSLA,S,1ESLN,S,0ESURF,ALLALLSELESEL,ALLESEL,S,TYPE,,2ESEL,A,TYPE,,3ESEL,R,REAL,,3/PSYMB,ESYS,1/PNUM,TYPE,1/NUM,1EPLOTESEL,ALLESEL,S,TYPE,,2ESEL,A,TYPE,,3ESEL,R,REAL,,3 CMSEL,A,_NODECM CMDEL,_NODECM CMSEL,A,_ELEMCM CMDEL,_ELEMCM CMSEL,S,_KPCM CMDEL,_KPCM CMSEL,S,_LINECMCMSEL,S,_VOLUCMCMDEL,_VOLUCM/GRES,cwz,gsavCMDEL,_TARGETCMDEL,_CONTACT/COM, CONTACT PAIR CREATION - END !⽤接触向导建⽴铜滑块与钢环间的接触/GO *DIM,pre,TABLE,2,1,1,time,press, !定义压⼒载荷随时间变化的TABLE表*SET,PRE(1,0,1) , 0*SET,PRE(1,1,1) , 20e6*SET,PRE(2,0,1) , 10*SET,PRE(2,1,1) , 20e6FLST,2,1,5,ORDE,1FITEM,2,10/GOSFA,P51X,1,PRES, %PRE% !在铜滑块上表⾯施加压⼒载荷TUNIF,20, !施加初始温度FLST,2,1,5,ORDE,1FITEM,2,3/GODA,P51X,UX, !在钢环的外表⾯施加UX向约束FLST,2,1,5,ORDE,1FITEM,2,3/GODA,P51X,UY, !在钢环的外表⾯施加UY向约束FLST,2,1,5,ORDE,1FITEM,2,3/GODA,P51X,UZ, !在钢环的外表⾯施加UZ向约束ESEL,S,MA T,,2 !选择铜滑块单元NSLE !选择铜滑块节点CSYS,1 !激活柱坐标系FLST,2,75,1,ORDE,2NROTA T,P51X !旋转节点坐标系FLST,2,75,1,ORDE,2FITEM,2,876FITEM,2,-950/GOD,P51X, ,0.01, , , ,UY, , , , , ! 给铜滑块加Y向位移载荷（即切向位移）ALLSEL !选择所有的节点，单元，⼏何体FINISH/SOL !进⼊求解器ANTYPE,4 !设置为瞬态分析TRNOPT,FULL !选择瞬态全分析⽅法LUMPM,0 !关闭集中质量矩阵计算格式NLGEOM,1 !定义⼤变形NSUBST,100,200,100 !定义载荷步OUTRES,ALL,20 !定义结果输出⼦步间隔AUTOTS,1 !打开时间开关KBC,0 !设置为斜坡载荷TIME,10 !定义求解时间SSTIF,0 !关闭应⼒刚化影响NROPT,FULL, , !选择瞬态全分析⽅法EQSLV, , ,0,PRECISION,0 !定义双精度求解MSA VE,0TOFFST,273, !定义温度偏移量SOLVE !求解/POST1 !进⼊后处理器SET,LAST !读取最后⼦步结果PLNSOL,TEMP !显⽰温度分布云图PLNSOL,S,EQV,2 !显⽰等效应⼒分布云图FINISH/EXIT,NOSAV !退出ANSYSFINISH/FILNAME,Exercise23 ! 定义隐式热分析⽂件名/PREP7 !进⼊前处理器ET,1,PLANE13,4 !选择单元类型MPTEMP,,,,,,,,MPTEMP,1,0MPDATA,C,1,,460 !定义材料1的⽐热MPTEMP,,,,,,,,MPTEMP,1,0MPDATA,KXX,1,,66.6 !定义材料1的热传导系数MPTEMP,,,,,,,,MPTEMP,1,0UIMP,1,REFT,,,30 !定义材料1的热膨胀系数的参考温度MPDATA,ALPX,1,,1.06e-5 !定义材料1的热膨胀系数MPTEMP,,,,,,,, MPTEMP,1,0MPDATA,EX,1,,206e9 !定义材料1的弹性模量MPDATA,PRXY,1,,0.3 !定义材料1的泊松⽐MPTEMP,,,,,,,,MPTEMP,1,0MPDATA,DENS,2,,8900 !定义材料2的密度MPTEMP,,,,,,,,MPTEMP,1,0MPDATA,C,2,,390 !定义材料2的⽐热MPTEMP,,,,,,,,MPTEMP,1,0MPDATA,KXX,2,,383 !定义材料2的热传导系数MPTEMP,,,,,,,,MPTEMP,1,0UIMP,2,REFT,,,30 !定义材料2的热膨胀系数的参考温度MPDATA,ALPX,2,,1.75e-5 !定义材料2的热膨胀系数MPTEMP,,,,,,,, MPTEMP,1,0MPDATA,EX,2,,103e9 !定义材料2的弹性模量MPDATA,PRXY,2,,0.3 !定义材料2的泊松⽐RECTNG,0,5,0,1.25, !建⽴钢固定块的⼏何模型RECTNG,0,1.25,1.25,2.5, !建⽴铜滑块的⼏何模型TYPE, 1MA T, 1REAL, !给钢固定块附材料属性ESYS, 0SECNUM,ESIZE,0.25,0, !设置单元划分尺⼨MSHAPE,0,2DMSHKEY,1CM,_Y,AREAASEL, , , , 1CM,_Y1,AREACHKMSH,'AREA'CMSEL,S,_YAMESH,_Y1CMDELE,_Y2 !给钢固定块划分单元TYPE, 1MA T, 2REAL, !给铜滑块附材料属性ESYS, 0SECNUM,CM,_Y,AREAASEL, , , , 2CM,_Y1,AREACHKMSH,'AREA'CMSEL,S,_YAMESH,_Y1CMDELE,_YCMDELE,_Y1CMDELE,_Y2 !给铜滑块划分单元CM,_NODECM,NODECM,_ELEMCM,ELEMCM,_KPCM,KPCM,_LINECM,LINECM,_AREACM,AREACM,_VOLUCM,VOLUMP,MU,1,0.2MA T,1MP,EMIS,1,R,3REAL,3ET,2,169ET,3,172R,3,,,1.0,0.1,0,RMORE,,,1.0E20,0.0,1.0,RMORE,0.0,0,1.0,,1.0,0.5RMORE,0,1.0,1.0,0.0,,1.0 KEYOPT,3,3,0 KEYOPT,3,4,0 KEYOPT,3,5,0 KEYOPT,3,7,0 KEYOPT,3,8,0 KEYOPT,3,9,1 KEYOPT,3,10,1 KEYOPT,3,11,0 KEYOPT,3,12,0 KEYOPT,3,2,0 KEYOPT,3,1,1LSEL,S,,,3CM,_TARGET,LINE TYPE,2NSLL,S,1ESLN,S,0CM,_CONTACT,LINE TYPE,3NSLL,S,1ESLN,S,0ESURF,ALLALLSELESEL,ALLESEL,S,TYPE,,2 ESEL,A,TYPE,,3 ESEL,R,REAL,,3 ESEL,ALLESEL,S,TYPE,,2 ESEL,A,TYPE,,3 ESEL,R,REAL,,3 CMSEL,A,_NODECM CMDEL,_NODECM CMSEL,A,_ELEMCM CMDEL,_ELEMCM CMSEL,S,_KPCM CMDEL,_KPCM CMSEL,S,_LINECM CMDEL,_LINECM CMSEL,S,_AREACM CMDEL,_AREACMCMSEL,S,_VOLUCMCMDEL,_VOLUCMCMDEL,_TARGETCMDEL,_CONTACTFLST,2,1,5,ORDE,1FITEM,2,1 !⽤接触向导建⽴滑块和固定块间的接触⾯/GODA,P51X,UX, !将下⾯的钢块施加X向约束FLST,2,1,5,ORDE,1FITEM,2,1/GODA,P51X,UY, !将下⾯的钢块施加Y向约束FLST,2,1,4,ORDE,1FITEM,2,6/GODL,P51X, ,UX,3.75 ! 给铜块右侧的线加X向位移载荷FINISH/SOL*DIM,pre,TABLE,2,1,1,time,press, !定义压⼒载荷随时间变化的TABLE表*SET,PRE(1,0,1) , 0*SET,PRE(1,1,1) , 10e6*SET,PRE(2,0,1) , 10*SET,PRE(2,1,1) , 10e6FLST,2,1,4,ORDE,1FITEM,2,7/GOSFL,P51X,PRES, %PRE% !在铜滑块上表⾯施加压⼒载荷TUNIF,30, !给铜块和钢块施加初始温度NLGEOM,ON !定义⼤变形TIME,3.75E-3 !定义求解时间TOFFSET,273 !定义温度偏移量AUTO,ON !打开时间开关NSUB,100,10000,100 !定义载荷步OUTRES,ALL,-10 !定义结果输出⼦步间隔SOLVE !求解/POST1 !进⼊后处理器SET,LAST !读取最后⼦步结果PLNSOL,TEMP ! 显⽰温度分布云图PLNSOL,S,EQV,2 !显⽰等效应⼒分布云图FINISH/EXIT,NOSAV !退出ANSYS2. 电磁三d命令流实例（论坛看到）/BATCH/COM,ANSYS RELEASE 12.0.1 UP20090415 11:10:14 03/29/2011 /input,menust,tmp,'',,,,,,,,,,,,,,,,1 /GRA,POWER/GST,ON/PLO,INFO,3/GRO,CURL,ON/CPLANE,1/REPLOT,RESIZEWPSTYLE,,,,,,,,0/TITLE, Test for EDGE-Formulation with SOLID117!! Parameters*SET,l,5e-6 ! length*SET,h,1e-6 ! height*SET,w,5e-6 ! width*SET,a,3e-6 ! 1st addition length of magnetic core*SET,b,5e-6 ! 2nd addition length of magnetic core!/PREP7 !!ET,1,117 !!MP,MURX,1,1MP,RSVX,1,17.2e-9 !MP,MURX,2,500MP,MURX,3,1BLOCK,-w,w,-h,h,-l,l !!BLOCK,-w-a,w+a,h+a,h+2*a,-l,lBLOCK,w+a,w+a+b,h+a,h+2*a,-l,lBLOCK,-w-a,-w-a-b,h+a,h+2*a,-l,lBLOCK,w+a,w+a+b,h,h+a,-l,lBLOCK,-w-a,-w-a-b,h,h+a,-l,l!BLOCK,-w-a,w+a,-(h+a),-(h+2*a),-l,lBLOCK,w+a,w+a+b,-(h+a),-(h+2*a),-l,lBLOCK,-w-a,-w-a-b,-(h+a),-(h+2*a),-l,lBLOCK,w+a,w+a+b,-h,-(h+a),-l,lBLOCK,-w-a,-w-a-b,-h,-(h+a),-l,l!BLOCK,-22e-6,22e-6,-15e-6,15e-6,-15e-6,15e-6 VOVLAP,ALL NUMCMP,VOLUVSEL,S,VOLU,,1V ATT,1VSEL,S,VOLU,,2,11,1V ATT,2VSEL,S,VOLU,,12V ATT,3VSEL,S,MA T,,1,2,1ASLV,SLSLA,SLESIZE,ALL,3e-6MSHKEY,1MSHAPE,0,3DVMESH,ALLVSEL,S,MA T,,3ASLV,SLSLA,S LESIZE,ALL,3e-6 MSHKEY,0 MSHAPE,1,3D VMESH,ALL ALLSEL,ALL FINISH/SOLUESEL,S,MA T,,1 BFE,ALL,JS,1,,,-1e9 ASEL,S,EXT DA,ALL,AZ,0 ALLSEL,ALL ANTYPE,0 NROPT,AUTO,, EQSLV,SPAR,,0, PRECISION,0 MSA VE,0 PIVCHECK,1/STA TUS,SOLU SOLVEFINISH/POST1! elements/VSCALE,1,1,0!。

ANSYS案例——20例ANSYS经典实例】

ANSYS案例——20例ANSYS经典实例】针对【典型例题】3.3.7(1)的模型，即如图3-19所示的框架结构，其顶端受均布力作用，用有限元方法分析该结构的位移。

结构中各个截面的参数都为：113.010PaE=，746.510mI-=，426.810mA-=，相应的有限元分析模型见图3-20。

在ANSYS平台上，完成相应的力学分析。

图3-19框架结构受一均布力作用图3-20单元划分、节点位移及节点上的外载解答对该问题进行有限元分析的过程如下。

1．基于图形界面的交互式操作(tepbytep)(1)进入ANSYS(设定工作目录和工作文件)程序→ANSYS→ANSYSInteractive→Workingdirectory(设置工作目录)→Initialjobname(设置工作文件名):beam3→Run→OK(2)设置计算类型(3)选择单元类型(4)定义材料参数ANSYSMainMenu:Preproceor→MaterialProp→MaterialModel→Struc tural→Linear→Elatic→Iotropic:E某:3e11(弹性模量)→OK→鼠标点击该窗口右上角的“”来关闭该窗口(5)定义实常数以确定平面问题的厚度ANSYSMainMenu:Preproceor→RealContant…→Add/Edit/Delete→Add→Type1Beam3→OK→RealContantSetNo:1(第1号实常数),Cro-ectionalarea:6.8e-4(梁的横截面积)→OK→Cloe(6)生成几何模型生成节点ANSYSMainMenu:Preproceor→Modeling→Creat→Node→InActiveCS→Nodenumber1→某:0,Y:0.96,Z:0→Apply→Nodenumber2→某:1.44,Y:0.96,Z:0→Apply→Nodenumber3→某:0,Y:0,Z:0→Apply→Nodenumber4→某:1.44,Y:0,Z:0→OK生成单元ANSYSMainMenu:Preproceor→Modeling→Create→Element→AutoNum bered→ThruNode→选择节点1，2(生成单元1)→apply→选择节点1，3(生成单元2)→apply→选择节点2，4(生成单元3)→OK(7)模型施加约束和外载左边加某方向的受力ANSYSMainMenu:Solution→DefineLoad→Apply→Structural→Force/Moment→OnNode→选择节点1→apply→Directionofforce:F某→VALUE：3000→OK→上方施加Y方向的均布载荷ANSYSMainMenu:Solution→DefineLoad→Apply→Structural→Preure→OnBeam→选取单元1(节点1和节点2之间)→apply→VALI：4167→VALJ：4167→OK左、右下角节点加约束(8)分析计算(9)结果显示(10)退出系统(11)计算结果的验证与MATLAB支反力计算结果一致。

悬索桥ansys命令流——很经典,很实用。

r,12,35.236,208.180,111.180,5.600,7.788,0.000
rmore,0.000,243.090,0.000,0.000,0.000,0.000
r,13,42.949,210.520,112.240,5.600,6.770,0.000
rmore,0.000,251.720,0.000,0.000,0.000,0.000
r,34,0.0196,0.0013,0.0002,0.4000,0.5980,0.000
rmore,0.000,0.0000,0.000,0.000,0.000,0.000 !ZF1
r,35,0.0177,0.0011,0.0001,0.3400,0.5980,0.000
!加劲梁截面实常数定义*********************************************************************
r,30,0.0445,0.0028,0.0028,0.6360,0.6360,0.000
rmore,0.000,0.0043,0.000,0.000,0.000,0.000 !ZS1 and ZX
ET,5,BEAM4 !中央扣单元
ET,6,BEAM4 !主塔单元
ET,7,BEAM4 !支撑桥面板小纵梁单元
r,10,19.616,164.530,85.363,5.600,8.260,0.000
rmore,0.000,181.050,0.000,0.000,0.000,0.000
r,11,19.145,149.890,82.626,5.600,7.966,0.000
rmore,0.000,171.430,0.000,0.000,0.000,0.000

ansys经典命令流

FINISH/CLEAR/PREP7esel,s,mat,,1 选择材料号为1的单元LCOMB,ALL!将两条线合并为一条线VOFFST,1,L!沿面1的法线方向偏移生成体LOCAL,12,1,L,,,,-90!定义圆柱局部坐标系LFILLT,2,3,R!两线之间生成圆角ARSYM,Y,ALL!镜像ARSYM,X,ALLnummrg,all合并所有项，合并相同位置处lgen,2,all,,,,sl复制项NUMCMP，Label（压缩编号消除空号）latt,1,3,1给线附属性材料号、实常数号、线单元类型号*do,i,3,7*enddo循环wpoff,,,h1-h2vsbw,all 移动工作面然后切割WPCSYS，WN，KCN将当前坐标系X-Y平面定义为工作平面/post1pldisp,1后处理，显示变形lesize,all,,,4lmesh,all定义网格划分长度为4，然后划分mshape,0mshape, key, dimension 指定单元形状key: 0 四边形（2D），六面体（3D）1 三角形(2D), 四面体(3D) Dimension: 2D 二维3D 三维mshkey,1mshkey, key 指定自由或映射网格方式0 自由网格划分1 映射网格划分2 如果可能的话使用映射，esll,sESLL, Type（选择已选线上的单元）pstres,on打开预应力分析sectype,1,beam,rectSECTYPE, ID, TYPE, SUBTYPE, NAME, REFINEKEY定义一个截面号，并初步定义截面类型ID: 截面号TYPE: BEAM:定义此截面用于梁SUBTYPE: RECT 矩形CSOLID:圆形实心截面CTUBE: 圆管I: 工字HREC: 矩形空管ASEC: 任意截MESH: 用户定义的划分网格NAME: 8字符的截面名称（字母和数字组成）REFINEKEY: 网格细化程度：0~5secdata,b,h描述梁截面lsymm,x,all线镜像nlgeom,1包括几何非线性lstr,8,116用两个点来定义一条直线ekill,all 杀死所有选择的单元NLIST,NODE1,NODE2,NINC,Lcoord,SORT1,SORT2,SORT3列出节点信息。

ansys命令流实例

ansys命令流实例1.结构静力分析－－梁分析/PREP7 !进入前处理器ET,1,BEAM3 !定义单元类型R,1,0.25,0.0052,0.5 !定义实常数MP,EX,1,210E6MP,PRXY,1,0.3 !定义材料属性N,1,0N,2,1N,3,5N,4,7N,5,7N,6,9N,7,11N,8,11N,9,13N,10,14 !定义节点E,1,2E,2,3E,3,4E,5,6E,6,7E,8,9E,9,10 !生成单元CP,1,UX,4,5CP,2,UY,4,5CP,3,UX,7,8CP,4,UY,7,8 !耦合节点FINISH/SOLU !进入求解器D,2,UXD,2,UYD,3,UYD,6,UYD,9,UY !施加位移约束F,10,FY,-4 !施加集中约束SFBEAM,1,1,PRES,4,4SFBEAM,2,1,PRES,4,4 !施加均布力SOLVE !求解/POST1 !进入后处理器PLDISP !绘制结构变形图PRDISP !列出各节点的位移ETABLE,IMOMENT,SMISC,6ETABLE,JMOMENT,SMISC,12ETABLE,ISHEAR,SMISC,2ETABLE,JSHEAR,SMISC,8 !将节点弯矩、剪力制表PRETAB !列表显示单元的弯矩、剪力/TITLE,SHEAR FORCE DISTRIBUTION !设置剪力分布图的标题PLLS,ISHEAR,JSHEAR !绘制剪力分布图/TITLE,BENDING MOMENT IDSTRIBUTION !设置弯矩分布图的标题PLLS,IMOMENT,JMOMENT !绘制弯矩分布图流固耦合实例RAD=0.8 !底面半径H=1G=9.8OMEGAR=2ROU=1000 !定义参数变量/PREP7 !进入前处理器ET,1,FLUID79 !选择单元类型KEYOPT,1,3,1 !设置单元关键字MP,EX,1,2E9 !设置杨氏模量MP,DENS,1,ROU !设置材料密度K,1K,2,RADK,3,RAD,HK,4,,H !生成关键点A,1,2,3,4 !连接关键点生成面积LESIZE,ALL,,,10 !设置网格划分精度AMESH,ALL !将面积划分网络/SOLU !进入求解器DL,2,,UXDL,1,,UYNSEL,S,LOC,XDSYM,SYMM,XD,ALL,UXNSEL,ALL !施加位移约束ACEL,,GOMEGA,,OMEGAR !施加惯性力SOLVE !求解/POST1SET,LAST !进入通用后处理器PLNSOL,U,X,0,1 !绘制应力云图UCENT=UY(22)UEDGE=UY(12)UELEV=UEDGE-UCENT !提取节点位移结构静力分析－－壳结构内力分析LENGTH=100YOUNG=200000THICKNESS=2FORCE=1000DENSITY=9E-6 !将材质、载荷、板的几何尺寸等参数化/PREP7 !进入前处理器MP,EX,1,YOUNGMP,NUXY,1,0.3MP,DENS,1,DENSITY !定义材质ET,1,SHELL63 !定义单元类型R,1,THICKNESS,THICKNESS,THICKNESS,THICKNESS !定义实常数!构建结构的几何模型K,1,0,0K,2,LENGTH,0K,3,LENGTH,LENGTHK,4,0,LENGTH !定义关键点A,1,2,3,4!通过关键点生成面LSEL,ALLLESIZE,ALL,,,16AMESH,ALL !设定网格划分参数，划分网格FINISH/SOLU !进入求解器NSEL,S,LOC,X,0,0D,ALL,ALL,0 !选择X＝0的节点将其固定NSEL,S,LOC,X,LENGTH,LENGTHD,ALL,ALL,0 !选择X＝LENGTH的节点将其固定NSEL,S,LOC,X,0.5*LENGTH,0.5*LENGTHNSEL,R,LOC,Y,0.5*LENGTH,0.5*LENGTHF,ALL,FZ,FORCEALLSEL !捕捉板的中心点并在中心点处施加集中力荷载SOLVE !求解FINISH/POST1 !进入后处理器/DSC,,10PLNSOL,U,Z,0,1 !绘图显示板的竖向变形NSEL,ALL !提取板的最大竖向变形NSORT,U,Z,1,1 !将节点的位移绝对值以升序排序*GET,MAXDEFLECTION,SORT,0,MAX !提取位移最大值并赋给变量NSEL,S,LOC,X,0NSEL,A,LOC,X,LENGTH,LENGTH !选择固定边节点NSORT,S,EQV,1,1 !将节点等效应力的绝对值以升序排序*GET,MAXSTRESS,SORT,0,MAX !提取等效应力最大值并赋给变量结构静力分析－－平面桁架分析/PREP7 !进入前处理器ET,1,LINK1 !选择单元R,1,0.1 !定义实常数MP,EX,1,30E6MP,PRXY,1,0.3 !定义材料属性N,1,0N,2,4N,3,8N,4,12N,5,0,3N,6,4,3N,7,8,3 !生成节点E,1,2E,2,3E,3,4E,4,7E,3,7E,2,7E,2,6E,2,5E,1,5E,5,6E,6,7 !生成单元FINISH/SOLU !进入求解器D,1,UXD,1,UYD,5,UX !实加位移约束F,2,FY,-15F,3,FY,-15F,4,FY,-15 !施加集中力SOLVE !求解/POST1 !进入能用后处理器PRESOL,FORC !列表显示反力。

ANSYS命令流实例

/PREP7 !进入前处理ANTYPE,STATIC !设置分析类型为静力结构分析PSTRES,ON !用于后面的模态分析中考虑预应力(该开关不影响静力分析)ET,1,LINK10 !选取单元类型1(单向杆单元)KEYOPT,1,3,0 !设置仅承受拉应力，KEYOPT(3)=0R,1,306796E-8,543248E-8 !设置实常数，包括绳索截面积(306796E-8)，初始应变(543248E-8) MP,EX,1,30E6 !定义材料的弹性模量(1号材料)MP,DENS,1,73E-5 !定义材料的密度(1号材料)N,1 ! 定义第1号节点N,14,100 ! 定义第14号节点FILL ! 均分填满第2号至第13号节点E,1,2 !由节点1及节点2生成单元EGEN,13,1,1 !依序复制生成13个单元D,ALL,ALL ! 对所有节点施加固定约束FINISH ! 前处理结束/SOLU ! 进入求解模块，求解预应力引起的应力状态SOLVE ! 求解FINISH ! 退出求解模块/POST1 ! 进入一般的后处理ETABLE,STRS,LS,1 !针对LINK10单元,建立单元列表STRS,通过LS及特征号1来获得单元的轴向应力*GET,STRSS,ELEM,13,ETAB,STRS !针对单元列表STRS, 提取13号单元的应力FINISH ! 后处理结束/POST26 ! 进入时间历程后处理，处理支反力RFORCE,2,1,F,X !将1号节点上的x方向支反力提取，并存储到2号变量中STORE ! 存储*GET,FORCE,V ARI,2,EXTREM,VMAX !将2号变量的最大值赋给参数FORCE/SOLU ! 再次进入求解模块，模态分析ANTYPE,MODAL ! 模态分析MODOPT,SUBSP,3 ! 选择子空间迭代法，求3阶模态MXPAND,3 ! 设定3阶模态扩展PSTRES,ON ! 用于在模态分析中考虑预应力(还需在前面的静力分析中也同时打开) DDELE,2,UX,13 ! 删除从2号节点到13号节点上的UX约束DDELE,2,UY,13 !删除从2号节点到13号节点上的UY约束SOLVE !求解*GET,FREQ1,MODE,1,FREQ ! 提取第1阶模态共振频率，并赋值给参数FREQ1*GET,FREQ2,MODE,2,FREQ ! 提取第2阶模态共振频率，并赋值给参数FREQ2*GET,FREQ3,MODE,3,FREQ ! 提取第3阶模态共振频率，并赋值给参数FREQ3*STATUS !列出所有参数的实际内容。

ansys经典命令流资料裂缝

发信人: wjc501 (阿昌), 信区: FEA标题: Re: wjc501请帮帮忙，有没有裂缝方面的ansys例题发信站: BBS 水木清华站 (Mon Jul 22 18:21:57 2002)供参考，不当之处请指正。

祝好运!(1) 如果你模拟结构体中裂缝扩展过程的模拟，在Ansys中可以用全解耦损伤分析方法来近似模拟裂缝扩展，我曾用Ansys软件中提供的可以定义10,000个材料参数和单元ekill/alive功能完成了层状路面体中表面裂缝和反射裂缝在变温作用下的扩展过程的模拟。

我模拟的过程相对来说比较简单，模拟过程中我们首先要知道裂缝的可能扩展方向，这样在裂缝可能扩展的带内进行网格加密处理，加密到什么程度依据计算的问题来确定。

(2) 如果采用断裂力学理论计算含裂缝结构体的应力强度因子，建模时只需在裂尖通过命令kscon生成奇异单元即可。

Ansys模块中存在的断裂力学模块可以计算I、II、III型应力强度因子(线弹性断裂力学)和J积分(弹塑性断裂力学)，在Ansys中verification里面有一个计算I型应力强度因子的例子vm143，参见该例子就可以了。

(3) 如果通过断裂力学模拟裂缝的扩展过程，需要采用动态网格划分，这方面我没有做，通过Ansys的宏命令流应该可以实现。

技术参考可参阅文献：杨庆生、杨卫.断裂过程的有限元模拟.计算力学学报，1997,14(4).(4) 我现在做动荷载作用下路面结构体中应力强度因子的分布规律，我是通过位移插值得到不同时间点处的应力强度因子。

如果想这样做，可参阅理论参考中关于应力强度因子计算说明。

【在 agz (agz) 的大作中提到: 】: 我在进行混凝土中裂缝模拟分析，急需这方面的例题，听说你在这方面是大希故而想请模: 谢谢啦。

--※来源:·BBS 水木清华站 ·[FROM: 210.32.156.200]发信人: agz (agz), 信区: FEA标题: Re: wjc501请帮帮忙，有没有裂缝方面的ansys例题发信站: BBS 水木清华站 (Mon Jul 22 21:01:07 2002)谢谢阿昌的热心指导，我的情况属于第二种，试了一个二维的问题，基本知道了一点儿关于裂缝分析的过程。

(完整版)ANSYS最常用命令流+中文注释(超级大全)

ANSYS最常用命令流+中文注释VSBV, NV1, NV2, SEPO, KEEP1, KEEP2 —Subtracts volumes from volumes，用于2个solid相减操作，最终目的是要nv1-nv2=?通过后面的参数设置，可以得到很多种情况：sepo项是2个体的边界情况，当缺省的时候，是表示2个体相减后，其边界是公用的，当为sepo的时候，表示相减后，2个体有各自的独立边界。

keep1与keep2是询问相减后，保留哪个体？当第一个为keep时，保留nv1,都缺省的时候，操作结果最终只有一个体，比如：vsbv,1,2,sepo,,keep,表示执行1-2的操作，结果是保留体2，体1被删除，还有一个1-2的结果体，现在一共是2个体（即1-2与2），且都各自有自己的边界。

如vsbv,1,2,,keep,,则为1-2后，剩下体1和体1-2，且2个体在边界处公用。

同理，将v换成a 及l是对面和线进行减操作！mp,lab, mat, co, c1,…….c4 定义材料号及特性lab: 待定义的特性项目（ex,alpx,reft,prxy,nuxy,gxy,mu,dens） ex: 弹性模量nuxy: 小泊松比alpx: 热膨胀系数reft: 参考温度reft: 参考温度prxy: 主泊松比gxy: 剪切模量mu: 摩擦系数dens: 质量密度mat: 材料编号（缺省为当前材料号）co: 材料特性值，或材料之特性，温度曲线中的常数项c1-c4: 材料的特性-温度曲线中1次项，2次项，3次项，4次项的系数定义DP材料：首先要定义EX和泊松比：MP，EX，MA T，……MP，NUXY，MAT，……定义DP材料单元表（这里不考虑温度）：TB，DP，MA T进入单元表并编辑添加单元表：TBDATA，1，CTBDATA，2，ψTBDATA，3，……如定义：EX=1E8，NUXY=0.3，C=27，ψ=45的命令如下：MP，EX，1，1E8MP，NUXY，1，0.3TB，DP，1TBDATA，1，27TBDATA，2，45这里要注意的是，在前处理的最初，要将角度单位转化到“度”，即命令：*afun,degVSEL, Type, Item, Comp, VMIN, VMAX, VINC, KSWP Type，是选择的方式，有选择（s）,补选（a），不选(u)，全选(all)、反选（inv）等,其余方式不常用Item, Comp 是选取的原则以及下面的子项如volu 就是根据实体编号选择，loc 就是根据坐标选取，它的comp就可以是实体的某方向坐标！其余还有材料类型、实常数等MIN, VMAX, VINC，这个就不必说了吧！,例：vsel,s,volu,,14vsel,a,volu,,17,23,2上面的命令选中了实体编号为14，17，19，21，23的五个实体VDELE, NV1, NV2, NINC, KSWP: 删除未分网格的体nv1:初始体号nv2:最终的体号ninc:体号之间的间隔kswp=0:只删除体kswp=1:删除体及组成关键点,线面如果nv1=all,则nv2,ninc不起作用其后面常常跟着一条显示命令VPLO,或aplo,nplo,这个湿没有参数的命令，输入后直接回车，就可以显示刚刚选择了的体、面或节点，很实用的哦！Nsel, type, item, comp, vmin, vmax, vinc, kabs 选择一组节点为下一步做准备Type: S: 选择一组新节点（缺省）R: 在当前组中再选择A: 再选一组附加于当前组U: 在当前组中不选一部分All: 恢复为选中所有None: 全不选Inve: 反向选择Stat: 显示当前选择状态Item: loc: 坐标node: 节点号Comp: 分量Vmin,vmax,vinc: ITEM范围Kabs: “0” 使用正负号“1”仅用绝对值下面是单元生死第一个载荷步中命令输入示例：!第一个载荷步TIME,... !设定时间值（静力分析选项）NLGEOM,ON !打开大位移效果NROPT,FULL !设定牛顿－拉夫森选项ESTIF,... !设定非缺省缩减因子（可选）ESEL,... !选择在本载荷步中将不激活的单元EKILL,... !不激活选择的单元ESEL,S,LIVE !选择所有活动单元NSLE,S !选择所有活动结点NSEL,INVE !选择所有非活动结点（不与活动单元相连的结点）D,ALL,ALL,0 !约束所有不活动的结点自由度（可选）NSEL,ALL !选择所有结点ESEL,ALL !选择所有单元D,... !施加合适的约束F,... !施加合适的活动结点自由度载荷SF,... !施加合适的单元载荷BF,... !施加合适的体载荷SA VESOLVE请参阅TIME,NLGEOM,NROPT,ESTIF,ESEL,EKILL,NSLE,NSEL,D, F,SF和BF命令得到更详细的解释。

ANSYS最常用命令流+中文注释(超级大全)

ANSYS最常用命令流+中文注释VSBV, NV1, NV2, SEPO, KEEP1, KEEP2 — Subtracts volumes from volumes，用于2个solid相减操作，最终目的是要nv1-nv2=?通过后面的参数设置，可以得到很多种情况：sepo项是2个体的边界情况，当缺省的时候，是表示2个体相减后，其边界是公用的，当为sepo的时候，表示相减后，2个体有各自的独立边界。

如vsbv,1,2,,keep,,则为1-2后，剩下体1和体1-2，且2个体在边界处公用。

同理，将v换成a及l是对面和线进行减操作！mp,lab, mat, co, c1,…….c4 定义材料号及特性lab: 待定义的特性项目（ex,alpx,reft,prxy,nuxy,gxy,mu,dens） ex: 弹性模量nuxy: 小泊松比alpx: 热膨胀系数reft: 参考温度reft: 参考温度prxy: 主泊松比gxy: 剪切模量mu: 摩擦系数dens: 质量密度mat: 材料编号（缺省为当前材料号）co: 材料特性值，或材料之特性，温度曲线中的常数项c1-c4: 材料的特性-温度曲线中1次项，2次项，3次项，4次项的系数定义DP材料：首先要定义EX和泊松比：MP，EX，MA T，……MP，NUXY，MAT，……定义DP材料单元表（这里不考虑温度）：TB，DP，MA T 进入单元表并编辑添加单元表：TBDATA，1，CTBDATA，2，ψTBDATA，3，……如定义：EX=1E8，NUXY=0.3，C=27，ψ=45的命令如下：MP，EX，1，1E8MP，NUXY，1，0.3TB，DP，1 TBDATA，1，27TBDATA，2，45这里要注意的是，在前处理的最初，要将角度单位转化到“度”，即命令：*afun,degVSEL, Type, Item, Comp, VMIN, VMAX, VINC, KSWP Type，是选择的方式，有选择（s）,补选（a），不选(u)，全选(all)、反选（inv）等,其余方式不常用Item, Comp 是选取的原则以及下面的子项如volu 就是根据实体编号选择，loc 就是根据坐标选取，它的comp就可以是实体的某方向坐标！其余还有材料类型、实常数等MIN, VMAX, VINC，这个就不必说了吧！,例：vsel,s,volu,,14vsel,a,volu,,17,23,2上面的命令选中了实体编号为14，17，19，21，23的五个实体VDELE, NV1, NV2, NINC, KSWP: 删除未分网格的体nv1:初始体号nv2:最终的体号ninc:体号之间的间隔kswp=0:只删除体kswp=1:删除体及组成关键点,线面如果nv1=all,则nv2,ninc不起作用其后面常常跟着一条显示命令VPLO,或aplo,nplo,这个湿没有参数的命令，输入后直接回车，就可以显示刚刚选择了的体、面或节点，很实用的哦！Nsel, type, item, comp, vmin, vmax, vinc, kabs 选择一组节点为下一步做准备Type: S: 选择一组新节点（缺省）R: 在当前组中再选择A: 再选一组附加于当前组U: 在当前组中不选一部分All: 恢复为选中所有None: 全不选Inve: 反向选择Stat: 显示当前选择状态Item: loc: 坐标node: 节点号Comp: 分量Vmin,vmax,vinc: ITEM范围Kabs: “0” 使用正负号“1”仅用绝对值下面是单元生死第一个载荷步中命令输入示例：!第一个载荷步TIME,... !设定时间值（静力分析选项）NLGEOM,ON !打开大位移效果NROPT,FULL !设定牛顿－拉夫森选项ESTIF,... !设定非缺省缩减因子（可选）ESEL,... !选择在本载荷步中将不激活的单元EKILL,... !不激活选择的单元ESEL,S,LIVE !选择所有活动单元NSLE,S !选择所有活动结点NSEL,INVE !选择所有非活动结点（不与活动单元相连的结点）D,ALL,ALL,0 !约束所有不活动的结点自由度（可选）NSEL,ALL !选择所有结点ESEL,ALL !选择所有单元D,... !施加合适的约束F,... !施加合适的活动结点自由度载荷SF,... !施加合适的单元载荷BF,... !施加合适的体载荷SA VESOLVE请参阅TIME,NLGEOM,NROPT,ESTIF,ESEL,EKILL,NSLE,NSEL, D,F,SF和BF命令得到更详细的解释。

ansys建模实例-混凝土箱梁(命令流)

/TITLE, BRIDGE XX/PREP7! 定义单元类型ET,1,SHELL63ET，2，SOLID45ET,3,LINK8!定义1,钢,2，混凝土，3,无限刚度MP，EX,1,2.1E11MP,DENS,1,7。

8E3MP,PRXY,1，0.3mp,alpx,1,1.0e—5MP，EX,2，3。

5E10MP，DENS，2,2.6E3MP,PRXY，2，0。

1667MP，EX，3，3。

5E15MP,DENS,3，7.8E3MP，PRXY，3,0。

3*set,nz，154!定义箱梁N1R,1,7，7，7,7!定义箱梁N2R,2,8,8,8，8!定义箱梁N3R,3，4，4,4，4!定义混凝土solid65的实常数R, 4!定义钢绞线的实常数R， 5， 140!降温=预加力/线膨胀系数*钢筋截面积*弹性模量tem1=—1190.4e6/（1。

0e-5*2。

11e11）tem2=-1339。

2e6/（1。

0e—5＊2.11e11）！定义纵向梯度参数表＊dim， z， array， nz*SET,z(1)，130,260，371,482，593，702,811,920,1029，1138，1247,1356，1465，1574，1683,1787.19,1891。

38，1995.57＊SET,z(19），2099。

76，2203。

95，2308。

14,2412.33,2516.52，2620.71,2724。

9,2829。

09，2933.28，3037。

47,3141.66,3245.85，＊SET，z(31），3350.04,3454。

23,3558.42,3662.61，3766。

8，3870。

99*SET，z（37),3975。

18,4079.37,4183.56,4287.75，*SET,z（41),4391.94,4496.13,4600。

32,4704。

51，4808.7，4912。

89,5017.08,5121.27，5225。

ANSYS土木工程经典实例命令流大全解析

1．ANSYS SOLID65环向布置钢筋的例子 (3)2.混凝土非线性计算实例（1）- MISO单压 (5)3.混凝土非线性计算实例（2）-MISO约束压 (6)4.混凝土非线性计算实例（3）- KINH滞回 (9)5.混凝土非线性计算实例（4）- KINH压-拉裂 (11)6.混凝土非线性计算实例（5） (12)7.混凝土非线性计算实例（6） (14)8.混凝土非线性计算实例（7）-MISO滞回 (16)9.混凝土非线性计算实例（8） (18)10.混凝土非线性计算实例（9）-梁平面应力 (20)11.四层弹簧－质点模型的地震分析 (22)12.悬臂梁地震分析 (48)13.用beam 54单元描述变截面梁的例子 (72)14.变截面梁实例 (73)15.拱桥浇筑过程分析-单元生死应用实例 (74)16.简支梁实体与预应力钢筋分析实例 (75)17. 简单的二维焊接分析－单元生死实例 (77)18.隧道开挖（三维）的命令流 (84)19.岩土接触分析实例 (101)20.钢筋混凝土管的动力响应特性分析实例 (109)21.隧道模拟开挖命令流（入门） (116)22.螺栓连接的模拟实现问题 (119)23.道路的基层、垫层模量与应力之间的关系 (129)23.滞回分析 (151)24.模拟某楼层浇注 (153)25.在面上施加移动的面力 (155)27.在任意面施加任意方向任意变化的压力 (159)28.预紧分析 (160)29.几何非线性+塑性+接触+蠕变 (162)30.埋设在地下的排水管道 (167)32.幕墙企业玻璃简化计算 (172)33.等截面杆单元生死应用实例 (188)34.梁板建模联系 (189)36.简单的例子－如何对结构的振动控制分析 (192)37.模态分析结果的输出实例 (194)38.火车过桥动态加载实例（部分） (196)39.悬索结构的找形和计算的例题 (213)40.陶瓷杆撞击铝板的例子 (218)41.求反作用力的APDL命令法 (221)42.LS-DYNA实例（部分） (222)43.路面分层填筑对路基的影响 (223)44.一个例子（含地震影响，求振兴与频率） (227)45.接触面上的压力总和 (231)46.施加位置函数荷载 (235)47.非线性分析考虑刚度退化 (236)48.一个圆形水池的静力分析 (237)49.ANSYS中混凝土模式预应力模拟的算例 (238)50.悬臂梁受重力作用发生大变形求其固有频率 (240)51.循环对称结构模态分析 (242)52.三角平台受谐波载荷作用的结构响应 (244)53.三角平台受一地震谱激励的应力分布和支反力 (246)54.三角平台受时程载荷作用的应力分布和变形过程 (248)55.经典层合板理论 (250)56．定易圆轨迹的例子 (257)57．模拟门式刚架施工－单元生死 (257)58．钢筋混凝土整体式模型例子 (260)59．在荷载步之间改变材料属性例子 (262)60．含预应力的特征值屈曲计算 (263)61．振型叠加计算及工况组合例子 (265)62．柱子稳定分析算（预应力，特征值屈曲，初始缺陷） (268)63．module MConcrete !混凝土模板 (271)64.混凝土开裂实例 (279)65．螺栓网格划分 (280)66．自由液面的土石坝平面渗流分析 (281)67．导出刚度矩阵 (285)68．某混凝土拱坝工程施工期及运行期温度场仿真分析 (286)69．移动温度荷载计算 (293)70．SHSD用于壳-实体装配实例An (295)71．ansys显示－隐式－回弹分析实例 (299)72．工况组合的经典例子 (314)1．ANSYS SOLID65环向布置钢筋的例子! 一个管道，环向配筋率为1%，纵向配筋率为0.5%，径向配筋率为0.1%!FINISH/CLEAR/PREP7!*! 单元属性ET,1,SOLID65!*KEYOPT,1,1,0KEYOPT,1,5,0KEYOPT,1,6,0KEYOPT,1,7,1!*!实参数1：不同方向配筋R,1,2,.001, , ,2, .01,RMORE, 90, ,2,.005 ,90 ,90 ,!Adds real constants to a set.!材料属性!混凝土基本材料属性MPTEMP,,,,,,,,MPTEMP,1,0MPDATA,EX,1,,30E2MPDATA,PRXY,1,,.2!屈服准则TB,MISO,1,1,5,TBTEMP,0TBPT,,0.0005,15TBPT,,0.001,21TBPT,,0.0015,24TBPT,,0.002,27TBPT,,0.003,24!破坏准则TB,CONC,1,1,9,TBTEMP,0TBDATA,,.5,.9,3,30,,!抗拉和抗压是10 倍的关系TBDATA,,,,1,,,MPTEMP,,,,,,,,MPTEMP,1,0!钢材基本属性MPDATA,EX,2,,200E3MPDATA,PRXY,2,,.27!屈服准则TB,BISO,2,1,2,TBTEMP,0TBDATA,,310,2E3,,,,! 管道内径和外径CYL4,0,0,3000, , , ,10000CYL4,0,0,2000, , , ,10000VSBV, 1, 2! 定义局部柱坐标CSWPLA,11,1,1,1,! Defines a local coordinate system at the origin of the working plane. KWPAVE, 11!move to keypoint11WPRO,,-90.000000,VSBW, 3WPCSYS,-1,0!Defines the working plane location based on a coordinate system. KWPAVE, 1WPRO,,,-90.000000VSBW,ALLESIZE,500,0,! 注意：设定单元局部坐标VATT, 1, 1, 1, 11!*VSWEEP,ALL/DEVICE,VECTOR,1/ESHAPE,1.0!Displays elements with shapes determined from the real constants or section/REPLO! 注意：红色代表最大配筋方向，绿色代表其次，蓝色表示最小配筋方向2.混凝土非线性计算实例（1）- MISO单压!MISO单压FINISH/CLEAR/PREP7ET,1,SOLID65R,1MP,EX,1,30E3MP,NUXY,1,0.2TB,CONCR,1,1TBTEMP,0TBDATA,1,0.3,0.5,2.5,-1TB,MISO,1,1,15TBPT,, 100E-6, 3.0TBPT,, 300E-6, 8.3TBPT,, 600E-6,14.6TBPT,, 900E-6,19.1TBPT,,1100E-6,21.0TBPT,,1250E-6,22.0TBPT,,1400E-6,22.6TBPT,,1550E-6,22.8TBPT,,1650E-6,22.7TBPT,,1800E-6,22.3TBPT,,2000E-6,21.4TBPT,,2800E-6,16.8TBPT,,3200E-6,14.7TBPT,,3800E-6,12.3TBPT,,4600E-6,9.9!TB,MKIN,1!TBTEMP,,STRAIN!TBDATA,,600E-6,1100E-6,1600E-6,3000E-6,4500E-6 !TBTEMP,0!TBDATA,,15,21.5,23,16,9BLOCK,0,50,0,50,0,50MSHAPE,0,3DMSHKEY,1ESIZE,10VMESH,ALLNSEL,S,LOC,X,0D,ALL,UX,0NSEL,S,LOC,Y,0D,ALL,UY,0NSEL,S,LOC,Z,0D,ALL,UZ,0NSEL,S,LOC,Y,50CP,1,UY,ALLNSEL,ALLFINISH/VIEW,1,1,1,1/REPLOT/SOLUANTYPE,STATICAUTOTS,OFFOUTRES,ALL,1TIME,50NSUBST,50D,1,UY,-.25NSEL,ALLSOLVEFINISH/POST26NSOL,2,1,U,Y,UYRFORCE,3,1,F,Y,FYADD,4,2,,,STRAIN,,,-1/50ADD,5,3,,,STRESS,,,-1/2500/AXLAB,X,STRAIN/AXLAB,Y,STRESSXVAR,4PLVAR,5FINISH3.混凝土非线性计算实例（2）-MISO约束压!MISO约束压FINISH/CLEAR/PREP7ANTYPE,STATICET,1,SOLID65R,1MP,EX,1,30E3MP,NUXY,1,0.2TB,CONCR,1,1TBTEMP,0TBDATA,1,0.3,0.5,2.5,-1TB,MISO,1,1,15TBPT,, 100E-6, 3.0TBPT,, 300E-6, 8.3TBPT,, 600E-6,14.6TBPT,, 900E-6,19.1TBPT,,1100E-6,21.0TBPT,,1250E-6,22.0TBPT,,1400E-6,22.6TBPT,,1550E-6,22.8TBPT,,1650E-6,22.7TBPT,,1800E-6,22.3TBPT,,2000E-6,21.4TBPT,,2800E-6,16.8TBPT,,3200E-6,14.7TBPT,,3800E-6,12.3TBPT,,4600E-6,9.9!TB,MKIN,1!TBTEMP,,STRAIN!TBDATA,,600E-6,1100E-6,1600E-6,3000E-6,4500E-6 !TBTEMP,0!TBDATA,,15,21.5,23,16,9BLOCK,0,50,0,50,0,50ESIZE,10VMESH,ALLMSHAPE,0,3DNSEL,S,LOC,X,0D,ALL,UX,0NSEL,S,LOC,Y,0D,ALL,UY,0NSEL,S,LOC,Z,0D,ALL,UZ,0NSEL,S,LOC,Y,50CP,1,UY,ALLD,ALL,UX,0D,ALL,UZ,0NSEL,ALLFINISH/VIEW,1,1,1,1/REPLOT/SOLUOUTRES,ALL,1TIME,20NSUBST,20,0,20D,1,UY,-.04LSWRITE,1TIME,30NSUBST,20,0,20D,1,UY,-.06LSWRITE,2TIME,40NSUBST,100D,1,UY,-.1LSWRITE,3LSSOLVE,1,3FINISH/POST26NSOL,2,1,U,Y,UY RFORCE,3,1,F,Y,FYADD,4,2,,,STRAIN,,,-1/50 ADD,5,3,,,STRESS,,,-1/2500 /AXLAB,X,STRAIN/AXLAB,Y,STRESSXVAR,4PLVAR,5FINISH4.混凝土非线性计算实例（3）- KINH滞回!KINH滞回FINISH/CLEAR/PREP7ANTYPE,STATICET,1,SOLID65R,1MP,EX,1,30E3MP,NUXY,1,0.2TB,CONCR,1,1TBTEMP,0TBDATA,1,0.3,0.5,2.5,-1TB,KINH,1,1,10TBPT,, 150E-6, 4.5TBPT,, 600E-6,14.8TBPT,,1000E-6,20.25TBPT,,1300E-6,22.3TBPT,,1480E-6,22.8TBPT,,1620E-6,22.8TBPT,,1800E-6,22.3TBPT,,2000E-6,21.4TBPT,,3500E-6,12.8TBPT,,5000E-6, 9.0BLOCK,0,50,0,50,0,50ESIZE,10VMESH,ALLMSHAPE,0,3DNSEL,S,LOC,X,0D,ALL,UX,0NSEL,S,LOC,Y,0D,ALL,UY,0NSEL,S,LOC,Z,0D,ALL,UZ,0NSEL,S,LOC,Y,50CP,1,UY,ALLNSEL,ALLFINISH/VIEW,1,1,1,1/REPLOT/SOLUOUTRES,ALL,ALLTIME,20NSUBST,20,0,20D,1,UY,-.04LSWRITE,1TIME,35NSUBST,15,0,15D,1,UY,-.01LSWRITE,2TIME,90NSUBST,55,0,55D,1,UY,-.12LSWRITE,3TIME,125NSUBST,35,0,35D,1,UY,-.05LSWRITE,4TIME,175NSUBST,50,0,50D,1,UY,-.15LSWRITE,5LSSOLVE,1,5FINISH/POST26NSOL,2,1,U,Y,UY RFORCE,3,1,F,Y,FYADD,4,2,,,STRAIN,,,-1/50 ADD,5,3,,,STRESS,,,-1/2500 /AXLAB,X,STRAIN/AXLAB,Y,STRESSXVAR,4PLVAR,5FINISH5.混凝土非线性计算实例（4）- KINH压-拉裂!KINH压——拉裂压到峰值，泄载、反向加载到拉裂FINISH/CLEAR/PREP7ANTYPE,STATICET,1,SOLID65R,1MP,EX,1,30E3MP,NUXY,1,0.2TB,CONCR,1,1TBTEMP,0TBDATA,1,0.3,0.5,2.5,-1TB,KINH,1,1,10TBPT,, 150E-6, 4.5TBPT,, 600E-6,14.8TBPT,,1000E-6,20.25TBPT,,1300E-6,22.3TBPT,,1480E-6,22.8TBPT,,1620E-6,22.8TBPT,,1800E-6,22.3TBPT,,2000E-6,21.4TBPT,,3500E-6,12.8TBPT,,5000E-6, 9.0BLOCK,0,50,0,50,0,50ESIZE,10VMESH,ALLMSHAPE,0,3DNSEL,S,LOC,X,0D,ALL,UX,0NSEL,S,LOC,Y,0D,ALL,UY,0NSEL,S,LOC,Z,0D,ALL,UZ,0NSEL,S,LOC,Y,50CP,1,UY,ALLNSEL,ALLFINISH/VIEW,1,1,1,1/REPLOT/SOLUOUTRES,ALL,ALLTIME,20NSUBST,20,0,20D,1,UY,-.04LSWRITE,1TIME,45NSUBST,25,0,25D,1,UY,.01LSWRITE,2LSSOLVE,1,2FINISH/POST26NSOL,2,1,U,Y,UYRFORCE,3,1,F,Y,FYADD,4,2,,,STRAIN,,,-1/50ADD,5,3,,,STRESS,,,-1/2500/AXLAB,X,STRAIN/AXLAB,Y,STRESSXVAR,4PLVAR,5FINISH6.混凝土非线性计算实例（5）FINISH/CLEAR/PREP7ANTYPE,STATICET,1,SOLID65R,1MP,EX,1,30E3MP,NUXY,1,0.2TB,CONCR,1,1 TBTEMP,0TBDATA,1,0.3,0.5,2.5,-1 TB,KINH,1,1,18 TBPT,, 100E-6, 3.0 TBPT,, 300E-6, 8.3 TBPT,, 600E-6,14.6 TBPT,, 900E-6,19.1 TBPT,,1100E-6,21.0 TBPT,,1250E-6,22.0 TBPT,,1400E-6,22.6 TBPT,,1550E-6,22.8 TBPT,,1650E-6,22.7 TBPT,,1800E-6,22.3 TBPT,,2000E-6,21.4 TBPT,,2400E-6,19.1 TBPT,,2800E-6,16.8 TBPT,,3200E-6,14.7 TBPT,,3600E-6,13.0 TBPT,,4100E-6,11.3 TBPT,,4600E-6,9.9BLOCK,0,50,0,50,0,50 ESIZE,10VMESH,ALLMSHAPE,0,3DNSEL,S,LOC,X,0D,ALL,UX,0NSEL,S,LOC,Y,0D,ALL,UY,0NSEL,S,LOC,Z,0D,ALL,UZ,0NSEL,S,LOC,Y,50CP,1,UY,ALLNSEL,ALLFINISH/VIEW,1,1,1,1/REPLOT/SOLUOUTRES,ALL,ALLTIME,50NSUBST,50,0,50D,1,UY,-.25SOLVEFINISH/POST26NSOL,2,1,U,Y,UYRFORCE,3,1,F,Y,FYADD,4,2,,,STRAIN,,,-1/50ADD,5,3,,,STRESS,,,-1/2500/AXLAB,X,STRAIN/AXLAB,Y,STRESSXVAR,4PLVAR,5FINISH7.混凝土非线性计算实例（6）!MISO约束压FINISH/CLEAR/PREP7ANTYPE,STATICET,1,SOLID65R,1MP,EX,1,30E3MP,NUXY,1,0.2TB,CONCR,1,1TBTEMP,0TBDATA,1,0.3,0.5,2.5,-1TB,KINH,1,1,15TBPT,, 100E-6, 3.0TBPT,, 300E-6, 8.3TBPT,, 600E-6,14.6TBPT,, 900E-6,19.1TBPT,,1100E-6,21.0TBPT,,1250E-6,22.0TBPT,,1400E-6,22.6TBPT,,1550E-6,22.8TBPT,,1650E-6,22.7TBPT,,1800E-6,22.3TBPT,,2000E-6,21.4TBPT,,2800E-6,16.8TBPT,,3200E-6,14.7TBPT,,3800E-6,12.3TBPT,,4600E-6,9.9!TB,MKIN,1!TBTEMP,,STRAIN!TBDATA,,600E-6,1100E-6,1600E-6,3000E-6,4500E-6 !TBTEMP,0!TBDATA,,15,21.5,23,16,9BLOCK,0,50,0,50,0,50ESIZE,10VMESH,ALLMSHAPE,0,3DNSEL,S,LOC,X,0D,ALL,UX,0NSEL,S,LOC,Y,0D,ALL,UY,0NSEL,S,LOC,Z,0D,ALL,UZ,0NSEL,S,LOC,Y,50CP,1,UY,ALLD,ALL,UX,0D,ALL,UZ,0NSEL,ALLFINISH/VIEW,1,1,1,1/REPLOT/SOLUOUTRES,ALL,1TIME,20NSUBST,20,0,20D,1,UY,-.04LSWRITE,1TIME,30NSUBST,20,0,20D,1,UY,-.06LSWRITE,2TIME,40NSUBST,100D,1,UY,-.1LSWRITE,3LSSOLVE,1,3FINISH/POST26NSOL,2,1,U,Y,UYRFORCE,3,1,F,Y,FYADD,4,2,,,STRAIN,,,-1/50ADD,5,3,,,STRESS,,,-1/2500/AXLAB,X,STRAIN/AXLAB,Y,STRESSXVAR,4PLVAR,5FINISH8.混凝土非线性计算实例（7）-MISO滞回!MISO滞回FINISH/CLEAR/PREP7ANTYPE,STATICET,1,SOLID65R,1MP,EX,1,30E3MP,NUXY,1,0.2TB,CONCR,1,1TBTEMP,0TBDATA,1,0.3,0.5,2.5,-1TB,MISO,1,1,15 TBPT,, 100E-6, 3.0 TBPT,, 300E-6, 8.3 TBPT,, 600E-6,14.6 TBPT,, 900E-6,19.1 TBPT,,1100E-6,21.0 TBPT,,1250E-6,22.0 TBPT,,1400E-6,22.6 TBPT,,1550E-6,22.8 TBPT,,1650E-6,22.7 TBPT,,1800E-6,22.3 TBPT,,2000E-6,21.4 TBPT,,2800E-6,16.8 TBPT,,3200E-6,14.7 TBPT,,3800E-6,12.3 TBPT,,4600E-6,9.9BLOCK,0,50,0,50,0,50 ESIZE,10VMESH,ALL MSHAPE,0,3DNSEL,S,LOC,X,0D,ALL,UX,0NSEL,S,LOC,Y,0D,ALL,UY,0NSEL,S,LOC,Z,0D,ALL,UZ,0NSEL,S,LOC,Y,50 CP,1,UY,ALLNSEL,ALLFINISH/VIEW,1,1,1,1/REPLOT/SOLUOUTRES,ALL,ALL TIME,20NSUBST,20,0,20D,1,UY,-.04 LSWRITE,1TIME,35NSUBST,15,0,15D,1,UY,-.01LSWRITE,2TIME,60NSUBST,25,0,25D,1,UY,-.06LSWRITE,3TIME,80NSUBST,20,0,20D,1,UY,-.02LSWRITE,4TIME,110NSUBST,30,0,30D,1,UY,-.08LSWRITE,5LSSOLVE,1,5FINISH/POST26NSOL,2,1,U,Y,UYRFORCE,3,1,F,Y,FYADD,4,2,,,STRAIN,,,-1/50ADD,5,3,,,STRESS,,,-1/2500/AXLAB,X,STRAIN/AXLAB,Y,STRESSXVAR,4PLVAR,5FINISH9.混凝土非线性计算实例（8）!MISO单压FINISH/CLEAR/PREP7ANTYPE,STATICET,1,SOLID65R,1MP,EX,1,30E3MP,NUXY,1,0.2TB,CONCR,1,1TBTEMP,0TBDATA,1,0.3,0.5,2.5,-1TB,MISO,1,1,10TBPT,, 150E-6, 4.5TBPT,, 600E-6,14.8TBPT,,1000E-6,20.25TBPT,,1300E-6,22.3TBPT,,1480E-6,22.8TBPT,,1620E-6,22.8TBPT,,1800E-6,22.3TBPT,,2000E-6,21.4TBPT,,3500E-6,12.8TBPT,,5000E-6, 9.0!TB,MKIN,1!TBTEMP,,STRAIN!TBDATA,,600E-6,1100E-6,1600E-6,3000E-6,4500E-6 !TBTEMP,0!TBDATA,,15,21.5,23,16,9BLOCK,0,50,0,50,0,50ESIZE,10VMESH,ALLMSHAPE,0,3DNSEL,S,LOC,X,0D,ALL,UX,0NSEL,S,LOC,Y,0D,ALL,UY,0NSEL,S,LOC,Z,0D,ALL,UZ,0NSEL,S,LOC,Y,50CP,1,UY,ALLNSEL,ALLFINISH/VIEW,1,1,1,1/REPLOT/SOLUOUTRES,ALL,1TIME,50NSUBST,100,0,100D,1,UY,-.25SOLVEFINISH/POST26NSOL,2,1,U,Y,UYRFORCE,3,1,F,Y,FYADD,4,2,,,STRAIN,,,-1/50ADD,5,3,,,STRESS,,,-1/2500/AXLAB,X,STRAIN/AXLAB,Y,STRESSXVAR,4PLVAR,5FINISH10.混凝土非线性计算实例（9）-梁平面应力!梁平面应力FINISH/CLEAR/PREP7ET,1,SOLID65R,1MP,EX,1,30E3MP,NUXY,1,0.2TB,CONCR,1,1TBTEMP,0TBDATA,1,0.3,0.5,2.5,-1TB,MISO,1,1,15TBPT,, 100E-6, 3.0TBPT,, 300E-6, 8.3TBPT,, 600E-6,14.6TBPT,, 900E-6,19.1TBPT,,1100E-6,21.0TBPT,,1250E-6,22.0TBPT,,1400E-6,22.6TBPT,,1550E-6,22.8TBPT,,1650E-6,22.7TBPT,,1800E-6,22.3TBPT,,2000E-6,21.4TBPT,,2800E-6,16.8 TBPT,,3200E-6,14.7 TBPT,,3800E-6,12.3 TBPT,,4600E-6, 9.9ET,2,LINK8R,2,30R,3,6MP,EX,2,200E3MP,NUXY,2,0.3TB,BKIN,2,1TBTEMP,0TBDATA,,300,0/ESHAPE,1/TRIAD,LBOT/VIEW,1,1,1,1N,1,0,0,0N,2,20,0,0NGEN,11,2,1,2,1,0,20 NGEN,43,22,1,22,1,0,0,20 E,1,2,4,3,23,24,26,25 EGEN,10,2,1TYPE,2REAL,2MAT,2E,3,25E,4,26REAL,3E,19,41E,20,42EGEN,42,22,1,14E,47,63E,48,64EGEN,10,88,589,590/NUMBER,1/PNUM,TYPE,1/REPLOTNSEL,S,LOC,Z,0D,ALL,UZ,0NSEL,S,LOC,Y,0NSEL,R,LOC,Z,740D,ALL,UY,0NSEL,ALLFINISH/SOLUANTYPE,STATICNSEL,S,LOC,Y,200NSEL,R,LOC,Z,200D,ALL,UY,-2TIME,200NSUBST,10LSWRITE,1D,ALL,UY,-0.21TIME,21NSUBST,100LSWRITE,2D,ALL,UY,-1TIME,100NSUBST,100LSWRITE,3NSEL,ALLAUTOTS,1OUTRES,ALL,1PRED,ONLSSOLVE,1,3FINISH11.四层弹簧－质点模型的地震分析包括频率分析，谱分析和时程分析作者：陆新征，清华大学土木工程系FINI/CLEAR/UNITS,SI/PREP7DAMPRATIO=0.02 !振型阻尼系数都为0.02ET,1,BEAM188KEYOPT,1,7,1KEYOPT,1,8,1ET,2,MASS21,,,4MP,EX,1,210E9MP,NUXY,1,0.3MP,DENS,1,7800MPDATA,DAMP,1,,DAMPRATIO SECTYPE,3,BEAM,ISECDATA,0.2,0.2,0.5,0.01,0.01,0.008 R,1,1.6E2,1.6E2R,2,1.2E2,1.2E2K,1,0,0,0K,2,0,3,0K,3,0,6,0K,4,0,9,0K,5,0,12,0K,6,0,0,10000L,1,2L,2,3L,3,4L,4,5LSEL,,,,1,4LATT,1,0,1,,6,,3ALLSLESIZE,ALL, , ,3 , , , , ,1 LMESH,ALLNUMMRG,ALLNUMCMP,ALLTYPE,2REAL,1E,2TYPE,2REAL,2E,8E,14E,20*SET,NT,1001*SET,DT,0.02*DIM,AC,,NT*VREAD,AC(1),RECORD,TXT(F8.3)/SOLUD,1,ALL,ANTYPE,0SOLVEFINISH/SOLU !模态分析ANTYPE,2MODOPT,SUBSP,8MXPAND,8, , ,1SOLVEFINI!得到自振频率1*GET,FREQ1,MODE,1,FREQ/ESHAPE,1.0/SOLUANTYPE,SPECTR !谱分析SPOPT,SPRSMDAMP,1,DAMPRATIO, , , , , ,SED,1SVTYPE,2,1FREQ,0.167,0.25,0.333,0.44444,0.5,0.667,1,1.25,1.667SV,,0.154350625,0.191590625,0.210210625,0.224175625,0.250716714,0.329 514922,0.484352764,0.598723486,0.786897371FREQ,2.222,10,100000SV,,1.034212766,1.034212766,0.3528!按7度多遇地震，地震影响系数0.08 第一组III场地TG＝0.45SRSS,0.01,DISPSOLVEFINISH/POST1/INP,,MCOMNSORT,U,SUM*GET,DMAX_SPE,SORT, ,MAX !反应谱最大变形!开始地震计算FINISH/CONFIG,NRES,20000/SOLUANTYPE,TRANSTRNOPT,FULLALPHAD,2*DAMPRATIO*FREQ1*2*3.1415926BETAD,2*DAMPRATIO/(FREQ1*2*3.1415926)*DO,I,1,1001ACEL,AC(I),0,0TIME,I*0.02OUTRES,ALL,ALLSOLVE*ENDDOFINISH/POST26FILENSOL,2,20,U,X,PLVAR,2, , , , , , , , , ,STORE,MERGEXVAR,1PLVAR,2,*GET,DMAX_TRANS,VARI,2,EXTREM,VMAX, , ! 时程分析最大位移数据文件record.txt-0.0014-0.0108-0.0101-0.0088-0.0095-0.012-0.0142-0.0128-0.011-0.0085-0.0085-0.0131-0.0176-0.0194-0.0162-0.0144-0.0108-0.0082-0.0042-0.0066-0.0131-0.019-0.0196-0.00660.0030.0141-0.0049-0.0144 -0.0203 -0.026 -0.0325 -0.0306 -0.0172 -0.0197 -0.0163 -0.0164 -0.0067 0.0025 0.015 0.0236 0.0252 0.0336 0.0463 0.0492 0.0419 0.0359 0.0271 0.0235 0.0339 0.0412 0.053 0.0639 0.0732 0.0652 0.0599 0.04 0.04 0.0063 -0.0515 -0.0787 -0.0603 -0.0484 -0.025 -0.0059 0.0134 0.0308 0.0499 0.071 0.0995 0.12190.1449 0.1155 0.0935 0.0892 0.0926 0.0839 0.0901 0.0993 0.1209 0.0328 -0.1475 -0.2066 -0.1989 -0.2034 -0.1816 -0.1725 -0.1752 -0.1753 -0.1805 -0.163 -0.1347 -0.1087 -0.0782 -0.0429 -0.0017 0.036 0.0785 0.1164 0.1598 0.196 0.2412 0.2729 0.3036 0.32 0.3417 0.2821 0.2324 -0.1198 -0.2373 -0.164 -0.1865 -0.1095 -0.07530.0113 0.0533 0.0895 0.1186 0.1757 0.0576 -0.2631 -0.1547 -0.1729 -0.1012 -0.0579 0.0237 -0.067 -0.198 -0.1641 -0.1685 -0.1481 -0.1231 -0.1001 -0.0751 -0.0523 -0.0271 -0.0044 0.0188 -0.0095 -0.0433 -0.0838 -0.0951 -0.0716 -0.0599 -0.0334 -0.0108 0.0185 0.042 0.0673 -0.0097 -0.0372 -0.004 0.0011 0.0344 0.0565 0.0883 0.1130.0219 0.0241 0.0683 0.0689 0.1318 0.1353 0.204 -0.0931 -0.1308 -0.0692 -0.0546 0.0072 0.0675 -0.1067 -0.1488 -0.1071 -0.1162 -0.0762 -0.0559 -0.0215 -0.0126 -0.0674 -0.0324 -0.0337 -0.0109 0.0017 0.0299 0.0488 0.0608 0.0222 -0.0032 -0.0245 0.0077 0.0211 0.0568 0.0826 0.1206 0.1478 0.1737 0.0421 0.0029 0.0259 0.0293-0.0147 0.0143 0.0206 0.0499 0.0645 0.0957 0.1128 0.1447 0.1629 0.1945 0.1856 0.1984 0.1769 0.125 -0.1207 -0.0542 -0.0384 -0.0311 -0.1118 -0.1661 -0.2464 -0.2025 -0.1835 -0.1317 -0.096 -0.0325 0.0154 0.0816 0.1319 0.1818 -0.0058 -0.0169 0.0285 0.0447 0.0983 0.1424 0.1853 0.2456 0.1685 -0.138 -0.0999 -0.1089 -0.0907-0.125 -0.2111 -0.1617 -0.1692 -0.1306 -0.1111 -0.0773 -0.051 -0.0544 -0.12 -0.1209 -0.1158 -0.1145 -0.0717 -0.0546 0.0064 -0.0804 -0.1634 -0.0859 -0.0961 -0.0396 -0.0147 0.0319 0.0648 0.0876 0.0472 0.0198 -0.0027 0.0292 0.0445 0.0785 0.1033 0.1352 0.1606 0.1861 0.1281 0.064 0.0204 0.0314 0.0373 0.0496 0.0235 -0.0084-0.0113 -0.0229 -0.0248 -0.0157 -0.0069 0.0147 0.0379 0.0579 0.0255 -0.0041 -0.0428 -0.0133 0.0095 0.023 -0.0129 -0.005 0.008 0.021 0.038 0.051 0.0157 -0.0032 -0.0111 0.0005 0.0076 0.0035 -0.0095 -0.0036 -0.0016 0.0038 0.0085 -0.0056 -0.0304 -0.0421 -0.0244 -0.0236 -0.0177 -0.0129 -0.0018 0.0203 -0.0108 -0.0091 -0.0034-0.0111 -0.0099 -0.0002 0.0073 0.0235 0.0355 0.0705 0.0779 0.0184 -0.0263 -0.0124 -0.0042 0.0159 0.0048 -0.0219 -0.0467 -0.0428 -0.0216 -0.0043 0.0159 0.032 0.0419 0.0123 -0.016 -0.0204 -0.0082 -0.0206 -0.0137 -0.0055 0.0053 0.0134 0.0266 0.0232 0.0079 -0.0008 0.02 0.0435 0.0492 0.0191 0.0092 -0.0022 -0.0021 0.00520.0255 0.0368 0.0525 0.0541 0.0425 0.0398 0.0559 0.0756 0.0365 0.0411 0.0098 -0.0204 -0.0249 -0.0405 -0.0413 -0.0471 -0.0433 -0.0458 -0.0057 0.0178 -0.0208 -0.0492 -0.053 -0.0362 -0.0405 -0.0308 -0.0316 -0.0265 -0.0265 -0.0269 -0.0345 -0.0309 -0.0217 -0.0078 0.0087 0.0281 0.031 0.0358 0.0341 0.0358 0.0287 0.0305 0.01120.0136 0.0384 -0.0861 -0.1349 -0.1342 -0.1354 -0.1193 -0.1042 -0.0829 -0.0651 -0.0444 -0.0258 -0.006 -0.0091 -0.0182 -0.0147 0.0085 0.0163 0.005 0.0264 0.0582 0.0867 0.12 0.1695 0.1111 -0.11 -0.0366 -0.0445 -0.0236 -0.096 -0.0656 -0.0597 -0.067 -0.0552 -0.0027 0.0378 0.1072 0.1669 0.0947 0.0408 0.0667 0.0132 -0.0095-0.0827 -0.1152 -0.115 -0.0803 -0.0369 0.0029 0.0545 0.1178 0.161 -0.027 0.0034 -0.0056 0.002 0.0146 0.0537 0.0798 -0.0205 -0.059 -0.0169 -0.0175 -0.0028 0.0074 0.0382 0.0567 0.0753 0.0801 0.0592 0.0304 0.0023 0.0064 -0.0406 -0.0451 -0.0079 0.0168 0.0567 0.0093 -0.0055 0.0044 -0.0123 -0.0282 -0.0437 -0.0352 -0.02550.0205 0.0519 0.0854 0.1144 0.0733 0.0237 -0.0368 -0.0271 -0.0217 -0.0873 -0.0973 -0.0589 -0.0336 0.0077 0.0259 0.0508 0.0361 0.0081 -0.0056 -0.0209 -0.0317 -0.0238 -0.0376 -0.055 -0.0722 -0.0803 -0.0523 -0.034 -0.0011 0.0065 -0.0037 -0.0005 -0.0168 -0.041 -0.008 0.0079 0.0374 0.0615 0.0665 0.0254 -0.0057 -0.0474 -0.0356-0.0048 0.0126 0.0379 0.0241 -0.0227 -0.0428 -0.0679 -0.0661 -0.059 -0.0513 -0.0408 -0.0309 -0.0266 -0.0541 -0.0628 -0.0908 -0.1107 -0.0881 -0.077 -0.0582 -0.0473 -0.0333 -0.0199 0.002 0.0211 0.0432 0.0613 0.0767 0.0933 0.1066 0.113 0.1187 0.1247 0.1334 0.1594 0.1797 0.2037 0.1236 0.0442 -0.014 -0.0666 -0.0555 -0.0693-0.1246 -0.1179 -0.105 -0.092 -0.0743 -0.0809 -0.085 -0.086 -0.0863 -0.0873 -0.0868 -0.0885 -0.0537 0.0052 0.0215 0.0245 0.058 0.0314 0.0236 0.0485 0.0589 0.0525 0.0355 0.0197 0.0199 0.0492 0.0343 0.0288 0.0432 0.0239 0.0088 0.0077 -0.0148 -0.0077 -0.0019 0.0075 0.0044 -0.0145 -0.0316 -0.0241 -0.0028 0.0182 0.04260.0512 0.0466 0.0479 0.0193 0.0222 0.0274 0.0393 0.0504 0.0577 0.0588 0.0822 0.0797 0.0949 0.0345 0.0045 -0.0123 -0.0347 -0.0426 -0.0416 -0.0275 -0.027 0.0074 0.0428 -0.0231 -0.0387 -0.0083 0.0139 0.0445 0.0027 -0.0697 -0.0796 -0.0251 -0.0135 0.0079 -0.0115 -0.0251 -0.0333 -0.0269 -0.0301 -0.02 -0.0067 -0.0038 0.01050.0344 0.0957 0.0898 0.0179 -0.0362 -0.0994 -0.0807 -0.0744 -0.0539 -0.033 -0.0128 0.0031 0.0148 0.0508 -0.0022 -0.0489 -0.0358 -0.0691 -0.0516 -0.0371 0.0088 0.0632 0.0841 0.1276 0.1388 0.1193 0.0751 0.0225 -0.0088 -0.0227 0.0074 0.0181 0.0544 0.0399 0.0045 -0.0082 -0.0185 -0.002 0.0006 -0.0117 -0.021 -0.0303 -0.0512-0.0579 -0.0266 -0.0178 0.004 0.0098 0.0137 0.0221 0.0437 0.0091 -0.0548 -0.0555 -0.0243 -0.0081 0.025 0.041 0.0182 -0.0027 -0.0243 -0.0015 0.0247 0.0482 0.0783 0.0622 0.0331 -0.0014 -0.0195 -0.0247 -0.0212 -0.011 0.005 0.0241 -0.0034 -0.0216 -0.0471 -0.0363 -0.0195 -0.0018 0.017 -0.008 0.0005 0.023 0.0374 0.06010.0432 0.0344 0.0505 0.0653 0.0683 0.0172 -0.017 -0.0527 -0.0664 -0.0387 -0.0222 -0.0033 0.0119 -0.0128 -0.0351 -0.0514 -0.0335 -0.0218 -0.0012 0.0142 0.007 -0.0063 -0.012 -0.0322 -0.0346 -0.0091 0.0073 0.0309 0.0472 0.0603 0.0576 0.033 -0.0073 -0.0777 -0.0608 -0.0438 -0.0209 0.0031 0.035 0.0293 0.0121 0.0338 0.03170.0206 0.0198 0.0174 0.0021 -0.0144 -0.0343 -0.0339 -0.0145 -0.0028 0.017 -0.0096 -0.0255 -0.0279 -0.0388 -0.0242 -0.0215 -0.0182 -0.0174 -0.0038 -0.0027 -0.0185 -0.0123 0.0087 0.0343 0.0695 0.091 0.0853 0.076 0.0513 0.0186 0.0015 -0.019 -0.0151 -0.0073 0.0021 0.0129 0.0215 0.0024 -0.0124 -0.0329 -0.0519 -0.0708 -0.0579-0.0307 -0.0145 -0.0009 -0.018 -0.0318 -0.0465 -0.0391 -0.0345 -0.0316 -0.0435 -0.0491 -0.0475 -0.042 -0.0361 -0.0277 -0.0258 -0.0139 -0.0068 0.0507 0.0722 0.0878 0.0782 0.0765 0.0439 0.008 0.0013 -0.0126 -0.0015 0.003 0.0104 0.0104 0.0193 0.0205 0.0074 -0.0056 -0.0072 0.007 0.0106 0.0147 -0.0009 -0.0159 -0.0187 -0.00070.0105 -0.0115 -0.0302 -0.0309 -0.0095 -0.0058 0.0004 0.002 0.005 0.0057 0.0097 0.0134 0.0177 0.0218 0.0261 0.0302 0.0346 0.0386 0.0474 0.0393 0.0238 0.0115 -0.0079 -0.0124 0.0054 0.0027 -0.025 -0.0566 -0.063 -0.0591 -0.0413 -0.0068 0.0272 0.0277 -0.0021 -0.006 -0.011 -0.0221 -0.0416 -0.0519 -0.0222 0.003 0.00790.0171 0.0253 0.0323 0.0391 0.0164 -0.0136 -0.0323 -0.0291 -0.0287 -0.0304 -0.0339 -0.0245 -0.0076 0.0125 0.0376 0.0402 0.0245 0.0156 -0.004 -0.0153 -0.0289 -0.0316 -0.0111 0.0094 0.0335 0.0576 0.0424 0.0143 -0.0007 -0.0135 -0.027 -0.0341 -0.0357 -0.0396 -0.0402 -0.0488 -0.048 -0.0406 -0.0407 -0.0351 -0.0187 -0.0057 0.0044-0.0072-0.0169-0.01150.0126finish/clear12.悬臂梁地震分析!定义单元类型/prep7ET,1,BEAM3R,1,0.09,0.000675,0.3, !0.3*0.3m截面MP,EX,1,2.06e11MP,NUXY,1,0.3MP,DAMP,1,0.02 !阻尼比0.02MP,dens,1,7.85e3!--------------------------------------------------!定义节点单元和位移约束*do,i,1,11n,i,(i-1)*0.5,0,0 !5m长,分10段*enddo*do,i,1,10e,i,i+1*enddod,1,ALL,0!--------------------------------------------------!定义和读入时程曲线NT=1000 !时程曲线有NT个点DT=0.01 !时间间隔*dim,ac,,NT !读入数据,这个数据文件可以用excel等软件来写/input,tianjin,txt !天津波东西向共10秒!--------------------------------------------------/SOLUNSUBST,1, , ,1 !1个子步OUTRES,ALL,1 !输出每个子步的结果ANTYPE,TRANS !时程分析*do,i,1,NTACEL,0,ac,0solve*enddo!--------------------------------------------------/POST26NSOL,2,11,U,Y,PLVAR,2,3, , ,下面是文件tianjin.txt，这个文件可以用excel等软件制作或编写很短的程序就能写出来先贴前500个点ac( 1)= -0.06334598ac( 2)= -0.04417088ac( 3)= -0.02188456ac( 4)= -0.00622243ac( 5)= 0.01599961ac( 6)= 0.03431334ac( 7)= 0.05332774ac( 8)= 0.07603510ac( 9)= 0.09322070ac( 10)= 0.11568701ac( 11)= 0.13706695ac( 12)= 0.15446468ac( 13)= 0.17966937ac( 14)= 0.16881546ac( 15)= 0.10252541ac( 16)= 0.03081323ac( 17)= -0.02835768ac( 18)= -0.04577148ac( 19)= 0.00371867ac( 20)= 0.05160179ac( 21)= 0.08542017ac( 22)= 0.14456214ac( 23)= 0.15453540ac( 24)= 0.07974091ac( 25)= -0.00708701ac( 26)= -0.08872422ac( 27)= -0.15415934ac( 28)= -0.16546641ac( 29)= -0.15415934ac( 30)= -0.15538390ac( 31)= -0.15533569ac( 32)= -0.15030567ac( 33)= -0.14931896ac( 34)= -0.14984284ac( 35)= -0.14375542 ac( 36)= -0.14322831 ac( 37)= -0.14130631 ac( 38)= -0.07721145 ac( 39)= 0.02666387 ac( 40)= 0.06092257 ac( 41)= 0.03247169 ac( 42)= 0.01391368 ac( 43)= -0.00049497 ac( 44)= -0.02224453 ac( 45)= -0.04070290 ac( 46)= -0.05431446 ac( 47)= -0.05096541 ac( 48)= -0.03585931 ac( 49)= -0.02533004 ac( 50)= -0.01427366 ac( 51)= -0.00180630 ac( 52)= 0.00995074 ac( 53)= 0.02199384 ac( 54)= 0.03315949 ac( 55)= 0.04807596 ac( 56)= 0.05553902 ac( 57)= 0.03931121 ac( 58)= -0.00056568 ac( 59)= -0.05451052 ac( 60)= -0.07338993 ac( 61)= -0.02214490 ac( 62)= 0.05202605 ac( 63)= 0.11738081 ac( 64)= 0.17217417 ac( 65)= 0.18113498 ac( 66)= 0.15603314 ac( 67)= 0.14236373 ac( 68)= 0.12819293 ac( 69)= 0.10743007 ac( 70)= 0.09416241 ac( 71)= 0.07652363 ac( 72)= 0.05719426 ac( 73)= 0.04319058 ac( 74)= 0.02332124 ac( 75)= 0.00492073 ac( 76)= -0.00940756 ac( 77)= -0.04927159 ac( 78)= -0.13574918。

ansys-土木工程应用实例命令流

Ansys-土木工程应用实例命令流/filname,huitan/units,si/prep7et,1,plane42keyopt,1,3,2mp,ex,1,2.160687961e11mp,nuxy,1,0.3TB,MISO,1,1,25,TBTEMP,0TBPT,,0.000814,175.88e6,TBPT,,0.001531,324.98e6,TBPT,,0.002819,591.26e6,TBPT,,0.003846,801.457e6,TBPT,,0.004239,873.22e6,TBPT,,0.004594,924.399e6,TBPT,,0.004955,956.88e6,TBPT,,0.00531,976.99e6,TBPT,,0.005699,990.2138e6,TBPT,,0.006261,1000e6,TBPT,,0.006957,1008.34e6,TBPT,,0.007804,1015.5e6,TBPT,,0.008779,1021.185e6,TBPT,,0.009827,1026.47e6,TBPT,,0.010929,1031.289e6,TBPT,,0.012066,1035.444e6,TBPT,,0.0132,1039.693e6,TBPT,,0.014354,1044.886e6,TBPT,,0.015494,1049.796e6,TBPT,,0.016631,1053.951e6,TBPT,,0.017785,1058.106e6,TBPT,,0.019039,1062.827e6,TBPT,,0.020276,1067.643e6,csys,1k,1,0.00512k,2,0.0064kgen,36,1,2,1,,10*do,j,1,69,2a,j,j+1,j+3,j+2*enddoa,71,72,2,1savelesize,all,,,4amesh,allnummrg,allnumcmp,allsavecp,1,uy,166,170,171,172,173finish/solunsel,s,loc,y,270d,all,allnplotallselnplotnsel,s,loc,y,90d,all,uxnplotallseloutres,all,alltime,1NSUBST,100nlgeom,ond,166,uy,-0.00128solve*get,def1,node,166,U,Ynsel,s,loc,y,90ddele,all,uyallself,166,fy,0solve*get,def2,node,166,U,Yfinish/post1 !进入后处理器set,list !显示所有计算的结果set,1,last,1 !读取第一个加载步最后子步的结果pldisp,1 !同时显示变形前与变形后结构的形状plnsol,s,eqv,1,1 !显示节点的等效应力set,2,last,1 !读取第二个加载步最后子步的结果pldisp,1 !同时显示变形前与变形后结构的形状。

ansys 温度荷载实例分析命令流

!XXXXXXXXXXXXXXXXXX!定义单元类型!XXXXXXXXXXXXXXXXXXFINI/CLE/TITLE,QSHang Pylon Test/UNITS,SICSYS,0!*afun,deg/PREP7/triad,on/view,1,1,1,1 !定义ISO查看/NERR,1000000et,1,plane55 !定义二维热单元et,2,surf153mp,kxx,1,3.3mp,kyy,1,3.3 !设定砼的导热系数mp,kzz,1,3.3!mp,dens,1,2400 !定义砼的密度mp,c,1,925 !定义砼的比热mp,prxy,1,0.2 !泊松比C50mp,ex,1,3.45e10 !混凝土的弹性模量C50 mp,alpx,1,1e-5!mptemp,1,0,2643,2750,2875 !定义铸钢的热性能!mpdata,kxx,2,1.44,1.54,1.22,1.22!mpdata,enth,2,0,128.1,163.8,174.2!mpplot,kxx,2!mpplot,enth,2!XXXXXXXXXXXXXXXXXXX!建立混凝土几何模型!XXXXXXXXXXXXXXXXXXX!************plane1K,1,0,0,0k,2,0,3.5,0k,3,6.885,3.5k,4,6.885,4k,5,21.2,4k,6,23.7,0l,1,2l,2,3l,3,4l,4,5local,11,1,21.2,0,0,,,,4/2.5,1 l,5,6CSYS,0k,7,18.9,0k,8,18.9,2k,9,19.4,2.5k,10,21,2.5k,11,22.2,0l,1,7l,7,8l,8,9l,9,10local,12,1,21,0,0,,,,2.5/1.2,1 l,10,11CSYS,0l,11,6al,all!************plane2 asel,nonelsel,nonek,12,0.8,0.25k,13,0.3,0.75k,14,0.3,2k,15,0.8,2.5k,16,5.2,2.5k,17,5.7,2k,18,5.7,0.75k,19,5.2,0.25a,12,13,14,15,16,17,18,19k,20,6.8,0.25k,21,6.3,0.75k,22,6.3,2k,23,6.8,2.5k,24,11.2,2.5k,25,11.7,2k,26,11.7,0.75k,27,11.2,0.25a,20,21,22,23,24,25,26,27k,28,12.8,0.25k,29,12.3,0.75k,30,12.3,2k,31,12.8,2.5k,32,15.3,2.5k,33,15.8,2k,34,15.8,0.75k,35,15.3,0.25a,28,29,30,31,32,33,34,35asel,scm,plane,areaallsasba,1,plane!************画网格aatt,1,1,AESIZE,all,0.02 !将所有面网格划分时的单元尺寸设置为0.02m MSHAPE,0,2D !利用四边形单元进行网格划分MSHKEY,2 !采用映射网格划分单元AMESH,ALL !对所有网格进行划分单元allscm,area1,areaarsym,x,area1,,,,0,0allsarsym,y,all,,,,0,0asel,allaglue,all!Nummrg,elem!Nummrg,node!Nummrg,kpnummrg,allnumcmp,all!**************定义组件allslsel,s,length,,6.885lsel,a,length,,0.5lsel,a,length,,14.32lsel,r,loc,y,0,100cm,Eligne,linensll,scm,Enode,nodeTYPE,2ESURFallslsel,s,length,,6.885lsel,a,length,,0.5lsel,a,length,,14.32lsel,r,loc,y,-100,0cm,Wligne,linensll,scm,Wnode,nodeTYPE,2ESURFallslsel,s,length,,5.173lsel,r,loc,x,0,100cm,Sligne,linensll,scm,Snode,nodeTYPE,2ESURFallslsel,s,length,,5.173lsel,r,loc,x,-100,0cm,Nligne,linensll,scm,Nnode,nodeALLSTYPE,2ESURF!**************计算日照时间内，即14个小时内的空气的温度值*DIM,t_outdoor,,13 !定义室外空气温度数组为t_outdoor*DO,t,6,18 !给室外空气温度数组t_outdoor赋值*SET,t_outdoor(t-5),31.5+5.5*sin(3.1415926535898*(t-9)/12)*ENDDO!**************日照强度*DIM,E_insolation,,13 !东边日照*DIM,S_insolation,,13 !东边日照*DIM,W_insolation,,13 !西边日照*DIM,N_insolation,,13 !北边日照E_insolation(1) =240 $S_insolation(1) =33 $ W_insolation(1) = 28 $ N_insolation(1) =89E_insolation(2) =521 $S_insolation(2) =80 $W_insolation(2) =69 $ N_insolation(2) = 160E_insolation(3) =629 $S_insolation(3) = 115 $W_insolation(3) =99 $ N_insolation(3) = 166E_insolation(4) =640 $S_insolation(4) =148 $W_insolation(4) = 127 $ N_insolation(4) = 152E_insolation(5) =536 $S_insolation(5) =184 $W_insolation(5) = 142 $ N_insolation(5) = 142E_insolation(6) =374 $S_insolation(6) =229 $W_insolation(6) = 153 $ N_insolation(6) = 153E_insolation(7) =178 $S_insolation(7) =240 $W_insolation(7) = 153 $ N_insolation(7) = 153E_insolation(8) =178 $S_insolation(8) =229 $W_insolation(8) = 322 $ N_insolation(8) = 153E_insolation(9) =165 $S_insolation(9) =184 $W_insolation(9) = 461 $ N_insolation(9) = 142E_insolation(10) =148 $S_insolation(10) =148 $W_insolation(10) =550 $ N_insolation(10) = 152E_insolation(11)=115 $S_insolation(11)=115 $W_insolation(11) =541 $ N_insolation(11) =166E_insolation(12)=80 $S_insolation(12)=80 $W_insolation(12)= 448 $ N_insolation(12) = 160E_insolation(13)=33 $S_insolation(13)=33 $W_insolation(13)= 206 $ N_insolation(13) = 89!**************综合换热系数h=12.47+3.33*1.5!**************综合大气温度*DIM,Et_synthetiser,,13*DIM,St_synthetiser,,13*DIM,Wt_synthetiser,,13*DIM,Nt_synthetiser,,13*DO,x,1,13,1*SET,Et_synthetiser(x),t_outdoor(x)+(E_insolation(x)/17.465)*0.65*SET,St_synthetiser(x),t_outdoor(x)+(S_insolation(x)/17.465)*0.65*SET,Wt_synthetiser(x),t_outdoor(x)+(W_insolation(x)/17.465)*0.65*SET,Nt_synthetiser(x),t_outdoor(x)+(N_insolation(x)/17.465)*0.65*ENDDO!**************求解/soluantype,trans !设置求解类型瞬态热分析trnopt,full ! 指定瞬态分析的求解方法为完全法timint,on ! 打开时间积分效应tunif,27.6 !设置构件的初始温度为20摄氏度!tref,20outres,all ! 将除SV AR和LOCI以外的所有计算数据写入数据库和文件中!antype,4!autots,on!trnopt,full!lnsrch,on!outres,basic,last!outpr,basic,lastl=0*do,m,1,13,1*do,r,3600,3600,3600time,(m-1)*3600+rl=l+1nsubst,1,100,1 !设置每个荷载的子步数数为1，最大值为100，最小值为1 autots,on ! 打开自动时间步长跟踪eqslv,JCG !指定方程求解器为JCGkbc,0 !使用递增方式加载线性荷载!施加对流荷载!alls!nsel,s,,,Enode!d,all,temp,Et_synthetiser(m)!alls!nsel,s,,,Wnode!d,ALL,temp,Wt_synthetiser(m)!alls!nsel,s,,,Snode!d,ALL,temp,St_synthetiser(m)!alls!nsel,s,,,Nnode!d,ALL,temp,Nt_synthetiser(m)!bf,Enode,hgen,E_insolation(m)!bf,Wnode,hgen,W_insolation(m)!bf,Snode,hgen,S_insolation(m)!bf,Nnode,hgen,N_insolation(m)allsSFl,Eligne,conv,17.456,,Et_synthetiser(m) SFL,Wligne,conv,17.456,,Wt_synthetiser(m) SFL,Sligne,conv,17.456,,St_synthetiser(m) SFL,Nligne,conv,17.456,,Nt_synthetiser(m)allslsel,s,,,Eligneesll,sesel,r,type,,1sfe,all,,HFLUX,,0.65*E_insolation(m)allslsel,s,,,Wligneesll,sesel,r,type,,1sfe,all,,HFLUX,,0.65*W_insolation(m)allslsel,s,,,Sligneesll,sesel,r,type,,1sfe,all,,HFLUX,,0.65*S_insolation(m)allslsel,s,,,Nligneesll,sesel,r,type,,1sfe,all,,HFLUX,,0.65*N_insolation(m)!alls!nsel,s,,,Enode!esln,s!esel,r,type,,1!sfe,all,,conv,,0.65*E_insolation(m)!alls!nsel,s,,,Wnode!esln,s!esel,r,type,,2!sfe,all,,conv,,0.65*W_insolation(m)!alls!nsel,s,,,Snode!esln,s!esel,r,type,,2!sfe,all,,conv,,0.65*S_insolation(m)!alls!nsel,s,,,Nnode!esln,s!esel,r,type,,2!sfe,all,,conv,,0.65*N_insolation(m)!sfa,5,,conv,10,t_outdoor(m) !给面5定义当前荷载步的对流换热系数及周围环境温度!sfa,21,,conv,10,t_outdoor(m) !给面21定义当前荷载步的对流换热系数及周围环境温度!sfa,18,,conv,10,t_outdoor(m) !给面18定义当前荷载步的对流换热系数及周围环境温度!sfa,10,,conv,10,t_outdoor(m) !给面10定义当前荷载步的对流换热系数及周围环境温度!施加太阳辐射荷载!bfv,1,HGEN,heat_eq(1,m)+long_wave(1) !给体1施加当前荷载步的等效生热率!bfv,4,HGEN,heat_eq(2,m)+long_wave(2) !给体4施加当前荷载步的等效生热率!bfv,3,HGEN,heat_eq(3,m)+long_wave(3) !给体3施加当前荷载步的等效生热率!bfv,2,HGEN,heat_eq(4,m)+long_wave(4) !给体2施加当前荷载步的等效生热率!bfe,all,hgen,1,1e10allssolve*enddo*enddo!*************************进入热应力求解/PREP7ETCHG,TTS !热到结构分析转换/SOLantype,4 !瞬态timint,1,struct !Turns on transient effectstimint,0,thermtimint,0,magtimint,0,elecautots,on !打开自动时间步长，有利于非线性收敛tref,20trnopt,full !完全瞬态nropt,full !设置牛顿-拉普森选项kbc,0 !渐变荷载nlgeom,on !打开大变形效应LUMPM,0PSTRES,ONNROPT,FULL, ,on!定义热应力计算参考温度TREF,20,!*************************边界条件DSYM,SYMM,Z,0csys,0LSEL,ALLDL,ALL,,UX,0DL,ALL,,UY,0DL,ALL,,UZ,0*DO,I,1,1355LDREAD,TEMP,,,1*i, ,'model_BD1','rth',' ' !读入热分析的计算结果OUTRES,ALL,ALL,TIME,1*iDELTIM,1,1,1SOLVE*ENDDO*DO,I,1,500LDREAD,TEMP,,,1355+10*i, ,'model_BD1','rth',' ' !读入热分析的计算结果OUTRES,ALL,ALL,TIME,1355+10*iDELTIM,10,10,10SOLVE*ENDDO*DO,I,1,300LDREAD,TEMP,,,6355+50*i, ,'model_BD1','rth',' ' !读入热分析的计算结果OUTRES,ALL,ALL,TIME,6355+50*iDELTIM,50,50,50SOLVE*ENDDO*DO,I,1,300LDREAD,TEMP,,,21355+100*i, ,'model_BD1','rth',' ' !读入热分析的计算结果OUTRES,ALL,ALL,TIME,21355+100*iDELTIM,100,100,100 SOLVE*ENDDOSET,LASTCSYS,0PATH,zdirec,2,50PPA TH,1,,0,0,0.17PPA TH,2,,0,0.7,0.17 PDEF,,S,xPLPATH,Sx/IMAGE,SA VE,zpath1,BMPSET,LASTCSYS,0PATH,zdirec,2,50PPA TH,1,,-1.72,0.4,0.17 PPA TH,2,,1,0.4,0.17 PDEF,,S,xPLPATH,Sx/IMAGE,SA VE,zpath2,BMPSET,LASTCSYS,0PATH,zdirec,2,50PPA TH,1,,0,0.7,0PPA TH,2,,0,0.7,0.35 PDEF,,S,zPLPATH,Sz/IMAGE,SA VE,zpath3,BMPSET,LASTCSYS,0PATH,zdirec,2,50PPA TH,1,,-1.72,0.7,0.17 PPA TH,2,,1,0.7,0.17 PDEF,,S,zPLPATH,Sz/IMAGE,SA VE,zpath4,BMPvsel,s,,,1 !选择编号为1的体元素eslv,r !选择当前所选体元素上的所有单元nsle,r !选择当前所选单元上的所有节点*get,node_total,node,0,count ! 提取当前所选节点的个数*get,node_num1(1),node,0,num,min ! 提取当前所选节点中节点的最小编号，并赋值给变量node_num1(1)*do,n,2,node_total,1*get,node_num1(n),node,node_num1(n-1),nxth ! 将当前所选择的节点的编号赋值给数组node_num1*enddo*do,f,1,node_total,1*get,node_temp1(f),node,node_num1(f),temp ! 提取当前所选节点的节点温度值，并赋值给数组node_temp1*enddototal=0*do,f,1,node_total,1total=total+node_temp1(f) ! 计算当前所选节点的温度值之和，并赋值给变量total*enddotarea=total/node_total+273 ! 计算当前所选节点温度值的平均值，并赋值给tareatsky=t_outdoor(m)-6+273 !计算当前时刻的天空温度值，并赋值给tsky*if,fuhao(1),eq,100,thenvalue_4=value_3*((tsky**4-tarea**4)*(1+cos_ph(1))/2) !计算与天空之间值长波辐射净强度*elsevalue_4=value_3*((t_earth(m)**4-tarea**4)*(1-cos_ph(1))/2) !计算与地面之间值长波辐射净强度*endif*set,long_wave(1,l),value_4/0.03 !计算长波辐射净强度的等效生热率allselvsel,s,,,4 !选择编号为4的体元素eslv,r !选择当前所选体元素上的所有单元nsle,r !选择当前所选单元上的所有节点*get,node_total,node,0,count ! 提取当前所选节点的个数*get,node_num1(1),node,0,num,min ! 提取当前所选节点中节点的最小编号，并赋值给变量node_num1(1)*do,n,2,node_total,1*get,node_num1(n),node,node_num1(n-1),nxth ! 将当前所选择的节点的编号赋值给数组node_num1*enddo*do,f,1,node_total,1*get,node_temp1(f),node,node_num1(f),temp ! 提取当前所选节点的节点温度值，并赋值给数组node_temp1*enddototal=0*do,f,1,node_total,1total=total+node_temp1(f) ! 计算当前所选节点的温度值之和，并赋值给变量total*enddotarea=total/node_total+273 ! 计算当前所选节点温度值的平均值，并赋值给tareatsky=t_outdoor(m)-6+273 !计算当前时刻的天空温度值，并赋值给tsky*if,fuhao(2),eq,100,thenvalue_4=value_3*((tsky**4-tarea**4)*(1+cos_ph(2))/2) !计算与天空之间值长波辐射净强度*elsevalue_4=value_3*((t_earth(m)**4-tarea**4)*(1-cos_ph(2))/2) !计算与地面之间值长波辐射净强度*endif*set,long_wave(2,l),value_4/0.03 !计算长波辐射净强度的等效生热率allselvsel,s,,,3 !选择编号为3的体元素eslv,r !选择当前所选体元素上的所有单元nsle,r !选择当前所选单元上的所有节点*get,node_total,node,0,count ! 提取当前所选节点的个数*get,node_num1(1),node,0,num,min ! 提取当前所选节点中节点的最小编号，并赋值给变量node_num1(1)*do,n,2,node_total,1*get,node_num1(n),node,node_num1(n-1),nxth ! 将当前所选择的节点的编号赋值给数组node_num1*enddo*do,f,1,node_total,1*get,node_temp1(f),node,node_num1(f),temp ! 提取当前所选节点的节点温度值，并赋值给数组node_temp1*enddototal=0*do,f,1,node_total,1total=total+node_temp1(f) ! 计算当前所选节点的温度值之和，并赋值给变量total*enddotarea=total/node_total+273 ! 计算当前所选节点温度值的平均值，并赋值给tareatsky=t_outdoor(m)-6+273 !计算当前时刻的天空温度值，并赋值给tsky*if,fuhao(3),eq,100,thenvalue_4=value_3*((tsky**4-tarea**4)*(1+cos_ph(3))/2) !计算与天空之间值长波辐射净强度*elsevalue_4=value_3*((t_earth(m)**4-tarea**4)*(1-cos_ph(3))/2) !计算与地面之间值长波辐射净强度*endif*set,long_wave(3,l),value_4/0.03 !计算长波辐射净强度的等效生热率allselvsel,s,,,2 !选择编号为2的体元素eslv,r !选择当前所选体元素上的所有单元nsle,r !选择当前所选单元上的所有节点*get,node_total,node,0,count ! 提取当前所选节点的个数*get,node_num1(1),node,0,num,min ! 提取当前所选节点中节点的最小编号，并赋值给变量node_num1(1)*do,n,2,node_total,1*get,node_num1(n),node,node_num1(n-1),nxth ! 将当前所选择的节点的编号赋值给数组node_num1*enddo*do,f,1,node_total,1*get,node_temp1(f),node,node_num1(f),temp ! 提取当前所选节点的节点温度值，并赋值给数组node_temp1*enddototal=0*do,f,1,node_total,1total=total+node_temp1(f) ! 计算当前所选节点的温度值之和，并赋值给变量total*enddotarea=total/node_total+273 ! 计算当前所选节点温度值的平均值，并赋值给tareatsky=t_outdoor(m)-6+273 !计算当前时刻的天空温度值，并赋值给tsky*if,fuhao(4),eq,100,thenvalue_4=value_3*((tsky**4-tarea**4)*(1+cos_ph(4))/2) !计算与天空之间值长波辐射净强度*elsevalue_4=value_3*((t_earth(m)**4-tarea**4)*(1-cos_ph(4))/2) !计算与地面之间值长波辐射净强度*endif*set,long_wave(4,l),value_4/0.03 !计算长波辐射净强度的等效生热率allsel*enddo*enddo。

一个经典的ansys热分析实例(流程序)

/PREP7/TITLE,Steady-state thermal analysis of pipe junction/UNITS,BIN ! 英制单位；Use U. S. Customary system of units (inches)! /SHOW, ! Specify graphics driver for interactive runET,1,90 ! Define 20-node, 3-D thermal solid elementMP,DENS,1,.285 ! Density = .285 lbf/in^3MPTEMP,,70,200,300,400,500 ! Create temperature tableMPDATA,KXX,1,,8.35/12,8.90/12,9.35/12,9.80/12,10.23/12! 指定与温度相对应的数据材料属性；导热系数；Define conductivity valuesMPDATA,C,1,,.113,.117,.119,.122,.125! Define specific heat values（比热）MPDATA,HF,2,,426/144,405/144,352/144,275/144,221/144! Define film coefficient;除144是单位问题，上面的除12也是单元问题! Define parameters for model generationRI1=1.3 ! Inside radius of cylindrical tankRO1=1.5 ! Outside radiusZ1=2 ! LengthRI2=.4 ! Inside radius of pipeRO2=.5 ! Outside pipe radiusZ2=2 ! Pipe lengthCYLIND,RI1,RO1,,Z1,,90 ! 90 degree cylindrical volume for tankWPROTA,0,-90 ! 旋转当前工作的平面；从Y到Z旋转-90度；；Rotate working plane to pipe axisCYLIND,RI2,RO2,,Z2,-90 ! 角度选择在了第四象限；90 degree cylindrical volume for pipeWPSTYL,DEFA ! 重新安排工作平面的设置；另外WPSTYL,STAT to list the status of the working plane；；Return working plane to default settingBOPT,NUMB,OFF ! 关掉布尔操作的数字警告信息；Turn off Boolean numbering warningVOVLAP,1,2 ! 交迭体；Overlap the two cylinders/PNUM,VOLU,1 ! 体编号打开；Turn volume numbers on/VIEW,,-3,-1,1/TYPE,,4 ! 精确面的显示；Precise hidden display/TITLE,Volumes used in building pipe/tank junctionVPLOTVDELE,3,4,,1 ! 修剪一些体与体相关的体的因素都删掉；Trim off excess volumes! Meshing 网格划分ASEL,,LOC,Z,Z1 ! Select area at remote Z edge of tankASEL,A,LOC,Y,0 ! Select area at remote Y edge of tankCM,AREMOTE,AREA ! 为面建立数组；Create area component called AREMOTE/PNUM,AREA,1/PNUM,LINE,1/TITLE,Lines showing the portion being modeledAPLOT/NOERASE ! 预防抹去LPLOT ! Overlay line plot on area plot/ERASEACCAT,ALL ! 连接面和线的准备映射；Concatenate areas and lines at remote tank edgesLCCAT,12,7LCCAT,10,5LESIZE,20,,,4 ! 4 divisions through pipe thicknessLESIZE,40,,,6 ! 6 divisions along pipe lengthLESIZE,6,,,4 ! 4 divisions through tank thicknessALLSEL ! Restore full set of entitiesESIZE,.4 ! Set default element size线的默认划分数MSHAPE,0,3D ! Choose mapped brick meshMSHKEY,1 ! 映射网格SAVE ! Save database before meshingVMESH,ALL ! Generate nodes and elements within volumes/PNUM,DEFA ! 重新安排数字规格/TITLE,Elements in portion being modeledEPLOTFINISH/COM, *** Obtain solution ***/SOLUANTYPE,STATIC ! Steady-state analysis typeNROPT,AUTO ! 自动选择牛顿-拉普森Program-chosenNewton-Raphson optionTUNIF,450 ! 给结点统一的温度；Uniform starting temperature at all nodesCSYS,1 ! 1 —Cylindrical with Z as the axis of rotation NSEL,S,LOC,X,RI1 ! Nodes on inner tank surfaceSF,ALL,CONV,250/144,450 ! 为结点指定表面载荷；对流；Convection（对流）；load at all nodesCMSEL,,AREMOTE ! 选择子集组合；Select AREMOTE component NSLA,,1 ! Nodes belonging to AREMOTED,ALL,TEMP,450 ! 设定边界温度条件Temperature constraints at those nodesWPROTA,0,-90 ! Rotate working plane to pipe axisCSWPLA,11,1 ! 在工作区声明本地的圆柱体系；Define local cylindrical c.s at working planeNSEL,S,LOC,X,RI2 ! Nodes on inner pipe surfaceSF,ALL,CONV,-2,100 ! 这里的-2表示材料2；；Temperature-dep. convection load at those nodesALLSEL/PBC,TEMP,,1 ! 边界符号的显示Temperature b.c. symbols on/PSF,CONV,,2 ! Convection symbols on 箭头显示/TITLE,Boundary conditionsNPLOTWPSTYL,DEFACSYS,0AUTOTS,ON ! Automatic time steppingNSUBST,50 ! Number of substepsKBC,0 ! Ramped loading (default)OUTPR,NSOL,LAST ! 显示最后一次的结点约束；Optional command for solution printoutSOLVEFINISH/COM, *** Review results ***/POST1/EDGE,,1 ! Displays only the "edges（刀口, 利刃, 锋, 优势, 边缘, 优势, 尖锐）" of an object；Edge display/PLOPTS,INFO,ON ! Legend column on/PLOPTS,LEG1,OFF ! Legend header off 圆柱数列的头部/WINDOW,1,SQUARE ! SQUA, form largest square window within the current graphics area；Redefine window size/TITLE,Temperature contours at pipe/tank junctionPLNSOL,TEMP ! Plot temperature contoursCSYS,11NSEL,,LOC,X,RO2 ! Nodes and elements at outer radius of pipeESLN ! 选择单元NSLE ! 选择结点/SHOW,,,1 ! 向量显示；Vector mode/TITLE,Thermal flux vectors at pipe/tank junctionPLVECT,TF ! Plot thermal flux（热通量）vectorsFINISH。