有限元分析软件ANSYS命令流中文说明
CommandVSBV, 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: 材料编号(缺省为当前材料号)
c 材料特性值,或材料之特性,温度曲线中的常数项
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
VSEL, Type, Item, Comp, VMIN, VMAX, VINC, KSWP
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 !设定牛顿-拉夫森选项
ESTIF,... !设定非缺省缩减因子(可选)
ESEL,... !选择在本载荷步中将不激活的单元
EKILL,... !不激活选择的单元
ESEL,S,LIVE !选择所有活动单元
NSLE,S !选择所有活动结点
NSEL,INVE !选择所有非活动结点(不与活动单元相连的结点)
D,ALL,ALL,0 !约束所有不活动的结点自由度(可 选)
NSEL,ALL !选择所有结点
ESEL,ALL !选择所有单元
D,... !施加合适的约束
F,... !施加合适的活动结点自由度载荷
SF,... !施加合适的单元载荷
BF,... !施加合适的体载荷
SAVE
SOLVE
请参阅TIME,NLGEOM,NROPT,ESTIF,ESEL,EKILL,NSLE,NSEL,D,F,SF和BF命令得到更详细的解释。
? 后继载荷步
在后继载荷步中,用户可以随意杀死或重新激活单元。象上面提到的,要正确的施加和删除约束和结点载荷。
用下列命令杀死单元:
Command:EKILL
GUI: Main Menu>Solution>-Load Step Opts-Other>Kill Elements
用下列命令重新激活单元:
Command: EALIVE
GUI: Main Menu>Solution>-Load Step Opts-Other>Activate Elem
!第二个(或后继)载荷步:
TIME,...
ESEL,...
EKILL,... !杀死选择的单元
ESEL,...
EALIVE,... !重新激活选择的单元
...
FDELE,... !删除不活动自由度的结点载荷
D,... !约束不活动自由度
...
F,... !在活动自由度上施加合适的结点载荷
DDELE,... !删除重新激活的自由度上的约束
SAVE
SOLVE
u /grid, key
key: “0” 或“off” 无网络
“1”或“on” xy网络
“2”或“x” 只有x线
“3”或“y” 只有y线
u xvar, n
n: “0”或“1” 将x轴作为时间轴
“n” 将x轴表示变量“n”
“-1” ?
u /axlab, axis, lab 定义轴线的标志
axis: “x”或“y”
lab: 标志,可长达30个字符
u plvar, nvar, nvar2, ……,nvar10 画出要显示的变量(作为纵坐标)
u rforce, nvar, node, item, comp, name 指定待存储的节点力数据
nvar: 变量号
node: 节点号
item comp
F x, y.z
M x, y,z
name: 给此变量一个名称,8个字符
u add, ir, ia,ib,ic,name,--,--,facta, factb, factc
将ia,ib,ic变量相加赋给ir变量
ir, ia,ib,ic:变量号
name: 变量的名称
NGEN,ITIME,INC,NODE1,NODE2,NINC,DX,DY,DZ,SPACE
是一个节点复制命令,
它是将一组节点在现有坐标系统下复制到其它位置。
ITIME: 复制的次数,包含自己本身。
INC: 每次复制节点时节点号码的增加量。
NODE1,NODE2,NINC: 选取要复制的节点,即要对哪些节点进行复制。
DX,DY,DZ: 每次复制时在现有坐标系统下,几何位置的改变量。
SPACE:间距比,是最后一个尺寸和第一个尺寸的比值。
Fini(退出四大模块,回到 BEGIN层)
/cle (清空内存,开始新的计算)
1. 定义参数、数组,并赋值.
2. /prep7(进入前处理)
定义几何图形:关键点、线、面、体
定义几个所关心的节点,以备后处理时调用节点号。
设材料线弹性、非线性特性
设置单元类型及相应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: 材料编号(缺省为当前材料号)
c 材料特性值,或材料之特性,温度曲线中的常数项
c1-c4: 材料的特性-温度曲线中1次项,2次项,3次项,4次项的系数
u Tb, lab, mat, ntemp,npts,tbopt,eosopt 定义非线性材料特性表
Lab: 材料特性表之种类
Bkin: 双线性随动强化
Bis 双线性等向强化
Mkin: 多线性随动强化(最多5个点)
Mis 多线性等向强化(最多100个点)
Dp: dp模型
Mat: 材料号
Ntemp: 数据的温度数
对于bkin: ntemp缺省为6
mis ntemp缺省为1,最多20
bis 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: 目录名(最多64个字符)缺省为当前
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….. 在非线性分析中是否打开预测器
sskey: off 不作预测(当有旋转自由度时或使用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: AUT程序选择
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为横向;3为切向
VALI,VALJ: I, J节点处压力值
VAL2I,VAL2J: 暂时无用
IOFFST, JOFFST: 线载距离I, J 节点距离
u lswrite, lsnum 将荷载与荷载选项写入荷载文件中
lsnum :荷载步文件名的后缀,即荷载步数
当 stat 列示当前步数
init 重设为“1”
缺省为当前步数加“1”
3.3.1 注意
1. 尽量加面载,不加集中力,以免奇异点
2. 面的切向荷载必须借助面单元
3.4 求解载荷步
u lssolve, lsmin, lsmax, lsinc 读入并求解多个荷载步
lsmin, lsmax, lsinc :荷载步文件范围
4 /post1(通用后处理)
u set, lstep, sbstep, fact, king, time, angle, nset 设定从结果文件读入的数据
lstep :荷载步数
sbstep:子步数,缺省为最后一步
time: 时间点(如果弧长法则不用)
nset: data set number
u dscale, wn, dmult 显示变形比例
wn: 窗口号(或all),缺省为1
dmult, 0或auto : 自动将最大变形图画为构件长的5%
u pldisp, kund 显示变形的结构
kund: 0 仅显示变形后的结构
1 显示变形前和变形后的结构
2 显示变形结构和未变形结构的边缘
u *get, par, node, n, u, x(y,z) 获得节点n的x(y,z)位移给参数par
等价于函数 ux(n),uy(n),uz(z)
node(x,y,z): 获得(x,y,z)节点号
arnode(x,y,z):获得和节点n相连的面
注意:此命令也可用于/solu模块
u fsum, lab, item 对单元之节点力和力矩求和
lab: 空 在整体迪卡尔坐标系下求和
rsys 在当前激活的rsys坐标系下求和
item: 空 对所有选中单元(不包括接触元)求和
cont: 仅对接触节点求和
u PRSSOL, ITEM, COMP 打印BEAM188、BEAM189截面结果
说明:只有刚计算完还未退出ANSYS时可用,重新进入ANSYS时不可用
item comp 截面数据及分量标志
S COMP X,XZ,YZ应力分量
PRIN S1,S2,S3主应力SINT应力强度,SEQV等效应力
EPTO COMP 总应变
PRIN 总主应变,应变强度,等效应变
EPPL COMP 塑性应变分量
PRIN 主塑性应变,塑性应变强度,等效塑性应变 u plnsol, item, comp, kund, fact 画节点结果为连续的轮廓线
item: 项目(见下表)
comp: 分量
kund: 0 不显示未变形的结构
1 变形和未变形重叠
2 变形轮廓和未变形边缘
fact: 对于接触的2D显示的比例系数,缺省为1
item comp discription
u x,y,z,sum 位移
rot x,y,z,sum 转角
s x,y,z,xy,yz,xz 应力分量
1,2,3 主应力
Int,eqv 应力intensity,等效应力
epeo x,y,z,xy,yz,xz 总位移分量
1,2,3 主应变
Int,eqv 应变intensity,等效应变
epel x,y,z,xy,yz,xz 弹性应变分量
1,2,3 弹性主应变
Int,eqv 弹性intensity,弹性等效应变
eppl x,y,z,xy,yz,xz 塑性应变分量
u PRNSOL, item, comp 打印选中节点结果
item: 项目(见上表)
comp: 分量
u PRETAB, LAB1, LAB2, ……LAB9 沿线单元长度方向绘单元表数据
LABn : 空: 所有ETABLE命令指定的列名
列名: 任何ETABLE命令指定的列名
u PLLS, LABI, LABJ, FACT, KUND 沿线单元长度方向绘单元表数据
LABI:节点I的单元表列名
LABJ:节点J的单元表列名
FACT: 显示比例,缺省为1
kund: 0 不显示未变形的结构
1 变形和未变形重叠
2 变形轮廓和未变形边缘
5 /post26 (时间历程后处理)
u nsol, nvar, node, item, comp,name
在时间历程后处理器中定义节点变量的序号
nvar:变量号(从2到nv(根据numvar定义))
node: 节点号
item comp
u x, y,z
rot x, y,z
u ESOL, NVAR, ELEM, NODE, ITEM, COMP, NAME 将结果存入变量
NVAR: 变量号,2以上
ELEM: 单元号
NODE: 该单元的节点号,决定存储该单元的哪个量,如果空,则给出平均值
ITEM:
COMP:
NAME: 8字符的变量名, 缺省为ITEM加COMP
u rforce, nvar, node, item, comp, name 指定待存储的节点力数据
nvar: 变量号
node: 节点号
item comp
F x, y.z
M x, y,z
name: 给此变量一个名称,8个字符
u add, ir, ia,ib,ic,name,--,--,facta, factb, factc
将ia,ib,ic变量相加赋给ir变量
ir, ia,ib,ic:变量号
name: 变量的名称
u /grid, key
key: “0” 或“off” 无网络
“1”或“on” xy网络
“2”或“x” 只有x线
“3”或“y” 只有y线
u xvar, n
n: “0”或“1” 将x轴作为时间轴
“n” 将x轴表示变量“n”
“-1” ?
u /axlab, axis, lab 定义轴线的标志
axis: “x”或“y”
lab: 标志,可长达30个字符
u plvar, nvar, nvar2, ……,nvar10 画出要显示的变量(作为纵坐标)
u prvar, nvar1, ……,nvar6 列出要显示的变量
6 PLOTCONTROL菜单命令
u pbc, ilem, ……,key, min, max, abs 在显示屏上显示符号及数值
item: u 所加的位移约束
rot 所加的转角约束
key: 0 不显示符号
1 显示符号
2 显示符号及数值
u /SHOW, FNAME, EXT, VECT, NCPL 确定图形显示的设备及其他参数
FNAME: X11:屏幕
文件名:各图形将生成一系列图形文件
JPEG: 各图形将生成一系列JPEG图形文件
说明:没必要用此命令,需要的图形文件可计算后再输出
7 参数化设计语言
u *do, par, ival, fval, inc 定义一个do循环的开始
par: 循环控制变量
ival, fval, inc:起始值,终值,步长(正,负)
u *enddo 定义一个do循环的结束
u *if,val1, oper, val2, base: 条件语句
val1, val2: 待比较的值(也可是字符,用引号括起来)
oper: 逻辑操作(当实数比较时,误差为1e-10)
eq, ne, lt, gt, le, ge, ablt, abgt
base: 当oper结果为逻辑真时的行为
lable: 用户定义的行标志
stop: 将跳出anasys
exit: 跳出当前的do循环
cycle: 跳至当前do循环的末尾
then: 构成if-then-else结构
注意:不允许跳出、跳进一个do,if循环至label句 ?
8 理论手册
1.方程组解法:(1)直接解法;(2)迭代解法
(1) 直接解法:a.稀疏矩阵法;b. 波前解法
a. 稀疏矩阵法:占内存大,但运算次数少;通过变换刚度矩阵的顺序使得非零元素最少
b. 波前解法: 占内存小
波前是指在还没有一个单元被解完的时候激活的方程数?
(2) 迭代解法:JCG法;PCG法;ICCG法
JCG法:可解实数、对称、非对称矩阵
PCG法:高效求解各种矩阵(包括病态),但仅解实、对称矩阵
ICCG法:类似JCG,但更强
2. 应变密度,等效应变,应力密度,等效应力
(1)应变密度(strain intensity)
应变密度
是三个主应变
(2)等效应变
有效泊松比 :用户由avprin 命令设定;0(如果不设定)
(3)应力密度(stress intensity)
应力密度
(4) 等效应力
等效应力
或
若
则有 (弹性状态下)
EGEN,ITIME, NINC, IEL1, IEL2, IEINC, MINC, TINC, RINC, CINC, SINC, DX, DY, DZ
单元复制命令是将一组单元在现有坐标下复制到其他位置,
但条件是必须先建立节点,节点之间的号码要有所关联。
ITIME:复制次数,包括自己本身。
NINC: 每次复制元素时,相对应节点号码的增加量。
IEL1,IEL2,IEINC:选取复制的元素,即哪些元素要复制。
MINC:每次复制元素时,相对应材料号码的增加量。
TINC:每次复制元素时,类型号的增加量。
RINC:每次复制元素时,实常数表号的增加量。
CINC:每次复制元素时,单元坐标号的增加量。
SINC:每次复制元素时,截面ID号的增加量。
DX, DY, DZ:每次复制时在现有坐标系统下,节点的几何位置的改变量。
mshape,key,dimension 指定网格化分时单元形状
key: 0 四边形(2D), 六面体(3D)
1 三角形(2D), 四面体(3D)
dimension: 2D 二维
3D 三维
定义局部坐标:
LOCAL,KCN,KCS,XC,YC,ZC,THXY,THYZ,THZX,PAR1,PAR2
KCN:坐标系统代号,大于10的任何一个号码都可以。
KCS:局部坐标系统的属性。
KCS=0 卡式坐标;KCS=1 圆柱坐标;KCS=2 球面坐标;KCS=3 自定义坐标;KCS=4 工作平面坐标;KCS=5 全局初始坐标。
XC,YC,ZC:局域坐标与整体坐标系统原点的关系。
THXY,THYZ,THZX:局域坐标与整体坐标系统X、Y、Z轴的关系。
声明单位:
/UNITS,LABEL
LABEL=SI (公制,米、千克、秒)
LABEL=CSG (公制,厘米、克、秒)
LABEL=BFT (英制,长度=ft英尺)
LABEL=BIN (英制,长度=in英寸)
定义节点:
N,NODE,X,Y,Z,THXY,THYZ,THZX
NODE:欲建立节点的号码;
X,Y,Z:节点在目前坐标系统下的坐标位置。
注意:若在圆柱坐标系统下x,y,z对应r,θ,z;在球面系统下对应 r,θ,?。
定义节点的集中力:
F,NODE,Lab,VALUE,VALUE2,NEND,NINC
NODE:节点号码。
Lab:外力的形式。
Lab=FX,FY,FZ,MX,MY,MZ(结构力学的方向、力矩方向)
=HEAT(热学的热流量)
=AMP,CHRG(电学的电流、载荷)
=FLUX(磁学的磁通量)
VALUE:外力的大小。
NODE,NEND,NINC:选取施力节点的范围和关联,故在建立节点时应先规划节点的号码,以方便整个程序的编辑。
定义作用于元素的分布力:
SFE,ELEM,LKEY,Lab,KVAL,VAL1,VAL2,VAL3,VAL4
ELEM:元素号码。
LKEY:建立元素后,依节点顺序,该分布力定义施加边或面的号码
Lab:力的形式。
Lab=PRES 结构压力
=CONV热学的对流
=HFLUX热学的热流率
VAL1~VAL4:相对应作用于元素边及面上节点的值。
例如:分布力位于编号为1的3d元素、第六个面,作用于此面的四个边上的力分别为:10,20,30,40。
SEF,1,6,PRES,,10,20,30,40
ansys中关于文件读取,保存,及退出程序的命令:
/Filname,fname,key 指定新的工作文件名
fname:文件名及路径,默认为先前设置的工作路径
key: 0 使用已有的log和error文件
1 使用新的log和error,但不删除旧的.
/Title,tile
指定一个标题
/Exit,slab,Fname,Ext,--, 退出程序
Slab: model, 仅保存模型数据文件(默认)
solu 保存模型及求解数据
all, 保存所有的数据文件
nosave, 不保存任何数据文件
/Input,Fname,Ext,--,LIne,log
读入数据文件
Fname,文件名及目录路径,默认为先前设置的工作目录
Ext, 文件扩展名
后面的几个参数一般可以不考虑.
(注): 用此命令时,文件名及目录路径都必须为英文,不能含有中文字符.
/Pbc,item,--,key,min,max,abs 在显示屏上显示符号及数值
item: u, 所加的位移约束
rot, 所加的转角约束
temp 所加的温度荷载
F 所加的集中力荷载
cp 耦合节点显示
ce 所加的约束方程
acel 所加的重力加速度
all 显示所有的符号及数值
key : 0 不显示符号
1 显示符号
2 显示符号及数值
[以上只列出了一些常用的item,详细的可参考帮助文档]
/plopts,vers,0 不在屏幕上显示ansys标记
1. wpoffs,xoff,yoff,zoff
移动工作平面
xoff-x方向移动的距离
yoff-y方向移动的距离
zoff-z方向移动的距离
2.csys,4 激活该局部坐标系
3.wprota,thxy,thyz,thzx
旋转工作平面
thxy-绕z轴旋转
thyz-绕x轴旋转
thzx-绕y轴旋转
4.改变划分网格后的单元
首先:esel,Type, Item, Comp, VMIN, VMAX, VINC, KABS
type中有
s-选择新的单元
r-在所选中的单元中再次选单元
a-再选别的单元
u-在所选的单元中除掉某些单元
all-选中所有单元
none-不选
inve-反选刚才没有被选中的所有单元
stat-显示当前单元的情况
其中
Item, Comp一般系统默认
VMIN-选中单元的最小号
VMAX-选中单元的最大号
VINC-单元号间的间隔
KABS:
0---核对号的选取
1----取绝对值
如:esel
其次:
emodif,IEL, STLOC, I1, I2, I3, I4, I5, I6, I7, I8
改变选中的单元类型为所需要的类型
5. 显示所有单元元素:
/eshape,SCALE
scale :
0--一般地显示面、体单元元素(系统默认)
1--显示所有的元素
如:/eshape,1
6. eplot,all
可以看到所有单元
7.lfillt,NL1, NL2, RAD, PCENT
对两相交的线进行倒圆
NL1-第一条线号
NL2-第二条线号
RAD-圆角半径
PCENT-是否生成关键点,一般为默认
如:lfillt,1,2,0.5
D, NODE, Lab, VALUE, VALUE2, NEND, NINC, Lab2, Lab3, Lab4, Lab5, Lab6 -- 定义节点的自由度约束.
NODE,节点编号,
Lab,自由度编号,如X向,Y向等
VALUE,约束点位移,实部,VALUE2,如果位移为复数,则为虚部
NEND, NINC, ,定义的终止节点编号和节点编号增量
Lab2, Lab3, Lab4, Lab5, 该部分节点的其他自由度编号。同lab
2.LATT, MAT, REAL, TYPE, --, KB, KE, SECNUM -- 定义线的属性,有限元划分用
MAT, REAL, TYPE分别为材料,实常数,单元类型编号,
KB,KE,定义截面的方向关键点,如beam18x系列,默认两关键点一致
SECNUM 定义的截面的编号
3.PLDISP, KUND 显示结构变形图
kund:0 只显示变形后的结构图
1 显示变形后的结构图+变形前的结构图
2 显示变形后的结构图+变形前的结构边界图
4.FLIST, NODE1, NODE2, NINC -- 列表节点力荷载
所列表的节点范围是:从节点编号node1到node2,以NINC的节点增加数
5.DLIST, NODE1, NODE2, NINC -- 列表节点约束.
所列表的节点范围是:从节点编号node1到node2,以NINC的节点增加数
6.LSEL, Type, Item, Comp, VMIN, VMAX, VINC, KSWP -- 选择一组线的子集
Type 定义选择集的类型 可以为
s-选择一个新的子集,默认如此
r-从当前选择子集中选择一部分作为新的子集
a-选择一个新的子集附加到当前选择集上
inve-觉得有时比较重要,对当前子集取数学上集合的逆操作
all-选择全部的线
还有u,none,stat等选项
Item ,comp
一般取item comp
line(材料mat 单元类型type 实常数R) 对应量的编号
loc坐标位置 x,y,z
VMIN, VMAX, VINC,根据Item ,comp取的量,而与之对应的量的数值范围;起始量的数值,终止量的数值,量的增加数值
KSWP
0 仅选择线
1 选择线外还将与线有关的属性比如关键点,单元,节点等一起选中 在ansys下的ls-dyna中编的程序里写入
edwrite,both
可生成d3plot文件,这样可在“独立”的ls-dyna中读入该文件。这是我的经验。
wpcsys,-1,0 将工作平面与总体笛卡尔系对齐
csys,1 将激活坐标系转到总体柱坐标系
antype,static 定义分析类型为静力分析
/post1中的几个命令:
set, lstep, sbstep, fact, king, time, angle, nset 设定从结果文件读入的数据
lstep :荷载步数
sbstep:子步数,缺省为最后一步
time: 时间点(如果弧长法则不用)
nset: data set number
dscale, wn, dmult 显示变形比例
wn: 窗口号(或all),缺省为1
dmult, 0或auto : 自动将最大变形图画为构件长的5%
pldisp, kund 显示变形的结构
kund: 0 仅显示变形后的结构
1 显示变形前和变形后的结构
2 显示变形结构和未变形结构的边缘
PRETAB, LAB1, LAB2, ……LAB9 沿线单元长度方向绘单元表数据
LABn : 空: 所有ETABLE命令指定的列名
列名: 任何ETABLE命令指定的列名
PLLS, LABI, LABJ, FACT, KUND 沿线单元长度方向绘单元表数据
LABI:节点I的单元表列名
LABJ:节点J的单元表列名
FACT: 显示比例,缺省为1
kund: 0 不显示未变形的结构
1 变形和未变形重叠
2 变形轮廓和未变形边缘
etable, lab,item,comp
将单元的某项结果制作成表格,以供pretable命令输出,
lab: 字段名称,自己指定
item: 结果的顶目名称,在每个单元的说明中有(在单元说明表中冒号左边的
comp, 结果项目名称的分量,在单元说明表中冒号右边的
比如将plane42单元的x应力分量制成表
etable,sx,x,x
LACAL,KCN,KCS,XC,YC,ZC,THXY,THYZ,THZX,PAR1,PAR2
定义区域坐标系统,该命令执行后,ANSYS坐标系统自动更改为新建立的坐标系统,故可以定义许多区域坐标系统,以辅助有限元模型的建立。
KCN:该区域坐标系统的确定代号,大于10的任何一个号码都可以。
KCS:该区域坐标系统的属性。0,1,2分别代表卡式坐标,圆柱坐标,球面坐标。
XC,YC,ZC:该区域坐标系统与整体坐标系统原点的关系。
THXY,THYZ,THZX:该区域坐标系统与整体系统X,Y,Z轴的关系
claer,nl1,nl2,lmesh
就是将后面的直线网格化之后的节点和元素都删除
但是共享节点依然存在
mshkey,key
声明是使用自由化网格(key=0)
对应网格(key=1)
或者是混合网格(key=2)
后面两种我因为是新手,所以不大会用,一般都用自由网格~~
关于工作平面:
KWPAVE, P1, P2, P3, P4, P5, P6, P7, P8, P9
把工作平面的中心移动到以上几点的平均点
最多9
如果只选一点,那么就是把工作平面的中心移动到此点
WPOFF, XOFF, YOFF, ZOFF
移动工作平面,注意xoff,yoff,zoff是相对当前点的移动量
而不是整体坐标
WPROT, THXY, THYZ, THZX
旋转工作平面
和上面的一样,是相对当前的工作平面选择一个角度,默认设置是角度为单位
wpstyl
关闭工作平面显示
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)
注意:可以先选择一部分项目,再执行合并。如果多次发生合并命令,一定要先合并节点,再合并关键点。合并节点后,实体荷载不能转化到单元,此时可合并关键点解决问题。
我也感觉和Glue效果一样,但是它有独到的好处的。
numcmp是压缩编号,对计算没有影响的。
一个条件命令
u *if,val1, oper, val2, base: 条件语句
val1, val2: 待比较的值(也可是字符,用引号括起来)
oper: 逻辑操作(当实数比较时,误差为1e-10)
eq, ne, lt, gt, le, ge, ablt, abgt
base: 当oper结果为逻辑真时的行为
lable: 用户定义的行标志
stop: 将跳出anasys
exit: 跳出当前的do循环
cycle: 跳至当前do循环的末尾
then: 构成if-then-else结构
一个循环命令
*do, par, ival, fval, inc 定义一个do循环的开始
par: 循环控制变量
ival, fval, inc:分别为起始值,终值,步长(可正可负)
……
*enddo 定义一个do循环的结束
一个网格划分命令
用desize定义单元尺寸时单元划分应遵守的级别
高:lesize
kesize
esize
desize
用smartzing定义单元尺寸时单元划分应遵守的级别
高:lesize
kesize
smartsize
定义表、数组等真的很好用哦
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
在ANSYS帮助系统中关于*SET 命令的注释下列出了ANSYS中可以使用的数学函数。所有这些数学函数均可以在ANSYS环境中使用,这些数学函数包括:
ABS(X) 求绝对值
ACOS(X) 反余弦
ASIN(X) 反正弦
ATAN(X) 反正切
ATAN2(X,Y) 反正切, ArcTangent of (Y/X) , 可以考虑变量X,Y 的符号
COS(X) 求余弦
COSH(X) 双曲余弦
EXP(X) 指数函数
GDIS(X,Y) 求以X为均值,Y为标准差的高斯分布,在使用蒙地卡罗法研究随机荷载和随机材料参数时,可以用该函数处理计算结果
LOG(X) 自然对数
LOG10(X) 常用对数(以10为基)
MOD(X,Y) 求 X/Y的余数. 如果 Y=0, 函数值为 0
NINT(X) 求最近的整数
RAND(X,Y) 取随机数,其中X 是下限, Y是上限
SIGN(X,Y) 取 X的绝对值并赋予Y的符号. Y>=0, 函数值为|X|, Y<0, 函数值为-|X|,.
SIN(X) 正弦
SINH(X) 双曲正弦
SQRT(X) 平方根
TAN(X) 正切
TANH(X) 双曲正切
esel,s,mat,,1 选择材料号为1的单元
*get,emin,elem,,num,min 获得最小的单元号
*get,emax,elem,,num,max 获得最大的单元号
*DO,I,emin,emax 作循环
*GET,V1,ELEM,I,VOLU 获得单元的体积存到V1的变量中
V=V+V1 求和获得材料1的总体积
*enddo
把一个矩阵的一列加起来的方法
提取当前选择集中的结点总数存入变量aaa1;
提取当前选择集中的结点的最小结点号存入变量aaa2;
定义aaa1×2数组aaa3;
开始循环:
aaa3数组的第一列存储结点号;
aaa3数组的第二列存储Sx;
下一个结点号存入变量aaa2;
循环结束。
/post1
*get,aaa1,node,0,count
*get,aaa2,node,0,num,min
*dim,aaa3,array,aaa1,2
*do,i,1,aaa1
aaa3(i,1)=aaa2
*get,aaa3(i,2),node,aaa2,s,x
aaa2=ndnext(aaa2)
*enddo
L, P1, P2, NDIV, SPACE, XV1, YV1, ZV1, XV2, YV2, ZV2:在两个关键点之间定义一条线。
功能:在当前激活坐标系统下,在两个指定关键点之间生成直线或曲线。
P1,P2:线的起点和终点。
NDIV:这条线的单元划分数。一般不用,指定单元划分数推荐用 LESIZE。这里需要说明一下:如果你的模型相对规则,为了得到高质量的网格,不妨在划线的时候指定单元划分数,这样,既方便又能按照自己的意愿来分网。
SPACE:间隔比。通常不用,指定间隔比推荐使用命令LESIZE。
说明: 线的形状由激活坐标系决定,直角坐标系中将产生一条直线,柱坐标系中,随关键的坐标不同可能产生直线,圆弧线或螺旋线。
KGEN, ITIME, NP1, NP2, NINC, DX, DY, DZ, KINC, NOELEM, IMOVE:通过一组关键点生成额外的关键点。
ITIME:生成操作总共执行的次数,如果要生成额外的点,该值必须大于 1。
NP1,NP2,NINC:被生成的那组关键点的编号为NP1至NP2,编号增量为NINC(缺省为 1)。
DX,DY,DZ:关键点在激活坐标系下的位置增量(柱坐标系和球面坐标系下要注意坐标的变换)。
KINC:生成的点集与原始点集之间的增量值。如果是0,则指定为最低可用关键点编号。
NOELEM:指定是否单元和节点也随之生成。0,生成;1,不生成。
IMOVE:指定关键点是否被移除或重新定义。0,按照ITIME要求生成额外关键点;1,移除原始关键点到新的位置,保持编号不变(ITIME,KINC,NOELEM被忽略)。
ANSYS的命令流里经常看到一些相对比较固定的代码,这些代码组合在一起构成 ANSYS 的一个操作。比如,通过旋转命令将面生成体的操作,命令流如下:
TYPE,2 !指定生成体的单元类型
EXTOPT,ESIZE,18,0, !指定单元划分数
EXTOPT,ACLEAR,1 !清除面网格
VROTAT,ALL,,,,,,15,16, !绕关键点15,16构成的轴将所划的面网格旋转360
这一组代码在GUI方式下通常由一个或几个对话框组成,还有的命令有时需要同时使用,否则就会出错。例如:想在自己定义的坐标系下编辑有限元模型,需要定义坐标系并激活,然后将所有节点移到当前坐标系中。命令流如下:
WPRO,,,90 !将工作平面绕Y轴旋转90度
CSWPLA,11,1,1,1, !在工作平面原点创建柱坐标系,并激活
NROTAT,ALL !将所有节点旋转到激活坐标系
为了能充分发挥命令流的优势,建议各位把自己常用的代码贡献出来,这样,不仅熟悉了命令,更重要的是,掌握了一种通过命令流来实现的 ANSYS操作。 ! DP材料参数
tb,dp,2
tbdata,,50.0e3,16.7,16.7!
!
mp,ex,2,40.0e6 !粉喷桩复合地基特性
mp,dens,2,1770
mp,nuxy,2,0.38
!
! DP材料参数 !基座粗砂垫层,20cm
tb,dp,3
tbdata,,0.0,25.0,25.0 !
mp,ex,3,46e6
mp,dens,3,1950
mp,nuxy,3,0.30
1.MP,Lab,MAT,C0,C1,C2,C3,C4
定义材料的属性(Material Property),材料属性为固定值时,其值为C0,当随温度变化时,由后四个参数控制。
MAT:对应ET所定义的号码(ITYPE),表示该组属性属于ITYPE。
Lab:材料属性类别,任何元素具备何种属性在元素属性表中均有说明。例如杨氏系数(Lab=EX,EY,EZ),密度(Lab=DENS),泊松比(Lab=NUXY,NUXYZ,NUZX),剪切模数(Lab=GXY,GYZ,GXZ),热膨胀系数(Lab=ALPX,ALPY,ALPZ)等。
2./ANTYPE,Antype,Status
声明分析类型,即欲进行哪种分析,系统默认为静力学分析。
Antype=STATIC or 0 静态分析(系统默认)
BUCKLE or 1屈曲分析
MODAL or 2 振动模态分析
HARMIC or 3 调和外力动和系统
TRANS or 4 瞬时动力系统分析
3.SFBEAM, ELEM, LKEY, Lab,VALI, VALJ, VAL2I, VAL2J, IOFFST, JOFFST
定义在梁元素上的分布力。
ELEM:元素号码。
LKEY:建立元素后,依节点顺序梁元素有四个面,该参为分力所施加的面号。
Lab RES(表示分布压力)。
VALI,VALJ:在I点及J点分布力的值。
4./pnum,label,key
!在有限元模块图形中显示号码。Label=欲显示对象的名称,node节点,elem元素,kp点,line线,area面积,volu体积;key=0为不显示号码(系统默认),=1为显示号码。
5.lesize,nl1,size,angsiz,ndiv,space,kforc,layer1,layer2 !定义所选择线段(nl1,nl1=all为目前所有的线段)进行元素网格化时元素的大小(size),元素的大小可用线段的长度(size)或该条线段要分割的元素数目(ndiv)来确定;space为间距比(最后一段长与最先一段长的比值,正值代表以线段方向为基准,负值以中央为基准,系统默认等间距)。
6.plnsol,item,comp !图标节点的解答。以连续的轮廓线表示。
Item为欲查看何种解答。 Item comp
S x,y,z,xy,yz,xz应力 S 1,2,3 主应力
S eqv,int 等效应力 F x,y,z 结构力
M x,y,z 结构力矩 u x,y,z,sum 位移分量及向量位移
rot x,y,z,sum 旋转位移分量及向量旋转位移 temp 温度
1./UNITS,LABEL
声明单位系统,表示分析时所用的单位,LABEL表示系统单位,如下所示
LABEL=SI (公制,公尺、公斤、秒)
LABEL=CSG (公制,公分、公克、秒)
LABEL=BFT (英制,长度=ft)
LABEL=BIN (英制,长度=in)
2.节点定义
有限元模型的建立是将机械结构转换为多节点和元素相连接,所以节点即为机械结构中一个点的坐标,指定一个号码和坐标位置。在 ANSYS中所建立的对象(坐标系、节点、点、线、面、体积等)都有编号。
相关命令
N,NODE,X,Y,Z,THXY,THYZ,THZX
定义节点,若在圆柱坐标系统下x,y,z对应r,θ,z,在球面系统下对应r,θ,?。
NODE:欲建立节点的号码
X,Y,Z:节点在目前坐标系统下的坐标位置
Menu Paths:Main Menu>Preprocessor>Create>Node>In Active CS
Menu Paths Main Menu>Preprocessor>Create>Node>On Working Plane
NDELE,NODE1,NODE2,NINC
删除在序号在NODE1号NODE2间隔为NINC的所有节点,但若节点已连成元素,要删除节点必先删除元素。例如:
NDELE,1,100,1 !删除从1到100的所有点
NDELE,1,100,99 !删除1和100两个点
Menu Paths:Main Menu>Preprocessor>Delete>Nodes
NPLOT,KNUM
节点显示,该命令是将现有卡式坐标系统下节点显示在图形窗口中,以供使用者参考及查看模块的建立。建构模块的显示为软件的重要功能之一,以检查建立的对象是否正确。有限元型的建立程中,经常会检查 各个对象的正确性及相关位置,包含对象视角、对象号码等,所以图形显示为有限元模型建立过程中不可缺少的步骤。KNUM=0不显示号码,为1显示同时显示节点号
Menu Paths:Utility Menu>plot>nodes
Menu Paths:Utility Menu>plot>Numbering…(选中NODE选项)
NLIST,NODE1,NODE2,NINC,Lcoord,SORT1,SORT2,SORT3
节点列式,该命令将现有卡式坐标系统下节点的资料列示于窗口中(会打开一个新的窗口),使用者可检查建立的坐标点是否正确,并可将资料保存为一个文件。如欲在其它坐标系统下显示节点资料,可以先行改变显示系统,例如圆柱坐标系统,执行命令DSYS,1。
Menu Paths:Utility Menu>List>Nodes
FILL,NODE1,NODE2,NFILL,NSTRT,NINC,ITIME,INC,SPACE
节点的填充命令是自动将两节点在现有的坐标系统下填充许多点,两节点间填充的节点个数及分布状态视其参数而定,系统的设定为均分填满。NODE1,NODE2为欲填充点的起始节点号码及终结节点号码,例如两节点号码为1(NODE1)和5(NODE2),则平均填充三个节点(2,3,4)介于节点1和5之间。
Menu Paths:Main Menu>Preprocessor>Create>Node>Fill between Nds
NGEN,ITIME,INC,NODE1,NODE2,NINC,DX,DY,DZ,SPACE
节点复制命令是将一组节点在现有坐标系统下复制到其它位置。
ITIME: 复制的次数,包含自己本身。
INC: 每次复制节点时节点号码的增加量。
NODE1,NODE2,NINC: 选取要复制的节点,即要对哪些节点进行复制。
DX,DY,DZ: 每次复制时在现有坐标系统下,几何位置的改变量。
Menu Paths:Main Menu>Preprocessor>(-Modeling-)Copy>(-Nodes-)Copy
ET,ITYPE,Ename,KOPT1,KOPT2,KOPT3,KOPT4,KOPT5,KOPT6,INOPR
元素类型(Element Type)为机械结构系统的含的元素类型种类,例如桌子可由桌面平面单元各桌脚梁单元构成,故有两个元素类型。ET命令是由ANSYS元素库中选择某个元素并定义该结构分析所使用的元素类型号码。
ITYPE:元素类型的号码
Ename:ANSYS元素库的名称,即使用者所选择的元素。
KOPT1~KOPT6:元素特性编码。
Menu Paths:Main Menu>Preprocessor Element Type>Add/Edit/Delete
MP,Lab,MAT,C0,C1,C2,C3,C4
定义材料的属性(Material Property),材料属性为固定值时,其值为C0,当随温度变化时,由后四个参数控制。
MAT:对应ET所定义的号码(ITYPE),表示该组属性属于ITYPE。
Lab:材料属性类别,任何元素具备何种属性在元素属性表中均有说明。例如杨氏系数(Lab=EX,EY,EZ),密度(Lab=DENS),泊松比(Lab=NUXY,NUXYZ,NUZX),剪切模数(Lab=GXY,GYZ,GXZ),热膨胀系数(Lab=ALPX,ALPY,ALPZ)等。
Menu paths:Main Menu>Preprocessor>Matial Props>Isotropic
R,NSET,R1,R2,R3,R4,R5,R6
定义”实常数”,即某一单元的补充几何特征,如梁单元的面积,壳单元的厚度。所带的的参数必须与元素表的顺序一致。
Menu paths:Main Menu>Preprocessor>Real Constants
E,I,J,K,L,M,N,O,P
SOLU
进入解题处理器,当有限元模型建立完以后,便可以进入/SOLU处理器,声明各种负载。但大部分负载的载声明也可在/PREP7中完成,建义全部负载在/SOLU处理中进行声明。
/ANTYPE,Antype,Status
声明分析类型,即欲进行哪种分析,系统默认为静力学分析。
Antype=STATIC or 0 静态分析(系统默认)
BUCKLE or 1屈曲分析
MODAL or 2 振动模态分析
HARMIC or 3 调和外力动和系统
TRANS or 4 瞬时动力系统分析
SFBEAM, ELEM, LKEY, Lab,VALI, VALJ, VAL2I, VAL2J, IOFFST, JOFFST
定义在梁元素上的分布力。
ELEM:元素号码。
LKEY:建立元素后,依节点顺序梁元素有四个面,该参为分力所施加的面号。
Lab RES(表示分布压力)。
VALI,VALJ:在I点及J点分布力的值。 前处理
/prep7
/pnum,label,key
!在有限元模块图形中显示号码。Label=欲显示对象的名称,node节点,elem元素,kp点,line线,area面积,volu体积;key=0为不显示号码(系统默认),=1为显示号码。
et,itype,ename,kopt1, kopt2, kopt3, kopt4, kopt5, kopt6,inopr
!元素类型定义。Itype为元素类型号码,通常由1开始;ename为ANSYS元素库的名称,如 beam3,plane42,solid45等;kopt1~kopt6为元素特性编码,如beam3的kopt6=1时,表示分析后的结果可输出节点的力及力矩,link1无需任何元素特性编码。
mp,lab,mat,c0,c1,c2,c3,c4
!定义材料特性。Lab为材料特性类别,如杨氏系数lab=ex、ey、ez,密度lab=dens,泊松比lab=nuxy、 nuyz、nuzx,剪力模数lab=gxy、gyz、gxz,热膨胀系数lab=alpx、alpy、alpz,热传导系数lab=kxx、kyy、 kzz,比热lab=c;mat对应前面定义的元素类型号码Itype;c0为材料特性类别的值。
r,nset,r1,r2,r3,r4,r5,r6
!元素几何特性。nset通常由1开始;r1~r6几何特性的值。
注:solid45元素不需要此命令,beam3单元有area截面积,惯性矩izz,高度height等。
例如:r,1,3e-4(截面积),2.5e-9(惯性矩),0.01(高度)
local,kcn,kcs,xc,yc,zc,thxy,thyz,thzx,par1,par2
!定义区域坐标系统。kcn区域坐标系统代号(大于10);kcs区域坐标系统属性(0为卡式坐标,1为圆柱坐标,2为球面坐标);xc,yc,zc(该区域坐标系统与整体坐标系统原点关系)。
csys,kcn !声明坐标系统,系统默认为卡式坐标(csys,0)。
k,npt,x,y,z !定义点。npt为点的号码;x,y,z为节点在目前坐标系统下的坐标位置。
kfill,np1,np2,nfill,nstrt,ninc,space
!点填充。np1和np2两点间,nfill为填充点的个数;nstrt,ninc,space为分布状态。
kgen,itime,np1,np2,ninc,dx,dy,dz,kinc,noelem,imove
!点复制。itime包含本身所复制的次数;knic为每次复制时点号码增加量;np1,np2,ninc点复制范围;dx,dy,dz每次复制在现有坐标下几何位置的改变量。
ksymm,ncomp,np1,np2,ninc, kinc,noelem,imove
!复制一组(np1,np2,ninc)点对称于某轴(ncomp);knic为每次复制时点号码增加量。
kl,nl1,ratio,nk1 !在已知线(nl1)上建立一个点(nk1),该点的位置由占全线段比例(radio)而定,比例为p1至nk1长度与p1至p2的长度。
kmodif,npt,x,y,z !修改现有点(npt)到新坐标(x,y,z)位置。
knode,npt,node !定义点(npt)于已知节点(node)上。
kdele,np1,np2,ninc !将一组点删除。
ksel,type,item,comp,vmin,vmax,vinc,kabs !选择有效点,type为选择方式。
Wpoffs,xoff,yoff,zoff !将工作平面中心点移到另外一点。
Wprota,thxy,thyz,thzx !将工作平面顺时针旋转一个角度。
l,p1,p2,ndiv,space,xv1,yv1,zv1,xv2,yv2,zv2 !由两点定义线段,此线段的形状可为直线(斜率)为0,或为曲线(以线段两端斜率xv1,yv1,zv1,xv2,yv2,zv2而定);ndiv为线段在进行网格化时欲分的元素数目。
Lstr,p1,p2 !用两个点来定义一条直线。
Lcomb,nl1,nl2,keep !将两条线合并为一条线,keep=0时原线段删除,keep=1时保留。
Ldiv,nl1,ratio,pdiv,ndiv,keep !将线分割为数条线,nl1为线段的号码;ndiv为线段欲分的段数(系统默认为两段),大于2时为均分;ratio为两段的比例(等于2时才作用);keep=0时原线段删除,keep=1时保留。
Lgen,itime,nl1,nl2,ninc,dx,dy,dz,kinc,noelem,imove !线段复制命令。itime包含本身所复制的次数;nl1,nl2,ninc为现有的坐标系统下复制到其他位置(dx,dy,dz);kinc为每次复制时线段号码的增加量。
Lfillt,nl1,nl2,rad,pcent !在两相交的线段nl1,nl2间产生一条半径等于rad的圆角曲线,同时自动产生三个点,其中两个点在nl1,nl2上,第三个点是新曲线的圆心定(若 pcent=0,则不产生该点)。
Larc,p1,p2,pc,rad !定义两点间的圆弧线,其半径为rad,pc为圆弧曲率中心部分的任何一点,不一定是圆心坐标。
Circle,pcent,rad,paxis,pzero,arc,nseg !产生圆弧线。Pcent为圆弧中心坐标点的号码;paxis 定义圆心轴正方向上任意点的号码;Pzero定义圆弧线起点轴上的任意点的号码,此点不一定在圆上;rad圆的半径;nseg为圆弧线欲划分的段数,完整为4。
Lang,nl1,p3,ang,phit,locat !产生一新的线段,此新的线段与已存在的线段nl1的夹角为ang,phit为新产生点的号码。
L2ang,nl1,nl2,angl,ang2,phit1,phit2 !产生新线段。此新线段与已存在的直线nl1夹角为ang1,与直线nl2的夹角为ang2。Phit1,Phit2为新产生两点的号码。
Ltan,nl1,P3,xv3,yv3,zv3 !产生三次曲线,该曲线方向为P2至P3,与已知曲线相切于P2。Xv3,y,v3,zv3为新线段在终点P3处的斜率。
L2tan,nl1,nl2 !建立新线段与已知两条相切的方式产生。若以负值输入,则相反。
Bspline,p1,p2,p3,p4,p5,p6,xv1,yv1,zv1,xv6,yv6,zv6 !通过6点曲线,并定义两端点的斜率。
spline,p1,p2,p3,p4,p5,p6,xv1,yv1,zv1,xv6,yv6,zv6 !通过6点曲线,每点之间形成一新线段,并可以定义两端点的斜率。
Ldele,nl1,nl2,ninc,kswp !kswp=0时只删除掉线段本身,=1时低单元点一并删除。
Lsel, type,item,comp,vmin,vmax,vinc,kabs !选择有效线段,type为选择方式。
A,p1,p2,p3,p4,p5,p6,p7,p8,p9 !由已知点定义面积
Al,l1,l2,l3,l4,l5,l6,l7,l8,l9,l10 !由已知线段定义面积
Agen, itime,na1,na2,ninc,dx,dy,dz,kinc,noelem,imove !面积复制命令。itime包含本身所复制的次数;na1,na2,ninc为现有的坐标系统下复制到其他位置(dx,dy,dz);kinc为每次复制时面积号码的增加量。
Arsym,ncomp,na1,na2,ninc,kinc,noelem,imove !复制一组面积na1,na2,ninc对称于轴ncomp;kinc为每次复制时面积号码的增加量。
Adrag, nl1,nl2,nl3,nl4,nl5,nl6, nlp1,nlp2,nlp3,nlp4,nlp5,nlp6 !面积的建立,沿某组线段路径,拉伸而成。
Arotat, nl1,nl2,nl3,nl4,nl5,nl6,pax1,pax2,arc,nseg !建立一组圆柱形面积,产生方式为绕着某轴(pax1,pax2,为轴上任意两点,并定义轴的方向)。Nseg为整个旋转角度方向中欲分段数目。
Aoffst,narea,dist,kinc !复制一块面积,产生方式为平移(offset)一块面积,以平面法线方向,平移距离为dist,kinc为面积号码增加量。
Afillt,na1,na2,rad !建立圆角面积,在两相交平面间产生曲面,rad为半径。
Askin,nl1,nl2,nl3,nl4,nl5,nl6 !沿已知线建立一个平滑薄层曲面。
Adele,na1,na2,ninc,kswp !kswp=0时只删除掉面积本身,=1时低单元点一并删除。
Asel, type,item,comp,vmin,vmax,vinc,kabs !选择有效面积,type为选择方式。
V,p1,p2,p3,p4,p5,p6,p7,p8 !由点定义体积。
Va,a1,a2,a3,a4,a5,a6,a7,a8,a9,a10 !由已知面定义体积
Vgen,itime,nv1,nv2,ninc,dx,dy,dz,kinc,noelem,imove !体积复制。
Vsymm,ncomp,nv1,nv2,ninc,kinc,noelem,imove!对称于轴(ncomp)复制一组体积
Vdrag,na1,na2,na3,na4,na5,na6,nlp1,…nlp6 !体积建立时将一组已知面积沿着某组线段路径,拉伸而成。
Vrotat,na1,na2,na3,na4,na5,na6,pax1,pax2,arc,nseg !建立一组圆柱形体积,产生方式为绕着某轴(pax1,pax2,为轴上任意两点,并定义轴的方向)。Nseg为整个旋转角度方向中欲分段数目。
Vdele,nv1,nv2,ninc,kswp !kswp=0时只删除掉体积本身,=1时低单元点一并删除。
Vsel, type,item,comp,vmin,vmax,vinc,kabs !选择有效体积,type为选择方式。
Rectng,x1,x2,y1,y2 !建立长方形面积
Pcirc,rad1,rad2,theta1,theta2 !建立平面圆面积,rad1,rad2为圆面积的内径及外径,theta1,theta2为圆面积的角度范围。系统默认值为0度到360度,每90度分段。
Rpoly,nsides,lside,majrad,minrad !建立一个以工作面中心点为基准的正多边形面积。边数为nsides,大小可由边长lside,或外接圆半径majard,或内切圆minrad。
Block,x1,x2,y1,y2,z1,z2 !建立一个长方体区块。
Blc4,xcorner,ycorner,width,height,depth !建立一个长方体区块。
Blc5,xcenter,ycenter,width,height,depth !建立一个长方体区块。区块体积中心点的x、y坐标。
Cylind,rad1,rad2,z1,z2,theta1,theta2 !建立一个圆柱体积,圆柱的方向为z方向并由z1,z2为z方向长度的范围;rad1,rad2为圆柱的内外半径;theat1,theta2为圆柱的起始、终结角度。
Cyl4,xcenter,ycenter, rad1, theta1, rad2,theta2,depth !建立一个圆柱体积。以圆柱体积中心点的x、y坐标为基准;rad1,rad2为圆柱的内外半径;theat1,theta2为圆柱的起始、终结角度。
Cyl5,xedge1,yedge1,xedge2,yedge2,depth !建立一个圆柱体积。xedge1,yedge1,xedge2,yedge2为圆柱上面或下面任一直径的x、y起点坐标与终点坐标。
Cone,rtop,rbot,z1,z2,theta1,theta2 !建立一个圆锥体积。Rtop,z1为圆锥上平面的半径与长度、rbot,z2为圆锥下平面的半径与长度;theat1,theta2为圆锥的起始、终结角度。
Rprism,z1,z2,nsides,lside,majrad,minrad !建立一个正多边形体积,z1,z2为z方向长度的范围,边数为nsides;边长lside;或外接圆半径majard;或内切圆minrad。
!声明元素大小、形状和网格种类
lesize,nl1,size,angsiz,ndiv,space,kforc,layer1,layer2 !定义所选择线段(nl1,nl1=all为目前所有的线段)进行元素网格化时元素的大小(size),元素的大小可用线段的长度(size)或该条线段要分割的元素数目(ndiv)来确定;space为间距比(最后一段长与最先一段长的比值,正值代表以线段方向为基准,负值以中央为基准,系统默认等间距)。
kesize,npt,size,fact1,fact2 !定义通过点(npt,npt=all为通过目前所有点的线段)的所有线段进行元素网格划分时元素的大小(size),不含lesize所定义的线段。元素的大小仅能用元素的长度(size)输入。该命令必须成对使用,因为线段基本上含两点。
esize,size,ndiv !定义元素网格化时元素的大小。该命令以目前所有对象为基准(不含lesize,kesize所定义的线段)。元素的大小可用元素的边长(size)或线段要分成元素数目(ndiv)来确定。
desize,minl,minh,mxel !系统默认元素大小(不含lesize,kesize,esize所定义)。
smrtsize,sizval,fac,expnd,trans,angl,angh,gratio,smhlc,smanc,mxitr,sprx !自由网格时,网格大小的高级控制(不含lesize,kesize,esize所定义)。一般由desize控制元素大小,desize及 smrtsize是相互独立的命令,仅能存在一个,执行smrtsize命令后desize自动无效。
mshkey,key !key=0自由网格(系统默认);key=1对应网格;key=2对应自由混合(仅适合2-D实体)。
mshape,key,dimension !声明网格化时元素的形状。2-D实体模型采用四边形(key=0)或全部为三角形(key=1)。
!进行网格化
xatt,mat,real,type,esys !type元素的形式号码,real元素的几何参数属性编号,mat元素的材料特性属性编号。Esys为建立元素时所在坐标系统号码。系统默认值为第一组及卡式坐标。
Xmesh !x对象网格化后,元素属性由xatt决定。
直接建模方法,求解器和后处理部分,自己整理部分:
三、 求解(solution)
/solu
antype,antype,status
!声明分析类型,系统默认为静力分析。antype=static or 0为静力分析(系统默认)。
f,node,lab,value,value2,nend,ninc !定义节点上的集中力。node为节点号码;lab为外力形式(结构力学中lab=fx,fy,fz,mx,my,mz力的方向、力矩方向;lab=heat热学中的热流量);value为外力的大小;node,nend,ninc为施力节点的范围。
d,node,lab,value,value2,nend,ninc,lab2,lab3,lab4,lab5,lab6
!定义节点自由度限制。node,nend,ninc自由度约束节点的范围;lab为自由度约束的形式(lab=ux、uy、uz、 rotx、roty、rotz;热学中lab=temp温度);结构力学中,lab=ux,lab2=uy。
sfbeam,elem,lkey,lab,vali,valj,val2i,val2j,ioffst,joffst
!定义分布力作用于梁。Elem元素号码;lkey定义分布力所施加面的号码(1,2,3,4);
lab=pres(表示分布压力);vali,valj为在I点及J点的分布力值。
sfe,elem,lkey,lab,kval,val1,val2,val3,val4
!定义分布力作用元素上。Elem元素号码;元素可分2D元素和3D元素,val1~val4的值为当初建元素时的节点顺序。 lab=pres(表示分布压力)。例如:sfe,4,2,pres,,20,60
sf,nlist,lab,value,value2
!定义节点间分布力。Nlist为分布力作用的边或面上的所有节点。通常用nsel命令选有效节点,然后设定 nlist=all;lab=pres结构力学的压力;value作用分布力的值。
outrp,item,freq,cname
!控制分析后的结果是否显示于输出窗口中。Item为欲选择结果的内容(item=all为所有结果,nsol为节点自由度结果,basic系统默认);freq为负载的次数,freq=all为最后负载。
outres,item,freq,cname
!控制分析结果存入数据库中的方式,通常使用其默认值。Item=all(系统默认)为选择结果的内容;freq=last(系统默认)为负载的次数。
solve!在解题过程中,质量矩阵、刚度矩阵、负载等资料都会保存在相关文件中。
lswrite,lsnum !将多重负载资料保存至文件中,所保存文件命名为jobname.sn,n=lsnum,n为2位数,第一负载n=01,第二负载n=02。
lssolve,slmin,lsmx,lsinc
!读取前所定义的多重负载,并求其解答。slmin,lsmx,lsinc为读取该阶段负载的范围。
ddele,node,lab,nend,ninc !将定义的约束条件删除。node,nend,ninc为欲删除约束条件节点的范围。Lab为欲删除约束条件的方向。
fdele,node,lab,nend,ninc !将已定义于节点上的集中力删除。node,nend,ninc为欲删除外力节点的范围。Lab为欲删除外力的方向。
sfdele,nlist,lab
!将定义的面负载删除。nlist为面负载所含节点。Lab=pres(结构力学)。
sfedele,elem,lkey,lab !将已定义的面负载从某元素上删除。Elem为元素号码;lkey为负载作用于元素边或面的号码;Lab=pres(结构力学)。 四、 后处理(postprocessing)
/post !一般后处理器,以便检查分析结果。
pldisp,kund !图标结构受外力的变形结果,kund=0为显示变形后的结构形状,kund=1为同时显示变形前和变形后的结构形状,kund=2为同时显示变形前和变形后的形状,但仅显示结构外观。
plesol,item,comp !图表元素的解答。以轮廓线方式表达,故会有不连续的状态,通常2-D及3-D元素才适用。Item为欲查看何种解答。
Item comp
S x,y,z,xy,yz,xz应力 S 1,2,3 主应力
S eqv,int 等效应力 F x,y,z 结构力
M x,y,z 结构力矩
plnsol,item,comp !图标节点的解答。以连续的轮廓线表示。
Item为欲查看何种解答。 Item comp
S x,y,z,xy,yz,xz应力 S 1,2,3 主应力
S eqv,int 等效应力 F x,y,z 结构力
M x,y,z 结构力矩 u x,y,z,sum 位移分量及向量位移
rot x,y,z,sum 旋转位移分量及向量旋转位移 temp 温度
prnsol,item,comp !打印节点的解答。Item为欲查看何种解答。
Item comp
U x,y,z 位移 U comp x,y,z方向及总向量方向的位移
S comp 应力 S prin 主应力,等效应力
etable,lab,item,comp
!将元素某项结果制作成表格形式。
pretab,lab1,lab2,lab3,lab4,lab5,lab6,lab7,lab8,lab9
!打印定义的表格资料。Lab1~lab9为前面所定义的表格字段名称。
pletab,itlab,avglab
!图标已定义的元素结果表格资料,图形的水平轴为元素号码,垂直轴为itlab值。Itlab为前面所定义的表格字段名称;avglab=noav不平均共同节点的值,avglab=avg平均共同节点的值。
plls,labi,labj,fact
!图标1-D线元素节点的结果。labi,labj为前面已定义I点及J点的结果。
set,lstep,sbstep,fact,kimg,time,ngle,nset
!当进行多重负载解题时,先行声明多重负载的号码lstep。例如,set,2表示欲检查第二个负载的结果。
save !保存目前所有的database资料。
resume, !回到最近save点重新开始。
/clear !清除目前所以的database资料,该命令在起始层才有效。
1 设置分析类型
ANTYPE,Antype,status,ldstep,action
其中antype表示分析类型
STATIC:静态分析
MODAL:模态分析
TRANS:瞬态分析
SPECTR:谱分析
2 KBC,KEY
制定载荷为阶跃载荷还是递增载荷
EKY=0 递增方式
KEY=1 阶跃方式
3 SOLVE 开始一个求解运算
4 LSSOLVE 读入并求解多个载荷步
5 TIME,time 设置求解时间
有时在分析中需要进入后处理,然后在保持进入后处理之前的状态的情况下接着算下去,可以使用以下的方法:
PARSAV,ALL,PAR,TXT
!PARSAV命令是储存ANSYS的参数,ALL代表所有参数,PAR是文件名,TXT是扩展名
/SOLU
ANTYPE,,REST,CruStep-1, ,CONTINUE
!ANTYPE是定义分析类型的命令,REST代表重启动,CruStep代表本载荷步的编号
PARRES,NEW,PAR,TXT
!PARRES是恢复参数的命令,NEW表示参数是以刷新状态恢复,PAR和TXT代表了储存了参数的文件名和扩展名
如果有单元生死的问题,可以这样处理:
ALLSEL,ALL
*GET,E_SUM_MAX,ELEM,,NUM,MAX !得到单元的最大编号,即单元的总数
ESEL,S,LIVE !选中“生”的单元
*GET,E_SUM_AL,ELEM,,COUNT
*DIM,E_POT_AL,,E_SUM_MAX !单元选择的指示
*DIM,E_NUM_AL,,E_SUM_AL !单元编号的数组
J=0
!读出所选单元号
*DO,I,1,E_SUM_MAX
*VGET,E_POT_AL(I),ELEM,I,ESEL
!对所有单元做循环,被选中的单元标志为“1”
*IF,E_POT_AL(I),EQ,1,THEN
J=J+1
E_NUM_AL(J)=I
*ENDIF
*ENDDO
ALLSEL,ALL
在重启动之后恢复单元生死状态
*if,E_SUM_AL,ne,0,then
*do,i,1,Num_Alive
esel,a,,,E_NUM_AL(i)
*enddo
ealive,all
allsel
*endif
/WINDOW, WN, XMIN, XMAX, YMIN, YMAX, NCOPY
注意x的坐标是 -1到1.67,y坐标是 -1到1
Xmin= off on, FULL, LEFT, RIGH, TOP, BOT, LTOP, LBOT, RTOP, RBOT
注意一个问题,除了1号窗口外,其他的不能用鼠标操作,只用先发/view和/dist,然后用/replot。NCOPY,指被拷贝的窗口
该命令可以比较两个窗口的不同点,从一个窗口拷贝到另外一个窗口,但是必须先试用命令/NOERASE ,然后再拷贝,使用 /ERASE,重新恢复
命令38:6.29
/DIST, WN, DVAL, KFACT
设定从观察人到焦点的距离
DVAL距离值
KFACT 0 代表用DVAl的实际值,1,代表DVAL为相对值,如0.5代表距离减少一半,也就是图像放大一倍
csys,kcn
kcn , 0 迪卡尔zuobiaosi
1 柱坐标
2 球
4 工作平面
5 柱坐标系(以Y轴为轴心)
n 已定义的局部坐标系
/grid, key
key: “0” 或“off” 无网络
“1”或“on” xy网络
“2”或“x” 只有x线
“3”或“y” 只有y线
xvar, n
n: “0”或“1” 将x轴作为时间轴
“n” 将x轴表示变量“n”
/axlab, axis, lab 定义轴线的标志
axis: “x”或“y”
lab: 标志,可长达30个字符
plvar, nvar, nvar2, ……,nvar10 画出要显示的变量(作为纵坐标)
Autots,key-----自动时间步长的设置
Key--自动时间步长的关键字,它有以下两种可选:
OFF -不打开自动时间步长
ON -- 打开自动时间步长.
再加各网格划分的吧:
MSHAPE,key,dimension---给单元进行网格划分
key-指定适合单元的类型,有以下几种:
0 -- 当是2D时,用四边形来给单元划分;当单元是3D时,用六边形给单元来划分。
1---当是2D时,用三角形来给单元划分;当单元是3D时,用四边形给单元来划分。
Dimension
指定所划分模型的维数。
2D -- 2维模型(面积划分)。
3D --3维模型 (体积划分)。
MSHKEY,key--给模型制定网格划分方法。
key-指定合适的类型,有以下几种:
0 -- free方式划分网格(缺省值)
1 -- 用mapped划分网格。
2 -- 如果可能的话用mapped划分网格,如果不行,就激活Smartsize采用free方式划分。
TYPE, itype — 指定元素类型指针。
itype--指定该单元的类型数。(缺省值为1)。
REAL, nset — 指定单元实常数指针。
nset--指定该值为后边定义单元的实常数值(缺省值为1)。
MAT, mat — 指定单元的材料属性指针。
mat--指定该值为后边定义单元的材料属性值。
1, 建立局部坐标系:
WPSTYL, SNAP, GRSPAC, GRMIN, GRMAX, WPTOL, WPCTYP, GRTYPE, WPVIS, SNAPANG: Controls the display
and style of the working plane.
snap:默认为0.05
grspac:默认为0.1
GRMIN, GRMAX:默认为-1,1
WPTOL:实体的精度值,默认为0.003
WPCTYP:坐标系类型,0,直角坐标系,1,柱面坐标系,2,球坐标系
GRTYPE:栅格类型,0,栅格和坐标都有,1仅有栅格,2 坐标(默认)
WPVIS:是否显示栅格,0,不显示GRTYPE(默认)1,显示GRTYPE
SNAPANG:角度的增量,只当wpcytp取1或2的时候使用,默认值是5度
2,CSYS,kcn:把先前定义的坐标系激活为下列坐标系:
0 笛卡儿整体坐标系
1,柱面坐标系,Z轴是旋转轴
2,球面坐标系
4或WP,工作坐标系
5,柱面坐标系,Y轴是旋转轴
11或更大的数,是你先前定义过的一些坐标系
3, D, node, lab, value, value2, nend, ninc, lab2, lab3, ……lab6 定义节点位移约束
Node : 预加位移约束的节点号,如果为all,则所有选中节点全加约束,要么就用
鼠标选取,或者用select选出要用的点再用all!
Lab: ux,uy,uz,rotx,roty,rotz,all,指的是某一个方向了!
Value,value2: 自由度的数值(缺省为0)
Nend, ninc: 节点范围为:node-nend,编号间隔为nincLab2-lab6: 将lab2-lab6以同样数值施加给所选节点。
注意:在节点坐标系中讨论
4,设置求解选项
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)求解
注:这个是转载的,觉得很有用,也写上!:)
等再总结再写上来,大家看看有用没?:)
/XFRM, LAB, X1, Y1, Z1, X2, Y2, Z2
定义旋转中心
LAB= NODE,KP ,LINE,AREA,VOLU,ELEM,XYZ,OFF
如果为实体,对应的X1 ,Y1 为实体的编号,如果为XYZ,对应的是两个点的坐标。可以只定义一个,然后该点即为旋转中心点。
控制进行动态缩转动时的中心点
*该命令用的不是很多,一般来讲焦点为默认的旋转中心,可以用该命令重新定义旋转中心。
/XRANGE, XMIN, XMAX
定义X轴显示的范围,一般要估计大小后确定。用/XRANGE,DEFAULT返回程序默认值,默认值为/GROPT中定义的值,程序自动标注对于对数标注通常显示的不准确。
XVAROPT, LAB
定义在X变量显示的参数,默认为SET NUMBER,
Y命令
/YRANGE, YMIN, YMAX, NUM
定义Y轴的范围,NUM为Y轴的数目
YMAX Y轴的最大值
YMIN Y轴的最小值
NUM Y轴的数目与命令/GRTYP设置有关,当/GRTYP,,2,数目为1-3,/GRTYP,2,数目为1-6
用/YRANGE,DEFAULT返回默认的程序自动选取标尺,整体的选项参照/GROPT命令
Z 命令
/ZOOM, WN, LAB, X1, Y1, X2, Y2
放大屏幕区域
WN 窗口号
LAB= OFF 重新返回最合适的状态
BACK 返回最后的状态
SCRN 屏幕 X1 Y1 为中心点 X2 Y2为角点
RECT 矩形 X1 Y1 ,X2 Y2对应的角点
GSUM
计算并显示实体模型的几何项目,(中心位置,惯性矩,长度面积,体积等),必须是被选择的点,线,面,体等,几何位置是整体坐标系中的位置,对于体和面,如果没有用AATT和VATT命令赋予材料号,则按单位密度来计算的,对于点和线,不管你使用了什么命令(LATT,KATT,MAT),都是按单位密度来计算。发出GSUM命令然后用*GET 和*VGET 命令来获得需要的数据,如果模型改变需要重新发出GSUM命令,该命令整合了KSUM,ASUM以及VSUM命令的功能。
/GRTYP, KAXIS
定义Y轴的数目
KAXIS= 1 单一轴,最多可以显示10条曲线
2 为每一条曲线定义一条Y轴,最多可以有三条曲线
3 同2,但是最多有6条曲线,而且是三维的可以采用等轴观看默认是 VIEW,1,1,2,3
1. /clear命令的应用:
Finish
/clear,start (nostart)
该命令删除所有的参数以及模型和结果,
要保留参数,可用一下命令,格式如下:
parsav,all,mypar !将参数保留起来
fini
/clear,start
parres,new,mypar !/clear后再释放保留的参数
注意:/clear命令并不删除*If 存储
/clear命令也不删除 *Do 存储 ,但将删除loop循环的所有参数
2.NUMCMP,KP
能将(相同实体或不同实体)上的kp按最初构建的顺序重新排序
如:k,1,
k,2,10,
k,4,20,
numcmp,kp
则线上的kp排序为1,2,3
3. NUMMRG,KP
能将统一类型实体上的kp连在一起,中间重复的kp只要一个,从新排序
k,1
k,2,10,
k,3,10,
k,4,20
nummrg,kp !此时kp的编号1,2,4 (k2,k3是同一个点)
numcmp,kp !此时kp的编号1,2,3
4.mshkey 定义网格的类型
mshkey,0 划分自由网格
mshkey,1 划分映射网格
mshkey,2 尽可能用映射,如若不能则用自由网格
5. ACLEAR,ALL
清除相关面上的节点和单元
6. 画线命令 L 和 LSTR
L :画线,与当前激活的坐标系的形式有关
LSTR:永远是 直角坐标系下的 直线形式,也就是不会随坐标系的变化而变化
7.Flst命令是GUI操作的拾取命令,总是与FITEM命令一起用,举例说明:
FLST,2,4,4,ORDE,2
!!第一个2表示拾取项作为后面命令的第一个条件,第一个4 表示拾取4项
!!第三个4 表示拾取直线号 最后一个2 表示有2项FITEM
FITEM,2,1
FITEM,2,-4 !负号表示与上面同类,即拾取1,2,3,4四条线
LCCAT,P51X !拾取的线作为LCCAT的第一个条件
8.circle 命令的使用
/PREP7
K,1,b/2,r
k,2,b/2,0
circle,1,r,,2,180,
!!第一个关键点指定圆心,第二个关键点确定零度点的位置,上述命令用来画一个-90~90的圆弧 1.*GET命令
*GET命令的使用格式为:
*GET, Par, Entity, ENTNUM, Item1, IT1NUM, Item2, IT2NUM
其中:
Par是存储提取项的参数名;
Entity是被提取项目的关键字,有效地关键字是NODE, ELEM, KP, LINE, AREA, VOLU, PDS等;
ENTNUM是实体的编号(若为0指全部实体);
Item1是指某个指定实体的项目名.例如,如果Entity是ELEM,那么Item1
要么是NUM(选择集中的最大或最小的单元编号),要么是COUNT (选择集中的单元数目).
可以把*GET命令看成是对一种树型结构从上至下的路径搜索,即从一般到特殊的确定.
2.删除参数,有两个方法!
(1)使用"="命令,右边为空,如aaa=表示删除参数aaa
(2)使用*set命令,不给参数赋值,如*set,aaa,来删除参数aaa
3.*REPEAT命令:最简单的循环命令,即按指定的循环次数执行上一条命令,而命令中的参数可以按固定的增量递增.
*REPEAT的用法为: NTOT, VINC1, VINC2, VINC3, VINC4, VINC5, VINC6, VINC7, VINC8, VINC9, VINC10, VINC11
NTOT表示当前命令被执行的次数(包括最初的一次)
VINC11~VINC11每执行一次第二个节点号加1.
注意:大多数以斜线(/)或星号(*)开头的命令,以及扩展名不是.mac的宏,都不可以重复调用.但是,以斜线(/)开头的图形命令可以重复调用.同时,要避免对交互式命令使用*REPEAT命令,诸如那些需要拾取或需要用户响应的命令!
4.*DOWHILE,parm 重复执行循环直到外部控制参数发生改变为止.
只要parm 为真,循环将不停的执行下去,如果parm为假,循环中止.
5.*cycle
当执行DO循环时,ANSYS程序如果需要绕过所有在*cycle和*ENDDO之间的命令,只需在下一次循环前执行它.
6.*ASK, Par, Query, DVAL:提示用户输入参数值
Par是数字字母名称,用于存储用户输入数据的标量参数的名称;
Query是文本串,向用户提示输入的信息,最多包含54个字符,不要使用具有特殊意义的字符,如"$"或"!";
DVAL是用户用空响应时赋给该参数的缺省值;该值可以是一个1-8个字符的字符串(括在单引号中),也可以是一个数值.如果没有赋缺省值,用户用空格响应时,该参数被删除.
7.*CFWRITE, Command:把ANSYS命令写到由*CFOPEN打开的文件中.
Command是将要写的命令或字符串.
8.*MSG, Lab, VAL1, VAL2, VAL3, VAL4, VAL5, VAL6, VAL7, VAL8
写输出信息通过ANSYS信号子程序.
该命令的VAL1到VAL8 参数均为字符参数.数据描述符%C用于在格式中指明字符数据(必须接在.*MSG命令后面).
9.PARRES, Lab, Fname, Ext, --从文件里面读参数,与PARSAV对应.
Lab:
NEW -- :用这些参数代替当前的参数
CHANGE -- :用这些参数扩展当前的参数,代替任意已经存在的
Fname:文件名和路径
Ext:扩展名
10.*VWRITE, Par1, Par2, Par3, Par4, Par5, Par6, Par7, Par8, Par9, Par10
通过该命令把数组中的数据写到格式化(表格式)的数据文件中.该命令最多可带有10个数组矢量作为参数,并把这些矢量中包含的数据写入当前打开的文件(*CFOPEN命令)中.
1.LMESH,NL1,NL2,NINC
对线划分网格的命令
参数说明:
NL1,NL2:划分网格的线的起止号
NINC: 线号的增量
【例如】Lmesh,1,3,1 !对线1,2,3划分网格
与定义截面有关的命令
2.SECWRITE, Fname, Ext, --, ELEM_TYPE
创建用户自定义截面,截面信息以ASCII形式存放
参数说明:
Fname:定义的截面名称
XT: 截面文件的扩展名,默认为 .sect
--: 空着不填
ELEM_TYPE:单元类型
【例如】SECWRITE,aa,,,,
用户自定义的截面,必须通过Plane82或Mesh200单元创建
3.SECREAD, Fname, Ext, --, Option
将用户自定义的截面读入Ansys中
参数说明:
Fname: 定义的截面名称,以及文件存放的路径
EXT: 截面文件的扩展名,默认为 .sect
--: 空着不填
Option:截面文件的来源
LIBRARY:来自截面库中,
MESH: 用户创建的截面文件
【例如】SECREAD,aa,,,MESH
4.SECOFFSET,Location,OFFSET1,OFFSET2,CG-Y,CG-Z,SH-Y,SH-Z
这个命令用来定义粱的节点与截面的位置位置关系
location:梁桥中节点的位置
ORIGIN:粱的节点置于截面的坐标原点
CENT: 粱的节点置于截面的形心
SHRC: 粱的节点置于截面的剪切中心
USER: 粱的节点与截面的位置关系由用户通过OFFSET1,OFFSET2指定
OFFSET1,OFFSET2只有在location为USER时起作用,
其值分别为相对截面的坐标原点的Y,Z轴的偏移量
1.*ABBR,Abbr,String--定义一个缩略语.
Abbr:用来表示字符串"String"的缩略语,长度不超过8个字符.
String:将由"Abbr"表示的字符串,长度不超过60个字符.
2.ABBRES,Lab,Fname,Ext-从一个编码文件中读出缩略语.
Lab:指定读操作的标题,
NEW:用这些读出的缩略语重新取代当前的缩略语(默认)
CHANGE:将读出的缩略语添加到当前缩略语阵列,并替代现存同名的缩略语.
Ext:如果"Fname"是空的,则缺省的扩展命是"ABBR".
3.ABBSAV,Lab,Fname,Ext-将当前的缩略语写入一个文本文件里
Lab:指定写操作的标题,若为ALL,表示将所有的缩略语都写入文件(默认)
4./UCMD,Cmd,SRNUM-给一个用户定义的命令名赋值.
Cmd:用户定义的命令名,只有前面的4个字符有意义.
SRNUM:对该命令来说,是编制好的用户子程序编号(1~10).
5.*AFUN,Lab-在参数表达式中,为角度函数指定单位.
Lab:指定将要使用的角度单位.有3个选项.
RAD:在角度函数的输入与输出中使用弧度单位(默认)
DEG:在角度函数的输入与输出中使用度单位.
STAT:显示该命令当前的设置(即是度还是弧度).
6.*DEL,Val1,Val2-删除一个或多个参数
Val1:有2个选项
ALL:删除所有用户定义的参数,或者是所有用户定义和系统定义的参数.
空:仅删除变量"Val2"指定的参数.
Val2:有下列选项!
Loc:若Val1=空,变量Val2可以指定参数在数组参数对话框中的位置他是按字母排列的结果:若VAl1=ALL时,这个选项无效
_PRM:若Val1=ALL时,表明要删除所有包含以下划线开头的参数(除了"_STATUS"和"_RETURN"),若 Val1为空, 表明仅删除以下划线开头的参数.
PRM_:若Val1=空,仅删除以下划线结尾的参数;若Val1=ALL,该选项无效.
空:若Val1=ALL,所有用户定义的参数都要删除.
7./INQUIRE,StrArray,FUNC-返回系统信息给一个参数.
StrArray:将接受返回值的字符数组参数名.
FUNC:指定系统信息返回的类型.
8.*CFCLOS — 关闭一个"命令"文件.
格式:*CFCLOS
9.*CFOPEN, Fname, Ext-打开一个"命令"文件
Ext:如果Fname为空,则其扩展名为"CMD"
10.*CREATE, Fname, Ext-打开或生成一个宏文件
Fname:若在宏里,使用命令"*USE"的Name选项读入文件时,不要使用路径名.
Ext:若在宏里,使用命令"*USE"的Name选项读入文件时,不要使用文件文件扩展名
1./PMACRO 指定宏的内容被写入到ANSYS的会话LOG文件中.
2./PSEARCH,Pname-为用户自定义的宏文件指定一个搜索目录.
Pname:将要搜索的中间目录路径名,长度不超过64个字符,最后必须是一个分界符.缺省时就是用户的根目录.
3./TEE,Lable,Fname,Ext-在命令被执行的同时,写一些列的命令到一个指定的文件中.
Lable:指导ANSYS软件对命令"/TEE"的处理方式.有下面选项:
NEW:将命令行的文本写入到文件Fname中,如果该文件Fname已经 存在,则将覆盖其内容.
APPEND:将命令行的文本添加到文件Fname中.
END:结束命领行文本写入或添加.
Ext:如果希望像执行ANSYS命令一样执行这个文件,则其扩展命为 ".mac" .
4.*ULIB,Fname,Ext-确定一个宏库文件.
5.*USE, Name, ARG1, ARG2, ARG3, ARG4, ARG5, ARG6, ARG7, ARG8, ARG9, AR10, AR11, AR12, AR13, AR14, AG15, AR16, AR17, AR18 -执行一个宏文件.
Name:用字母开头且长度不超过32个字符的名称,它可以是一个宏文件名,或者是一个宏库文件里的宏块名.
ARG1~ AR18 :将值传递给宏文件或宏块中ARG1~ AR18参数被引用的地方.
6./WAIT,DTIME-在读下一个命令时引起的一个延时.
DTIME:延时时间,单位为秒,最大的延时时间为59秒.
7.*GO,Base-在输入文件里,程序执行指定行.
Base:将要"进行"的动作.选项有:
:lable是一个用户定义的标题,必须以":"开头,后面的字符最多 不超过8个.命令读入器会跳到与":lable"相匹配的那行.
STOP:它会引起ANSYS程序从当前位置退出.
1.*MFOURI, Oper, COEFF, MODE, ISYM, THETA, CURVE -计算一个傅立叶的系数或者求出其值.
Oper:傅立叶运算的类型.有下面的选项:
FIT:根据 MODE, ISYM, THETA, CURVE 求出傅立叶的系数COEFF.
EVAL:根据COEFF, MODE, ISYM, THETA计算傅立叶曲线的CURVE
COEFF:包含傅立叶系数的数组参数名.
MODE:包含着预期傅立叶项模态数的数组参数名.
ISYM:包含着相应傅立叶级数项对称字的数组参数名.
THETA, CURVE :分别包含着θ和 CURVE 描述的数组参数名.
2.*MFUN, ParR, Func, Par1 -对一个数组参数矩阵进行复制或转置.
ParR:结果数组参数名,这个参数必须是一个具有维数大小的数组.
Func:复制或转置函数.若Func=COPY,Par1 被复制到ParR里,若
Func=TRAN,Par1 被转置到ParR里,其中矩阵Par1 中的行号(m)和列号(n)被转置为矩阵中的列号和行号.
Par1: 输入将要复制或转置的数组参数矩阵!
3.*TOPER, ParR, Par1, Oper, Par2, FACT1, FACT2, CON1 — — 对表格参数进行操作.
ParR:结果表格参数.
Par1:第一个表格参数的名称.
Oper:将要完成的操作.如ADD表示:
ParR(i,j,k)=FACT1*Par1(i,j,k)+FACT2*Par2(i,j,k)+CON1
Par2:第2个表格参数的名称.
FACT1:与第1个表格参数相乘的因子,缺省为1.0;
FACT2:与第2个表格参数相乘的因子,缺省为1.0;
CON1 :偏移的常数增量,缺省为0.
4.*VABS, KABSR, KABS1, KABS2, KABS3 —给函数或数组参数施加绝对值.
KABSR:结果参数的绝对值.若为0,不取绝对值,若为1,取绝对值.
KABS1, KABS2, KABS3 :分别对1、2、3个参数取绝对值的控制键,若为0,不取绝对值,若为1,取绝对值.
绝对值施加到操作进行之前的每个输入参数上和操作完成之后的结果上。
5.*VCOL, NCOL1, NCOL2 — 在矩阵运算中指定列标号。
NCOL1, NCOL2 :与命令“*MXX”运算中,分别对Par1、Par2所使用的列标号。默认值就是填充数组结果的值。
注意:在数组参数矩阵运算中,指定列标号。子矩阵的大小将由从运算命令中定义的左上角数组元素的开始处到右下角的元素来确定,右下角元素的列标号将由本命令来指定,右下角元素的行标号将由“*VLEN”命令来指定. {:{39}:} 好东西,支持一下 好东西,谢谢楼主 好东西,谢谢楼主!收藏了,慢慢学习 比较全面的汇总资料 谢谢楼主分享。{:{23}:} 很全面很喜欢啊