加载求解及后处理技术

土木工程中的有限元辽宁科技大学资源与土木工程学院于新第4章加载求解及后处理技术4.1荷载及其施加4.2荷载步选项及设置4.3分析类型与求解控制选项4.4求解代价估计4.5重启动分析4.6时变结构的多荷载步求解4.7通用后处理4.8时间历程后处理4.1荷载及其施加---荷载在ANSYS中荷载包括边界条件和作用力，对结构分析如下：位移、力、压力、温度、重力一般可将荷载分为六类，如表4-1所示。4.1荷载及其施加---荷载荷载即可施加在几何模型（关键点、硬点、线、面、体）上，也可施加在有限元模型（节点、单元）上，或者二者混合使用。施加在几何模型上的荷载独立于有限元网格，不必为修改网格而重新加载；施加在有限元模型上且要修改网格，则必须先删除荷载再修改网格，然后重新施加荷载。不管施加到何种模型上，在求解时荷载全部转换（自动或人工）到有限元模型上。4.1荷载及其施加---施加自由度约束在结构分析中自由度共有7个，自由度的方向均依从节点坐标系。约束可施加在节点、关键点、线和面上。4.1荷载及其施加---施加自由度约束1.节点自由度约束及相关命令⑴对节点施加自由度约束命令：D,NODE,Lab,VALUE,VALUE2,NEND,NINC,Lab2,Lab3,Lab4,Lab5,Lab6NODE---拟施加约束的节点号，其值可取ALL、组件名。Lab---自由度标识符，如UX、ROTZ等。如为ALL，则为所有适宜的自由度。VALUE---自由度约束位移值或表式边界条件的表格名称。VALUE2---约束位移值的第二个数，如为复数输入时，VALUE为实部，而VALUE2为虚部。NEND,NINC---节点编号范围和编号增量，缺省时NEND=NODE，NINC=1。Lab2,Lab3,Lab4,Lab5,Lab6---其它自由度标识符，VALUE对这些自由度也有效。各自由度的方向用节点坐标系确定，转角约束位移用弧度输入例如：D,ALL,ALL!对所选节点的全部自由度施加约束D,18,UX,,,,,UY,UZ!对节点18的3个平动自由度全部施加约束D,20,UX,1.0e-4!对节点20的UX施加约束，且约束位移值为1.0e-4D,22,UX,0.1,,25,,UY,ROTY!对节点22～25的UX,UY,ROTY施加约束，且位移值均为0.14.1荷载及其施加---施加自由度约束⑷在节点上施加对称和反对称约束命令：DSYM,Lab,Normal,KCNLab---对称标识，如为SYMM则生成对称约束，如为ASYM则生成反对称约束。Normal---约束的表面方向标识，一般垂直于参数KCN坐标系中的坐标方向。其值有：=X（缺省）：表面垂直于X方向，非直角坐标系为R方向；=Y：表面垂直于Y方向，非直角坐标系为θ方向；=Z：表面垂直于Z方向，球和环坐标系为Φ方向；KCN---用于定义表面方向的整体或局部坐标系的参考号。4.1荷载及其施加---施加自由度约束2.关键点自由度约束及相关命令命令：DK,KPOI,Lab,VALUE,VALUE2,KEXPND,Lab2,Lab3,Lab4,Lab5,Lab6KPOI---关键点编号，也可取ALL或元件名。KEXPND---扩展控制参数。如为0则仅施加约束到关键点上的节点；如为1则扩展到关键点之间（两关键点所连线）的所有节点上，且包括关键点上的节点，当然约束位移值相同。其余参数同D命令中的参数。列表和删除关键点自由度约束的命令分别为：列表：DKLIST,KPOI删除：DKDELE,KPOI,Lab例如：DK,ALL,ALL!约束所选择全部关键点的全部自由度DK,1,UY!对关键点1施加UY自由度约束DK,2,UX,0.01,,,UY,ROTZ!对关键点2的UX,UY,ROTZ施加约束，且位移值均为0.014.1荷载及其施加---施加自由度约束3.对线施加自由度约束命令：DL,LINE,AREA,Lab,Value1,Value2LINE---线编号，也可为ALL（缺省）或元件名。AREA---包含该线的面编号，并假定对称与反对称面垂直于该面，且线位于对称或反对称面内，缺省为包含该线的所选择面中的最小编号。如不是对称或反对称约束，则此面号无意义。Lab---自由度标识符，其值可取：=SYMM：对称约束，按DSYM命令的方式生成；=ASYM：反对称约束，按DSYM命令的方式生成；=UX,UY,UZ,ROTX,ROTY,ROTZ,WRAP：各自由度约束；=ALL：所有适宜的自由度约束（与单元相关）。Value1---自由度约束位移值或表格边界条件的表格名称。表格边界条件仅对UX、UY、UZ、ROTX、ROTY、ROTZ有效，且Value1=%tabname%。Value2---仅对FLOTRAN分析时有用，对结构分析无意义。该命令对线上的所有节点施加自由度约束。而列表和删除线上自由度约束的命令分别为：列表：DLLIST,LINE删除：DLDELE,LINE,Lab4.1荷载及其施加---施加自由度约束!EX4.2对线施加约束并转换finish$/clear$/prep7et,1,95$blc4,,,10,10,10!定义单元类型、创建长方体dl,7,,ux,0.1!线7施加UX自由度约束，位移值为0.1dl,5,,all!线5施加全部自由度约束dl,11,6,symm!线11施加对称约束，面号为6dl,10,6,asym!线10施加反对称约束，面号为6dl,6,,symm!线6施加对称约束，面号缺省DLLIST!列表显式线约束信息esize,2$vmesh,all!划分单元dtran$DLIST!转换约束到有限元模型，并列表显示4.1荷载及其其施加---施加自由由度约束束4.对面施加加自由度度约束命令：DA,AREA,Lab,Value1,Value2其中AREA为拟施加加约束的的面号，，也可为为ALL或元件名名，其余余同DL命令中的的参数。。该命令对对面上的的所有节节点施加加自由度度约束。。列表和删删除面上上自由度度约束的的命令分分别为：：列表：DALIST,AREA删除：DADELE,AREA,Lab5.约束转换换命令仅转换约约束自由由度命令令：DTRAN边界条件件和荷载载转换命命令：SBCTRAN这两命令令将几何何模型施施加的约约束和荷荷载转换换到有限限元模型型上。也也可不执执行这两两个命令令而在求求解时由由系统自自动转换换。4.1荷载及其其施加---施加自由由度约束束6.自由度约约束的冲冲突对于DK、DL和DA命令施加加的自由由度约束束参数可可能会发发生冲突突，例如如：DL指定会与与相邻线线（有公公共关键键点）上上的DL指定冲突突；DL指定会与与任一关关键点上上的DK指定冲突突；DA指定会与与相邻面面（有公公共关键键点和公公共线））上的DA指定冲突突；DA指定会与与任一线线上的DL指定冲突突；DA指定会与与任一关关键点上上的DK指定冲突突。按下列顺顺序将施施加到几几何模型型上的自自由度约约束转换换到有限限元模型型上：①按面号号增加的的顺序，，将DA的自由度度约束转转换到面面上的所所有节点点；②按面号号增加的的顺序，，将DA约束的SYMM和ASYM转换到面面上的所所有节点点；③按线号号增加的的顺序，，将DL自由度约约束转换换到线上上的所有有节点；；④按线号号增加的的顺序，，将DL的SYMM和ASYM约束转换换到线上上的所有有节点；；⑤将DK自由度约约束转换换到关键键点上的的所有节节点。所以，对对冲突的的约束，，DK命令改写写DL命令，DL命令改写写DA命令，施施加在较较大编号号图素上上的约束束改写较较低编号号上的约约束。这种冲突突的处理理与命令令执行的的前后顺顺序没有有关系，，但当发发生冲突突时，系系统会发发出警告告信息。。4.1荷载及其其施加---施加集中中荷载结构分析析中的集集中荷载载及其标标识符为为力FX,FY,FZ及力矩MX,MY,MZ。4.1荷载及其其施加---施加集中中荷载1.施加节点点集中荷荷载命令：F,NODE,Lab,VALUE,VALUE2,NEND,NINCNODE---节点编号号，也可可为ALL或元件名名。Lab---集中荷载载标识符符，如FX,FY,FZ,MX,MY,MZ其中任一一。VALUE---集中荷载载值或表表式边界界条件的的表格名名称。VALUE2---集中荷载载值的第第二个数数，如为为复数输输入时，，VALUE为实部部，而而VALUE2为虚部部。NEND,NINC---节点编编号范范围和和编号号增量量。节点集集中荷荷载列列表：：FLIST删除节节点集集中荷荷载：：FDELE4.1荷载及及其施施加---施加集集中荷荷载2.施加关关键点点集中中荷载载命令：：FK,KPOI,Lab,VALUE,VALUE2其中KPOI为关键键点号号，也也可取取ALL或元件件名。。其余余参数数同F命令。。FKLIST命令和和FKDELE命令。。转换命命令FTRAN仅将集集中荷荷载转转换到到有限限元模模型的的节点点上。。★不管管在何何种模模型上上施加加集中中荷载载，都都与节节点坐坐标系系相关关。★如果果尚没没有生生成有有限元元模型型，因因无节节点存存在，，对节节点坐坐标系系操作作无效效，所所施加加的荷荷载仅仅与总总体坐坐标系系相关关。★如果果几何何模型型和有有限元元模型型同时时存在在，则则节点点坐标标系的的设置置就有有效。。不管管是在在何时时何模模型上上施加加的荷荷载，，如果果节点点坐标标系重重新设设置了了，则则荷载载也跟跟着一一并改改变。。所以以在改改变节节点坐坐标系系时应应慎重重，以以避免免出现现错误误。4.1荷载及及其施施加---施加集集中荷荷载finish$/clear$/prep7et,1,beam4!定义单单元类类型k,1$k,2,5$k,3,10!创建3个关键键点l,1,2$l,2,3!创建2条线local,12,0,,,,90!设置12号局部部坐标标系，，其X12轴与总总体直直角坐坐标系系的Y轴相同同，!而其Y12轴与总总体坐坐标系系的X轴平行行，但但方向向相反反。nrotat,all!此时对对节点点坐标标系的的操作作无效效dk,1,all!关键点点1自由度度全部部约束束fk,2,fy,-1000!在当前前节点点坐标标系（（与总总体坐坐标系系相同同）于于关键键点2施加FY=-1000!其力的的作用用方向向与总总体直直角坐坐标系系的Y轴平行行。esize,1$lmesh,all!划分网网格，，生成成有限限元模模型nrotat,all!设置所所有节节点的的节点点坐标标系与与当前前激活活坐标标系相相同（（12号坐标标系））LPLOT!关键点点2上的FY=-1000方向与与Y12轴平行行，而而与总总体坐坐!标系的的X轴平行行了（（节点点坐标标系改改变了了，荷荷载跟跟着改改变））fk,3,fy,1000!在关键键点3施加FY=1000，方向向与Y12轴平行行f,6,fx,-1000!在节点点6施加FX=-1000，其方方向与与X12轴平行行sbctran!转换所所有边边界条条件到到有限限元模模型EPLOT!显示单单元与与边界界条件件4.1荷载及及其施施加---施加面面荷载载结构分分析中中的面面荷载载为压压力，，其标标识符符为PRES。虽然线线分布布荷载载和面面分布布荷载载都称称为压压力，，但对对不同同的单单元类类型，，其荷荷载单单位不不尽不不同。。对于2D面单元元，无无论面面荷载载施加加在单单元边边或边边界线线（LINE），其其荷载载单位位都是是“力力/面积””。对于SHELL单元，，施加加中面面法向向的面面荷载载单位位为““力/面积””，而而单元元边或或单元元边界界线上上的面面荷载载单位位为““力/长度””。对于梁梁单元元，其其分布布荷载载单位位为““力/长度””，单单元端端部荷荷载单单位为为“力力”。。对于3D实体单单元，，其面面荷载载的单单位为为“力力/面积””。4.1荷载及及其施施加---施加面面荷载载4.1荷载及及其施施加---施加面面荷载载1.施加节节点面面荷载载⑴对对节点点群施施加面面荷载载命令：：SF,Nlist,Lab,VALUE,VALUE2Nlist---节点群群，可可取ALL或元件件名。。Lab---面荷载载标识识符，，结构构分析析为PRES。VALUE---面荷载载值或或表格格型面面荷载载的表表格名名称。。VALUE2---复数输输入时时面荷荷载值值的第第二个个值。。★对于于单个个节点点不能能使用用该命命令。。★对于于3D体单元元面，，由Nlist节点群群能够够确定定多少少个单单元面面就施施加多多少单单元面面（与与几何何面无无关）），与与单元元是否否被单单独选选择无无关。。利用用该命命令可可以解解决大大面上上局部部加载载的问问题。。★对于于2D面单元元，当当在单单元外外部边边界（（不是是单元元边））上加加载时时，可可仅选选择外外部边边界上上的节节点群群即可可加载载；当当节点点群不不在单单元外外部边边界时时，尚尚须单单独选选择包包含这这些节节点的的单元元，否否则不不予施施加。。面荷荷载的的方向向与单单元面面平行行，且且指向向单元元面边边界。。该特特点对对于单单元周周边施施加相相同面面荷载载时比比较简简单，，当然然也可可施加加单元元任一一边的的面荷荷载，，但稍稍稍麻麻烦些些。4.1荷载及及其施施加---施加面面荷载载!EX4.4A3D单元SF加载示示例finish$/clear$/prep7et,1,95$blc4,,,10,10,20!定义单单元类类型，，创建建长方方体esize,,4$vmesh,all!定义单单元网网格数数目，，划分分单元元网格格asel,s,loc,y,10!选择Y=10的几何何面nsla,s!选择与与面相相关的的节点点，但但不包包含面面边界界节点点sf,all,pres,1000!施加节节点群群压力力荷载载（力力/面积）），仅仅4个单元元面asel,s,loc,z,20!选择Z=20的面nsla,s,1!选择与面相相关的所有有节点sf,all,pres,1000!施加节点群群压力荷载载（力/面积），所所有单元面面4.1荷载及其施施加---施加面荷载载4.1荷载及其施施加---施加面荷载载!EX4.4B2D单元SF加载示例FINISH$/CLEAR$/PREP7!①定义单元，，创建带孔孔面ET,1,82$BLC4,,,100,200blc4,30,60,40,80$asba,1,2wprota,,-90$wpoff,,,60$asbw,all!②切分面，以以便划分网网格wpoff,,,80$asbw,all$wprota,,,90wpoff,,,30$asbw,all$wpoff,,,40asbw,allwpcsys,-1$ESIZE,5$AMESH,ALL/psf,PRES,NORM,24.1荷载及其施施加---施加面荷载载SF,ALL,PRES,100!对所有单元元施加面荷荷载，即外外部边界加加载sfdele,all,pres!删除上述面面荷载4.1荷载及其施施加---施加面荷载载nsel,s,loc,x,0!选择X=0的节点群sf,all,pres,100!对上述节点点群施加面面荷载nsel,s,loc,x,15,20!选择X=15~20的节点esln,s,1!选择上述节节点能够确确定的全部部单元nsel,r,loc,x,15!从中选择X=15的节点群sf,all,pres,110!对上述节点点群施加面面荷载（内内部单元的的一边上））4.1荷载及其施施加---施加面荷载载nsel,s,loc,x,40,60!选择X=40~60的节点nsel,r,loc,y,10,30!从中选择Y=10~30的节点esln,s,1!选择上述节节点能够确确定的全部部单元sf,all,pres,100!对上述节点点群施加面面荷载（内内部单元的的周边上））lsel,s,loc,x,100!选择X=100的线nsll,s,1$esln,s!选择与线相相关的全部部节点，再再选择与节节点相关的的全部单元元nsel,s,loc,x,95!重新选择节节点群（在在上述单元元范围内））sf,all,pres,-100!对上述节点点群施加面面荷载（内内部单元的的一边上））allsel$eplot4.1荷载及其施施加---施加面荷载载4.1荷载及其施施加---施加面荷载载⑵定义节节点号与面面荷载的函函数关系命令：SFFUN,Lab,Par,Par2Lab---面荷载标识识符，结构构分析为PRES。Par---储存面荷载载值的参数数名（数组组参数）。。Par2---用于复数输输入时的第第二个值。。该命令定义义节点号与与面荷载的的函数关系系，数组中中值的位置置（数组下下标）表示示节点号，，数组值表表示面荷载载的大小。。该命令对对于施加由由其它软件件计算出的的节点面荷荷载时比较较有用，但但对于ANSYS自动生成的的有限元模模型，其节节点编号由由系统自动动确定，显显然要直接接应用这种种函数关系系并不容易易。该命令令所定义的的函数关系系，可用于于SF和SFE命令。4.1荷载及其施施加---施加面荷载载!EX4.5节点号及其其荷载函数数finish$/clear$/prep7et,1,45$blc4,,,10,10,20!定义单元类类型，创建建长方体esize,5$vmesh,all!定义单元尺尺寸，划分分网格*dim,mypres,,100!定义数组mppres*do,i,1,100$mypres(i)=i*10.0$*enddo!为数组赋值值sffun,pres,mypres(1)!定义节点号号与面荷载载函数关系系nsel,s,loc,y,10$sf,all,pres,10!选择节点群群，施加面面荷载SFLIST!该面荷载的的节点上的的值为10+i*10*DO,I,1,100$MYPRES(I)=I*50.0$*ENDDO!为数组重新新赋值，定定义另组关关系NSEL,S,LOC,z,20$SF,ALL,PRES,0!选择节点群群，并施加加面荷载ALLSEL$SFLIST!列表显示所所有面荷载载的值4.1荷载及其施施加---施加面荷载载⑶定义面面荷载梯度度命令：SFGRAD,Lab,SLKCN,Sldir,SLZER,SLOPELab---面荷载标识识符，结构构分析为PRES。SLKCN---斜率坐标系系的参考号号，缺省为为0（总体直角角坐标系））。Sldir---在SLKCN坐标系中梯梯度（或斜斜率）的方方向，其值值可取：=X（缺省）：：沿X方向的斜率率，对非直直角坐标系系为R方向；=Y：沿Y方向的斜率率，对非直直角坐标系系为θ方向；=Z：沿Z方向的斜率率，对球或或环坐标系系为φ方向；SLZER---斜率基值为为0的坐标位置置。如为角角度则单位位为度，如如果奇点在在180°则SLZER在±180°°之间，如果果奇点在0°，则SLZER在0°～360°之间。SLOPE---斜率值，即即单位长度度或单位角角度的荷载载值，沿Sldir正方向递增增为正，递递减为负。。该命令所定定义的梯度度（斜率））可为随后后的SF、SFE、SFL和SFA命令使用，，每个节点处处的荷载按按下式计算算：CVALUE=VALUE+(SLOPE×(COORD-SLZER))其中VALUE是命令SF、SFE、SFL和SFA中的参数值值，COORD为节点坐标标。4.1荷载及其施施加---施加面荷载载★定义的梯梯度仅在当当前被激活活，后面定定义的梯度度将替代前前面的。★一旦设定定了荷载梯梯度，则对对随后的荷荷载施加命命令都有效效。★取消荷载载梯度，无无参数的SFGRAD命令★命令SFGRAD,STAT可显示当前前的状态。。该命令不不能对PIPE系列单元施施加梯度荷荷载，且该该命令不能能采用表格格型边界条条件。其余命令如如SFSCALE、SFCUM、SFLIST和SFDELE等使用方法法与前面同同类命令类类似。但SFSUM仅对节点群群荷载有效效（SF命令施加的的荷载），，对于SFE、SFL及SFA无效。4.1荷载及其施施加---施加面荷载载2.施加单元荷荷载⑴在单元元上施加面面荷载命令：SFE,ELEM,LKEY,Lab,KVAL,VAL1,VAL2,VAL3,VAL4ELEM---拟施加面荷荷载的单元元号，也可可为ALL或元件名。。LKEY---与面荷载相相关的荷载载控制参数数，缺省为为1，在每个单单元的帮助助中有说明明。Lab---面荷载标识识符，结构构分析为PRES。KVAL---当Lab=PRES时，KVAL=0或1表示VAL1~VAL4为压力的实实部KVAL=2表示VAL1~VAL4为压力的虚虚部。VAL1---第一个面荷荷载值或表表格边界条条件名称，，比较典型型的是在面面上的第1个节点上，，节点的顺顺序在单元元中明确地地给定。VAL2~VAL4---为面上节点点的第2、3、4个面荷载值值，如果为为空，则与与VAL1相等；如果果为0或其它空值值则均为0；4.1荷载及其施施加---施加面荷载载对于2D平面单元，，可对单元元的任一面面（实为单单元边界））施加面荷荷载，荷载载施加到该该单元面的的角节点上上（高次单单元的中间间节点荷载载由系统自自动处理）），相邻角角节点的数数值可以不不等。对于3D体单元，用用SFE施加面荷载载时也要确确定面号及及方向才能能保证正确确（可根据据单元节点点列表确定定单元面号号），同样样也可施加加不同的荷荷载值使得得该面上各各节点荷载载不同。对于SHELL单元，其①①和②面为为底面和顶顶面，其余余为侧面（（侧边）。。SF和SFE比较而言，，对2D平面单元，，SF施加单元周周边面荷载载较为方便便，而SFE则施加单元元任一边面面荷载较为为方便；对对于3D体单元，SF施加的面荷荷载对各节节点是等值值的（除非非使用SFFUN定义），而而SFE可施加各节节点不等值值和等值两两种面荷载载；对于SHELL单元，SFE较SF方便。一般般而言，对对于通过几几何模型生生成的有限限元模型，，通过SFL和SFA命令施加荷荷载更加便便捷，且不不易出错。。4.1荷载及其施施加---施加面荷载载⑵在梁单元元施加面荷荷载命令：SFBEAM,ELEM,LKEY,Lab,VALI,VALJ,VAL2I,VAL2J,IOFFST,JOFFSTELEM---拟施加面荷荷载的单元元号，也可可为ALL或元件名。。LKEY---荷载面号(缺省为1)，在每个梁梁单元的帮帮助中有说说明。Lab---面荷载标识识符，结构构分析为PRES。VALI,VALJ---节点I和J附近的荷载载数值。如如VALJ为空则与VALI相同，否则则为其输入入值。VAL2I,VAL2J---当前未启用用。IOFFSET---VALI荷载值的作作用点离开开I节点的距离离。JOFFSET---VALJ荷载值的作作用点离开开J节点的距离离。该命令是对对梁单元（（BEAM系列）施加加单元荷载载的惟一命命令，施加加到梁单元元线（LINE）上的荷载载不能转换换到有限元元模型。梁梁单元荷载载有线性分分布荷载、、局部线性性分布荷载载、跨间集集中力三种种。对于梁单元元的垂直和和切向分布布荷载其单单位为“力力/长度”，而而对于端部部荷载则为为“力”。。4.1荷载及其施施加---施加面荷载载★线性分布荷载载：如节点I和节点J的横向分布集度度分别为q1和q2，则命令为：sfbeam,elem,1,pres,q1,q2★局部线性分布荷荷载，q1到节点I的距离为a1，q2到节点J的距离为a2,则命令为：sfbeam,elem,1,pres,q1,q2,,,a1,a2★跨间集中力：设设集中力为p1,到节点I的距离为a1，则命令为sfbeam,elem,1,pres,p1,,,,a1,-1!注意JOFFSET必须设为-1。所有荷载均相对对于单元而言，，对每个单元可可施加多个LKEY不同的荷载，但但对于同一LKEY值，只能施加一一种。如BEAM3单元LKEY=1为垂直单元轴线线的荷载，LKEY=2为平行单元轴线线的分布荷载，，而LKEY=3或4时为单元端部面面荷载（力）；；同时可利用keyopt(10)设置长度或长度度比确定IOFFSET或JOFFSET。4.1荷载及其施加---施加面荷载!EX4.7在梁单元上施加加荷载finish$/clear$/prep7et,1,beam3!定义单元类型k,1$k,2,10$l,1,2!创建关键点和线线esize,,10$lmesh,all!定义单元数目，，划分单元/pnum,elem,1!设置单元号显示示sfbeam,3,1,pres,50,100!单元3施加垂直线性分分布荷载，值分分别为50和100sfbeam,5,1,pres,100!单元5施加垂直均布荷荷载，值为100sfbeam,7,1,pres,50,100,,,0.2,0.1!单元7施加垂直局部线线性分布荷载，，值为50和100!50距离I节点0.2，100距离J节点为0.1sfbeam,9,1,pres,100,,,,0.4,-1!单元9施加集中荷载100，距离I节点0.4sfbeam,3,2,pres,50,100,,,0.2,0.1!单元3施加切向局部线线性分布荷载4.1荷载及其施加---施加面荷载3.在线上施加面荷荷载命令：SFL,LINE,Lab,VALI,VALJ,VAL2I,VAL2JLINE---拟施加荷载的线线号，也可为ALL或元件名。Lab---面荷载标识符，，结构分析为PRES。VALI---线始端关键点处处的面荷载值，，也可为表格型型边界条件的表表格名。VALJ---线末端关键点处处的面荷载值，，也可为表格型型边界条件的表格名。如为空空（缺省）与VALI相等，否则采用用输入数据。VAL2I,VAL2J---为复数输入时的的虚部，而VALI和VALJ则为实部。该命令仅对2D面单元的边界（（线）、轴对称称单元本身、壳壳单元边界（线线）有效，对3D实体单元的线无无效。对于2D面单元，其输入入的面荷载值为为“力/面积”；而对壳壳单元，其输入入的面荷载值为为“力/长度”，这点需需要特别注意。。4.1荷载及其施加---施加面荷载4.在面上施加面荷荷载命令：SFA,AREA,LKEY,Lab,VALUE,VALUE2AREA---拟施加面荷载的的面号，也可为为ALL或元件名。LKEY---荷载施加的面号号(缺省为1)。如果面为体单单元的表面，则则LKEY将被忽略；对壳壳单元LKEY可取1或2，而其它值无效效，单元帮助中中有详细说明。。Lab---面荷载标识符，，结构分析为PRES。VALUE---面荷载值，也可可为表格名称。。VALUE2---对结构分析无意意义。该命令对壳单元元和3D体单元的面施加加法向面荷载，，对2D面单元无效。4.1荷载及其施加---施加面荷载5.面荷载梯度及其其加载定义面荷载梯度度后，可在SF、SFE、SFL和SFA命令中使用。如如前所述，SFE、SFL及SFBEAM命令可以直接施施加线性分布荷荷载，采用SFGRAD命令定义荷载梯梯度后，使用SF和SFA命令施加线性分分布荷载比较方方便，如静水压压力、圆柱体分分布压力等。!EX4.9利用荷载梯度在在直角坐标系下下的施加方法FINISH$/CLEAR$/PREP7et,1,82$blc4,,,10,60!定义单元类型，，创建面esize,2$amesh,all!定义单元尺寸，，划分网格/PSF,PRES,NORM,2!设置荷载显示方方式sfgrad,pres,,y,0,-5!定义荷载梯度，，SLZER=0，沿Y正方向递减5单位/长度nsel,s,loc,x,0!选择X=0，且Y=0~40的节点群sel,r,loc,y,0,40sf,all,pres,600!对节点群施加面面荷载，基值（（Y=SLZER=0处）为600!上述荷载施加结结果：Y=0处为600，Y=40处为600+(40-0)×(-5)=400sfgrad,pres,,y,30,-20!再重新定义梯度度荷载，SLZER=30，斜率为-20nsel,s,loc,x,10!选择X=10的节点群sf,all,pres,0!对节点群施加面面荷载，基值（（Y=SLZER=30处）为0allsel$eplot!如图4-1所示4.1荷载及其施加---施加面荷载4.1荷载及其施加---施加面荷载finish$/clear$/prep7et,1,solid95$et,2,surf154!定义SOLID95单元和表面效应应单元surf154blc4,,,10,10,40!创建长方体esize,5$vmesh,all!定义单元尺寸，，划分网格/psf,pres,tany,2!设置压力显式方方式（单元坐标标系Y切向）nsel,s,loc,y,10!选择Y=10的所有节点type,2$esurf!设置单元类型2,生成表面效应单单元esel,s,type,,2!选择单元类型为为2的单元（表面效效应单元）sfe,all,3,pres,,100!施加LKEY=3的面荷载（切向向）allsel!可查看单元荷载载（求解过程与与SURF154无关）4.1荷载及其施加---施加面荷载4.1荷载及其施加---施加体荷载在结构分析中，，ANSYS的体荷载只有温温度，其标识符符为TEMP。4.1荷载及其施加---施加体荷载几个主要的体荷荷载施加命令如如下：BF,NODE,Lab,VAL1BFE,ELEM,Lab,STLOC,VAL1,VAL2,VAL3,VAL4BFK,KPOI,Lab,VAL1BFL,LINE,Lab,VAL1BFA,AREA,Lab,VAL1BFV,VOLU,Lab,VAL1其使用方法与面面荷载施加命令令类似，例如第第1个参数均为图素素编号，也可为为ALL或元件名；第2个参数Lab=TEMP或FLUE；VAL1～VAL4为体荷载值，其其中VAL2~VAL4为单元不同位置置上的体荷载值值；STLOC为VAL1指定一个对应的的起始位置。4.1荷载及其施加---施加惯性荷载惯性荷载有加速速度、角速度和和角加速度。4.1荷载及其施加---施加惯性荷载惯性荷载没有删删除命令，要删删除惯性荷载，，需将荷载值设设为0，且为斜坡荷载载。ACEL、OMEGA和DOMEGA命令分别用于施施加在总体直角角坐标系中的加加速度、角速度度和角加速度。。需要注意的是是ACEL命令施加的是加加速度不是重力力场，因此要施施加一个-Y方向的重力场，，必须施加一个个＋Y方向的加速度。。使用CGOMGA和DCGOMG命令定义一转动动物体的加速度度和角加速度，，但为相对于参参考坐标系转动动时的物理量（（该物体绕参考考坐标系转动））。CGLOC命令用于指定参参考坐标系相对对于整个笛卡尔尔坐标系的位置置。CMOMEGA和CMDOMGA命令在单元元件件上施加参考坐坐标系下的角速速度和角加速度度。ANSYS定义的三种类型型转动为：①整个结构绕总总体直角坐标系系转动（OMEGA和DOMEGA命令输入）；②单元元件绕参参考坐标系轴的的转动（CMOMEGA和CMDOMEGA命令）；③整体直角坐标标系绕加速度原原点的转动（CGOMGA、DCGOMG和CGLOC命令）。以上三种类型转转动中，可两两两组合同时施加加到结构上。此处仅介绍ACEL命令及使用方法法，命令如下：：命令：ACEL,ACELX,ACELY,ACELZ其中ACELX,ACELY,ACELZ分别为总体直角角坐标系X轴、Y轴和Z轴的结构线加速速度值。4.1荷载及其施加---施加体荷载在耦合场分析中中，通常将包含含一个分析中的的结果施加在第第二个分析中作作为荷载，例如如可将热分析中中计算得到的节节点温度，作为为体积荷载施加加到结构分析中中，形成耦合场场荷载。施加耦合场荷载载的命令为LDREAD命令，该命令是是从一个结果文文件读出数据然然后作为荷载施施加到模型上。。因此该命令不不仅仅在施加耦耦合场荷载时使使用，也可用于于其它分析目的的，例如可用于于结构分析中读读入反作用力作作为进一步分析析的荷载等。命令：LDREAD,Lab,LSTEP,SBSTEP,TIME,KIMG,Fname,Ext4.1荷载及其施加---初应力荷载及施施加初应力（InitialStress）可以指定为一一种“荷载”进进行施加，但仅仅在静态分析和和全瞬态分析中中可以使用，可可以用于线性分分析或非线性分分析。初应力荷荷载只能在第一一个荷载步中施施加，且ANSYS支持初应力荷载载的单元类型有有：PLANE2、PLANE42、PLANE82、PLANE182、PLANE183、SOLID45、SOLID92、SOLID95、SOLID185、SOLID186、SOLID187、SHELL181、SHELL208、SHELL209、LINK180、BEAM188、BEAM189。初应力荷载是单单元坐标系下的的值，如果单元元坐标系与总体体坐标系不同应应谨慎。初应力力荷载只能在求求解层施加.初应力荷载的施施加采用覆盖方方式，即多次施施加时后面命令令结果覆盖前面面的命令结果。。初应力荷载施加加在被选择的单单元上，如果单单元选择集为空空或不选择某些些单元，则不施施加初应力荷载载到这些单元上上。4.1荷载及其施加---初应力荷载及施施加4.1荷载及其施加---初应力荷载及施施加1.施加初始常应力力荷载命令：ISTRESS,Sx,Sy,Sz,Sxy,Syz,Sxz,MAT1,MAT2,MAT3,MAT4,MAT5,MAT6,MAT7,MAT8,MAT9,MAT10Sx,Sy,Sz,Sxy,Syz,Sxz---初始的常应力值值。MAT1～MAT10---初应力拟施加到到的材料编号，，如没有指定，，则施加到所有有材料上。该命令对所选择择的单元施加一一组初始常应力力值。4.1荷载及其施加---初应力荷载及施施加2.从文件施加初应应力荷载命令：ISFILE,Option,Fname,Ext,--,LOC,MAT1,MAT2,MAT3,MAT4,MAT5,MAT6,MAT7,MAT8,MAT9,MAT10Option---初应力荷载操作作控制参数，其其值可取：=READ（缺省）：从文文件读入初应力力数据；=LIST：列出已经读入入的初应力；=DELE：删除已经读入入的初应力。Fname---当Option=READ时，Fname为一目录和文件件名。当Option=LIST或DELE时，Fname为列表或删除单单元编号上的初初应力。Ext---文件扩展名或层层号，当Fname为空时，Ext缺省为“IST”。如果Option=LIST或DELE则Ext为层壳单元的层层号。4.1荷载及其施加---初应力荷载及施施加LOC---总体位置标志，，确定每个单元元内初应力要施施加的位置，其其值可取：=0（缺省）：在单单元质心上施加加初应力；=1：单元积分点上上施加初应力；；=2：在单元指定位位置上施加初应应力。即由初应应力文件确定将将初应力荷载施施加到什么位置置，此时各个单单元施加的位置置可以不相同。。=3：常应力状态。。用初应力文件件中的第一个应应力数据将所有有单元初始化为为一个常应力。。MAT1～MAT10---初应力拟施加到到的材料编号。。该命令对所选择择的单元施加初初应力荷载，初初应力的单元号号与所选择的单单元号相对应。。4.1荷载及其施加---初应力荷载及施施加3.生成初应力文件件命令：ISWRITE,Switch其中Switch参数控制初应力力文件是否生成成文件，其中可可取：ON：以工作文件及及扩展名IST生成初应力文件件，并写入数据据；OFF：不生成初应力力文件。该命令仅在求解解层有效，如果果已有同名文件件存在则覆盖之之。该命令不支支持CDWRITE命令。用ISWRITE命令写出的应力力为单元积分点点应力，对于非非线性分析，写写入的应力数据据为收敛后应力力；对于线性分分析，为求解完完成后的应力。。因此其初应力力文件标志区数数据总为eis,elemno,1，其中elemno为单元号，而1表示为积分点应应力的位置标识识。在用ISFILE命令读入时，如如果位置标志为为0，则采用各单元元的第一个应力力记录；如果位位置标志为2，则采用初应力力文件中的位置置标志（即1）；如果位置标标志为3，则采用应力文文件的第一个应应力数据。4.1荷载及其施加---初应力荷载及施施加!EX4.12初应力荷载finish$/clear$/filname,colu1$/prep7!定义工作文件名名为colu1et,1,plane82$mp,ex,1,2e5$mp,nuxy,1,0.3!定义单元类型和和材料属性blc4,,,1,10$esize,2$amesh,all!创建面，定义网网格尺寸，划分分网格nsel,s,loc,y,0$d,all,uy$d,1,ux!施加约束条件nsel,s,loc,y,10$sf,all,pres,-10!施加节点面荷载载allsel$finish/solu$iswrite,on!进入求解层，打打开初应力文件件生成开关solve$finish!求解生成初应应力文件（在在当前工作目目录中）4.1荷载及其施加加---初应力荷载及及施加Colu1.ist文件格式及内内容!*******INITIALSTRESSFILEcolu1.ist!INITIALSTRESSRECORDFORELEMENT1!SXSYSZSXYSYZSXZeis,1,1-0.283773E-1210.00000.00000-0.189735E-120.000000.00000-0.112355E-1210.00000.00000-0.200682E-120.000000.000000.195843E-1210.00000.00000-0.162630E-120.000000.000000.200728E-1210.00000.00000-0.196512E-120.000000.00000!!INITIALSTRESSRECORDFORELEMENT2!SXSYSZSXYSYZSXZeis,2,1-0.205613E-1210.00000.00000-0.562951E-130.000000.00000-0.259348E-1210.00000.00000-0.316921E-120.000000.000000.501821E-1210.00000.000000.263753E-120.000000.000000.525358E-1210.00000.00000-0.577546E-120.000000.00000!!INITIALSTRESSRECORDFORELEMENT5!SXSYSZSXYSYZSXZeis,5,1-0.774758E-1110.00000.000000.397402E-110.000000.00000-0.771516E-1110.00000.00000-0.398236E-110.000000.000000.100182E-1010.00000.000000.604651E-120.000000.000000.100338E-1010.00000.00000-0.587971E-120.000000.00000单元号位置标识4.1荷载及其施加加---初应力荷载及及施加!为说明问题，，这里重新建建模finish$/clear$/filname,colu2$/prep7!定义工作文件件名为colu2et,1,plane82$mp,ex,1,2e5$mp,nuxy,1,0.3!定义单元类型型和材料属性性blc4,,,1,10$esize,2$amesh,all!创建面，定义义网格尺寸，，划分网格nsel,s,loc,y,0$d,all,uy$d,1,ux$alls!施加约束条件件/solu!进入求解层loc=2!定义位置参数数，改变此参参数可得到不不同的加载效效果isfile,read,colu1,ist,,loc!从文件colu1.ist中读入初应力力并作为荷载载施加isfile,list!查看施加的初初应力荷载solve!求解并可查看看结果4.1荷载及其施加加---初应力荷载及及施加初应力荷载不不是施加“应应力历史”而而是一种“荷荷载”。因此对于用ISWRITE命令生成的初初应力文件，，再用ISFILE命令读入后，，当仅有初应应力荷载时，，其效果是模型中应力为为零而位移与与原荷载产生生的位移反向向。要消除由于初初应力荷载引引起的位移且且保持模型中中应力不变，可可将原荷载一一并施加，此此时模型中应应力与原荷载产生的的应力相同，，但位移场为为零（位移很很小，可认为是零））。一悬臂梁在端端部受集中力力作用，先生生成初应力文文件；然后再施施加初应力荷荷载和集中力力；计算后可可得到荷载作用下的的应力场但无无位移场。4.1荷载及其施加加---初应力荷载及及施加!EX4.13悬臂梁荷载作作用下的应力力场但无位移移场finish$/clear$/filname,cant1$/prep7!定义工作文件件名等et,1,PLANE42$mp,ex,1,2e5$mp,nuxy,1,0.3!定义单元类型型、材料特性性等blc4,,,10,1$esize,0.5$amesh,all!创建面、划分分网格等nsel,s,loc,x,0$d,all,all$allsel!施加约束等f,2,fy,-10$finish!施加端部集中中荷载等/solu$iswrite,on$solve!求解，生成初初应力文件等等/post1$pldisp$plnsol,s,x!查看求解结果果等/solu!再次进入求解解层。也可重重新开始一个个工作和恢复复模型等isfile,read,cant1,ist,,2!读入初应力文文件，作为荷荷载施加!fdele,all,all!如果删除原荷荷载则仅有初初应力荷载，，否则为二者者共同作用solve!求解（此时荷荷载为原荷载载和初应力荷荷载）/post1$pldisp$plnsol,s,x!查看求解结果果4.1荷载及其施加加---初应力荷载及及施加4.1荷载及其施加加---初应力荷载及及施加4.1荷载及其施加加---荷载步及相关关概念与荷载有关的的几个术语或或概念为：荷载步（LoadSteps）荷载子步（Substeps）斜坡荷载（RampedLoads）阶跃荷载（SteppedLoads）时间（Time）及时间步（（Timestep）平衡迭代（EquilibriumIterations）。与土木工程相相同的概念如如荷载工况和和荷载组合等等，将在后处处理中予以介介绍。4.1荷载及其施加加---荷载步及相关关概念1.荷载步、荷载载子步和平衡衡迭代荷载步是为求求解而定义的的荷载配置，，可根据荷载载历程（时间间和空间）在在不同的荷载载步内施加不不同的荷载。。例如在结构构线性静态分分析中，可将将结构自重和和外荷载分两两步施加到结结构上，第一一个荷载步可可施加自重，，第二个荷载载步可施加外外荷载等。荷载子步是在在某个荷载步步之内的求解解点（由程序序定义荷载增增量），不同同分析中荷载载子步有不同同的目的。例例如在线性静静态或稳态分分析中，使用用子步逐渐增增加荷载可获获得精确解；；在瞬态分析析中使用子步步可得到较小小的积分步长长，以满足瞬瞬态时间积累累法则；在谐谐分析中，使使用子步可获获得不同频率率下的解。平衡迭代是在在给定子步下下为了收敛而而进行的附加加计算。在非非线性分析中中，平衡迭代代作为一种迭迭代修正具有有重要作用，，迭代计算多多次收敛后得得到该荷载子子步的解。4.1荷载及其施加加---荷载步及相关关概念4.1荷载及其施加加---荷载步及相关关概念2.斜坡荷载和阶阶跃荷载当在一个荷载载步中设置一一个以上子步步时，就必须须定义荷载是是斜坡荷载或或是阶跃荷载载。阶跃荷载指荷荷载全值施加加在第一个荷荷载子步，其其余荷载子步步内荷载保持持不变。对于于荷载步2按要求是由荷荷载步1的全值荷载突突然卸载，而而程序实际上上是从荷载步步1的终点到荷载载步2的第一个子步步内完成的，，所以可增加加荷载步2的子步数（减减小时间增量量）以模拟突突然卸载过程程。斜坡荷载指在在每个荷载子子步，荷载逐逐渐增加，在在该荷载步结结束时达到荷荷载全值。载载步内子步的的荷载采用线线性内插。4.1荷载及其施加加---荷载步及相关关概念3.时间及时间步步在所有静态和和稳态分析中中，不管是否否与时间“真真实”相关，，ANSYS都使用时间作作为跟踪参数数。在瞬态分析或或与速率有关关的静态分析析（如蠕变或或粘塑性）中中，时间代表表实际的按年年月顺序的时时间，可用小小时、分、秒秒等表示。在在指定荷载历历程的同时，，在每个荷载载步终点给时时间赋值。对于与速率无无关的静态分分析，时间仅仅仅成为识别别荷载步和子子步的计数器器，每一个荷荷载步和子步步都与惟一的的时间点对应应，故子步也也称时间步。。因此这种情情况下，“time”可用任意单位位和数值。在在缺省情况下下，程序自动动对time赋值，例如在在荷载步1结束时time=1，在荷载步2结束时time=2等等。该time在后处理的时时间（荷载））---变形曲线中非非常有用。当采用弧长法法求解时，时时间不必单调调增加，可以以为负值。荷载步和子步步都与时间点点对应，即荷荷载步或子步步是一定时间间间隔内的系系列荷载。两两个连续子步步之间的时间间差称为时间间步长或时间间增量。平衡衡迭代就是为为收敛在给定定时间点上进进行迭代求解解的方法。4.2荷载步选项及及设置---输出选项1.控制写入数据据库和结果文文件的结果数数据命令：OUTRES,Item,FREQ,CnameItem---写入数据库和和结果文件的的解项（结果果）控制参数数，其值可取取：=ALL（缺省）：写写入所有解项项；=ERASE：将当前设置置恢复到ANSYS缺省状态；=STAT：当前设置状状态列表；=BASIC：仅写入NSOL,RSOL,NLOAD,STRS,FGRAD,FFLUX；=NSOL：仅写入节点点DOF结果，如UX,UY,UZ,ROTX,ROTY,ROTZ等；=RSOL：仅写入节点点反力结果；；=V：仅写入瞬态态分析时的节节点速度；=A：仅写入瞬态态分析时的节节点加速度；；=ESOL：仅写单元结结果，包括单单元节点力NLOAD，单元节点应应力STRS，单元弹性应应变EPEL，单元热、初初始和膨胀应应变EPTH，单元塑性应应变EPPL，单元蠕变应应变EPCR，单元节点梯梯度FGRAD，单元节点流流量FFLUX，积分点位置置LOCI，状态变量SVAR（仅USERMAT时），单元表表数据MISC等。4.2荷载步选项及及设置---输出选项FREQ---写入内容的频频率（即写入入哪个子步的的结果），其其值可取：=NONE：禁止写入所所有子步的内内容；=ALL：写入每个子子步的内容，，是谐分析或或EXPASS打开时的缺省省状态；=LAST：写入每个荷荷载步的最后后子步内容，，是静态或瞬瞬态分析的缺缺省；=n：写入荷载步步中每隔n个子步的内容容（包括最后后子步）；=-n：写入荷载步步中按子步均均匀分割的n个子步的内容容（仅为自动动时间步打开开时）。=%array%：存有N个时间值的数数组，程序根根据这些值写写入数据，时时间值为升序序，且数值介介于荷载步的的开始和结束束时间之间。。多荷载步时时，必须改变变时间值以保保证在荷载步步开始和结束束时间之间（（可重新定义义数组和时间间）。Cname---为CM命令创建的存存放单元或节节点的元件名名。如Item=ALL或BASIC或RSOL等就不能使用用元件名。该命令控制写写入数据库和和结果文件的的数据，当分分析生成的结结果文件特别别大时，可采采用该命令有有选择的写入入数据，以控控制RST文件的大小并并节约资源。。该命令如重复复执行，则采采用后执行的的设置，即可可利用此特性性先后设置不不同的参数，，写入不同荷荷载步下的不不同数据内容容。4.2荷载步选项及及设置---输出选项2.结果输出控制制命令：OUTPR,Item,FREQ,Cname其中参数意义义同OUTRES命令的参数，，但该命令是是控制向输出出文件.OUT写入的内容。。利用OUTRES和OUTPR组合可严格控控制输出的内内容和数据。。3.图形求解追踪踪器命令：/GST,Lab其中Lab为打开或关闭闭图形求解追追踪器控制参参数，如Lab=ON则打开，如Lab=OFF则关闭。GST将保存以GST为扩展名的文文件中（GST文件以ANSYS的DISPLAYW程序阅读）。。GST方式仅适用于于非线性结构构分析等模拟拟，如求解时时平衡迭代与与收敛过程的的图形。4.2荷载步选项及及设置---输出选项4.惯性释放计算算控制命令：IRLF,KEY其中KEY为是否考虑惯惯性释放控制制参数，其值值可取：=0（缺省）：无无惯性释放计计算；=1：用惯性释放放力平衡荷载载；=-1：仅为输出预预计算总质量量，不考虑惯惯性释放。可使用IRLIST命令显示惯性性释放计算结结果。惯性释放（inertiarelief）就是通过计计算加速度施施加惯性力来来平衡外荷载载，利用惯性性释放可求解解“全自由””结构，如水水中的船舶、、空中的飞机机、航天器等等结构的静力力分析。当然然也可利用惯惯性释放求解解外荷载自平平衡结构的内内力和变形（（加速度为0）。4.2荷载步选项及及设置---输出选项在求解全自由由结构或自平平衡结构的内内力和变形时时，可不施加加任何约束，，也可施加且且仅可施加一一个节点的““虚约束”，，该节点约束束仅仅为防止止刚体位移或或刚体转动所所必须的约束束（如对2D结构，可施加加3个自由度约束束；对3D结构可施加6个自由度约束束）。在求解解时，程序先先计算在外荷荷载作用下结结构各节点的的加速度，然然后将加速度度转化为惯性性力反向施加加到每个节点点上，由此构构造一个平衡衡的力系（支支座反力为0）,从而求解得到到所有节点相相对“虚支座座”的位移及及结构的荷载载响应。各节节点加速度通通过结构质量量矩阵和外荷荷载计算得到到，当然包括括平移加速度度和转动加速速度。惯性释释放仅仅用于于静态态分析析，且且非线线性、、子结结构、、点单单元及及轴对对称单单元等等不支支持惯惯性释释放，，梁单单元和和层壳壳单元元的质质心偏偏置和和变截截面梁梁单元元也不不支持持惯性性释放放；同同时使使用2D和3D单元的的结构构模态态分析析也不不推荐荐使用用。如如果在在第二二及后后续荷荷载步步使用用惯性性释放放，则则必须须在第第一荷荷载步步打开开EMATWRITE命令，，以便便使用用单元元刚度度矩阵阵。4.2荷载步步选项项及设设置---输出出选项项例如一一开口口框架架受自自平衡衡的外外荷载载作用用，求求其结结构内内力和和变形形。finish$/clear$/prep7et,1,beam3$mp,ex,1,2.1e11!定义单单元类类型及及材料料特性性mp,prxy,1,0.3$mp,dens,1,7800$r,1,0.2,0.05,0.2!定义质质量密密度和和实常常数k,1,10,10$k,2,15,10$k,3,15,15$k,4,10,15!创建关关键点点l,1,2$l,2,3$l,3,4$lesize,all,,,10$lmesh,all!创建线线，并并划分分网格格等/solu!进入求求解层层antype,0$irlf,1!定义分分析类类型，，打开开惯性性释放放d,19,all!对节点点19施加约约束（（可对对任意意一个个节点点施加加约束束）f,22,fy,1000$f,1,fy,-1000!施加一一对集集中力力（自自平衡衡）solve$finish$/post1!求解并并进入入后处处理层层pldisp,1!显示变变形图图（约约束仅仅对刚刚体位位移有有影响响，对对变形形形状状无影影响））etable,mi,smisc,6!定义单单元弯弯矩表表etable,mj,smisc,12plls,mi,mj,-1!显示弯弯矩图图（与与静力力计算算结果果相等等）irlist!显示惯惯性释释放加加速度度等数数据PRRSOL!显示支支座反反力（（极小小，属属于计计算误误差））4.2荷载步步选项项及设设置---输出出选项项4.2荷载步步选项项及设设置---输出出选项项4.2荷载步步选项项及设设置---其它它选项项1.定义参参考温温度命令：：TREF,TREFV其中TREFV为热膨膨胀的的参考考温度度，缺缺省值值为0.0℃。如果果没有有使用用命令令TUNIF定义均均布温温度，，也可可用此此值作作为均均布温温度使使用。。该命令令为结结构和和显式式动力力分析析中的的热应应变计计算设设定参参考温温度，，热应应变将将按α×(T-TREFV)计算，，其中中α为热膨膨胀系系数，，T是单元元温度度。如如果α为随温温度变变化时时，TREFV也应有有一定定的温温度范范围。。参考温温度也也可在在MP命令中中使用用REFT参数定定义，，如MP,REFT,MAT,C0；