ansys命令流汇总

VSEL, Type, Item, Comp, VMIN, VMAX, VINC, KSWP

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，MAT，……

MP，NUXY，MAT，……

定义 DP 材料单元表（这里不考虑温度）：TB，DP，MAT 进入单元表并编辑添加单元表：TBDATA，1，C

TBDATA，2，ψ

TBDATA，3，……

如定义：EX=1E8，NUXY=0.3，C=27，ψ=45 的命令如下： MP，EX，1，1E8

MP，NUXY，1，0.3

TB，DP，1

TBDATA，1，27

TBDATA，2，45  这里要注意的是，在前处理的᳔初，要将角度单 位转化到“度”，即命令：*afun,deg

Type，是选择的方式，有选择（s）,补选（a），不选(u)，全选(all)、

反选（inv）等,其余方式不常用

Item, Comp 是选取的原则以及下面的子项

如 volu 就是根据实体编号选择，

    loc 就是根据坐标选取，它的 comp 就可以是 实体的某方向坐标！ 其余还有 材料类型、实常数等

MIN, VMAX, VINC，这个就不必说了吧！

,例：vsel,s,volu,,14

vsel,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 !设定牛顿－拉夫森选项    u /grid, key

ESTIF,... !设定非缺省缩减因子（可选）    key: “0” 或“off” 无网络

ESEL,... !选择在本载荷步中将不激活的单元    “1”或“on” xy 网络

EKILL,... !不激活选择的单元    “2”或“x” 只有 x 线

ESEL,S,LIVE !选择所有活动单元    “3”或“y” 只有 y 线

NSLE,S !选择所有活动结点    u xvar, n

NSEL,INVE !选择所有非活动结点（不与活动单元相连的结点）    n: “0”或“1” 将 x 轴作为时间轴

D,ALL,ALL,0 !约束所有不活动的结点自由度（可     选）    “n” 将 x 轴表示变量“n”

NSEL,ALL !选择所有结点    “-1”  ？

ESEL,ALL !选择所有单元    u /axlab, axis, lab 定义轴线的标志

D,... !施加合适的约束    axis: “x”或“y”

F,... !施加合适的活动结点自由度载荷    lab: 标志，可长达 30 个字符

SF,... !施加合适的单元载荷    u plvar, nvar, nvar2, ……,nvar10 画出要显示的变量（作为纵坐标）

BF,... !施加合适的体载荷

SAVE

SOLVE    u rforce, nvar, node, item, comp, name 指定待存储的节点力数据

请    参    阅    nvar: 变量号

TIME,NLGEOM,NROPT,ESTIF,ESEL,EKILL,NSLE,NSEL,D,F,SF 和    node: 节点号

BF 命令得到更详细的解释。    item comp

? 后继载荷步    F x, y.z

在后继载荷步中，用户可以随意杀死或重新激活单元。象上面提到    M x, y,z

的，要正确的施加和删除约束和结点载荷。    name: 给此变量一个名称，8 个字符

用下列命令杀死单元：    u add, ir, ia,ib,ic,name,--,--,facta, factb, factc

Command:EKILL    将 ia,ib,ic 变量相加赋给 ir 变量

GUI: Main Menu>Solution>-Load Step Opts-Other>Kill Elements    ir, ia,ib,ic：变量号

用下列命令重新激活单元：    name: 变量的名称

Command: EALIVE

GUI: Main Menu>Solution>-Load Step Opts-Other>Activate Elem    NGEN,ITIME,INC,NODE1,NODE2,NINC,DX,DY,DZ,SPACE

!第二个（或后继）载荷步：    是一个节点复制命令,

TIME,...    它是将一组节点在现有坐标系统下复制到其它位置。

ESEL,...    ITIME: 复制的次数，包含自己本身。

EKILL,... !杀死选择的单元    INC: 每次复制节点时节点号码的增加量。

ESEL,...    NODE1,NODE2,NINC: 选取要复制的节点，即要对哪些节点进行复

EALIVE,... !重新激活选择的单元    制。

    DX,DY,DZ: 每次复制时在现有坐标系统下，几何位置的改变量。

FDELE,... !删除不活动自由度的结点载荷    SPACE:间距比,是᳔后一个尺寸和第一个尺寸的比值。

D,... !约束不活动自由度

    Fini(退出四大模块，回到 BEGIN 层)

F,... !在活动自由度上施加合适的结点载荷    /cle (清空内存，开始新的计算)

DDELE,... !删除重新激活的自由度上的约束    1．  定义参数、数组，并赋值.

SAVE    2．    /prep7(进入前处理)

SOLVE    定义几何图形：关键点、线、面、体

定义几个所关心的节点，以备后处理时调用节点号。

设材料线弹性、非线性特性

设置单元类型及相应 KEYOPT

设置实常数

设置网格划分，划分网格

根据需要耦合某些节点自由度

定义单元表

存盘

3．/solu

加边界条件

设置求解选项

定义载荷步

求解载荷步

4./post1（通用后处理）

5./post26  （时间历程后处理）

6.PLOTCONTROL 菜单命令

7.参数化设计语言

8.理论手册

Fini(退出四大模块，回到 BEGIN 层)

/cle (清空内存，开始新的计算)

1  定义参数、数组，并赋值.

u  dim, par, type, imax, jmax, kmax, var1, vae2, var3 定义数组

par: 数组名

    type：   array 数组，如同 fortran,下标᳔小号为 1，可以多达三维 （缺省）

char 字符串组（每个元素᳔多 8 个字符）

table

imax,jmax, kmax 各维的᳔大下标号

    var1,var2,var3 各维变量名，缺省为 row,column,plane(当 type 为 table 时)

2    /prep7(进入前处理)

2.1  定义几何图形：关键点、线、面、体

u  csys,kcn

kcn , 0 迪卡尔 zuobiaosi

1 柱坐标

2 球

4 工作平面

5 柱坐标系（以 Y 轴为轴心）

n 已定义的局部坐标系

u  numstr, label, value 设置以下项目编号的开始

node

elem

kp

line

area

volu

注意：vclear, aclear, lclear, kclear 将自动设置节点、单元开始号为᳔ 高号，这时 如需要自定义起始号，重发 numstr

u  K, npt, x,y,z, 定义关键点

Npt：关键点号，如果赋 0，则分配给᳔小号

u  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 u  A, P1, P2, ……… P18 由关键点生成面

u  AL, L1,L2, ……,L10 由线生成面

面的法向由 L1 按右手法则决定，如果 L1 为负号，则反向。（线需 在某一平面内坐标值固定的面内）

u  vsba, nv, na, sep0,keep1,keep2 用面分体

u  vdele, nv1, nv2, ninc, kswp 删除体

kswp: 0 只删除体

1 删除体及面、关键点（非公用）

u  vgen, itime, nv1, nv2, ninc, dx, dy, dz, kinc, noelem, imove 移动或 拷贝体

itime: 份数

nv1, nv2, ninc：拷贝对象编号

dx, dy, dz  ：位移增量

kinc: 对应关键点号增量

noelem,：0：同时拷贝节点及单元

1：不拷贝节点及单元

imove：  0：拷贝体

1：移动体

u  cm, cname, entity 定义组元，将几何元素分组形成组元

cname: 由字母数字组成的组元名

    entity: 组元的类型（volu, area, line, kp, elem, node） u  cmgrp, aname, cname1, ……,cname8 将组元分组形成组元集合

aname: 组元集名称

    cname1……cname8: 已定义的组元或组元集名称 u  cmlist,name

u  cmdele,name 

u  cmplot, label1

2.2  定义几个所关心的节点，以备后处理时调用节点号。 u  n,node,x,y,z,thxy, thyz, thzx 根据坐标定义节点号

如果已有此节点，则原节点被重新定义，一般为᳔大节点号。 2.3  设材料线弹性、非线性特性

u  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 次项的 系数

u  Tb, lab, mat, ntemp,npts,tbopt,eosopt 定义非线性材料特性表 Lab: 材料特性表之种类

Bkin: 双线性随动强化

Biso: 双线性等向强化

Mkin: 多线性随动强化(᳔多 5 个点)

Miso: 多线性等向强化（᳔多 100 个点）

Dp: dp 模型

Mat: 材料号

Ntemp: 数据的温度数

对于 bkin: ntemp 缺省为 6

miso: ntemp 缺省为 1，᳔多 20

biso: ntemp 缺省为 6，᳔多为 6

dp: ntemp, npts, tbopt 全用不上

Npts: 对某一给定温度数据的点数

u  TBTEMP,temp,kmod 为材料表定义温度值

temp: 温度值

kmod: 缺省为定义一个新温度值

如果是某一整数，则重新定义材料表中的温度值

注意：此命令一发生，则后面的 TBDATA 和 TBPT 均指此温度，应 该按升序

    若  Kmod  为  crit,  且  temp  为空 ，则其后的  tbdata  数据为 solid46,shell99,solid191 中所述破坏准则

    如果 kmod 为 strain,且 temp 为空，则其后 tbdata 数据为 mkin 中 特性。

u  TBDATA, stloc, c1,c2,c3,c4,c5,c6

给当前数据表定义数据（配合 tbtemp,及 tb 使用）

stloc: 所要输入数据在数据表中的初始位置，缺省为上一次的位置 加 1

    每重新发生一次 tb 或 tbtemp 命令上一次位置重设为 1， （发生 tb 后第一次用空闲此项，则 c1 赋给第一个常数） u  tbpt, oper, x,y 在应力-应变曲线上定义一个点

oper: defi 定义一个点

dele 删除一个点

x,y：坐标

2.4  设置单元类型及相应 KEYOPT

u  ET, itype, ename, kop1……kop6, inopr 设定当前单元类型 Itype：单元号

Ename：单元名设置实常数

u  Keyopt, itype, knum, value

itype: 已定义的单元类型号

knum: 单元的关键字号

value: 数值

注意：如果  ,则必须使用 keyopt 命令，否则也可在 ET 命令中输入 2.5  设置网格划分，划分网格

2.5.1  映射网格划分

1.面映射网格划分

条件：a. 3 或 4 条边

b.面的对边必须划分为相同的单元或其划分与一个过渡形网格的划 分相匹配

    c. 该面如有 3 条边，则划分的单元不必须为偶数，并且各边单 元数相等

d. mahkey

e. mshpattern

* 如果多于四条边，可将线合并成 Lcomb

可用 amap 命令，先选面，再选 4 个关键点即可

* 指定面的对边的分割数，以生成过渡映射四边形网格，只适用于 有四条边的面？

2. 体映射网格划分

（1）若将体划分为六面体单元，必须满足以下条件     a. 该体的外形为块状（六面体）、楔形或棱形（五面体）、四面 体

    b. 对边必须划分为相同的单元数，或分割符合过渡网格形式     c. 如果体是棱形或四面体，三角形面上的单元分割数必须是偶 数

（2） 当需要减少围成体的面数以进行映射网格划分时，可以对面 相加或连接。如果连接而有边界线，线也必须连接在一起。 （3）体扫掠生成网格 

步骤：

    a. 确定体的拓扑是否能够进行扫掠。侧面不能有孔；体内不能 有封闭腔;源面与目标面必须相对

b. 定义合适的单元类型

c. 确定扫掠操作中如何控制生成单元层的数目 lesize

d. 确定体的哪一个边界面作为源面、目标面

e. 有选择地对源面、目标面和边界面划分网格

3. 关于连接线和面的一些说明

连接仅是映射网格划分的辅助工具

4. 用 desize 定义单元尺寸时单元划分应遵守的级别

高：lesize

kesize

esize

desize

用 smartzing 定义单元尺寸时单元划分应遵守的级别 高：lesize

kesize

smartsize

u  LESIZE,NL1,Size, Angsiz,ndiv,space,kforc,layer1,layer2,kyndiv

为线指定网格尺寸

NL1: 线号，如果为 all,则指定所有选中线的网格。

Size: 单元边长，（程序据 size 计算分割份数，自动取整到下一个整 数）？

Angsiz: 弧线时每单元跨过的度数？

Ndiv: 分割份数

Space: “+”:  ᳔后尺寸比᳔先尺寸

“-“: 中间尺寸比两端尺寸

free: 由其他项控制尺寸

kforc 0: 仅设置未定义的线，

1：设置所有选定线，

2：仅改设置份数少的，

3：仅改设置份数多的

kyndiv: 0，No,off 表示不可改变指定尺寸

1，yes,on 表示可改变

u  ESIZE,size,ndiv 指定线的缺省划分份数

    （已直接定义的线，关键点网格划分设置不受影响） u  desize, minl, minh,…… 控制缺省的单元尺寸

minl: n 每根线上低阶单元数（缺省为 3）

defa 缺省值

stat 列出当前设置

off 关闭缺省单元尺寸

minh: n 每根线上（高阶）单元数（缺省为 2）

u  mshape, key, dimension 指定单元形状

key: 0 四边形（2D），六面体（3D）

1 三角形  (2D), 四面体(3D)

Dimension: 2D 二维

3D 三维

u  smart,off 关闭智能网格

u  mshkey, key 指定自由或映射网格方式

key: 0 自由网格划分

1 映射网格划分

    2 如果可能的话使用映射，否则自由（即使自由 smartsizing 也不管用了）

u  Amesh, nA1,nA2,ninc 划分面单元网格

nA1,nA2,ninc 待划分的面号，nA1 如果是 All,则对所有选中面划分 u  SECTYPE, ID, TYPE, SUBTYPE, NAME, REFINEKEY

定义一个截面号，并初步定义截面类型

ID: 截面号

TYPE: BEAM:定义此截面用于梁

SUBTYPE: RECT 矩形

CSOLID:圆形实心截面

CTUBE: 圆管

I: 工字形

HREC: 矩形空管

ASEC: 任意截面

MESH: 用户定义的划分网格

NAME: 8 字符的截面名称（字母和数字组成）

REFINEKEY: 网格细化程度：0~5（对于薄壁构件用此控制，对于 实心截面用 SECDATA 控制）

u  SECDATA, VAL1, VAL2, …….VAL10 描述梁截面 说明：对于  SUBTYPE=MESH,  所需数据由  SECWRITE  产生， SECREAD 读入

u  SECNUM,SECID 设定随后梁单元划分将要使用的截面编号 u  LATT, MAT, REAL, TYPE, --, KB, KE, SECNUM

为准备划分的线定义一系列特性

MAT: 材料号

REAL: 实常数号

TYPE: 线单元类型号

KB、KE: 待划分线的定向关键点起始、终止号 SECNUM: 截面类型号

u  SECPLOT,SECID,MESHKEY 画梁截面的几何形状及网格划分

SECID:由 SECTYPE 命令分配的截面编号

MESHKEY：0：不显示网格划分

1：显示网格划分

u   /ESHAPE, SCALE 按看似固体化分的形式显示线、面单元 SCALE: 0:简单显示线、面单元

1：使用实常数显示单元形状 

u  esurf, xnode, tlab, shape 在已存在的选中单元的自由表面覆盖产

生单元

xnode: 仅为产生 surf151 或 surf152 单元时使用

tlab: 仅用来生成接触元或目标元

    top 产生单元且法线方向与所覆盖的单元相同，仅对梁或 壳有效，对实体单元无效

    Bottom 产生单元且法线方向与所覆盖的单元相反，仅对 梁或壳有效，对实体单元无效

Reverse 将已产生单元反向

Shape: 空 与所覆盖单元形状相同

Tri 产生三角形表面的目标元

注意：选中的单元是由所选节点决定的，而不是选单元，如同将压 力加在节点上而不是单元上

u  Nummrg,label,toler, Gtoler,action,switch 合并相同位置的 item label: 要合并的项目

node: 节点，  Elem,单元，kp: 关键点（也合并线，面及点） mat: 材料，type: 单元类型，Real: 实常数

    cp：耦合项，CE：约束项，CE: 约束方程，All：所有项 toler: 公差

Gtoler：实体公差

Action: sele 仅选择不合并

空 合并

switch: 较低号还是较高号被保留（low, high）

注意：可以先选择一部分项目，再执行合并。如果多次发生合并命 令，一定要先合并节点，再合并关键点。合并节点后，实体荷载不 能转化到单元，此时可合并关键点解决问题。

u  Lsel, type, item, comp, vmin, vmax, vinc, kswp 选择线

type: s 从全部线中选一组线

r 从当前选中线中选一组线

a 再选一部线附加给当前选中组

au

none

u(unselect)

inve: 反向选择

item: line 线号

loc 坐标

length 线长

comp: x,y,z

kswp: 0 只选线

1 选择线及相关关键点、节点和单元

u  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”仅用绝对值

u  NSLL,type, nkey 选择与所选线相联系的节点

u  nsla, type, nkey: 选择与选中面相关的节点

type：s 选一套新节点

r 从已选节点中再选

a 附加一部分节点到已选节点

u 从已选节点中去除一部分

nkey: 0 仅选面内的节点

    1 选所有和面相联系的节点（如面内线，关键点处的节点） u  esel, type, item, comp, vmin, vmax, vinc, kabs 选择一组单元 Type: S: 选择一组单元（缺省）

R: 在当前组中再选一部分作为一组

A: 为当前组附加单元

U: 在当前组中不选一部分单元

All: 选所有单元

None: 全不选

Inve: 反向选择当前组（？）

Stat: 显示当前选择状态

Item： Elem: 单元号

Type: 单元类型号

Mat: 材料号

Real: 实常数号

Esys: 单元坐标系号

u  ALLSEL, LABT, ENTITY 选中所有项目

LABT: ALL: 选所有项目及其低级项目

BELOW: 选指定项目的直接下属及更低级项目

ENTITY: ALL: 所有项目（缺省）

VOLU:体 高级

AREA:面

LINE :线

KP:关键点

ELEM:单元 

NODE:节点 低级

u  Tshap,shape 定义接触目标面为 2D、3D 的简单图形

Shape: line:直线

Arc:顺时针弧

Tria:3 点三角形

Quad:4 点四边形

………….

2.6  根据需要耦合某些节点自由度

u  cp, nset, lab,,node1,node2,……node17

nset: 耦合组编号

lab: ux,uy,uz,rotx,roty,rotz

node1-node17: 待耦合的节点号。如果某一节点号为负，则此节点从 该耦合组中删去。如果 node1=all,则所有选中节点加入该耦合组。 注意：1，不同自由度类型将生成不同编号

2，不可将同一自由度用于多套耦合组

u  CPINTF, LAB, TOLER 将相邻节点的指定自由度定义为耦合自 由度

LAB：UX,UY,UZ,ROTX,ROTY,ROTZ,ALL

TOLER: 公差，缺省为 0.0001

说明：先选中欲耦合节点，再执行此命令

2.7  定义单元表

说明：1，单元表仅对选中单元起作用，使用单元表之前务必选择一 种类型的单元

    2，单元表各行为选中各单元，各列为每单元的不同数据 u  ETABLE, LAB, ITEM, COMP 定义单元表，添加、删除单元表某 列

LAB:用户指定的列名（REFL, STAT, ERAS 为预定名称） ITEM: 数据标志（查各单元可输出项目）

COMP: 数据分量标志

2.8  存盘

u  save, fname, ext,dir, slab 存盘

fname : 文件名（᳔多 32 个字符）缺省为工作名

ext: 扩展名（᳔多 32 个字符）缺省为 db

dir: 目录名（᳔多  个字符）缺省为当前

slab: “all” 存所有信息

“model” 存模型信息

“solv” 存模型信息和求解信息

3   /solu

u   /solu 进入求解器

3.1  加边界条件

u  D, node, lab, value, value2, nend, ninc, lab2, lab3, ……lab6 定义节 点位移约束

Node  : 预加位移约束的节点号，如果为 all,则所有选中节点全加约 束，此时忽略 nend 和 ninc.

Lab: ux,uy,uz,rotx,roty,rotz,all

Value,value2: 自由度的数值（缺省为 0）

Nend, ninc: 节点范围为：node-nend，编号间隔为 ninc Lab2-lab6: 将 lab2-lab6 以同样数值施加给所选节点。 注意：在节点坐标系中讨论

3.2  设置求解选项

u  antype, status, ldstep, substep, action

antype: static or 1 静力分析

buckle or 2 屈曲分析

modal or 3 模态分析

trans or 4 瞬态分析

status: new 重新分析（缺省），以后各项将忽略

rest 再分析，仅对 static,full transion 有效

    ldstep: 指定从哪个荷载步开始继续分析，缺省为᳔大的，runn 数（指分析点的᳔后一步）

    substep: 指定从哪个子步开始继续分析。缺省为本目录中， runn 文件中᳔高的子步数

action, continue: 继续分析指定的 ldstep,substep

说明：继续以前的分析（因某种原因中断）有两种类型 singleframe restart: 从停止点继续

需要文件：jobname.db 必须在初始求解后马上存盘

jobname.emat 单元矩阵

jobname.esav  或 .osav :  如果.esav  坏了，将.osav  改

为.esav

    results file: 不必要，但如果有，后继分析的结果也将很 好地附加到它后面

注意：如果初始分析生成了.rdb, .ldhi, 或 rnnn 文件。必须删除再做 后继分析

步骤：  （1）进入 anasys 以同样工作名

（2）进入求解器，并恢复数据库

（3）antype, rest

（4）指定附加的荷载

（5）指定是否使用现有的矩阵（jobname.trl）（缺省重新生成）

kuse: 1 用现有矩阵

（6）求解

multiframe restart:从以有结果的任一步继续（用不着） u  pred,sskey, --,lskey….. 在非线性分析中是否打开预测器

（当有旋转自由度时或使用 solid65 时缺省为 off）

    on 第一个子步后作预测（除非有旋转自由度时或使用 solid65 时缺省为 on）

--    ： 未使用变量区

lskey: off 跨越荷载步时不作预测（缺省）

    on 跨越荷载步时作预测（此时 sskey 必须同时 on） 注意：此命令的缺省值假定 solcontrol 为 on 

u  autots, key 是否使用自动时间步长

key:on: 当 solcontrol 为 on 时缺省为 on

off: 当 solcontrol 为 off 时缺省为 off

    1: 由程序选择（当 solcontrol 为 on 且不发生 autots 命令时 在  .log 文件中纪录“1”

注意：当使用自动时间步长时，也会使用步长预测器和二分步长 u  NROPT, option,--,adptky 指定牛顿拉夫逊法求解的选项 OPTION: AUTO:程序选择

FULL:完全牛顿拉夫逊法

MODI:修正的牛顿拉夫逊法

INIT：使用初始刚阵

UNSYM：完全牛顿拉夫逊法，且允许非对称刚阵

ADPTKY:ON: 使用自适应下降因子

OFF：不使用自适应下降因子

u  NLGEOM，KEY

KEY: OFF:不包括几何非线性（缺省）

ON：包括几何非线性

u  ncnv, kstop, dlim, itlim, etlim, cplim 终止分析选项

kstop: 0 如果求解不收敛，也不终止分析

1 如果求解不收敛，终止分析和程序（缺省）

2 如果求解不收敛，终止分析,但不终止程序

dlim：᳔大位移，缺省为 1.0e6

itlim: 累积迭代次数，缺省为无穷多

etlim：程序执行时间（秒），缺省为无穷

cplim：cpu 时间（秒），缺省为无穷

u  solcontrol  ,key1, key2,key3,vtol 指定是否使用一些非线性求解缺 省值

key1: on 激活一些优化缺省值（缺省）

CNVTOL  Toler=0.5%Minref=0.01(对力和弯矩)

NEQIT   ᳔大迭代次数根据模型设定在 15~26 之间

ARCLEN  如用弧长法则用较 ansys5.3 更先进的方法 PRED  除非有 rotx,y,z 或 solid65，否则打开

LNSRCH  当有接触时自动打开

CUTCONTROL  Plslimit=15%, npoint=13

SSTIF  当 NLGEOM,on 时则打开

NROPT,adaptkey  关闭（除非：摩擦接触存在；单元 12,26,48,49,52 存在；当塑性存在且有单元 20,23,24,60 存在）

AUTOS  由程序选择

off 不使用这些缺省值

key2: on 检查接触状态（此时 key1 为 on）

    此时时间步会以单元的接触状态（据 keyopt(7)的假定）为基 础

当 keyopt(2)=on 时，保证时间步足够小

key3: 应力荷载刚化控制，尽量使用缺省值

空：缺省，对某些单元包括应力荷载刚化，对某些不包括（查）

nopl:对任何单元不包括应力刚化

incp:对某些单元包括应力荷载刚化（查）

vtol：

u  outres, item, freq, cname 规定写入数据库的求解信息     item: all 所有求解项

basic 只写 nsol, rsol, nload, strs

nsol 节点自由度

rsol 节点作用荷载

nload 节点荷载和输入的应变荷载（？）

strs 节点应力

freq: 如果为 n,则每 n 步（包括᳔后一步）写入一次

none: 则在此荷载步中不写次项

all: 每一步都写

last: 只写᳔后一步（静力或瞬态时为缺省）

3.3  定义载荷步

u  nsubst, nsbstp, nsbmx, nsbmn, carry 指定此荷载步的子步数     nsbstp: 此荷载步的子步数

    如果自动时间步长使用 autots,则此数定义第一子步的长度； 如果 solcontrol 打开，且 3D 面-面接触单元使用，则缺省为 1-20 步； 如果 solcontrol 打开，并无 3D 接触单元，则缺省为 1 子步；如果 solcontrol 关闭，则缺省为以前指定值；如以前未指定，则缺省为 1） nsbmx, nsbmn：᳔多，᳔少子步数（如果自动时间步长打开）？ u  time, time 指定荷载步结束时间

注意：第一步结束时间不可为“0”

u  f, node, lab, value, value2, nend, ninc 在指定节点加集中荷载     node:节点号

lab: Fx,Fy,Fz,Mx,My,Mz

value: 力大小

value2: 力的第二个大小（如果有复数荷载）

    nend,ninc：在从 node 到 nend 的节点（增量为 ninc）上施加同 样的力

注意：（1）节点力在节点坐标系中定义，其正负与节点坐标轴正向 一致

u  sfa, area, lkey, lab, value, value2 在指定面上加荷载

area: n 面号

all 所有选中号

lkey: 如果是体的面，忽略此项

lab: pres

value: 压力值

u  SFBEAM,  ELEM,  LKEY,  LAB,  VALI,  VALJ,  VAL2I,  VAL2J, IOFFST, JOFFST

对梁单元施加线荷载

ELEM: 单元号，可以为 ALL,即选中单元 

LKEY: 面载类型号，见单元介绍。对于 BEAM188，1  为竖向；2    item  comp  截面数据及分量标志

为横向；3 为切向    S  COMP  X,XZ,YZ 应力分量

VALI,VALJ: I, J 节点处压力值    PRIN  S1,S2,S3 主应力 SINT 应力强度，SEQV 等效应力

VAL2I,VAL2J: 暂时无用    EPTO  COMP  总应变

IOFFST, JOFFST: 线载距离 I, J 节点距离    PRIN  总主应变，应变强度，等效应变

u  lswrite, lsnum 将荷载与荷载选项写入荷载文件中    EPPL  COMP  塑性应变分量

lsnum  ：荷载步文件名的后缀，即荷载步数    PRIN  主塑性应变，塑性应变强度，等效塑性应变

当 stat 列示当前步数

    init 重设为“1”    u  plnsol, item, comp, kund, fact 画节点结果为连续的轮廓线

缺省为当前步数加“1”    item: 项目（见下表）

3.3.1  注意    comp: 分量

1. 尽量加面载，不加集中力，以免奇异点    kund: 0 不显示未变形的结构

2. 面的切向荷载必须借助面单元    1 变形和未变形重叠

3.4  求解载荷步    2 变形轮廓和未变形边缘

u  lssolve, lsmin, lsmax, lsinc 读入并求解多个荷载步    fact: 对于接触的 2D 显示的比例系数，缺省为    1

    lsmin, lsmax, lsinc  ：荷载步文件范围    item  comp  discription

4   /post1（通用后处理）    u  x,y,z,sum  位移

u  set, lstep, sbstep, fact, king, time, angle, nset 设定从结果文件读入    rot  x,y,z,sum  转角

的数据    s  x,y,z,xy,yz,xz  应力分量

    lstep :荷载步数    1，2，3  主应力

    sbstep：子步数，缺省为᳔后一步    Int,eqv  应力 intensity,等效应力

    time： 时间点（如果弧长法则不用）    epeo  x,y,z,xy,yz,xz  总位移分量

    nset： data set number    1,2,3  主应变

u  dscale, wn, dmult 显示变形比例    Int,eqv  应变 intensity,等效应变

    wn: 窗口号（或 all）,缺省为 1    epel  x,y,z,xy,yz,xz  弹性应变分量

    dmult, 0 或 auto : 自动将᳔大变形图画为构件长的 5%    1，2，3  弹性主应变

u  pldisp, kund 显示变形的结构    Int,eqv  弹性 intensity,弹性等效应变

    kund：  0 仅显示变形后的结构    eppl  x,y,z,xy,yz,xz  塑性应变分量

    1 显示变形前和变形后的结构    u  PRNSOL, item, comp 打印选中节点结果

    2 显示变形结构和未变形结构的边缘    item: 项目（见上表）

u   *get, par, node, n, u, x(y,z) 获得节点 n 的 x(y,z)位移给参数 par    comp: 分量

等价于函数 ux(n),uy(n),uz(z)    u  PRETAB, LAB1, LAB2,  ……LAB9 沿线单元长度方向绘单元表

    node(x,y,z): 获得(x,y,z)节点号    数据

    arnode(x,y,z)：获得和节点 n 相连的面    LABn : 空： 所有 ETABLE 命令指定的列名

注意：此命令也可用于/solu 模块    列名： 任何 ETABLE 命令指定的列名

u  fsum, lab, item 对单元之节点力和力矩求和    u  PLLS, LABI, LABJ, FACT, KUND 沿线单元长度方向绘单元表数

    lab: 空 在整体迪卡尔坐标系下求和    据

    rsys 在当前激活的 rsys 坐标系下求和    LABI:节点 I 的单元表列名

    item: 空 对所有选中单元（不包括接触元）求和    LABJ:节点 J 的单元表列名

    cont: 仅对接触节点求和    FACT: 显示比例，缺省为 1

u  PRSSOL, ITEM, COMP 打印 BEAM188、BEAM1 截面结果    kund: 0 不显示未变形的结构

说明：只有刚计算完还未退出 ANSYS 时可用，重新进入 ANSYS    1 变形和未变形重叠

时不可用    2 变形轮廓和未变形边缘

5   /post26  （时间历程后处理） 

u  nsol, nvar, node, item, comp,name    rot 所加的转角约束

在时间历程后处理器中定义节点变量的序号    key: 0 不显示符号

    nvar：变量号（从 2 到 nv（根据 numvar 定义））    1 显示符号

    node: 节点号    2  显示符号及数值

item  comp    u   /SHOW, FNAME, EXT, VECT, NCPL 确定图形显示的设备及其

u  x, y,z    他参数

rot  x, y,z    FNAME: X11:屏幕

u  ESOL, NVAR, ELEM, NODE, ITEM, COMP, NAME 将结果存入    文件名：各图形将生成一系列图形文件

变量    JPEG: 各图形将生成一系列 JPEG 图形文件

    NVAR: 变量号，2 以上    说明：没必要用此命令，需要的图形文件可计算后再输出

ELEM: 单元号    7  参数化设计语言

NODE: 该单元的节点号，决定存储该单元的哪个量，如果空，则    u   *do, par, ival, fval, inc 定义一个 do 循环的开始

给出平均值    par: 循环控制变量

ITEM:    ival, fval, inc：起始值，终值，步长（正，负）

COMP:    u   *enddo 定义一个 do 循环的结束

NAME: 8 字符的变量名， 缺省为 ITEM 加 COMP    u   *if,val1, oper, val2, base: 条件语句

u  rforce, nvar, node, item, comp, name 指定待存储的节点力数据    val1, val2: 待比较的值（也可是字符，用引号括起来）

nvar: 变量号    oper: 逻辑操作（当实数比较时，误差为 1e-10）

node: 节点号    eq, ne, lt, gt, le, ge, ablt, abgt

item  comp    base: 当 oper 结果为逻辑真时的行为

F  x, y.z    lable: 用户定义的行标志

M  x, y,z    stop: 将跳出 anasys

name: 给此变量一个名称，8 个字符    exit: 跳出当前的 do 循环

u  add, ir, ia,ib,ic,name,--,--,facta, factb, factc    cycle: 跳至当前 do 循环的末尾

将 ia,ib,ic 变量相加赋给 ir 变量    then: 构成 if-then-else 结构

    ir, ia,ib,ic：变量号    注意：不允许跳出、跳进一个 do,if 循环至 label 句 ？

    name: 变量的名称    8 理论手册

u   /grid, key    1.方程组解法：(1)直接解法；(2)迭代解法

    key: “0” 或“off” 无网络    (1) 直接解法：a.稀疏矩阵法；b. 波前解法

    “1”或“on” xy 网络    a. 稀疏矩阵法：占内存大，但运算次数少；通过变换刚度矩阵的顺

    “2”或“x” 只有 x 线    序使得非零元素᳔少

“3”或“y” 只有 y 线    b. 波前解法： 占内存小

u  xvar, n    波前是指在还没有一个单元被解完的时候激活的方程数？

    n: “0”或“1” 将 x 轴作为时间轴    (2) 迭代解法：JCG 法；PCG 法；ICCG 法

    “n” 将 x 轴表示变量“n”    JCG 法：可解实数、对称、非对称矩阵

    “-1”  ？    PCG 法：高效求解各种矩阵（包括病态），但仅解实、对称

u   /axlab, axis, lab 定义轴线的标志    矩阵

    axis: “x”或“y”    ICCG 法：类似 JCG,但更强

    lab: 标志，可长达 30 个字符    2. 应变密度，等效应变，应力密度，等效应力

u  plvar, nvar, nvar2, ……,nvar10 画出要显示的变量（作为纵坐标）    （1）应变密度（strain intensity）

u  prvar, nvar1, ……,nvar6 列出要显示的变量    应变密度

6  PLOTCONTROL 菜单命令    是三个主应变

u  pbc, ilem, ……,key, min, max, abs 在显示屏上显示符号及数值    （2）等效应变

    item: u 所加的位移约束    有效泊松比 ：用户由 avprin 命令设定；0（如果不设定）

（3）应力密度(stress intensity)    LABEL=CSG  （公制，厘米、克、秒）

应力密度    LABEL=BFT  （英制，长度=ft 英尺）

(4) 等效应力    LABEL=BIN  （英制，长度=in 英寸）

等效应力

或

若    定义节点：

则有 （弹性状态下）    N,NODE,X,Y,Z,THXY,THYZ,THZX

    NODE：欲建立节点的号码；

    X,Y,Z：节点在目前坐标系统下的坐标位置。

EGEN,ITIME, NINC, IEL1, IEL2, IEINC, MINC, TINC, RINC,    注意：若在圆柱坐标系统下 x,y,z 对应 r,θ,z；在球面系统下对应 r,θ,