ANSYS模块简介
APDL换行与续行-
APDL规定每行72个字符
如果要写表达式A=C1+C2 (C1与C2都为表达式
可以用
B=C1
A=B+C2
将一行拆成两行来做但是如果不是表达式,而是输入一个命令参数过多的话,可以用续行命令RMORE,格式如下:
RMORE, R7, R8, R9, R10, R11, R12
这个命令每次也只能输入6个参数,如果多于6个,可以重复使用RMORE就可以输入13-18,19-24等等。另外,于上面续行相应的是换行,一行命令太短可以使用多个命令共一行
$”,没有双引号。这样就可以将一行变成多行使。:) 换行符是“
ANSYS常见术语/命令对照表
Utility Menu 实用菜单
SAVE_DB 存储数据库 RESUME_DB 恢复数据库
Select Entity 选择实体 Comp/Assembly 组元/集合
Plot/Replot 画图/重新画图 Pan,Zoom,Rotate… 平移,缩放,旋转…
WorkPlane(WP) 工作平面 Coordinate System(CS) 坐标系
Macro 宏 Preference… 优先设置…
Preprocessor 前处理 General Postproc 通用后处理
TimeHist Postproc 时间历程后处理 APDL ANSYS参数化设计语言
Line Fillet 在两条线的过渡生成线 Arbitrary 任意形状
Cylinder 圆柱体 Prism 棱柱体 Cone 圆锥形
Sphere 球形 Polygon 多边形 Stress 应力
Strain 应变 Displacement 位移 DOF 自由度
Von Mises(Stress) 平均应力 Contour 等高线(图)
Deformed/Undeformed shape 变形后/未变形的形状 Results Summary 结果摘要
Radiation Matrix 辐射矩阵 Modeling 建模
Meshing 划分网格 Attribute 属性 LS (Load Step) 载荷步
ansys的常用命令介绍
对ANSYS学习也有一个来月的时间了,可是还是什么都不会!郁闷!整理了一些ANSYS常用的命令;但深知自己的水平,还不敢保证完全正确;给大家一些参考,望指正:
A,P1,P2,…,P17,P18(以点定义面)
AADD,NA1,NA2,…NA8,NA9(面相加)
AATT,MAT,REAL,TYPE,ESYS,SECN(指定面的单元属性)
【注】ESYS为坐标系统号、SECN为截面类型号。
*ABBR,Abbr,String(定义一个缩略词)
ABBRES,Lab,Fname,Ext(从文件中读取缩略词)
ABBSAVE,Lab,Fname,Ext(将当前定义的缩略词写入文件)
ABS,IR,IA,–,--,Name,–,--,FACTA(取绝对值)
ACCAT,NA1,NA2(连接面)
ACEL,ACEX,ACEY,ACEZ(定义结构的线性加速度)
ACLEAR,NA1,NA2,NINC(清除面单元网格)
ADAMS,NMODES,KSTRESS,KSHELL
ADAPT, NSOLN, STARGT, TTARGT, FACMN, FACMX, KYKPS, KYMAC
ADD,IR, IA, IB, IC, Name, --,-- , FACTA, FACTB, FACTC(变量加运算)
ADELE,NA1,NA2,NINC,KSWP(删除面)
【注】KSWP =0删除面但保留面上关键点、1删除面及面上关键点。
ADRAG,NL1,NL2,…,NL6,NLP1,NLP2,…,NLP6(将既有线沿一定路径拖拉成面)
AESIZE,ANUM,SIZE(指定面上划分单元大小)
AFILLT,NA1,NA1,RAD(两面之间生成倒角面)
AFSURF,SAREA,TLINE(在既有面单元上生成重叠的表面单元)
*AFUN, Lab(指定参数表达式中角度单位)
AGEN, ITIME, NA1, NA2, NINC, DX, DY, DZ, KINC, NOELEM, IMOVE(复制面)
AGLUE,NA1,NA2,…,NA8,NA9(面间相互粘接)
AINA,NA1,NA2,…,NA8,NA9(被选面的交集)
AINP,NA1,NA2,…,NA8,NA9(面集两两相交)
AINV,NA,NV(面体相交)
AL,L1,L2,…,L9,L10(以线定义面)
ALIST,NA1,NA2,NINC,Lab(列表显示面的信息)
【注】Lab=HPT时,显示面上硬点信息,默认为空。
【注】LabT=ALL(指定实体及其所有下层实体)、BELOW(指定实体及其下一层实体);
Entity=ALL、VOLU、AREA、LINE、KP、ELEM、NODE。
AMESH,AREA,KP1,KP2,KP3,KP4(通过点划分面单元)
/AN3D,Kywrd,KEY(三维注释)
ANCNTR,NFRAM,DELAY,NCYCL(在POST1中生成结构变形梯度线的动画)
ANCUT,NFRAM,DELAY,NCYCL,QOFF,KTOP,TOPOFF,NODE1,NODE2,NODE3(在POST1中生成等势切面云图动画)
ANDATA,DELAY,NCYCL,RSLTDAT,MIN,MAX,INCR,FRCLST,AUTOCNTRKY(生成某一范围内的结果数据的顺序梯度线动画)
ANDSCL,NFRAM,DELAY,NCYCL(在POST1中生成结构变形的动画)
ANFLOW,NFRAM,DELAY,NCYCL,TIME,SPACING,SIZE,LENGTH(生成粒子流或带电粒子运动的动画)
/ANGLE,WN,THETA,Axis,KINCR(绕指定轴旋转视图)
ANHARM,NFRAM,DELAY,NCYCL(生成谐波分析的动画)
ANIM,NCYCL,KCYCL,DELAY(动画显示图形序列)
ANISOS,NFRAM,DELAY,NCYCL(在POST1中生成等势面云图动画)
ANMODE,NFRAM,DELAY,NCYCL,KACCEL(在POST1中生成结构变形模态的动画)
/ANNOT,Lab,VAL1,VAL2(激活图形显示注释)
【注】Lab=OFF、ON、DELE、SAVE、SCALE、XORIG、YORIG、SNAP、STAT、DEFA、REFR、TMODE。
【注】NOEFLIP=0、1。
ANTIME,NFRAM,DELAY,NCYCL,AUTOCNTRKY,RSLTDAT,MIN,MAX(在指定时间段内生成动画)
ANTYPE,Antype,Status,LDSTEP, SUBSTEP, Action(定义分析类型)
【注】Label =STATIC、BUCKLE、MODAL、HARMIC、TRANS…;
Status=NEW、REST;
Action= CONTINUE、ENDSTEP、RSTCREATE。
/ANUM,NUM,TYPE,XHOT,YHOT(指定注释的数目、类型和热点位置)
AOFFST,NAREA,DIST,KNIC(偏移生成面)
AOVLAP,NA1,NA2,…,NA8,NA9(面搭接)
APLOT,NA1,NA2,NINC,DEGEN,SCALE(显示所选面)
APPEND,Lstep,SBSTEP,FACT,KIMG,TIME,ANGLE,NSET(读入载荷结果数据)
APTN,NA1,NA2,…,NA8,NA9(面间相互分割)
AREFINE,NA1,NA2,NINC,LEVEL,DEPTH,POST,RETAIN(将面附近的单元网格细化)
【注】LEVEL:指定细化的程度(1、2、3、4、5);DEPTH:指定细化的深度;
POST=CLEAN(对细化区域进行光滑和清理工作)、SMOOTH(只作光滑工作)、OFF;
RETAIN=ON(对于全是四边形的网格,细化不会将三角形引入)、OFF(可能将三角形引入)
AREVERSE,ANUM,NOEFLIP(将面的法线方向反向)
AROTAT,NL1,NL2,NL3,NL4,NL5,NL6,PAX1,PAX2,ARC,NSEG(绕轴旋转生成面)
【注】PAX1,PAX2为定义轴的关键点;ARC为旋转角度。
ARSCALE,NA1,NA2,NINC,RX,RY,RZ,KINC,NOELEM,IMOVE(面缩放)
ASBA,NA1,NA2,SEPO,KEEP1,KEEP2(面减面)
ASBL,NA,NL,SEPO,KEEP1,KEEP2(面减线)
ASBV,NA,NV,SEPO,KEEP1,KEEP2(面减体)
ASBW,NA,SEPO,KEEP(工作平面分离面)
ASEL,Type,Item,Comp,VMIN,VMAX,VINC,KSWP(选择面)
【注】Item =HPT时,选择包含硬点的面。
ASKIN,NL1,NL2,…,NL8,NL9(通过引导线由蒙皮生成光滑曲面)
ASUB,NA1,P1,P2,P3,P4(选择面的一部分生成新面)
ARSYM,Ncomp,NA1,NA2,NINC,KINC,NOELEM,IMOVE(面镜像)
ATRAN,KCNTO,NA1,NA2,NINC,KINC,NOELEM,IMOVE(将面转化到另一坐标系)
/AUTO,WN(启动自动调整模式)
AUTOTS,Key(设定自动时间步长)
AVPRIN,KEY,EFFNU(指定在同一节点处先计算主应力或矢量和,再进行平均)
/BATCH,Lab(进入批处理模式)
【注】Lab=LIST(批处理的输出包括输入文件列表)、black。
BF,NODE,Lab,VAL1,VAL2,VAL3,PHASE(在节点上施加体载荷)
BFA,AREA,Lab,VAL1,VAL2,VAL3,PHASE(在面上施加体载荷)
BFADELE,AREA,Lab(删除面上的体载荷)
BFALIST,AREA,Lab(列表显示面上的体载荷)
BFCUM,Lab,Oper,FACT,TBASE(设置节点上体载荷的施加方式)
【注】Oper=REPL(后定义的值替换原值)、ADD(后定义的值与原值相加)、IGNO(忽略后值);
BFDELE,NODE,Lab(删除节点上的体载荷)
BFE,ELEM,Lab,STLOC,VAL1,VAL2,VAL3,VAL4(在单元上施加体载荷)
BFECUM,Lab,Oper,FACT,TBASE(设置单元上体载荷的施加方式)
BFEDELE,ELEM,Lab(删除单元上的体载荷)
BFELIST,ELEM,Lab(列表显示单元上的体载荷)
BFESCALE,Lab,FACT,TBASE(按比例缩放节单元上的体载荷)
BFK,KPOI,Lab,VAL1,VAL2,VAL3,PHASE(在关键上施加体载荷)
BFKDELE,KPOI,Lab(删除关键点上的体载荷)
BFKLIST,KPOI,Lab(列表显示线关键点上的体载荷)
BFL,LINE,Lab,VAL1,VAL2,VAL3,PHASE(在线上施加体载荷)
BFLDELE,Line,Lab(删除线上的体载荷)
BFLIST,NODE,Lab(列表显示节点上的体载荷)
BFLLIST,LINE,Lab(列表显示线上的体载荷)
BFSCALE,Lab,FACT,TBASE(按比例缩放节点上的体载荷)
BFTRAN(将实体模型中的体载荷转换到有限元模型)
BFUNIF,Lab,VALUE(在所有节点施加均匀的体载荷)
BFV,VOLU,Lab,VAL1,VAL2,VAL3,PHASE(在体上施加体载荷)
BFVDELE,VOLU,Lab(删除体上的体载荷)
BFVLIST,VOLU,Lab(列表显示体上的体载荷)
BLC4,XCORNER,YCORNER,WIDTH,HEIGHT,DEPTH(指定角点位置生成矩形或长方体)
BLC5,XCENTER,YCENTER,WIDTH,HEIGHT,DEPTH(指定中心位置生成矩形或长方体)
BLOCK,X1,X2,Y1,Y2,Z1,Z2(根据两角点生成长方体)
BOPTN,Label,Value(对布尔运算进行设置)
【注】Label=KEEP,设定是否保留源图元,Value(YES、NO)
Label=PTOL,设定警告或错误信息,Value(0、2、-1)
Label=NWARN,设定布尔运算程序版本,Value(RV52、RV51)
Label=VERSION,设定运算公差,Value
BSPLIN,P1,P2,P3,P4,P5,P6,XV1,YV1,ZV1,XV6,YV6,ZV6(拟合样条曲线)
CE,NEQN,CONST,NODE1,Lab1,C1,NODE2,Lab2,C2,NODE3,Lab3,C3(生成约束方程)
CECMOD,NEQN,CONST(求解过程中修改约束方程得常数项)
CEDELE,NEQN1,NEQN2,NINC,Nsel(删除约束方程)
CEINTF,TOLER,DOF1,DOF2,DOF3,DOF4,DOF5,DOF6,MoveTol(在接触面生成约束方程)
CELIST,NEQN1,NEQN2,NINC,Nsel(列表显示约束方程)
CENTER,NODE,NODE1,NODE2,NODE3,RADIUS(将弧线的曲率中心定义为节点)
CERIG,MASTE,SLAVE,Ldof,Ldof2,Ldof3,Ldof4,Ldof5(生成刚性区域)
CESGEN,ITIME,INC,NSET1,NSET2,NINC(从既有约束方程生成新的约束方程)
CGLOC,XLOC,YLOC,ZLOC(定义加速度坐标系相对于整体直角坐标系的位置)
CGOMGA,CGOMX,CGOMY,CGOMZ(指定旋转物体的角速度)
CHECK,Sele,Lev1(检查当前数据库数据的完整性)
【注】Sele=blank(检查所有数据库数据)、ESEL(检查所选单元数据);
当Sele= ESEL时,Lev1=WARN(警告信息单元数据)、ERR(仅错误信息单元数据)。
CHKMSH,Comp(检查面或体的分网结果)
CIRCLE,PCENT,RAD,PAXIS,PZERO,ARC,NSEG(生成圆或弧)
/CLABEL,WN,KEY(定义梯度线的标签)
CLOCAL,KCN,KCS,XL,YL,ZL,THXY,THYZ,THZX,PAR1,PAR2(根据激活的坐标系定义局部坐标系)
CM, cname, entity(定义组元,将几何元素分组形成组元)
【注】cname: 由字母数字组成的组元名
entity: 组元的类型(volu, area, line, kp, elem, node)
【注】aname: 组元集名称
cname1……cname8: 已定义的组元或组元集名称
/COLOR,PBAK,Key_ On_Off,KEY_TYPE,KEY_INDEX(为背景添加纹理)
CON4,XCENTER,YCENTER,RAD1,RAD2,DEPTH(任意位置生成锥体或圆台)
CONE,RBOT,RTOP,Z1,Z2,THETA1,THETA2(以原点为中心生成锥体或圆台)
/CONTOUR,WN,NCONT,VMIN,VINC,VMAX(自定义等间隔梯度线)
CP,NSET,Lab,NODE1,NODE2,NODE3,…,NODE16,NODE17(定义耦合集)
CPDELE,NSET1,NSET2,NINC,Nsel(删除耦合)
CPINTF,Lab,TOLER(耦合重合节点自由度)
/CPLANE,KEY(定义切平面)
CPLGEN,NSETF,Lab1,Lab2,Lab3,Lab4,Lab5(取与一既有耦合集相同的节点生成不同自由度的耦合集)
CPLIST,NSET1,NSET2,NINC,Nsel(列表显示耦合)
CPSGEN,ITIME,INC,NSET1,NSET2,NINC(在选定节点生成与一既有耦合集具有相同自由度的新耦合集)
CS,KCN,KCS,NORIG,NXAX,NXYPL,PAR1,PAR2(根据已有的三个节点定义局部坐标系)
CSDELE,KCN1,KCN2,KCINC(删除局部坐标系。【注】ALL为全部)
CSKP,KCN,KCS,PORIG,PXAXS,PXYPL,PAR1,PAR2(根据三个关键点定义局部坐标系)
CSLIST,KCN1,KCN2,KCINC(列表显示局部坐标系的信息)
CSWPLA,KCN,KCS,PAR1,PAR2(根据当前定义的工作平面定义局部坐标系)
CSYS,KCN(声明坐标系统)
【注】CSYS,WP(强迫激活的坐标系在建模时与工作平面一起移动)
/CTYPE,KEY,DOTD,DOTS,DSHP,TLEN(定义梯度线显示类型)
CURR2D(计算电磁场中二维导体中的电流)
/CVAL,WN,V1,V2,V3,V4,V5,V6,V7,V8(自定义不等间隔梯度线)
CYCLIC,NSECTOR,ANGLE,KCN,NAME,USRCOMP(指定一个循环对称分析定义扇区)
【注】NSECTOR=扇区数目或者STATUS、OFF、UNDOUBLE。
CYL4,XCENTER,YCENTER,RAD1,THETA1,RAD2,THETA2,DEPTH(任意点生成圆柱或扇环柱体)
CYL5,XEDGE1,YEDGE1,XEDGE2,YEDGE2,DEPTH(以两点为底面直径端点生成圆柱体)
CYLIND,RAD1,RAD2,Z1,Z2,THETA1,THETA2(以工作平面圆点为中心生成圆柱体)
D,NODE,Lab,VALUE,VALUE2,NEND,NINC,Lab2,Lab3,Lab4,Lab5,Lab6(施加约束)
DCUM,Oper,RFACT,IFACT,TBASE(重新设置约束)
DELTIM,DTIME,DTMIN,DTMAX,Carry(定义时间步长)
DESOL,ELEM,Item,Comp,V1,V2,V3,V4,V5,V6(修改单元解数据)
DETAB,ELEM,Lab,V1,V2,V3,V4,V5,V6(修改单元表数据)
/DEVICE,Lab1,KEY(定义梯度线图显示方式)
Lab1=VECTOR:按等值线图显示;
=DITHER:按颤动云图显示。
/DEVICE,FONT,1,Val1、Val2、 Val3、 Val4、 Val5、 Val6(设置图例字体)
/DEVICE,FONT,2,Val1、Val2、 Val3、 Val4、 Val5、 Val6(设置实体编号字体)
/DEVICE,FONT,3,Val1、Val2、 Val3、 Val4、 Val5、 Val6(设置标注/图形字体)
【注】Par: 数组名
Type: array 数组,如同fortran,下标最小号为1,可以多达三维(缺省)
char 字符串组(每个元素最多8个字符)
table
IMAX,JMAX,KMAX 各维的最大下标号
Var1,Var2,Var3 各维变量名,缺省为row,column,plane(当type为table时)
/DIST,WN,DVAL,KFACT(对视图进行缩放)
DNSOL,NODE,Item,Comp,V1,V2,V3,V4,V5,V6(修改节点解数据)
DOFSEL,Type,Dof1,Dof2,Dof3,Dof4,Dof5,Dof6(选择集中载荷标识)
DOMEGA,DOMGX,DOMGY,DOMGZ(定义结构在整体指教坐标系中的角加速度)
DSCALE,RFACT,IFACT,TBASE(A按比例缩放约束)
/DSCALE,WN,DMULT(改变显示位移时所采用的比例因子)
DSYM,Lab,Normal,KCN(定义节点的约束条件对称于某轴)
【注】Lab为对称的方式:正对称(Lab=SYMM)或反对称(Lab=ASYM),Normal为对称面在目前坐标系统(KCN)的法线方向Normal=(X、Y、Z)
DSYS,KCN(定义显示坐标系)
DTRAN(将实体模型中的约束转换到有限元模型中)
E,I,J,K,L,M,N,O,P(定义元素的连接方式)
EDELE,IEL1,IEL2,INC(元素消除)
/EDGE,WN,KEY,ANGLE(定义单元边界显示方式)
/EFACET,NUM(定义单元边界分段数目)
EGEN,ITIME,NINC,IEL1,IEL2,IEINC,MINC,TINC,RINC,CINC,SINC,DX,DY,DZ(元素复制:自动编号)
EINTF,TOLER,K,TLAB,KCN,DX,DY,DZ,KNONROT(用二维线单元连接重合的节点)
ELIST,IEL1,IEL2,INC,NNKEY,RKEY(元素列表)
EMAGERR(计算电磁场分析中的相对误差)
EMF(电磁场分析中计算沿路径的电动势和电压降)
EMID,Key,Edges(增加或删除中间节点)
EMODIF,IEL,STLOC,I1,I2,I3,I4,I5,I6,I7,I8(调整单元坐标系方向)
EMORE,Q,R,S,T,U,V,W,X(单元节点超过个时,在E命令后使用)
EMUNIT, Lab, VALUE(定义磁场单位)
EN,IEL,IJ,K,L,M,N,O,P(通过节点生成指定单元)
ENGEN,IINC,ITIME,NINC,IEL1,IEL2,IEINC,MINC,TINC,RINC,CINC,SINC,DX,DY,DZ(元素复制:用户自己进行编号)
ENORM,ENUM(重新定义壳单元的法线方向)
169.ENSYM,IINC,–,NINC,IEL1,IEL2,IEINC(镜像生成新单元:用户自己进行编号)
EPLOT(元素显示)
ERASE(擦除当前图形窗口显示的内容)
EREFINE,NE1,NE2,NINC,LEVEL,DEPTH,POST,RETAIN(将单元附近的单元网格细化)
ERESX,Key(控制单元积分点解的外推方式)
Key=DEFA(线形材料单元节点解由积分点解外推得到)
YES(节点解由积分点解外推得到)
NO(节点解由积分点解拷贝得到)
ERNORM,Key(定义是否进行误差估计)
ERRANG,EMIN,EMAX,EINC(从文件读入单元数据)
ESEL, Type, Item, Comp, VMIN, VMAX, VINC, KABS(选择单元子集)
/ESHAPE,SCALE(显示单元形状)
ESIZE,SIZE,NDIV(指定线划分单元的默认数目)
ESLA, Type(选择已选面上的单元)
ESLL, Type(选择已选线上的单元)
ESLN, Type, EKEY, NodeType(选择已选节点上的单元)
ESORT,Item,Lab,ORDER,KABS,NUMB(对单元数据指定新的排序方式)
ESURF,XNODE,Tlab,Shape(在既有单元表面生成表面单元)
ESYM,–,NINC,IEL1,IEL2,IEINC(镜像生成新单元:自动编号)
ESYS,KCN(定义单元坐标系。【注】只能通过局部坐标系定义)
ET,ITYPE,Ename,KOPT1,KOPT2,KOPT3,KOPT4,KOPT5,KOPT6,INOPR(定义单元)
【注】KOPT1~KOPT6为元素特性编码,BEAM3的KOPT6=1时,表示分析后的结果可输出节点的力或力矩。
【注】Lab为字段名,最多8个字符;Item,Comp分别为单元输出表中的名称和分量。
ETABLE,ERAS(删除单元表)
ETABLE,Lab,ERAS(删除单元表中数据项)
ETLIST,ITYPE1,ITYPE2,INC(列表显示模型中使用的单元类型)
EUSORT(恢复最初的单元数据排序方式
EWRITE,Fname,Ext,–,KAPPND,Format(将单元数据写入文件)
/EXIT,Slab,Fname,Ext,–(退出)
EXTOPT,Lab,Val1,Val2,Val3(扫掠体生成控制选项)
EXTOPT,ACLEAR,Val1(指定在体扫掠完成后是否删除源面上的网格)
EXTOPT,ATTR,Val1,Val2,Val3(定义扫掠体单元属性)
EXTOPT,ESIZE,Val1,Val2(定义扫掠方向的单元尺寸)
EXTOPT,VSWE,AUTO,Val2(指定在体扫掠中是否自动选择源面和目标面)
EXTOPT,VSWE,TETS,Val2(指定在体扫掠中对无法扫掠的体用四面体划分网格)
F,NODE,Lab,VALUE,VALUE2,NEND,NINC(定义节点上的集中力)
/FACET,Lab(定义实体显示的面的表示法)
FCUM,Oper,RFACT,IFACT(重新设置集中载荷)
FDELE,NODE,Lab,NEND,NINC(删除集中载荷)
/FILNAM,Fname,Key(更改文件名称)
【注】Key=0/OFF(应用原先的log和err文件)、1/ON(新建log和err文件,但不删去原文件)。
FILL,NODE1,NODE2,NFILL,NSTRT,NINC,ITIME,INC,SPACE(节点填充)
FK,KPOL,Lab,VALUE,VALUE2(在关键点上施加集中载荷)
FKDELE,KPOI,Lab(删除关键点的集中载荷)
FKLIST,KPOI,Lab(列表显示关键点集中载荷的信息)
FLIST,NODE1,NODE2,NINC(列表显示节点集中载荷的信息)
FLST, NFIELD, NARG, TYPE, Otype, LENG
FLUXV(电磁场分析中计算中通过一闭合环路的电通量)
FMAGSUM,Cnam1,Cnam2,…,Cnam8, Cnam9(电磁场分析中单元分量上电磁力的和)
/FOCUS,WN,XF,UF,ZF,KTRANS(平移视图)
FOR2D(计算体上的电磁力)
/FORMAT,NDIGIT,Ftype,NWIDTH,DSIGNF,LINE,CHAR(定义列表格式)
/FORMAT,DEFA(恢复使用默认格式)
FSCALE,RFACT,IFACT(按比例缩放集中载荷)
FSUM,LAB,ITEM(对所选单元节点力和节点弯矩进行求和)
FTRAN(将实体模型中的集中载荷转换到有限元模型中)
FVMESH,KEEP(从分离的面单元出发,生成四面体单元)
GCGEN,Ccomp,Tcomp,NUMC,RADC,Tlab,Shape(生成接触单元)
*GET, Par, NODE, N, Item1, IT1NUM, Item2, IT2NUM(得到数值并将其存储为标量参数或者数组中元素)
/GLINE,WN,STYLE(定义单元轮廓线线型)
GPLOT(多窗口显示实体)
/GRAPHICS,Key(定义图形显示模式)
【注】Key =POWER(激活PowerGraphics显示模式)、FULL(激活全模式显示)。
/GRESUME,Fname,Ext,–(从文件中读取图形显示设置)
/GSAVE,Fname,Ext,–(将图形显示设置保存到文件中)
/GST,Lab(控制求解跟踪的图形显示) 【注】Lab=ON(跟踪求解过程中不进行图形显示)、 OFF(跟踪求解过程中不进行图形显示)。
/GTYPE,WN,Label,KEY(为各窗口选择显示内容)
/HEADER,Header,Stitle,Idstmp,Notes,Colhed,Minmax(定义列表表头显示的内容)
/HEADER,STAT(恢复使用默认的表头格式)
HPTCREAT,LINE,ENTITY,NHP,RATIO,VAL1(线上通过比率生成硬点)
HPTCREAT,LINE,ENTITY,NHP,COORD,VAL1,VAL2,VAL3(线上指定坐标生成硬点)
HPTCREAT,AREA,ENTITY,NHP,COORD,VAL1,VAL2,VAL3(面上指定坐标生成硬点)
【注】ENTITY为线/面名。
HPTDELETE,NP1,NP2,NINC(删除硬点)
IMPD,Vpath,Ipath,Vsymm,Isymm(计算特定参考平面上装置的阻抗)
INRES,Item1,Item2,Item3,…,Item7,Item8(从数据项中选择要恢复的数据)
INTSRF,Lab(合成表面上的节点结果)
【注】Lab=PRES(压力)、TAUW(剪力)、FLOW(压力及剪力)。
IRLF,Key(定义执行惯性释放计算)
IRLIST(输出惯性释放计算结果)
ISFILE,Option,Fname,Ext,–,LOC,MAT1,…,MAT10(对从文件中读取的初应力操作)
ISFILE,READ,Fname,Ext,Dir,LOC,MAT1,…,MAT10(从文件中读取初应力)
ISFILE,LIST,Fname,Ext,Dir,LOC,MAT1,…,MAT10(列表显示读取的初应力)
ISFILE,DELE,Fname,Ext,Dir,LOC,MAT1,…,MAT10(删除读取的初应力)
ISTRESS,Sx,Sy,Sz,Sxy,Syz,Sxz,MAT1,…,MAT10(施加恒定初应力)
ISWRITE,Switch(生成包含初应力的ASCII文件)
K,NPT,X,Y,Z(定义关键点)
【注】NPT:关键点号,如果赋0,则分配给最小号
KATT,MAT,REAL,TYPE,ESYS(指定关键点的单元属性)
KBC,KEY(指定载荷增加方式)
KBETW,KP1,KP2,KPNEW,Type,VALVE(在两个关键点之间生成关键点)
【注】Type=RATIO、DIST。
KCALC,KPLAN,MAT,KCSYM,KLOCPR(断裂力学分析中计算应力强度因子)
KCENTER,Type,VAL1,VAL2,VAL3,VAL4,KPNEW(在三点定义圆中心生成关键点)
【注】Type=KP、LINE。
KCLEAR,NP1,NP2,NINC(清除点单元网格)
KDELE,NP1,NP2,NINC(删除关键点)
KDIST,KP1,KP2(计算关键点间距)
KESIZE,NPT,SIZE,FACT1,FAC2(指定关键点附近单元尺寸)
KEYOPT, ITYPE, KNUM, VALUE(单元主要选项设置)
KFILL,NP1,NP2,NFILL,NSTRT,NINC,SPACE(在两个关键点之间生成一批关键点)
KGEN,itime,Np1,Np2,Ninc,Dx,Dy,Dz,kinc,noelem,imove
【注】Itime:拷贝份数
Np1,Np2,Ninc:所选关键点
Dx,Dy,Dz:偏移坐标
Kinc:每份之间节点号增量
noelem: “0” 如果附有节点及单元,则一起拷贝。
“1”不拷贝节点和单元
imove: “0” 生成拷贝
“1”移动原关键点至新位置,并保持号码,此时(itime,kinc,noelem)被忽略
注意:MAT,REAL,TYPE 将一起拷贝,不是当前的MAT,REAL,TYPE
KL,NL1,RATIO,NK1(在线上定义关键点)
KLIST,NP1,NP2,NINC,Lab(列表显示关键点信息)
【注】Lab=HPT时,显示硬点信息。
KMESH,NP1,NP2,NINC(在关键点处生成点单元)
KMODIF,NPT,X,Y,Z(修改关键点信息)
KMOVE,NPT,KC1,X1,Y1,Z1,KC2,X2,Y2,Z2(移动关键点到一交点处)
KNODE,NPT,NODE(在既有节点位置生成关键点。【注】NPT为关键点号)
KPLOT,NP1,NP2,NINC,Lab(显示关键点信息)
【注】Lab=HPT时,显示硬点信息。
KPSCALE,NP1,NP2,NINC,RX,RY,RZ,KINC,NOELEM,IMOVE(点缩放)
KREFINE,NP1,NP2,NINC,LEVEL,DEPTH,POST,RETAIN(将关键点附近的单元网格细化)
KSEL,Type,HPT,Comp,VMIN,VMAX,VINC,KABS(选择硬点)
KSYMM,Ncomp,NP1,NP2,NINC,KINC,NOELEM,IMOVE(点镜像)
KUSE,KEY(指定是否使用新的三角矩阵)
KWPAVE,P1,P2,…,P8,P9(将一组关键点中心位置作为工作平面原点)
KWPLAN,,KORIG,KXAX,KPLAN(通过三个关键点定义工作平面)
L,P1,P2,NDIV,SPACE,XV1,YV1,ZV1,XV2,YV2,ZV2(定义线)
L2ANG,NL1,NL2,ANG1,ANG2,PHIT1,PHIT2(生成直线与两直线均成一定角度)
L2TAN,NL1,NL2(生成直线与两直线均相切)
LANG,NL1,P3,ANG,PHIT,LOCAT(生成直线与已知直线成一定角度)
LARC,P1,P2,PC,RAD(生成弧线)
/LARC,XCENTER,YCENTER,XLRAD,ANGLE1,ANGLE2(使用弧线注释)
LAREA,P1,P2,NAREA(在面上两关键点之间生成一条最短的线)
LATT,MAT,REAL,TYPE,–,KB,KE,SECNUM(指定线的单元属性)
LCABS,LCNO,KABS(指定是否对载荷工况取绝对值)
LCASE,LCNO(将载荷工况读入)
LCDEF,LCNO,LSTEP,SBSTEP,KIMG(从结果文件中定义载荷工况)
LCDEF,LCNO,ERASE(删除一载荷工况)
LCFACT,LCNO,FACT(指定载荷工况的比例因子)
LCFILE,LCNO,Fname,Ext,–(从载荷工况文件中定义载荷工况)
LCLEAR,NL1,NL2,NINC(清除线单元网格)
LCOMB,NL1,NL2,KEEP(线线合并)
LCOPER,Oper1,LCASE,Oper2,LCASE2(载荷工况的组合运算)
【注】Oper1=ADD(加)、SUB(减)、SQUA(平方)、SQRT(平方根)、SRSS(平方和求平方根)、MIN(比较存储最小值)、MAX(比较存储最大值)。
LCOPER,LPRIN(重新计算线单元的主应力)
LCSEL,Type,LCMIN,LCMAX,LCINC(选择载荷工况)
LCWRITE,LCNO,Fname,Ext,–(将当前载荷工况写入载荷工况文件中)
LCZERO(清空数据库中以前的数据)
LDELE,NL1,NL2,NINC,KSWP(删除线)
【注】KSWP=0删除线但保留线上关键点、1删除线及线上关键点。
LDIV,NL1,RATIO,PDIV,NDIV,KEEP(将线分割)
LDRAG,NK1,NK2,…,NK6,NL1,NL2,…,NL6(将一组既有关键点按一定路径拖拉成线)
LDREAD,Lab,LSTEP,SBSTEP,TIME,KIMG,Fname,Ext,–(施加耦合场载荷)
LESIZE,NL1,SIZE,ANGSIZ,NDIV,SPACE,KFORC,LAYER1,LAYER2,KYNDIV(指定所选线上单元数)
LEXTND,NL1,NK1,DIST,KEEP(将线延伸)
LFILLT,NL1,NL2,RAD,PCENT(两条相交线生成倒角)
LFSURF,SLINE,TLINE(在既有面单元上生成重叠的表面单元)
LGLUE,NL1,NL2,…,NL8,NL9(线间相互粘接)
/LIGHT,WN,NUM,INT,XV,YV,ZV,REFL(为模型添加光源)
LINA,NL,NA(线面相交)
/LINE,X1,Y1,X2,Y2(使用线注释)
LINL,NL1,NL2,…NL8,NL9(被选线的交集)
LINP,NL1,NL2,…NL8,NL9(线集两两相交)
LINV,NL,NV(线体相交)
LLIST,NL1,NL2,NINC,Lab(列表显示线信息)
【注】Lab=HPT时,显示线上硬点信息。
LMATRIX,Symfac,Coilname,Curname,Indname(计算任一导体的自感应矩阵)
LMESH,NL1,NL2,NINC(划分线生成线单元)
LOCAL,KCN,KCS,XC,YC,ZC,THXY,THYZ,THZX,PAR1,PAR2(定义局部坐标系统)
LOVLAP,NL1,NL2,…,NL8,NL9(线搭接)
LPLOT,NL1,NL2,NINC(显示线的信息)
LPSCALE,NL1,NL2,NINC,RX,RY,RZ,KINC,NOELEM,IMOVE(线缩放)
LPTN,NL1,NL2,…,NL8,NL9(线间相互分割)
LREFINE,NL1,NL2,NINC,LEVEL,DEPTH,POST,RETAIN(将选定线附近的单元网格细化)
LREVERSE,LNUM,NOEFLIP(将线的法线方向反向)
LROTAT,NK1,NK2,NK3,NK4,NK5,NK6,PAX1,PAX2,ARC,NSEG(绕轴旋转生成线)
LSBA,NL,NA,SEPO,KEEPL,KEEPA(线减面)
LSBL,NL1,NL2,SEPO,KEEP1,KEEP2(线减线)
LSBV,NL,NV,SEPO,KEEPL,KEEPV(线减体)
LSBW,NL,SEPO,KEEP(工作平面分离线)
LSCLEAR,Lab(删除所有得载荷)
LSDELE,LSMIN,LSMAX,LSINC(删除载荷步文件)
LSEL,Type,Item,Comp,VMIN,VMAX,VINC,KSWP(选择线)
【注】Item =HPT时,选择包含硬点的线。
/LSPEC,LCOLOR,LINSTL,XLNWID(定义线注释属性)
LSREAD,LSNUM(读入载荷步文件)
LSSOLVE,LSMIN,LSMAX,LSINC(顺序读取载荷步文件求解)
LSTR,P1,P2(两点间生成直线)
LSWRITE,LSNUM(将载荷步写入文件)
/LSYMBOL,X,Y,SYMANG,SYMTYP,SYMSIZ,KEYBMP(使用符号注释)
LSYMM,Ncomp,NL1,NL2,NINC,KINC,NOELEM,IMOVE(线镜像)
LTAN,NL1,P3,XV3,YV3,ZV3(生成直线与已知直线相切)
LWPLAN,WN,NL1,RATIO(通过指定线上一点并垂直于此线定义工作平面)
/MPLIB,R-W_opt,PATH(定义材料库路径)
MAGSOLV,OPT,NRAMP,CNVCSG,CNVFLUX,NEQIT,BIOT(定义磁场求解选项)
MAT,MAT(指定材料特性)
MCHECK,Lab(检查分网连接性)
MMF(计算沿某路径的磁动力)
MODMSH,Lab (实体或有限元模型的选项设定)
【注】Lab= STAT(给出当前网格划分状态)、NOCHECK(禁止模型相互对照检查)、CHECK(恢复模型相互对照检查)、DETACH(分离有限元模型和几何实体模型)。
【注】模态提取方法Method=LANB、SUBSP、REDUC、UNSYM、DAMP、QRDAMP。
Lab =EXPND(控制网格扩展模式)、PYRA(控制生成过度金字塔形单元)、TIMP(控制四面体单元改进程度)、TRANS(控制网格过度模式)等。
MOVE,NODE,KC1,X1,Y1,Z1,KC2,X2,Y2,Z2(移动一节点到坐标系平面的交点)
MP,Lab,MAT,C0,C1,C2,C3,C4(定义线性材料特性)
MPCHG,MAT,ELEM(改变单元材料)
MPDATA,Lab,MAT,STLOC,C1,C2,C3,C4,C5,C6
MPLIST,MAT1,MAT2,INC,L,TEVL(列表显示材料线性特性)
MPPLOT,Lab,MAT,TMIN,TMAX,PMIN,PMAX(图形显示线性材料关于温度的变化曲线)
MPREAD,Fname,Ext,Dir,LIB(从文件中读入材料特性)
MPTEMP,STLOC,T1,T2,T3,T4,T5,T6(在材料特性中定义温度)
MPWRITE,Fname,Ext,–,LIB,MAT(将材料特性保存到文件)
MSHAPE,KEY,Dimension(指定单元划分的形状)
【注】单元类型标志KEY=0(2D时划分为四边形、3D时划分为六面体)、1(2D时划分为三角形、3D时划分为四面体);维数Dimension=2D、3D。
【注】 KEY=0(自由分网)、1(映射分网)、2(自主选择)。
【注】 KEY= 0:中间节点在单元曲线边界上;
1:中间节点在单元直线边界上;
2:不生成中间节点。
N,NODE,X,Y,Z,THXY,THYZ,THZX(通过坐标定义节点)
NDELE,NODE1,NODE2,NINC(删除节点)
NDIST,ND1,ND2(计算两节点之间距离)
NDSURF,Snode,Telem,DIMN(在既有单元表面叠加生成表面单元)
NFORCE,ITEM(计算每个所选节点上的节点力和节点弯矩,并对与所选节点相关单元进行求和)
NGEN,ITIME,INC,NODE1,NODE2,NINC,DX,DY,DZ,SPACE(节点复制)
NKPT,NODE,NPT(将既有关键点定义为节点)
NLIST,NODE1,NODE2,NINC,Lcoord,SORT1,SORT2,SORT3(节点列表)
NMODIF ,NODE,X,Y,Z,THXY,THYZ,THZX(修改节点信息)
NOORDER,Lab(单元重新排序)
NPLOT,KNUM (节点显示:KNUM=0不显示号码;KNUM=1显示号码)
NROTAT,NODE1,NODE2,NINC(使节点坐标系方向与当前坐标系方向一致)
NREAD,Fname,Ext,–(从节点文件中读入节点数据)
NREFINE,NN1,NN2,NINC,LEVEL,DEPTH,POST,RETAIN(将节点附近的单元网格细化)
NRRANG,NMIN,NMAX,NINC(指定从节点文件中读取节点的范围)
NSCALE,INC,NODE1,NODE2,NINC,RX,RY,RZ(比例缩放原有节点生成新节点)
NSEL,Type,Item,Comp,VMIN,VMAX,VINC,KABS(节点的选择)
【注】 选择方式Type=S(选择一组节点为Active的节点)、R(在现有的Active节点中,再选择某些节点为Active
节点)、A(再选择某些节点,加入现有Active节点中)、U(在现有Active节点中,排除某些节点)、ALL(全部);
资料卷标Item=NODE(用节点号码选取)、LOC(用节点坐标选取);
Comp=无(Item=NODE)、X(Y、Z)(表示节点X(Y、Z)坐标为准,当Item=LOC)。
【注】NKEY=0(仅选择所选线内部的节点)、1(选择所有节点包括关键点上的节点)。
NSORT,Item,Comp,ORDER,KABS,NUMB,SEL(为节点数据指定新的排序方式)
NSUBST,NSBSTP,NSBMX,NSBMN,Carry(直接指定子步数)
NSYM,Ncomp,INC,NODE1,NODE2,NINC(通过镜像既有节点生成新节点)
/NUMBER,NKEY(设置显示实体编号时,编号与颜色是否一起显示)
NUMCMP,Label(压缩编号消除空号)
NUMMEG,Label,TOLER,GTOLER,Action,Switch(合并重复项)
NUMOFF,Label,VALUE(设置编号偏差值)
NUMSTR,Label,VALUE(设置起始编号)
NUSORT(恢复最初的节点数据排序方式)
NWPAVE,N1,N2,N3,…,N8,N9(将一组节点中间位置作为工作平面原点)
NWPLAN,WN,NORIG,NXAX,NPLAN(通过三个节点定义工作平面)
NWRITE,Fname,Ext,–,KAPPND(将节点数据写入到一文件中)
OMEGA,OMEGAX,OMEGAY,OMEGAZ,KSPIN(定义结构等速旋转负荷效应)
OUTPR,Item,FREQ,Cname(控制输出文件的记录内容)
【注】在SOLVE之前,欲选择结果的内容Item=ALL(全部)、NOSL(节点自由度结果、BASIC(系统默认);负载的次数FREQ=ALL(最后负载)。
/OUTPUT,Fname,Ext,–,Loc(求解之前将求解信息写入文件或输出到屏幕)
OUTRES, Item,FREQ,Cname(控制数据库和结果文件中的记录内容)
【注】欲选择结果的内容Item=ALL(系统默认);负载的次数FREQ=LAST(最后负载)。
PADELE,DELOPT(删除路径)
/PAGE,ILINE,ICHAR,BLINE,BCHAR(定义输出时列表和屏幕页的尺寸)
PAGET,PARRAY,POPT(将路径数据保存到预先定义好的数组中)
PAPUT,PARRAY,POPT(从数组中恢复路径及路径中的数据)
PARESU,Lab,Fname,Ext,–(从文件中恢复路径及路径中的数据)
PASAVE,Lab,Fname,Ext,–(将路径上的数据保存到文件中)
PATH,NAME,nPts,nSets,nDiv(定义路径属性)
PATH,STATUS(查看所定义的路径)
/PBC,Item,–,Key,MIN,MAX,ABS(显示边界条件符号和数值)
/PBF,Item,–,KEY(在图形上显示体载荷)
PCALC,Oper,LabR,Lab1,Lab2,FACT1,FACT2,CONST(路径上数据的数学运算)
【注】Oper=ADD(加)、MULT(乘)、DIV(除0、EXP(幂)、DERI(微分)、INTG(积分)、(自然对数)、SIN(正玄)、COS(余玄)、ASIN(反正玄)、ACOS(反余玄)
PCIRC,RAD1,RAD2,THETA1,THETA2(以工作平面圆点为中心生成圆或环行面)
/PCIRCLE,XCENTER,YCENTER,XLRAD(使用圆形注释)
PCROSS,LabXR,LabYR,LabZR,LabX1,LabY1,LabZ1,LabX2,LabY2,LabZ2(求两路径矢量的点积)
PDEF,Lab,Item,Comp,Avglab(将数据影射到路径上)
PDEF,STAT(查看路径所包含的数据项)
PDEF,CLEAR(清除所有用户定义的数据项)
389.PDOT,LabR,LabX1,LabY1,LabZ1,LabX2,LabY2,LabZ2(求两路径矢量的点积)
PLCRACK,LOC,NUM(显示裂纹及破碎区)
PLDISP,KUND(显示变形结果)
【注】KUND=0(显示变形后的结构形状)、1(同时显示变形前及变形后的结构形状)、2(同时显示变形前及变形后的结构形状,但仅显示结构外观)。
PLESOL,Item,Comp(显示单元解答)
【注】Item——欲查看何种解答;Comp——Item所定义的分量。
Item: S S S F M
Comp: X,Y,Z,XY,YZ,XZ:应力 1,2,3:主应力 EQV,INT:等效应力 X,Y,Z:结构力 X,Y,Z:结构力矩
PLESOL,SERR,,KUND,Fact(图表显示单元能量误差)
PLESOL,SERR(列表显示单元能量误差)
【注】图形横轴为单元号码,纵轴为Itlab值。
Itlab为定义的表格字段名称;Avglab=NOAV(不平均共同节点的值)、AVG(平均共同节点的值)
PLF2D,NCONT,OLAY,ANUM,WIN(生成等势轮廓线图)
PLLS,LabI,LabJ,Fact,KUND(梯度线图显示线单元解)
PLLS,LabI,LabJ,Fact(显示1-D线单元节点结果)
【注】LabI,LabJ为前面定义的I、J点结果表。
【注】以连续轮廓线方式表达,会有连续的状态,比PLESOL好。
Item(+31): U ROT TEMP
Comp(+31): X,Y,Z,SUM:位移分量及向量位移 X,Y,Z,SUM:旋转位移分量及向量旋转位移 无:温度
/PLOPTS,Label,KEY(设置窗口显示选项)
PLPAGM,Item,Gscale,Nopt(沿路径以梯度线图显示数据项的变化)
PLPATH,Lab1,Lab2,Lab3,Lab4,Lab5,Lab6(绘制路径上一数据项关于另一数据项的变化曲线图)
PLSECT,Item,Comp,RHO,KBR(沿某路径显示薄膜应力或薄膜应力加弯曲应力薄膜变化)
PLTRAC,Analopt,Item,Comp,TRPNum,Name,MXLOOP,TOLER,OPTION(图形显示运动轨迹
PLVECT,Item,Lab2,Lab3,LabP,MODE,Loc,Edge(定义矢量图显示方式)
PMAP,FORM,DISCON(定义路径点间分段点插补方式)
PMGTRAN,Fname,FREQ,Fcnam1,Fcnam2,Pcnam1,Pcnam2,Ecnam1,Ccnam1(对瞬态电磁分析的结果求和)
/PMORE,–,X5,Y5,X6,Y6,X7,Y7,X8,Y8(使用多边形注释)
/PNUM,Label,KEY(设置图形显示风格)
【注】Label=NODE、ELEM、KP、LINE、AREA、VOLU、DOMAIN、SEC、MAT、TYPE、REAL、LOC、SVAL、TABNAM、STAT、ESYS、DEFA。
/POLYGON,NVERT,X1,Y1,X2,Y2,X3,Y3,X4,Y4(使用多边形注释)
/POST1(进入通用后处理器)
/POST26(进入时间历程后处理器)
POWERH(计算导体的RMS功率损失)
PPATH,POINT,NODE,X,Y,Z,CS(定义路径上的点)
PRANGE,LINC,VMIN,VMAX,XVAR(定义曲线图X轴数据变量及范围)
PRERR(列表显示所选单元的结构能量(热能)百分比误差)
PRESOL,Item,Comp(列表显示单元解数据)
【注】以单元为单位,将Item所声明项次的结果显示在窗口中,使用者可以保存该资料。
Item: S F M
Comp: X,Y,Z,XY,YZ,XZ:应力 X,Y,Z:结构力 X,Y,Z:结构力矩
PRETAB,Lab1,Lab2,Lab3,Lab4,Lab5,Lab6,Lab7,Lab8,Lab9(列表显示单元表数据)
PRITER(列表求解概要包括载荷步信息、迭代信息以及收敛信息)
PRNLD,Lab,TOL,Item(列表显示单元节点力)
PRNSOL,Item,Comp(列表显示节点解数据)
【注】以节点为单位将Item所声明项次的结果显示在窗口中,可保存。
Item: U U S S
Comp:X,(Y,Z):位移 COMP:X,Y,Z方向及总向量方向的位移 COMP:应力 PRIN:主应力,等效应力
PRPATH,Lab1,Lab2,Lab3,Lab4,Lab5,Lab6(沿路径列表显示数据项)
PRRSOL,Lab(列表显示反作用力数据)
PRSECT,RHO,KBR(沿某路径将应力分解为薄膜应力或薄膜应力加弯曲应力,列表输出)
PRSSOL,Item,Comp(列表显示截面上的结果(仅对BEAM188和BEAM189单元有效))
PRVECT,Item,Lab2,Lab3,LabP(列表显示矢量数据项的大小及方向余玄)
/PSF,Item,Comp,KEY,KSHELL,Color(图形显示面载荷)
/PSPEC,PCOLOR,KFILL,KBORDR(定义多边形注释属性)
/PSTATUS,WN(查看各窗口的图形设置信息)
PSTRES,Key(是否计入预应力)
【注】Key=OFF(不计入)、ON(计入)。
/PSYMB,Label,KEY(符号显示控制)
PVECT,Oper,LabXR,LabYR,LabZR(定义沿路径的法向量、切向量或位置向量)
/PWEDGE,XCENTER,YCENTER,XLRAD,ANGLE1,ANGLE2(使用饼形注释)
QFACT(从模型频率解中计算高频电磁振荡的品质因子)
QUAD,NODE1,NINTR,NODE2,NFILL,NSTRT,NINC,PKFAC(从三节点的二次线上生成节点)
R,NSET,R1,R2,R3,R4,R5,R6(定义实常数)
【注】实常数是与单元类型一起使用的,实常数中R1、R2等值对不同的单元类型有不同的意义。
RALL(给出模型所有统计信息)
/RATIO,WN,RATOX,RATOY(改变图形显示的纵横比例)
REAL,NSET(指定实常数)
RECTNG,X1,X2,Y1,Y2(两角(左下、右上)生成矩形)
REFLCOEF,Portin,Pvolt,Pang,Pdist,Vpathy(计算电压反射系数、驻波比(VSWR)及在同轴装置中的回流损失)
*REPEAT, NTOT, VINC1, VINC2, VINC3, …VINC11(重复前一命令)
/REPLOT,Label(重新生成当前显示)
/RESET(恢复缺省的图形设置)
RFILSZ(估算生成文件大小)
RITER,NITER(获得分析的迭代次数)
RLIST,NSET1,NSET2,NINC(列表显示模型中使用的实常数)
PLF2D, NCONT, OLAY, ANUM, WIN
RMEMRY(给出内存统计量)
RPOLY,NSIDES,LSIDE,MAJRAD,MINRAD(以原点为中心生成正多边形面)
【注】NSIDES为边数;LSIDE为边长; MAJRAD为外接圆半径;MINRAD为内接圆半径。
【注】RADIUS为外接圆半径。
【注】Z1为底面高度;Z2为顶面高度;NSIDES为边数;LSIDE为边长; MAJRAD为外接圆半径;MINRAD为内接圆半径。
RSPEED,MIPS,SMFLOP,VMFLOP(获得计算机的性能)
RSTAT(给出模型单元、节点的信息)
RSYS,KCN(定义结果坐标系)
RTHICK,Par,ILOC,JLOC,KLOC,LLOC(对某些板壳单元定义不同的厚度)
RTIMST(获得求解时间估算值)
RWFRNT(估计内存要求)
SABS,Key(单元表中数据是否取绝对值)
SADD,LabR,Lab,Lab2,FACT1,FACT2,CONST(对单元表中两列数据项进行加运算)
SALLOW,STRS1,STRS2,STRS3,STRS4,STRS5,STRS6(定义允许应力值)
SBCTRAN(将实体模型中的所有边界条件转换到有限元模型)
SDELETE,SFRIST,SLAST,SINC,KNOCLEAN(删除截面类型)
SECDATA,VAL1,VAL2,…,VAL9,VAL10(定义截面几何数据)
SECNUM,SECID(指定横截面)
SECOFFSET,Location,OFFSET1,OFFSET2,CG-Y,CG-Z,SH-Y,SH-Z(指定单元节点在截面上的位置)
SECPLOT,SECID,VAL1,VAL2(图形显示截面)
SECREAD,Fname,Ext,–,Option(读入截面文件)
SECTYPE,SECID,Type,Subtype,Name,REFINEKEY(选择截面类型:共11种)
467.SECWRITE,Fname,Ext,–,ELEM_TYPE(将截面节点及单元信息保存到.sec文件中)
SEDLIST,Sename,KOPT(列表显示超单元的自由度解)
/SEG,Lvabel,Aviname,DELAY(将图象数据存储在本地)
/SEG,MULTI,Aviname,DELAY(定义保存图形序列到独立的图形段)
/SEG,OFF(关闭捕捉图片函数)
/SEG,STAT(查看图形段状态)
/SEG,DELE(删除图形段)
SENERGY,OPT,ANTYPE(确定存储的磁能或共轭能量)
SET,Lstep,SBSTEP,FACT,KIMG,TIME,ANGLE,NSET(读入载荷结果数据)
SEXP,LabR,Lab1,Lab2,EXP1,EXP2(对单元表中两列数据项进行幂运算及乘法运算)
SF,Nlist,Lab,VALUE,VALUE2(在节点上施加面载荷)
SFA,AREA,LKEY,Lab,VALUE,VALUE2(在面上施加面载荷)
SFACT,TYPE(对安全系数种类进行定义)
SFADELE,AREA,LKEY,Lab(删除面上的面载荷)
SFALIST,AREA,Lab(列表显示面上的面载荷)
SFBEAM,ELEM,LKEY,Lab,VALI,VALJ,VAL2I,VAL2J,IOFFS,JOFFS(在梁上定义分布载荷) 【注】ELEM为分布载荷的施加单元名,如果为P,则在图形界面中选取;LKEY默认值为1。
SFCALC,LabR,LabS,LabT,TYPE(计算安全系数)
SFCUM,Oper,FACT,FACT2(设置面载荷施加方式)
SFDELE,Nlist,Lab(删除节点上的面载荷)
SFE,ELEM,LKEY,Lab,KVAL,VAL1,VAL2,VAL3,VAL4(分布力作用单元上)
SFEDELE,NODE,LKEY,Lab(删除单元上的面载荷)
SFELIST,ELEM,Lab(列表显示单元上的面载荷)
SFGRAD,Lab,SLKCN,Sldir,SLZER,SLOPE(指定面载荷的梯度)
SFL,LINE,Lab,VALI,VALJ,VAL2I,VAL2J(在线上施加面载荷)
SFLDELE,Line,Lab(删除线上的面载荷)
SFLIST,NODE,Lab(列表显示节点上的面载荷)
SFLLIST,LINE,Lab(列表显示线上的面载荷)
SFSCALE,Lab,FACT,FACT2(按比例缩放面载荷)
SFTRAN(将实体模型中的面载荷转化到有限元模型中)
/SHADE,WN,TYPE(指定表面明暗类型)
SHPP,Lab,VALUE1,VALUE2(单元形状检查)
SHPP,ON,VALUE1(设置个别形状检查)
SHPP,SUMMARY(查看形状检查结果)
SHPP,STATUS(查看形状检查参数设置)
SHPP,MODIF,VALUE1,VALUE2(改变形状参数限制)
SHRINK,RATIO(收缩显示)
SLIST,SFRIST,SLAST,SINC,Details,Type(列表显示截面)
SMAX,LabR,Lab,Lab2,FACT1,FACT2(比较并存储两列数据项中较大值)
SMIN,LabR,Lab,Lab2,FACT1,FACT2(比较并存储两列数据项中较小值)
SMRTSIZE,SIZLVL,FAC,EXPND,TRANS,ANGL,ANGH,GRATIO,SMHLC,SMANC,MXITR,SPRX(定义自动分网参数)
SMULT,LabR,Lab1,Lab2,FACT1,FACT2(对单元表中两列数据项进行乘法运算)
/SOLU(进入求解器单元)
SOLVE(开始求解)
SPARM,Porti,Portj(计算同轴波导或TE10模式激励矩形波导的两端口间的散射系数)
503.SPH4,XCENTER,YCENTER,AD1,RAD2(工作平面任意位置生成球体或空心球体)
SPH5,XEDGE1,YEDGE1,XEDGE2,YEDGE2(以两点确定的直径生成球体)
SPHERE,RAD1,RAD2,THETA1,THETA2(以原点为中心生成球体或空心球体)
SPLINE,P1,P2,P3,P4,P5,P6,XV1,YV1,ZV1,XV6,YV6,ZV6(生成分段样条曲线)
SPOINT,NODE,X,Y,Z(计算某点的力矩和)
/SSCALE,WN,SMULT(定义等值线的外形)
SSUM(计算并打印单元表中每列之和)
STAT(显示数据库设置状态)
/STITLE,NLINE,Title(定义列表的副标题)
SUBSET,Lstep,SBSTEP,FACT,KIMG,TIME,ANGLE,NSET(读入载荷结果数据)
SUMTYPE,Label(指定载荷工况组合时的应力类型)
TALLOW,TEMP1,TEMP2,TEMP3,TEMP4,TEMP5,TEMP6(定义允许应力值相应的温度)
TB,Lab,MAT,NTEMP,NPTS,TBOPT,EOSOPT(在非线性材料属性或特殊单元输入中激活一单元表格)
TBDATA,STLOC,C1,C2,C3,C4,C5,C6(定义单元表格中的数据)
TBLIST,Lab,MAT(列表显示材料非线性特性)
518.TBPLOT,Lab,MAT,TBOPT,TEMP,SEGN(图形显示非线性材料的应力-应变曲线)
【注】oper: defi 定义一个点
dele 删除一个点
x,y:坐标
TCHG,ELEM1,ELEM2,ETYPE2(将四面体退化单元转化为非退化单元)
TIME,TIME(通过时间定义载荷步)
TIMP,ELEM,CHGBND,IMPLEVEL(对不附属于体的四面体单元进行改进)
/TLABEL,XLOC,YLOC,Text(使用文字注释)
TOFFST,VALUE(选择温度的单位)
TORQ2D(计算磁场中物体上的扭矩)
TORQC2D,RAD,NUMN,LCSYS(计算磁场中物体上环行路径的扭矩)
TORQSUM,Cnam1,Cnam2,…,Cnam8, Cnam9(对2-D平面问题中单元上的电磁麦克斯韦和虚功扭矩求和)
TORUS,RAD1,RAD2,RAD3,THETA1,THETA2(生成环体)
【注】RAD1,RAD2,RAD3中最大直径为主半径,最小为内半径,中间值为外半径。
529.TRANSFER,KCNTO,INC,NODE1,NODE2,NINC(将节点模式转换到另一坐标系中)
TREF,TREF(定义参考温度)
/TRIAD,Lab(控制是否显示整体坐标系标志,并对其位置进行定义) 【注】Lab=ORIG(在原点显示坐标系)、OFF(关闭显示)、LBOT(在左下角显示坐标系)、RBOT(在右下角显示坐标系)、LTOP(在左上角显示坐标系)、RTOP(在右上角显示坐标系)。
/TRLCY,Lab,TLEVEL,N1,N2,NINC(透明显示)
TRPDEL,NTRP1,NTRP2,TRPINC(删除轨迹点)
TRPLIS,NTRP1,NTRP2,TRPINC(列表显示轨迹点信息)
TRPOIN,X,Y,Z,VX,VY,VZ,CHRG,MASS(定义粒子流轨迹上的点)
TRTIME,TIME,SPACING,OFFSET,SIZE,LENGTH(定义流动轨迹时间间隔)
/TSPEC,TCOLOR,TSIZE,TXTHIC,PANGLE,IANGLE(定义文字标注属性)
TUNIF,TEMP(定义结构中所有节点的温度)。
【注】适用于均匀温度负载时使用)
/TXTRE,VOLU,NUM,N1,N2,NINC(为体选择纹理)
/TXTRE,ON(激活纹理显示)
/TYPE,WN,Type(定义显示类型)
TYPE,ITYPE(指定单元类型)
/UDOC,Wind,Class,Key(指定图例栏中图例和文本在窗口中的位置)
UIMP,MAT,Lab1,Lab2,Lab3,VAL1,VAL2,VAL3(求解过程中修改材料特性)
/UNITS,Label,LENFACT,MASSFACT,TIMEFACT,TEMPFACT,TOFFSET,CHARGEFACT,FORCEFACT, HEATFACT(选择单位制)
/USER,WN(冻结自动调整模式产生的视距和焦点)
V,P1,P2,P3,P4,P5,P6,P7,P8(由关键点生成体)
【注】P1~P4为底面四点,P5~P8为顶面四点
VA,A1,A2,…,A9,A10(由面生成体)
VADD,NV1,,NV2,…NV8,NV9(体相加)
VATT,MAT,REAL,TYPE,ESYS(指定体的单元属性)
VCLEAR,NV1,NV2,NINC(清除体单元网格)
VCROSS,LabXR,LabYR,LabZR,LabX1,LabY1,LabZ1,LabX2,LabY2,LabZ2(计算并存储单元表中两矢量的叉积)
VDELE,NV1,NV2,NINC,KSWP(删除体)
【注】KSWP =0删除体但保留体上关键点、1删除体及体上关键点
VDOT,LabXR,LabX1,LabY1,LabZ1,LabX2,LabY2,LabZ2(计算并存储单元表中两矢量的点积)
VDRAG,NA1,NA2,…,NA6,NLP1,NLP2,…,NLP6(将既有面沿一定路径拖拉成体)
VEXT,NA1,NA2,NINC,DX,DY,DZ,RX,RY,RZ(将面按比例延伸生成体)
VGEN, itime, nv1, nv2, ninc, dx, dy, dz, kinc, noelem, imove (移动或拷贝体)
VGLUE,NV1,NV2,…,NV8,NV9(体间相互粘接)
/VIEW,WN,XV,YV,ZV(改变观察方向)
VIMP,VOL,CHGBND,IMPLEVEL(对选定体内的四面体单元进行改进)
VINP,NV1,,NV2,…NV8,NV9(体集两两相交)
VINV,NV1,,NV2,…NV8,NV9(被选体的交集)
VLIST,NV1,NV2,NINC,Label(列表显示体的信息)
【注】Lab=HPT时,显示体上硬点信息。
VLSCALE,NV1,NV2,NINC,RX,RY,RZ,KINC,NOELEM,IMOVE(体缩放)
VMESH,NV1,NV2,NINC(划分体生成体单元)
VOFFST,NRREA,DIST,KNIC(将面偏移生成体)
VOVLAP,NV1,NV2,…,NV8,NV9(体搭接)
VPLOT,NV1,NV2,NINC,DEGEN,SCALE(显示所选体)
VPTN,NV1,NV2,…,NV8,NV9(体间相互分割)
VROTAT,NA1,NA2,NA3,NA4,NA5,NA6,PAX1,PAX2,ARC,NSEG(绕轴旋转生成体)
VSBA,NV,NA,SEPO,KEEPV,KEEPA(体减面)
VSBV,NV1,NV2,SEPO,KEEP1,KEEP2(体减体)
VSBW,NV,SEPO,KEEP(工作平面分离体)
VSCALE,WN,VRATIO,KEY(设置矢量箭头显示比例)
VSWEEP,VNUM,SRCA,TRGA,LSMO(体扫掠)
VSYMM,Ncomp,NV1,NV2,NINC,KINC,NOELEM,IMOVE(体镜像)
VTRAN,KCNTO,NV1,NV2,NINC,KINC,NOELEM,IMOVE(将体转化到另一坐标系)
/VUP,WN,Label(指定整体坐标系参考方向)
WAVES,Wopt,OLDMAX,OLDRMS(进行重排序)
WERASE(清除所有重新排序波表)
WFRONT,KPRNT,KCALC(按当前编号对模型中最大波前数进行预估)
/WINDOW,WN,XMIN,XMAX,YMIN,YMAX,NCOPY(定义窗口大小位置)
/WINDOW,WN,ON(OFF)(激活或关闭窗口)
/WINDOW,WN,DELE(删除窗口)
WMORE,NODE1,NODE2,NINC,ITIME,INC(将更多的节点加入初始波表)
WPAVE,X1,Y1,Z1,X2,Y2,Z2,X3,Y3,Z3(将一组特殊点中心作为工作平面原点)
WPCSYS,WN,KCN(将当前坐标系X-Y平面定义为工作平面)
WPLANE,,XORIG,YORIG,ZORIG,XXAX,YXAX,ZXAX,XPLAN,YPLAN,ZPLAN(通过三个坐标点定义工作平面)
WPOFFS,XOFF,YOFF,ZOFF(偏移工作平面)
WPROTA,THXY,THYZ,THZX(旋转工作平面)
WPSTYL,SNAP,GRSPAC,GRMIN,GRMAX,WPTOL,WPCTYP,GRTYPE,WPOVIS,SNAPANG
(工作平面设置)
【注】SNAP:捕捉增量;GRSPAC:栅格之间距离;GRMIN,GRMAX:栅格区大小;WPTOL:工作平面容差;WPCTYP:坐标系类型(0、1、2);GRTYPE:栅格显示类型(0栅格与坐标系标志、1仅栅格、2仅坐标系);WPOVIS栅格显示(0显示、1不显示);SNAPANG:捕捉角度增量。
WPSTYL,STAT(获得工作平面状态。【注】可用WPSTYL,DEFA恢复缺省状态)
589.WSORT,Lab,KORD,–,Wopt,OLDMAX,OLDRMS(按几何性质对单元重新排序)
WSTART,NODE1,NODE2,NINC,ITIME,INC(定义初始波表)
/XFRM,LAB,X1,Y1,Z1,X2,Y2,Z2(指定动态旋转中心)
/ZOOM,WN,Lab,X1,Y1,X2,Y2(对图形显示窗口的某一区域进行缩放)
有限元基本概念和原理
有限元分析(FEA,Finite Element Analysis)的基本概念是用较简单的问题代替复杂问题后再求解。它将求解域看成是由许多称为有限元的小的互连子域组成,对每一单元假定一个合适的(较简单的)近似解,然后推导求解这个域总的满足条件(如结构的平衡条件),从而得到问题的解。这个解不是准确解,而是近似解,因为实际问题被较简单的问题所代替。由于大多数实际问题难以得到准确解,而有限元不仅计算精度高,而且能适应各种复杂形状,因而成为行之有效的工程分析手段。
有限元是那些集合在一起能够表示实际连续域的离散单元。有限元的概念早在几个世纪前就已产生并得到了应用,例如用多边形(有限个直线单元)逼近圆来求得圆的周长,但作为一种方法而被提出,则是最近的事。有限元法最初被称为矩阵近似方法,应用于航空器的结构强度计算,并由于其方便性、实用性和有效性而引起从事力学研究的科学家的浓厚兴趣。经过短短数十年的努力,随着计算机技术的快速发展和普及,有限元方法迅速从结构工程强度分析计算扩展到几乎所有的科学技术领域,成为一种丰富多彩、应用广泛并且实用高效的数值分析方法。
有限元方法与其他求解边值问题近似方法的根本区别在于它的近似性仅限于相对小的子域中。20世纪60年代初首次提出结构力学计算有限元概念的克拉夫(Clough)教授形象地将其描绘为:“有限元法=Rayleigh Ritz法+分片函数”,即有限元法是Rayleigh Ritz法的一种局部化情况。不同于求解(往往是困难的)满足整个定义域边界条件的允许函数的Rayleigh Ritz法,有限元法将函数定义在简单几何形状(如二维问题中的三角形或任意四边形)的单元域上(分片函数),且不考虑整个定义域的复杂边界条件,这是有限元法优于其他近似方法的原因之一。
对于不同物理性质和数学模型的问题,有限元求解法的基本步骤是相同的,只是具体公式推导和运算求解不同。有限元求解问题的基本步骤通常为:
第一步:问题及求解域定义:根据实际问题近似确定求解域的物理性质和几何区域。
第二步:求解域离散化:将求解域近似为具有不同有限大小和形状且彼此相连的有限个单元组成的离散域,习惯上称为有限元网络划分。显然单元越小(网络越细)则离散域的近似程度越好,计算结果也越精确,但计算量及误差都将增大,因此求解域的离散化是有限元法的核心技术之一。
第三步:确定状态变量及控制方法:一个具体的物理问题通常可以用一组包含问题状态变量边界条件的微分方程式表示,为适合有限元求解,通常将微分方程化为等价的泛函形式。
第四步:单元推导:对单元构造一个适合的近似解,即推导有限单元的列式,其中包括选择合理的单元坐标系,建立单元试函数,以某种方法给出单元各状态变量的离散关系,从而形成单元矩阵(结构力学中称刚度阵或柔度阵)。
为保证问题求解的收敛性,单元推导有许多原则要遵循。 对工程应用而言,重要的是应注意每一种单元的解题性能与约束。例如,单元形状应以规则为好,畸形时不仅精度低,而且有缺秩的危险,将导致无法求解。
第五步:总装求解:将单元总装形成离散域的总矩阵方程(联合方程组),反映对近似求解域的离散域的要求,即单元函数的连续性要满足一定的连续条件。总装是在相邻单元结点进行,状态变量及其导数(可能的话)连续性建立在结点处。
第六步:联立方程组求解和结果解释:有限元法最终导致联立方程组。联立方程组的求解可用直接法、选代法和随机法。求解结果是单元结点处状态变量的近似值。对于计算结果的质量,将通过与设计准则提供的允许值比较来评价并确定是否需要重复计算。
简言之,有限元分析可分成三个阶段,前处理、处理和后处理。前处理是建立有限元模型,完成单元网格划分;后处理则是采集处理分析结果,使用户能简便提取信息,了解计算结果。
对称边界条件
有些结构由於具有某些对称性,我们可以在对称面上施加适当的对称边界条件,这样只要建立部分的模型,既省时又省力。
这里讲的对称,不只是几何形状的对称,还包括边界条件、外力施加、材料性质的对称。如果仅有形状对称,其他条件只要有一样不对称,就不能用对称模型求解。
举一个简单的例子。图一是一个含圆孔的平板,边界条件如图所示。我们知道图一的模型上下、左右对称,因此可以简化成图二。
为什麽图二的边界条件可以设为如此?为回答这个问题,我们先将图一分为上下两部分。假设A点与B点是上下两部分相对称的点(图三),则A点与B点的「y」方向位移方向相反,大小相等。想像A与B点越来越靠近对称面,直到几乎重合在一起(图四),此时A、B两点仍可视为对称於对称面,两点的y方向位移仍为大小相等,方向相反。然而此时两点几乎重合,因此A、B两点的位移应该相同。为了满足上面的条件,唯一的可能就是在对称面上y方向位移等於零。图二中的另一边x方向位移等於零,可用相同方法解释。
事实上,图二的两个对称边界,还有一个条件,那就是剪应力xy为零。这个道理的解释方法与位移类似。参考图三,A、B两点由於对称,剪应力也要一模一样;当A、B两点越来越靠近对称面(图四),两点的剪应力就必须满足牛顿第三运动定律(作用力等於反作用力,方向相反),此时唯一的可能就是剪应力为零。因此在对称面上,剪应力为零
总之,在对称面上,垂直於对称面的位移以及作用於对称面上的剪应力皆为
在破坏力学的应用
破坏力学是固体力学的一个分支,这门学科主要在探讨四个部分:
含裂纹结构的受到外力时的应力分布;
含裂纹结构受到多大的外力,裂纹会成长;
结构中,裂纹一旦成长,会往那个方向成长;
各种工程结构抵抗裂纹成长的能力,这部分通常由实验决定。
首先先介绍破坏力学的基本概念。图一为含裂纹结构的破坏模式,共分为第一型(mode I or opening mode)、第二型(mode II or sliding mode)以及第三型(mode III or tearing mode)。
裂纹尖端为奇异点(应力正比於根号 r分之一(stress ~ 1/r^0.5),r为结构中某一点与裂纹尖端的距离),去探讨裂纹尖端(或附近)的应力有多大是没有意义的。有鉴於此,我们引入了一个参数:应力强度因子(stress intensity factor)。应力强度因子共有三个,一般写做KI、KII以及KIII,分别对应到三种不同的破坏模式。
在线弹性破坏力学(linear elastic fracture mechanics,简称LEFM)中,外力的大小与应力强度因子成正比关系。此外,应力强度因子也和几何参数(例如裂纹长度、外力与裂纹的距离等)有关。对於同一结构,应力强度因子可视为含裂纹结构所受外力大小的一个指标;换言之,应力强度因子越大,结构越危险。当应力强度因子超过「破坏韧性(fracture toughness),记为Kc」,裂纹开始成长。破坏韧性通常由制作标准试片由实验求得,不同的材料有不同的破坏韧性,所以破坏韧性可视为材料性质。
至於裂缝会往那个方向成长?有好几个理论可以预测,例如能量释放率理论(energy release rate)、最大周向应力理论(maximum circumferential stress)、J积分理论(J-integral)、应变能密度理论(strain energy density theory)。以上几个理论,除了最大周向应力之外,基本上都和应变能有关,这个部分较为深入,不详述了。
接着我们要探讨有限元素法在破坏力学的应用。前面讲过,探讨裂缝尖端附近的应力有多大是没有意义的。因此我们的重点在於,如何运用有限元素法求得裂纹尖端的应力强度因子。最常见的就是利用「四分之一节点(quarter point)」元素来模拟裂纹尖端,如图二。严格来说,四分之一节点元素不是「三角形元素」,它是由四边形元素退化而成,退化方式如图三。一般我们利用节点上的位移,求得裂纹尖端的开口位移(crack tip opening displacement)来反推得到应力强度因子。MARC及ANSYS皆用四分之一节点元素来算应力强度因子。
以上所说的四分之一节点元素,其最原始的统御方程式和弹性力学并没有两样,只是将其中几个节点变了位置,因此形状函数也跟着改变,最後得到一个结果:这种元素内部的应力呈现根号r分之一的奇异性,因此可用来模拟裂纹。
另外有些软体,例如NASTRAN,则利用另一种特殊元素,共有18个节点,如图四。这种元素最原始的统御方程式,和一般的弹性力学稍有不同,它一开始就假设应力正比於根号 r分之一。
不管如何,裂纹尖端的元素与一般的元素不同,裂纹尖端网格的大小,由算出来的应力强度因子准不准来决定。如何知道准不准?可建立一个简单模型,有理论解的题目,即可比较,例如图四。以平面问题为例,若采用四分之一节点元素,元素边长大约在裂纹长度的百分之六左右;若采用NASTRAN的特殊元素,其长度更只要裂纹长度的百分之十至二十即可。当然以上的准则并非万用,不同的题目最好还是多试几个网格粗细,再做决定。
方向
元素有方向性。有些元素如果方向搞错,则跑不出来;有些虽然跑得出来,但结果却有问题。 以平面元素(包括平面应变、平面应力以及轴对称)为例,四个节点的编号必须是「逆时钟方向」。如果是顺时钟方向,则在有限元素的定义中,这个元素的面积小於零,就跑不出来了。
板(或壳)元素也有方向性。板元素的方向由右手定则决定,也就是元素的正向方向,由元素编号顺序的方向决定,如此也定义了板元素的上表面或下表面。当板受到bending时,上表面和下表面分别会受到张应力以及压应力(当然也可能颠倒,视bending的方向而定)。图一为两个相邻的板元素受到bending作用,元素1的方向和元素2的方向不同,因此它们的上、下表面也不一致(参考图一)。在受到图一所示的bending时,元素1的上表面为张应力,下表面为压应力;元素2的上表面则为压应力,下表面为张应力。有些有限元素软体(如MARC)在计算板元素的「上表面」节点应力时,是根据两个元素的「上表面」高斯点上应力经外差到节点上再平均求得,但注意图一中两个元素的上下表面定义颠倒,所以平均後的上表面节点应力就很不准了。
再解释清楚一点,我们先将两个元素分开,同时画上厚度,比较清楚,如图二,图一中的A点在图二则可分为A1至A4四点。如果单独看元素1,A1点应该受到张应力,A2为压应力;单独看元素2,A3为张应力,A4为压应力。当有限元素法在计算A点「上表面」应力时,是将元素1以及元素2的「上表面」应力平均(即将A1及A4的应力平均),本来应该上表面是张应力的,被这麽一搞,应力下降,误差就来了。
一般有限元素法软体,在做automesh时,元素的方向与surface(或area)的方向一致,因此同一个surface(或area)做出来的网格,方向都相同。但如果两个surface(或area)相邻,就要特别注意方向有没有一致。如果是「手动」建立网格,更要特别注意建立出来的元素的方向,以免造成困扰。
形状函数
形状函数在有限元素法中是非常重要的一个概念,它定义了元素内部位移的分布。以一维线性元素为例(图一) 此种元素有两个节点,分别以节点1及2表示。节点1的座标为xi = -1,节点2则为xi = -2。元素中有几个节点,就有几个形状函数。因此我们有两个形状函数。形状函数有个特点:考虑第n个形状函数,若代入第n个节点的座标,其函数值为1;若代入其他节点的座标,则函数值为0。由於图一为两节点元素,所以形状函数为线性,其表示式为:其他种类的元素的形状函数,一般有限元素法书籍(或商用有限元素法软体的使用手册)皆有提到,这里不多写了。
元素的特性主要由形状函数所主宰。例如三节点之三角形元素,沿着元素内部任一方向,位移皆为线性分布。为什麽我们知道它是线性分布?只要查查该元素的形状函数立刻就明白了。由於应变为位移的对空间的一次微分,所以在三节点三角形元素内部,应变是保持不变的,因此该元素又称为等应变元素(constant strain triangle)。应变既然不变,应力也不变。同样道理,四节点金字塔3D元素,应力及应变也是保持不变的。而四边形四节点元素,沿着边长,位移为线性,但若沿着其他方向,位移则为二次曲线分布。
当我们充分掌握各种元素的特性後,利用有限元法解题,就比较能够掌握下列几个问题:
该用哪种元素?
网格大小?
预期会有何种结果?
如何判读结果?
当然要回答以上的问题,不止要对元素特性有充分了解,还要对材料力学、弹性力学等基础理论有一定的熟悉。以2D元素模拟工程梁问题为例。选用哪种元素比较好?我们知道梁的问题,基本上是受到bending作用,因此轴向应力沿着横向的分布为线性。因此我们不会选用三节点三角形元素以及四节点四边形元素,采用八节点四边形元素是比
较恰当的选择。当然我们可以选择线性元素,然後网格密一点。这对简单几何形状的模型而言,虽然行得通;但若针对较复杂的模型,在前处理建网格时,会比较花时间。一般而言,要看使用者的情形而定。
针对同一个题目,使用相同大小的网格,但一个采用低阶元素,另一个则采用高阶元素。一般而言,低阶元素的模型位移较小,高阶元素的位移较大。因为高阶元素采用较高阶的形状函数,元素在变形方面,显得较为自由,容易变形。这个特性也会反映在模态分析上。一般而言,越硬的东西,自然频率较高,因此,高阶元素的模型,算出来的自然频率较高。以上两个模型,若元素数量够多,基本上结果不会有太大差异(即元素够多,计算结果收敛)。
除此之外,形状函数还会影响到分布力、体力(body force)在有限元素法中的输入、高斯积分等
Shear Locking 误差
Shear Locking 这种误差, 许多有限元素的初学者可能没有听过,
但是, 若你的 FEM model 没做好, 或元素使用不当,
这种误差就会出现, 而且答案差了十万八千里.
最糟的是, 可能算出了错误答案还不知道
shear locking
shear locking 是 FEM 造成的数值误差, 发生於细长结构的分析(尤其bending),
现象: 算出 之 shear strain energy 过大, 大到不合理.
细长结构bending 之正常现象应是: shear strain energy << bending strain energy.
例: 使用low order linear element (如4-node plane element, 8-node solid element)於 bending 分析, 会发生以上现象.
这是因为 low order element 变形会 overstiffness.
相对的, 使用high order element (有 mid-side node) 较不会发生 shear locking.
但若 high order element 的 aspect ratio 过大, 仍有可能 shear locking.
解决对策:
reduced integration . (用於 low order 或 high order element )
incompatible mode (extra shape function, ANSYS 有使用)
用於 low order element .
例如 ansys 的 plane 42 & solid 45 element,
内定为 incompatible mode , 以防止 shear locking.
所以, 於 细长梁厚度方向, 使用 少量 的 plane 42来分析bending cantilever beam , 也能得到好的结果. reduced integration 用於 low order element , 会有 hourglass 困扰,
所以 incompatible mode 用於 low order element 是较好的方案, 据书上的说法, 以板元素来说, 薄板理论 (K-L plate) 因为没考虑 shear stress, 所以无 shear locking 问题.
但厚板理论 (Mindlin plate) 有考虑 shear stress进来, 所以可能会有 shear locking 发生, 解决方法(可能)同上
ANSYS小知识(经典)
1.绘制等值线
最近俺写文章,期刊上大都不用彩色,所以打出的云图一片模糊,无法识别,这时候可以选择出等值线图,但是等值线图也是彩色的,如何把它转成黑白的呢?开始是抓图后用Photoshop处理,太麻烦,ansys自己行不行呢?小弟琢磨了一阵,终于弄出来了,现贴出来给大家参考一下。
1 用命令jpgprf,500,100,1将背景变为白色;
2 plotctrls>device option中,把vector mode改为on,画出等值线图;
3 plotctrls>style>contour>contour labeling, 将key vector mode contour labels设为on every Nth ele,对N输入一个数值,值越大,图中的label越少;
4 plotctrls>style>colors>contour colors,将所有的系列都改为黑色;
5 如果不喜欢ansys给出的MX,MN标志,可以用plotctrls>window controls>window options把它们去掉,将MINM 后的Mix-Min Symbols改为off就可以了。
这时候,一幅清晰的黑白等值线图就出来了。
2.改彩色为灰度显示
可以自己合成颜色,我自己调配的一个颜色文件给你看看,在默认的9段表示时,可以达到灰度由白到黑,当然你还可以调整使它由黑到白(文件中的4到12颠倒即可)。plotctrls>style>color>load color map加载我上面提供的附件即可
*** COLOR MAP CREATED FOR THE WIN32 DRIVER ***
0 0 0 0
1 7 0 60
2 100 0 100
3 62 0 100
4 95 95 95
5 88 88 88
6 80 80 80
7 70 70 70
8 60 60 60
9 51 50 50
10 41 41 41
11 30 30 30
12 20 20 20
13 62 62 62
14 78 78 78
15 100 100 100
3.黑白等值线
对体和面来说,ANSYS默认的结果输出格式是云图格式,而这种彩色云图打印为黑白图像时对比很不明显,无法表达清楚,这对于发表文章来说是非常不便的。发文章所用的结果图最好是等值线图,并且最好是黑白的等值线图。笔者原来进行这项工作时一般借用photoshop等第三方软件,很麻烦,并且效果不好。现通过摸索,发现通过灵活运用ansys本身也能实现这项功能。现将步骤写给大家,感谢simwe对我的帮助。
(1)将要输出的结果调出,这时为彩色云图;
(2)将云图转换为等值线图的形式
GUI:plotCtrls—>Device Options—>[/DEVI]中的vector mode 选为on
命令:/DEVICE,VECTOR,1
这时结果为彩色等值线,若直接输出,打印为黑白图像时仍然不清晰,为此需进行以下几步将图像转换为黑白形式;
(3) 将背景变为白色
命令:jpgprf,500,100,1
/rep
(4)对等值线中的等值线符号(图中为A,B,C等)的疏密进行调整
GUI:plotCtrls—>Style—>Contours—> Contours Labeling 在Key Vector mode contour label 中选中on every Nth elem,然后在N= 输入框中输入合适的数值,例如5,多试几次,直到疏密合适
命令:/clabel,1,5
(5)将彩色等值线变为黑色
GUI:plotCtrls—>Style—>Colors—>Contours Colors 将Items Numbered 1,Items Numbered 2等复选框中的颜色均选为黑色,图像即可变为黑白等值线图像
命令:/color,cntr,whit,1 等等
(6)最后一步:出图
GUI:plotCtrls—>Capture Image
希望对大家能有所帮助。