第二章FLAC3D原理及入门指南
2-1定义一个FISH函数
new
def abc
abc = 25 * 3 + 5
End
print abc
2-2使用一个变量
new
def abc
hh = 25
abc = hh * 3 + 5
End
Print hh
Print abc
2-3对变量和函数的理解
new
def abc
hh = 25
abc = hh * 3 + 5
End
set abc=0 hh=0
print hh
print abc
print hh
new
def abc
abc = hh * 3 + 5
end
set hh=25
print abc
set abc=0 hh=0
print hh
print abc
print hh
2-4获取变量的历史记录
new
gen zone brick size 1 2 1
model mohr
prop shear=1e8 bulk=2e8 cohes=1e5 tens=1e10
fix x y z range y -0.1 0.1
apply yvel -1e-5 range y 1.9 2.1
plot set rotation 0 0 45
plot block group
def get_ad
ad1 = gp_near(0,2,0)
ad2 = gp_near(1,2,0)
ad3 = gp_near(0,2,1)
ad4 = gp_near(1,2,1)
end
get_ad
def load
load=gp_yfunbal(ad1)+gp_yfunbal(ad2)+gp_yfunbal(ad3)+gp_yfunbal(ad4) end
hist load
hist gp ydis 0,2,0
step 1000
plot his 1 vs -2
2-5用FISH函数计算体积模量和剪砌模量
new
def derive
s_mod = y_mod / (2.0 * (1.0 + p_ratio))
b_mod = y_mod / (3.0 * (1.0 - 2.0 * p_ratio))
end
set y_mod = 5e8 p_ratio = 0.25
derive
print b_mod
print s_mod
2-6 在FLAC输入中使用符号变量
New
def derive
s_mod = y_mod / (2.0 * (1.0 + p_ratio))
b_mod = y_mod / (3.0 * (1.0 - 2.0 * p_ratio))
end
set y_mod = 5e8 p_ratio = 0.25
derive
gen zone brick size 2,2,2
model elastic
prop bulk=b_mod shear=s_mod
print zone prop bulk
print zone prop shear
2-7 控制循环
New
def xxx
sum = 0
prod = 1
loop n (1,10)
sum = sum + n
prod = prod * n
end_loop
end
xxx
print sum, prod
new
gen zone brick p0 (0,0,0) p1 (-10,0,0) p2 (0,10,0) p3 (0,0,-10) model elas
plot set rotation 0 0 45
plot block group
def install
pnt = zone_head
loop while pnt #null
z_depth = -z_zcen(pnt)
y_mod = y_zero + cc * sqrt(z_depth)
z_prop(pnt, ’shear’) = y_mod / (2.0*(1.0+p_ratio))
z_prop(pnt, ’bulk’) = y_mod / (3.0*(1.0-2.0*p_ratio))
pnt = zone_next(pnt)
end_loop
end
set p_ratio=0.25 y_zero=1e7 cc=1e8
install
2-8 拆分命令行
new ;example of a sum of many things
def long_sum
temp = v1 + v2 + v3 + v4 + v5 + v6 + v7 + v8 + v9 + v10
long_sum = temp + v11 + v12 + v13 + v14 + v15
end
2-9 变量类型
new
def haveone
aa = 2
bb = 3.4
cc = ’Have a nice day’
dd = aa * bb
ee = cc + ’, old chap’
end
haveone
print fish
2-10 IF条件语句
new
def abc
if xx > 0 then
abc = 1000
else
abc = -1000
end_if
end
set xx = 10
print abc
set xx = 0
print abc
2-11 索单元自动生成
new
gen zone brick size 10 3 5
plot set rotation 0 0 45
plot block group
def place_cables
loop n (1,5)
z_d = float(n) - 0.5
command
sel cable beg 0.0,1.5,z_d end 7.0,1.5,z_d nseg 7 end_command
end_loop
end
place_cables
plot grid sel geom red
new
gen zone brick size 10 3 5
plot set rotation 15 0 60
plot block group
mod mohr
prop bulk 1e8 shear .3e8 fric 35
prop coh 1e3 tens 1e3
ini dens 1000
set grav 0,0,-10
fix x y z range z -.1 .1
fix y range y -.1 .1
fix y range y 2.9 3.1
fix x range x -.1 .1
fix x range x 9.9 10.1
set large
hist unbal
solve
save cab_str.sav
ini xdis 0 ydis 0 zdis 0
hist gp xdisp 0,1,5
def place_cables
loop n (1,5)
z_d = 5.5 - float(n)
z_t = z_d + 0.5
z_b = z_d - 0.5
command
free x range x -.1,.1 z z_b z_t
solve
sel cable beg 0.0,0.5,z_d end 7.0,0.5,z_d nseg 7
sel cable beg 0.0,1.5,z_d end 7.0,1.5,z_d nseg 7
sel cable beg 0.0,2.5,z_d end 7.0,2.5,z_d nseg 7
sel cable prop emod 2e10 ytension 1e8 xcarea 1.0 & gr_k 2e10 gr_coh 1e10 gr_per 1.0
end_command
end_loop
end
place_cables
save cab_end.sav
plot sketch sel cable force red
2-12圆形隧道开挖模拟计算
;建立模型
gen zon radcyl p0 0 0 0 p1 6 0 0 p2 0 1 0 p3 0 0 6 &
size 4 2 8 4 dim 3 3 3 3 rat 1 1 1 1.2 group outsiderock
gen zone cshell p0 0 0 0 p1 3 0 0 p2 0 1 0 p3 0 0 3 &
size 1 2 8 4 dim 2.7 2.7 2.7 2.7 rat 1 1 1 1 group concretliner fill group insiderock
gen zon reflect dip 90 dd 90 orig 0 0 0
gen zon reflect dip 0 dd 0 ori 0 0 0
gen zon brick p0 0 0 6 p1 6 0 6 p2 0 1 6 p3 0 0 13 size 4 2 6 group outsiderock1
gen zon brick p0 0 0 -12 p1 6 0 -12 p2 0 1 -12 p3 0 0 -6 size 4 2 5 group outsiderock2
gen zon brick p0 6 0 0 p1 21 0 0 p2 6 1 0 p3 6 0 6 size 10 2 4 group outsiderock3
gen zon reflect dip 0 dd 0 orig 0 0 0 range group outsiderock3
gen zon brick p0 6 0 6 p1 21 0 6 p2 6 1 6 p3 6 0 13 size 10 2 6 group outsiderock4 gen zon brick p0 6 0 -12 p1 21 0 -12 p2 6 1 -12 p3 6 0 -6 size 10 2 5 group outsiderock5
gen zon reflect dip 90 dd 90 orig 0 0 0 range x -0.1 6.1 z 6.1 13.1
gen zon reflect dip 90 dd 90 orig 0 0 0 range x -0.1 6.1 z -6.1 -12.1
gen zon reflect dip 90 dd 90 orig 0 0 0 range x 6.1 21.1 z -12.1 13.1
;绘制模型图
plot block group
plot add axes red
;plot set rotation 0 0 45 用于显示三维模型
;设置重力
set gravity 0 0 -10
;给定边界条件
fix z range z -12.01,-11.99
fix x range x -21.01,-20.99
fix x range x 20.99,21.01
fix y range y -0.01 0.01
fix y range y 0.99,1.01
;求解自重应力场
model mohr
ini density 1800 ;围岩的密度
prop bulk=1.47e8 shear=5.6e7 fric=20 coh=5.0e4 tension=1.0e4 ;体积、剪切、摩擦角、凝聚力、抗拉强度
set mech ratio=1e-4
solve
save Gravsol.sav
plot cont zdisp outl on
plot cont szz
;毛洞开挖计算
initial xdisp=0 ydisp=0 zdisp=0
model null range group insiderock any group concretliner any
plot block group
plot add axes red
set mech ratio=5e-4
solve
save Kaiwsol.sav
plot cont zdisp
plot cont sdisp
plot cont szz
plot cont xzz
;模筑衬砌计算
model elas range group concretliner any
plot block group
plot add axes red
ini density 2500 range group concretliner any ;衬砌混凝土的密度
prop bulk=16.67e9,shear=12.5e9 range group concretliner any ;衬砌混凝土的体积弹模、剪切弹模
set mech ratio=1e-4
solve
save zhihusol.sav
plot cont zdisp
plot cont sdisp
plot cont szz
plot cont xzz
;完成计算分析
第四章 FISH语言及建模技术
4-1数组
new
def afill ;fill matrix with random numbers
array var(4,3)
loop m (1,4)
loop n (1,3)
var(m,n) = urand
endloop
endloop
end
def ashow ;display contents of matrix
loop m (1,4)
hed = ’ ’
msg = ’ ’+string(m)
loop n (1,3)
hed = hed + ’ ’+string(n)
msg = msg + ’ ’+string(var(m,n))
endloop
if m = 1
dum = out(hed)
endif
dum = out(msg)
endloop
end
afill
Ashow
4-2函数操作
new
def xxx
aa = 2 * 3
xxx = aa + bb
end
new
def stress_sum
stress_sum = 0.0
pnt = zone head
loop while pnt # null
Stress_sum = stress_sum + z_sxx(pnt)
pnt=z_next(pnt)
end_loop
end
new
def stress_sum
sum = 0.0
pnt = zone head
loop while pnt # null
sum = sum + z_sxx(pnt)
pnt=z next(pnt)
end loop
stress_sum = sum
end
new
def h_var_1
ipz = z_near(1,2,3)
H_var_1 = sxx(ipz) + fac
H_var_2 = syy(ipz) + fac
H_var_3 = szz(ipz) + fac
H_var_4 = sxy(ipz) + fac
H_var_5 = syz(ipz) + fac
H_var_6 = sxz(ipz) + fac
end
4-3函数删除与重定义
new
def joe
ii=out(‘This is A Function’)
end
def fred
joe
end
fred
def joe
ii=(‘This is A New Function’)
end
fred
4-4字符串
new
def in_def
xx = in(msg+’(’+’default:’+string(default)+’):’) if type(xx) = 3
in_def = default
else
in_def = xx
endif
end
;
def moduli_data
default = 1.0e9
msg=’Input Young‘s modulus ’
Y_mod = in_def
default = 0.25
msg=’Input Poisson‘s ratio ’
p_ratio = in_def
if p_ratio = 0.5 then
ii = out(’ Bulk mod is undefined at Poisson‘s ratio = 0.5’)
ii = out(’ Select a different value --’)
p_ratio = in_def
endif
s_mod = y_mod / (2.0 * (1.0 + p_ratio))
b_mod = y_mod / (3.0 * (1.0 - 2.0 * p_ratio))
end
;
moduli_data
gen zone brick size 2,2,2
model elastic
prop bulk = b_mod shear = s_mod
print p_ratio y_mod b_mod s_mod
pause
print zone prop bulk
pause
print zone prop shear
4-5马蹄形隧道网格
new
gen zone radcyl size 5 10 6 12 rat 1 1 1 1.2 &
p0 0,0,0 p1 100,0,0 p2 0,200,0 p3 0,0,100
gen zone radtun size 5 10 5 12 rat 1 1 1 1.2 &
p0 0,0,0 p1 0,0,-100 p2 0,200,0 p3 100,0,0
gen zone reflect dip 90 dd 270 origin 0,0,0
plot set rotation 30 0 30
plot block group
4-6复杂形状网格生成
new
; 主隧道网格
gen zon radcyl p0 15 0 0 p1 23 0 0 p2 15 50 0 p3 15 0 8 &
size 4 10 6 4 dim 4 4 4 4 rat 1 1 1 1 fill
gen zon reflect dip 90 dd 90 orig 15 0 0
gen zon reflect dip 0 ori 0 0 0
; 辅助隧道网格
gen point id 1 (2.969848,0.0,-0.575736)
gen point id 2 (2.969848,50.0,-0.575736)
gen zon radcyl p0 0 0 -1 p1 7 0 0 p2 0 50 -1 p3 0 0 8 p4 7 50 0 & p5 0 50 8 p6 7 0 8 p7 7 50 8 p8 point 1 p10 point 2 &
八天学会Ansys命令流 八天学会Ansys命令流
为方便大家的交流和学习,特推出"跟我学命令流"课程本课程分为三部分:前处理,加载求解,后处理 每部分的学习时间:10天,共计30天 每天学习大约10个命令 希望本课程对大家能有所帮助 第一天 目标:熟悉ANSYS基本关键字的含义 k --> Keypoints 关键点 l --> Lines 线 a --> Area 面 v --> Volumes 体 e --> Elements 单元 n --> Nodes 节点 cm --> component 组元 et --> element type 单元类型 mp --> material property 材料属性 r --> real constant 实常数 d --> DOF constraint 约束 f --> Force Load 集中力 sf --> Surface load on nodes 表面载荷
bf --> Body Force on Nodes 体载荷 ic --> Initial Conditions 初始条件 第二天 目标:了解命令流的整体结构,掌握每个模块的标识 !文件说明段 /BATCH /TITILE,test analysis !定义工作标题 /FILENAME,test !定义工作文件名 /PREP7 !进入前处理模块标识 !定义单元,材料属性,实常数段 ET,1,SHELL63 !指定单元类型 ET,2,SOLID45 !指定体单元 MP,EX,1,2E8 !指定弹性模量 MP,PRXY,1,0.3 !输入泊松比 MP,DENS,1,7.8E3 !输入材料密度 R,1,0.001 !指定壳单元实常数-厚度...... !建立模型 K,1,0,0,, !定义关键点 K,2,50,0,,
A N S Y S命令流解释大 全 Document serial number【LGGKGB-LGG98YT-LGGT8CB-LGUT-
一、定义材料号及特性 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=，C=27，ψ=45的命令如下：
MP，EX，1，1E8 MP，NUXY，1， TB，DP，1 TBDATA，1，27 TBDATA，2，45这里要注意的是，在前处理的最初，要将角度单位转化到“度”，即命令：*afun,deg 三、单元生死载荷步 !第一个载荷步 TIME,... !设定时间值（静力分析选项） NLGEOM,ON !打开大位移效果 NROPT,FULL !设定牛顿－拉夫森选项 ESTIF,... !设定非缺省缩减因子（可选） ESEL,... !选择在本载荷步中将不激活的单元 EKILL,... !不激活选择的单元 ESEL,S,LIVE !选择所有活动单元 NSLE,S !选择所有活动结点 NSEL,INVE !选择所有非活动结点（不与活动单 元相连的结点） D,ALL,ALL,0 !约束所有不活动的结点自由度（可 选） NSEL,ALL !选择所有结点 ESEL,ALL !选择所有单元
在ANSYS中，命令流是由一条条ANSYS的命令组成的一个命令组合，这些命令按照一定顺序排布，能够完成一定的ANSYS功能，这些功能一般来说通过菜单操作也能够实现（而那些命令流能够实现，菜单操作实现不了的单个命令比较少见）。以下命令是结合我自身经验，和前辈们的一些经验而总结出来的，希望对大家有帮助。 （1）．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 选择线及相关关键点、节点和单元 （2）．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”仅用绝对值 （3）．Esel, type, item, comp, vmin, vmax, vinc, kabs选择一组单元 type: S: 选择一组单元（缺省） R: 在当前组中再选一部分作为一组 A: 为当前组附加单元 U: 在当前组中不选一部分单元 All: 选所有单元 None: 全不选 Inve: 反向选择当前组 Stat: 显示当前选择状态 Item：Elem: 单元号 Type: 单元类型号
1、A,P1,P2,P3,P4,P5,P6,P7,P8,P9 此命令用已知的一组关键点点（P1~P9 ）来定义面（Area）， 最少使用三个点才能围成面，同时产生转围绕些面的线。 点要依次序输入，输入的顺序会决定面的法线方向。 如果超过四个点，则这些点必须在同一个平面上。Men uPaths:Ma inMenu >Preprocessor>Cre ate>Arbitrary>ThroughKPs 2、G ABBR，Abbr,String ――定义一个缩略语. Abbr:用来表示字符串"String "的缩略语，长度不超过8个字符. String :将由"Abbr "表示的字符串，长度不超过6 0个字符. 3、A BBRES，Lab，Fname，EGt —从一个编码 文件中读出缩略语. Lab :指定读操作的标题， NEW :用这些读出的缩略语重新取代当前的缩略语（默认） CHANGE :将读出的缩略语添加到当前缩略语阵列，并替代现存同名的缩略语. EGt:如果"Fname "是空的，则缺省的扩展命是"ABBR". 4、ABBSAV , Lab , Fname , EGt —将当前的 缩略语写入一个文本文件里 Lab :指定写操作的标题，若为ALL，表示将所有 的缩略语都写入文件（默认） 5、 add,ir,ia,ib,ic,name,--,--,facta,factb,f actc 将ia,ib,ic变量相加赋给ir变量 ir,ia,ib,ic :变量号 name:变量的名称 6、 Adele,na1,na2,ninc,kswp ! kswp=O 时只 删除掉面积本身，=1时低单元点一并删除。 7、 Adrag ,n 11, nl2, nl3, nl4, nl5, nl6, nlp1, nlp2, nlp3 ,n lp4, nlp5, nlp6 !面积的建立，沿某组线段路径，拉伸而成。 8、Afillt,na1,na2,rad !建立圆角面积，在两相 交平面间产生曲面，rad为半径。 9、GAFUN，Lab 在参数表达式中，为角度函数指定单位.
一、定义材料号及特性 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 三、单元生死载荷步 !第一个载荷步 TIME,... !设定时间值（静力分析选项） NLGEOM,ON !打开大位移效果 NROPT,FULL !设定牛顿－拉夫森选项 ESTIF,... !设定非缺省缩减因子（可选） ESEL,... !选择在本载荷步中将不激活的单元 EKILL,... !不激活选择的单元 ESEL,S,LIVE !选择所有活动单元 NSLE,S !选择所有活动结点 NSEL,INVE !选择所有非活动结点（不与活动单 元相连的结点） D,ALL,ALL,0 !约束所有不活动的结点自由度（可 选） NSEL,ALL !选择所有结点
!ANSYS命令流学习笔记5 workbench中命令流的一些应用 学习重点： 1. 定义单元类型 2. 使用各向异性材料时，定义其单元为圆柱坐标系 3. 有角度吊装时，定义吊装约束 4. workbench、APDL的联合仿真 案例如下： 如下图模型，四个顶点通过杆件连接，进行吊装时的有限元分析。 1. 建立模型，设定必要的坐标系。 分成两个solid，内部的圆柱，剩余的矩形部分。建立两个坐标系，分别用于指定各向异性材料的属性、吊装的固定点。 下图，建立圆柱坐标系，编号100，用于指定各向异性材料。 下图，建立直角坐标系，编号12，用于指定吊装固定点。 2. 建立named selection，方便在命令流中选择必要的元素。 下图，将四个吊装点，中间的圆柱，分别定义为任何名称，必须是英文才能用于APDL命令中。 3. 定义边界条件，施加重力加速度，在static structural 下插入command(APDL) ，内容如下 /prep7 alls cmsel,s,c1,elem !选择c1单元所有节点,既圆柱体的所有单元 emodif,all,esys,100 !其坐标系转换为100坐标系，因为缠绕的各向异性材料必须在圆柱坐标系下定义单元的坐标。 !完成对各向异性材料的坐标系设定。 et,10,10 !定义编号为10的，link10单元 r,10,0.01 !定义编号为10的实常数0.01，用于定义link10单元的截面积0.01mm^2 *get,nmax,node,,num,max !获取node的最大数值，储存在nmax名称的变量里 csys,12 !调用csys12坐标系 n,nmax+1, !csys坐标原点建立node，后面会将其固定，既吊装的固定点 mat,1 type,10 !选取编号10的单元类型 real,10 !选取编号10的实常数 cmsel,s,kk1,node !选择kk1点，kk1已经定义为named selection *get,k1,node,,num,max !获取已选节点中的节点数最大值，既kk1的节点编号，取值为k1 cmsel,s,kk2,node ! k2,k3,k4方法类型k1 *get,k2,node,,num,max cmsel,s,kk3,node *get,k3,node,,num,max cmsel,s,kk4,node *get,k4,node,,num,max
ansys——ANSYS命令流（Ⅰ） 1. A，P1，P2，…，P17，P18（以点定义面） 2. AADD，NA1，NA2，…NA8，NA9（面相加） 3. AATT，MAT，REAL，TYPE，ESYS，SECN（指定面的单元属性） 【注】ESYS为坐标系统号、SECN为截面类型号。 4. *ABBR，Abbr，String（定义一个缩略词） 5. ABBRES，Lab，Fname，Ext（从文件中读取缩略词） 6. ABBSAVE，Lab，Fname，Ext（将当前定义的缩略词写入文件） 7. ABS，IR，IA，--，--，Name，--，--，FACTA（取绝对值） 【注】************* 8. ACCAT，NA1，NA2（连接面） 9. ACEL，ACEX，ACEY，ACEZ（定义结构的线性加速度） 10. ACLEAR，NA1，NA2，NINC（清除面单元网格） 11. ADAMS，NMODES，KSTRESS，KSHELL 【注】************* 12. ADAPT, NSOLN, STARGT, TTARGT, FACMN, FACMX, KYKPS, KYMAC 【注】************* 13. ADD，IR, IA, IB, IC, Name, --,-- , FACTA, FACTB, FACTC（变量加运算） 14. ADELE，NA1，NA2，NINC，KSWP（删除面） 【注】KSWP =0删除面但保留面上关键点、1删除面及面上关键点。 15. ADRAG，NL1，NL2，…，NL6，NLP1，NLP2，…，NLP6（将既有线沿一定路径拖拉成面） 16. AESIZE，ANUM，SIZE（指定面上划分单元大小） 17. AFILLT，NA1，NA1，RAD（两面之间生成倒角面） 18. AFSURF，SAREA，TLINE（在既有面单元上生成重叠的表面单元） 19. *AFUN, Lab（指定参数表达式中角度单位） 20. AGEN, ITIME, NA1, NA2, NINC, DX, DY, DZ, KINC, NOELEM, IMOVE（复制面） 21. AGLUE，NA1，NA2，…，NA8，NA9（面间相互粘接） 22. AINA，NA1，NA2，…，NA8，NA9（被选面的交集） 23. AINP，NA1，NA2，…，NA8，NA9（面集两两相交） 24. AINV，NA，NV（面体相交） 25. AL，L1，L2，…，L9，L10（以线定义面） 26. ALIST，NA1，NA2，NINC，Lab（列表显示面的信息） 【注】Lab=HPT时，显示面上硬点信息，默认为空。 27. ALLSEL，LabT，Entity（选择所有实体） 【注】LabT=ALL（指定实体及其所有下层实体）、BELOW（指定实体及其下一层实体）； Entity=ALL、VOLU、AREA、LINE、KP、ELEM、NODE。 28. AMESH，NA1，NA2，NINC（划分面生成面单元） AMESH，AREA，KP1，KP2，KP3，KP4（通过点划分面单元） 29. /AN3D，Kywrd，KEY（三维注释） 30. ANCNTR，NFRAM，DELAY，NCYCL（在POST1中生成结构变形梯度线的动画） 31. ANCUT，NFRAM，DELAY，NCYCL，QOFF，KTOP，TOPOFF，NODE1，NODE2，NODE3（在P OST1中生成等势切面云图动画）
FINISH /CLEAR /PREP7 CSYS,0 !单元类型定义 ********************************************************************* ******** ET,1,BEAM4 !车道单元 ET,2,LINK10 !主缆单元 ET,3,BEAM4 !加劲梁单元 ET,4,LINK10 !吊索单元 ET,5,BEAM4 !中央扣单元 ET,6,BEAM4 !主塔单元 ET,7,BEAM4 !支撑桥面板小纵梁单元 ET,8,BEAM4 !桥面板单元 !材料性质定义 ********************************************************************* ******** MP,EX,1,3.5E10 !车道常数 MP,DENS,1,7850
MP,NUXY,1,0.3 MP,EX,2,2E11 !主缆材料常数MP,DENS,2,7850 MP,NUXY,2,0.3 KEYOPT,2,2,2 KEYOPT,2,3,0 MP,EX,3,2.1E11 !加劲梁材料常数MP,DENS,3,10000 MP,NUXY,3,0.3 MP,EX,4,2E11 !吊索材料常数MP,DENS,4,7850 MP,NUXY,4,0.3 KEYOPT,4,2,0 KEYOPT,4,3,0 MP,EX,5,2E11 !中央扣材料常数MP,DENS,5,7850 MP,NUXY,5,0.3 MP,EX,6,3.5E10 !主塔材料常数MP,DENS,6,2600 MP,NUXY,6,0.167
ANSYS软件APDL命令流建模的体会ANSYS软件APDL命令流建模的体会首先申明，本人学习ANSYS基本上是靠自己一点一点琢磨出来的，由于本人喜欢用APDL命令流，故总结出来的几点经验也就比较适合用APDL命令的朋友。 1、多看help，ANSYS的help为我们提供了很强大的功能，我最喜欢的是其中对各个命令有关参数的说明和解释部分，不管是建模、加载、后处理等，都可以通过apdl命令来实现。只要你知道命令，如“aatt ”，在help搜索栏输入“aatt”，回车，弹出aatt的有关页码，一般其中有一个只有“aatt”的一项，确认，即可看到你要查询的aatt命令的有关参数意义，本人常用的命令有: et---定义单元类型 mp---定义材料属性 k----建关键点， l----建线条 a---由关键点建立面 al---由线建立面 v----由关键点建立体 vl---由线建立体 va--由面建立体 lsel---在很多很多线中选择你需要的目标线，数量可以无限多…… asel---在很多很多面中选择你需要的目标面，数量也可以无限多…… vsel---在很多很多体中选择你需要的目标体，数量也可以无限多…… latt----给选中的线按材料编号赋属性(前提是首先已定义好材料) aatt---给选中的面按材料编号赋属性
vatt-----给选中的体按材料编号赋属性 acel---按坐标轴赋体积力， lmesh，amesh，vmesh---对线、面、体进行剖分 d---在节点上加约束边界 dl---在线上加载约束边界 da----在面上加载约束边界 2、以上只是列出了常见的几个命令，但是ansys提供的命令是很多的，我们不可能都记得，计算记得，也不知道其有关参数是如何定义的，那不要紧，我们可以与界面操作结合起来学习。我们先利用界面操作实现，然后在保存路径里面找到文件“file.log”，在该文件里有该操作等价的apdl命令，那以后我们就可以使用了。 3、复合命令，很多命令是复合命令，通过几个命令的组合以实现一定的目标，如FITEM、FLST等。这里不予以详述，大家可在学习中慢慢体会。 4、ansys提供的apdl语言可像fortain、c语言一样，可以编程，有条件语句、逻辑语句、文件读写等，但是这些语句语法有个特点，就是在相应的语句前要加“*”，以示其与以上apdl命令的区别。 以上只是一点小小的总结，希望对大家有帮助。 K, NPT, X, Y, Z Defines a keypoint. Npt: Reference number for keypoint. If zero, the lowest available number is assigned X,y,z: Keypoint location in the active coordinate system (may be R, θ, Z or R, θ, Φ). If X = P, gra phical picking is enabled and all other fields (including NPT) are ignored (valid only in the GUI).
Q235属性：弹性模量E=2.1e5N/mm2 密度=7.85e-6kg/mm3 泊松比=0.3 mp,ex,1,2.1e5 mp,prxy,1,0.3 mp,dens,1,7.85e-6 1,ksymm镜像点 2,arsym镜像面 3,kgen复制点 4.adele删除面 6，kdist,k1,k2测量两关键点的距离 7,adele,a,,,1删除area and below 8,创建圆柱面： circle创建圆 然后创建直线 然(轴线) 利用拉伸命令创建圆柱面creat__areas__by Lines adrag线拉伸成面modeling>operate>extrude>lines>>along lines VDRAG面拉伸成体modeling>operate>extrude>areas>>along lines ！创建空心圆柱体 CYLIND,RAD1,RAD2,Z1,Z2,THETA1,THETA2这个命令 Main Menu>Preprocessor>Modeling>Create>Volumes>Cylinder>By Dimensions Main Menu>Preprocessor>Trefftz Domain>TZ Geometry>Create>Volume>Cylinder>By Dimensions 9,aptn分割面 10,asbw用工作平面切割面 11.wpoffs 12.wprota
https://www.doczj.com/doc/774181303.html,ng过圆外一点做圆的切线(0°或180°) 14，nummrg将重复的点消除 15，asba面减去面 16,两个圆柱面的相贯线作法：做出两个相穿的圆柱面，利用APTN命令 17，选择面，不选择一部分面 asel,u,loc,z,kz(735) 18.在工作平面上生成一个矩形面 RECTING,X1,X2,Y1,Y2 X1，X2——矩形在工作平面X方向坐标值的变化范围 Y1,Y2——矩形在工作平面Y方向坐标值的变化范围 18,圆阵列 建立工作平面与圆柱的横截面平行，在工作平面情况下建立局部坐标系（柱坐标系），然后利用agen命令复制。 19，转换成局部柱坐标系 20,kfill在两个关键点之间生成一个或多个关键点 21.网格划分 aatt,1,14,1,!aatt,mat,real,type,esys,secn aesize,all,1000!aesize,anum,size,单元尺寸 mshape,0,2d!mshape,key,dimension指定划分单元形状amesh,all k,1,24000,33000,2230 k,2,24000,33000,-2230 k,3,-24000,33000,-2230 k,4,-24000,33000,2230 kfill,2,3,23,5,1,1 kfill,1,4,23,28,1,1 *do,i,5,26 l,i,i+1 *enddo
ANSYS最常用命令流+中文注释 VSBV, NV1, NV2, SEPO, KEEP1, KEEP2 —Subtracts volumes from volumes，用于2个solid相减操作，最终目的是要nv1-nv2=?通过后面的参数设置，可以得到很多种情况：sepo项是2个体的边界情况，当缺省的时候，是表示2个体相减后，其边界是公用的，当为sepo的时候，表示相减后，2个体有各自的独立边界。keep1与keep2是询问相减后，保留哪个体？当第一个为keep时，保留nv1,都缺省的时候，操作结果最终只有一个体，比如：vsbv,1,2,sepo,,keep,表示执行1-2的操作，结果是保留体2，体1被删除，还有一个1-2的结果体，现在一共是2个体（即1-2与2），且都各自有自己的边界。如vsbv,1,2,,keep,,则为1-2后，剩下体1和体1-2，且2个体在边界处公用。同理，将v换成a 及l是对面和线进行减操作！ mp,lab, mat, co, c1,…….c4 定义材料号及特性 lab: 待定义的特性项目（ex,alpx,reft,prxy,nuxy,gxy,mu,dens） ex: 弹性模量 nuxy: 小泊松比 alpx: 热膨胀系数 reft: 参考温度 reft: 参考温度 prxy: 主泊松比 gxy: 剪切模量 mu: 摩擦系数 dens: 质量密度 mat: 材料编号（缺省为当前材料号） co: 材料特性值，或材料之特性，温度曲线中的常数项 c1-c4: 材料的特性-温度曲线中1次项，2次项，3次项，4次项的系数 定义DP材料： 首先要定义EX和泊松比：MP，EX，MA T，…… MP，NUXY，MAT，…… 定义DP材料单元表（这里不考虑温度）：TB，DP，MA T 进入单元表并编辑添加单元表：TBDATA，1，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: 节点号
《史上最全的ANSYS命令流查询与解释》 【1】************************************************************************************* 对ansys主要命令的解释 1，/PREP7 ! 加载前处理模块 2，/CLEAR,NOSTART ! 清除已有的数据, 不读入启动文件的设置(不加载初始化文件)初始化文件是用于记录用户和系统选项设置的文本文件 /CLEAR, START !清除系统中的所有数据,读入启动文件的设置 /FILENAME, EX10.5 ! 定义工程文件名称 /TITLE, EX10.5 SOLID MODEL OF AN AXIAL BEARING ! 指定标题 4，F,2,FY,-1000 ! 在2号节点上施加沿着-Y方向大小为1000N的集中力 6，FINISH ! 退出模块命令 7，/POST1 ! 加载后处理模块 8，PLDISP,2 ! 显示结构变形图,参数“2”表示用虚线绘制出原来结构的轮廓 9，ETABLE,STRS,LS,1 ! 用轴向应力SAXL的编号”LS,1”定义单元表STRS ETABLE, MFORX,SMISC,1 ! 以杆单元的轴力为内容, 建立单元表MFORX ETABLE, SAXL, LS, 1 ! 以杆单元的轴向应力为内容, 建立单元表SAXL ETABLE, EPELAXL, LEPEL, 1 ! 以杆单元的轴向应变为内容, 建立单元表EPELAXL ETABLE,STRS_ST,LS,1 !以杆件的轴向应力“LS,1”为内容定义单元表STRS_ST ETABLE, STRS_CO, LS,1 !以杆件的轴向应力“LS,1”定义单元表STRS_CO ETABLE,STRSX,S,X ! 定义X方向的应力为单元表STRSX ETABLE,STRSY,S,Y ! 定义Y方向的应力为单元表STRSY *GET,STRSS_ST,ELEM,STEEL_E, ETAB, STRS_ST !从单元表STRS_ST中提取STEEL_E单元的应力结果，存入变量STRSS_ST; *GET, STRSS_CO,ELEM,COPPER_E,ETAB,STRS_CO”从单元表STRS_CO中提取COPPER_E单元的应力结果，存入变量STRSS_CO 10 FINISH !退出以前的模块 11, /CLEAR, START ! 清除系统中的所有数据,读入启动文件的设置 12 /UNITS, SI !申明采用国际单位制 14 /NUMBER, 2 !只显示编号, 不使用彩色 /NUMBER, 0 ! 显示编号, 并使用彩色 15 /SOLU ! 进入求解模块:定义力和位移边界条件,并求解 ANTYPE, STATIC ! 申明分析类型是静力分析(STA TIC或者0) OUTPR, BASIC, ALL ! 在输出结果中, 列出所有荷载步的基本计算结果 OUTPR,BASIC,ALL !指定输出所有节点的基本数据 OUTPR,BASIC,LAST ! 选择基本输出选项,直到最后一个荷载步 OUTPR,,1 ! 输出第1个荷载步的基本计算结果 OUTPR,BASIC,1 ! 选择第1荷载步的基本输出项目 OUTPR,NLOAD,1 ! 指定输出第1荷载步的内容 OUTRES,ALL,0 !设置将所有数据不记录到数据库。 NSUBST,1 ! 指定当前求解的荷载步 16 /AUTO,1 ! 设置模型显示的最佳比例 17 /VUP,1,X ! 设置X轴向上 /ANGLE,1,0 ! 水平轴夹角0度 18 SMRT,OFF ! 关闭智能化网格功能 18 MP, EX, 1, 207E9 ! 定义第1类材料弹性模量EX=207GPa
ANSYS APDL命令流悬臂梁分析教程 本文通过分析悬臂梁介绍了ANSYS APDL相关命令流方法。 考虑悬臂梁如图2-2，求x=L变形量。已知条件：杨氏系数E=200E9;截面参数：t=0.01m, w=0.03m, A=3E-4,I=2.5E-9;几何参数：L=4m, a=2m, b=2m;边界外力F=2N,q=0.05N/m. 使用ANSYS解决该问题的命令如下： /FILNAM,EX2-1? ! 定义文件名 /TITLE,CANTILEVER BEAM DEFLECTION? !定义分析的标题 /UNITS,SI !定义单位制(注意观察输出窗口的单位) /PREP7 !进入前置处理 ET,1,3? !定义元素类型为beam3 MP,EX,1,200E9 ! 定义杨氏模量 R,1,3E-4,2.5E-9,0.01 !定义实常数(要严格根据该元素类型的说明文档所给出的实常数格式) N,1,0,0!定义第1号节点X坐标为0，Y坐标为0 N,2,1,0!定义第2号节点X坐标为1，Y坐标为0 N,3,2,0 !定义第3号节点X坐标为2，Y坐标为0 N,4,3,0 !定义第4号节点X坐标为3，Y坐标为0 N,5,4,0!定义第5号节点X坐标为4，Y坐标为0 E,1,2!把1、2号节点相连构成单元，系统将自定义为1号单元 E,2,3!把2、3号节点相连构成单元，系统将自定义为2号单元 E,3,4!把3、4号节点相连构成单元，系统将自定义为3号单元 E,4,5!把4、5号节点相连构成单元，系统将自定义为4号单元 FINISH? !退出该处理层 /SOLU!进入求解处理器 D,1,ALL,0 !对1节点施加约束使它X，Y向位移都为0
第二章FLAC3D原理及入门指南 2-1定义一个FISH函数 new def abc abc = 25 * 3 + 5 End print abc 2-2使用一个变量 new def abc hh = 25 abc = hh * 3 + 5 End Print hh Print abc 2-3对变量和函数的理解 new def abc hh = 25 abc = hh * 3 + 5 End set abc=0 hh=0 print hh print abc print hh new def abc abc = hh * 3 + 5 end set hh=25 print abc set abc=0 hh=0 print hh print abc print hh
2-4获取变量的历史记录 new gen zone brick size 1 2 1 model mohr prop shear=1e8 bulk=2e8 cohes=1e5 tens=1e10 fix x y z range y -0.1 0.1 apply yvel -1e-5 range y 1.9 2.1 plot set rotation 0 0 45 plot block group def get_ad ad1 = gp_near(0,2,0) ad2 = gp_near(1,2,0) ad3 = gp_near(0,2,1) ad4 = gp_near(1,2,1) end get_ad def load load=gp_yfunbal(ad1)+gp_yfunbal(ad2)+gp_yfunbal(ad3)+gp_yfunbal(ad4) end hist load hist gp ydis 0,2,0 step 1000 plot his 1 vs -2 2-5用FISH函数计算体积模量和剪砌模量 new def derive s_mod = y_mod / (2.0 * (1.0 + p_ratio)) b_mod = y_mod / (3.0 * (1.0 - 2.0 * p_ratio)) end set y_mod = 5e8 p_ratio = 0.25 derive print b_mod print s_mod 2-6 在FLAC输入中使用符号变量 New def derive s_mod = y_mod / (2.0 * (1.0 + p_ratio)) b_mod = y_mod / (3.0 * (1.0 - 2.0 * p_ratio)) end set y_mod = 5e8 p_ratio = 0.25
ADAMS常用快捷方式 T+左键平动模型 R+左键旋转模型 Z+左键动态缩放 F或Ctrl+F 以最大比例全面显示模型 S+左键沿着垂直于屏幕的轴线旋转 W+左键将屏幕的一部分放大 C+左键定制旋转中心 E+左键将某构件的XY平面作为观察面 G 切换工作格栅的隐藏与显示 V 切换图标的隐藏与显示 M 打开信息窗口 Esc 结束当前得操作 F1 根据当前的状态，打开相应的帮助 F2 打开读取命令文件的对话框 F3 打开命令输入窗口 F4 打开坐标窗口 F8 进入后处理模块 Ctrl+N 新建数据库 Ctrl+O 打开数据库 Ctrl+S 保存数据库 Ctrl+P 打印 Ctrl+Q 退出ADAMS/View Ctrl+Z 取消上一步操作 Ctrl+Shift+z 恢复上一步的撤销操作 Ctrl+E 编辑一个元素 Ctrl+C 复制一个元素 Ctrl+X 删除一个元素 ADAMS宏命令大全 1. A，P1，P2，…，P17，P18（以点定义面） 2. AADD，NA1，NA2，…NA8，NA9（面相加） 3. AATT，MAT，REAL，TYPE，ESYS，SECN（指定面的单元属性）【注】ESYS为坐标系统号、SECN为截面类型号。 4. *ABBR，Abbr，String（定义一个缩略词） 5. ABBRES，Lab，Fname，Ext（从文件中读取缩略词） 6. ABBSA VE，Lab，Fname，Ext（将当前定义的缩略词写入文件）
7. ABS，IR，IA，--，--，Name，--，--，FACTA（取绝对值） 【注】************* 8. ACCAT，NA1，NA2（连接面） 9. ACEL，ACEX，ACEY，ACEZ（定义结构的线性加速度） 10. ACLEAR，NA1，NA2，NINC（清除面单元网格） 11. ADAMS，NMODES，KSTRESS，KSHELL 【注】************* 12. ADAPT, NSOLN, STARGT, TTARGT, FACMN, FACMX, KYKPS, KYMAC 【注】************* 13. ADD，IR, IA, IB, IC, Name, --,-- , FACTA, FACTB, FACTC（变量加运算） 14. ADELE，NA1，NA2，NINC，KSWP（删除面） 【注】KSWP =0删除面但保留面上关键点、1删除面及面上关键点。 15. ADRAG，NL1，NL2，…，NL6，NLP1，NLP2，…，NLP6（将既有线沿一定路径拖拉成面） 16. AESIZE，ANUM，SIZE（指定面上划分单元大小） 17. AFILLT，NA1，NA1，RAD（两面之间生成倒角面） 18. AFSURF，SAREA，TLINE（在既有面单元上生成重叠的表面单元） 19. *AFUN, Lab（指定参数表达式中角度单位） 20. AGEN, ITIME, NA1, NA2, NINC, DX, DY, DZ, KINC, NOELEM, IMOVE（复制面） 21. AGLUE，NA1，NA2，…，NA8，NA9（面间相互粘接） 22. AINA，NA1，NA2，…，NA8，NA9（被选面的交集） 23. AINP，NA1，NA2，…，NA8，NA9（面集两两相交） 24. AINV，NA，NV（面体相交） 25. AL，L1，L2，…，L9，L10（以线定义面） 26. ALIST，NA1，NA2，NINC，Lab（列表显示面的信息） 【注】Lab=HPT时，显示面上硬点信息，默认为空。 27. ALLSEL，LabT，Entity（选择所有实体） 【注】LabT=ALL（指定实体及其所有下层实体）、BELOW（指定实体及其下一层实体）； Entity=ALL、VOLU、AREA、LINE、KP、ELEM、NODE。 28. AMESH，NA1，NA2，NINC（划分面生成面单元） AMESH，AREA，KP1，KP2，KP3，KP4（通过点划分面单元） 29. /AN3D，Kywrd，KEY（三维注释） 30. ANCNTR，NFRAM，DELAY，NCYCL（在POST1中生成结构变形梯度线的动画） 31. ANCUT，NFRAM，DELAY，NCYCL，QOFF，KTOP，TOPOFF，NODE1，NODE2，NODE3（在POST1中生成等势切面云图动画） 32. ANDATA，DELAY，NCYCL，RSLTDAT，MIN，MAX，INCR，FRCLST，AUTOCNTRKY（生成某一范围内的结果数据的顺序梯度线动画）
ANSYS命令流使用技巧分享（收录汇总） 谈到Ansys使用技巧，不得不说APDL二次开发，针对二次开发并结合本人多年使用经验，有以下几点经验与各位分享。技巧毕竟很多，也欢迎各位专家留言补充，我们也可以整理汇总以待分享。 技巧一：ansys apdl语言高亮编辑器 命令流在编写时时常会把命令记错，如果写错了未察觉到，在计算时就会非常麻烦，因此一个帮助修正错误命令的编辑器必不可少。这种工具很多，我一直用的是UE，成功掌握二次开发必不可少。 blob.png 技巧二：建模画网格按照Number来区分各部件. 在ansys可以通过建立component来区分每个部件，有利于查看和编辑。采用下面的命令在建立模型和划分网格时，所有编号都从设定的值开始。 vsel,none asel,none lsel,none ksel,none ! NSS= NSS=150001 !设定值 NUMSTR,KP,NSS,
NUMSTR,LINE,NSS, NUMSTR,AREA,NSS, NUMSTR,VOLU,NSS, 技巧三：巧用循环语句*do和*enddo 可以采用较少的命令选择或者建立类似的模型，具体格式如下： *do,i,1,6 cmsel,a,JIECHU_pinghengmx_xia_%i% *enddo 技巧四:在计算结果中，实体云图中切面显示结果 命令流 /cplane /type 图形界面操作 <1.设置工作面为切面 <2.PlotCtrls-->Style-->Hidden line Options 将[/TYPE]选项选为section 将[/CPLANE]选项选为working plane 技巧五:某结点或单元的应力-应变关系曲线 由于手头没有相关的东西资料可以操作，待到有资料和时间的时候会补充上来图片。 1，定义变量： 拾取主菜单：Main Menu>Time Hist postproc>Define Variables>在随之弹出的对话框中点击Add键，定义第一个变量序号为2，选取第一个变量stress,确定与之对应的下一级选项（如Y-direction SY等）；返回定义变量对话框，再点击add键，定义第二个变量序号为3，选取第二个变量strain-elastic及以及对赢得下一级选项（如Y-dir'n EPEL Y等，在应力-应变图中，其向量的取向应相同）。同理再定义变量4，选取变量strain-plastic及与之对应的下一级选项如Y-dir'n EPEL Y等），在应力-应变图中，应变是弹性应变和塑性应变累加的总应变。为使其实现相加，还需进行以下操作：拾取主菜单：Main Menu>Time Hist