《隧道及地下工程FLAC解析方法》全部命令流汇总
第二章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 val ue --’)
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 & size 3 10 6 4 dim 3 3 3 3 rat 1 1 1 1
gen zon radcyl p0 0 0 -1 p1 0 0 -8 p2 0 50 -1 p3 7 0 0 p4 0 50 -8 & p5 7 50 0 p6 7 0 -8 p7 7 50 -8 p9 point 1 p11 point 2 & size 3 10 6 4 dim 3 3 3 3 rat 1 1 1 1
;衬砌网格
sel shell range cyl end1 0 0 -1 end2 0 50 -1 rad 3
; 隧道外围边界网格
gen zone radtun p0 7 0 0 p1 50 0 0 p2 7 50 0 p3 15 0 50 p4 50 50 0 & p5 15 50 50 p6 50 0 50 p7 50 50 50 &
p8 23 0 0 p9 7 0 8 p10 23 50 0 p11 7 50 8 &
size 6 10 3 10 rat 1 1 1 1.1
gen zone brick p0 0 0 8 p1 7 0 8 p2 0 50 8 p3 0 0 50 &
p4 7 50 8 p5 0 50 50 p6 15 0 50 p7 15 50 50 &
size 3 10 10 rat 1 1 1.1
;对称得出1/2模型
gen zon reflect dip 0 ori 0 0 0 range x 0 23 y 0 50 z 8 50
gen zon reflect dip 0 ori 0 0 0 range x 23 50 y 0 50 z 0 50
;建立主隧道和辅助隧道块名称
group service range cyl end1 0 0 -1 end2 0 50 -1 rad 3
group main range cyl end1 15 0 0 end2 15 50 0 rad 4
;对称得出完整模型
gen zon reflect dip 90 dd 270 ori 0 0 0
;网格显示
plot set rotation 30 0 30
plot block group
4-7网格连接
new
gen zone brick size 4 4 4 p0 0,0,0 p1 4,0,0 p2 0,4,0 p3 0,0,2
gen zone brick size 8 8 4 p0 0,0,2 p1 4,0,2 p2 0,4,2 p3 0,0,4
attach face range z 1.9 2.1
model elas
prop bulk 8e9 shear 5e9
fix z range z -.1 .1
fix x range x -.1 .1
fix x range x 3.9 4.1
fix y range y -.1 .1
fix y range y 3.9 4.1
apply szz -1e6 range z 3.9 4.1 x 0,2 y 0,2
hist unbal
;网格显示
plot set rotation 30 0 30
plot block group
solve
save att.sav
plot cont zdisp outl on
单一网格(小)程序
new
gen zone brick size 8 8 8 p0 0,0,0 p1 4,0,0 p2 0,4,0 p3 0,0,4 model elas
prop bulk 8e9 shear 5e9
fix z range z -.1 .1
fix x range x -.1 .1
fix x range x 3.9 4.1
fix y range y -.1 .1
fix y range y 3.9 4.1
apply szz -1e6 range z 3.9 4.1 x 0,2 y 0,2
hist unbal
;网格显示
plot set rotation 30 0 30
plot block group
solve
save noattx.sav
plot cont zdisp outl on
单一网格(大)程序
new
gen zone brick size 4 4 4 p0 0,0,0 p1 4,0,0 p2 0,4,0 p3 0,0,4 model elas
prop bulk 8e9 shear 5e9
fix z range z -.1 .1
fix x range x -.1 .1
fix x range x 3.9 4.1
fix y range y -.1 .1
fix y range y 3.9 4.1
apply szz -1e6 range z 3.9 4.1 x 0,2 y 0,2
hist unbal
;网格显示
plot set rotation 30 0 30
plot block group
solve
save noattx.sav
plot cont zdisp outl on
4-8立方体洞穴网格生成
new
def parm
rad=4.0
len=10.0
in_size=6
rad_size=10
end
parm
gen zone radbrick edge len size in_size in_size in_size rad_size &
rat 1.0 1.0 1.0 1.2 dim rad rad rad
;对称得出完整模型
gen zon reflect dip 0 dd 0 ori 0 0 0
gen zon reflect dip 90 dd 270 ori 0 0 0
;网格显示
plot set rotation 30 0 30
plot block group
4-9球体洞穴网格生成
new
def parm
rad=4.0
len=10.0
in_size=6
rad_size=10
end
parm
def make_sphere
p_gp=gp_head
loop while p_gp#null
px=gp_xpos(p_gp)
py=gp_ypos(p_gp)
pz=gp_zpos(p_gp)
dist=sqrt(px*px+py*py+pz*pz)
if dist>0 then
k=rad/dist
ax=px*k
ay=py*k
az=pz*k
maxp=max(px,max(py,pz))
k=len/maxp
bx=px*k
by=py*k
bz=pz*k
u=(maxp-rad)/(len-rad)
gp_xpos(p_gp)=ax+u*(bx-ax)
gp_ypos(p_gp)=ay+u*(by-ay)
gp_zpos(p_gp)=az+u*(bz-az)
end_if
p_gp=gp_next(p_gp)
end_loop
end
gen zone radbrick edge len size in_size in_size in_size rad_size & rat 1.0 1.0 1.0 1.2 dim rad rad rad
make_sphere
;对称得出完整模型
gen zon reflect dip 0 dd 0 ori 0 0 0
gen zon reflect dip 90 dd 270 ori 0 0 0
;网格显示
plot set rotation 15 0 30
plot block group
4-10应力边界
new
gen zone brick size (4,4,4) p0 (0,0,0) p1 (4,0,0) &
p2 (0,4,0) p3 (2,0,3.464)
model elastic
prop bulk 1e8 shear .3e8
apply nstress -1e6 range plane dip 60 dd 270 origin 0.1,0,0 above step 1
;网格显示
plot set rotation 30 0 30
plot block group
plot grid red fap green
4-11改变应力边界条件
new
gen zone brick size 6 6 6
model elas
prop bulk 1e8 shear 7e7
fix x range x -0.1 0.1
def superstep
loop ns (1,n_steps)
x_stress = stress_inc
command
apply sxx add x_stress range x 5.9,6.1 y 0,6 z 0,2
step 100
end_command
end_loop
end
set n_steps=100 stress_inc=-1e3
hist zone sxx 6,0,0
plot create sxx_hist
plot show
Superstep
;网格显示
plot set rotation 30 0 30
plot block group
plot grid red fap green
new
gen zone brick size 6 6 6
model elas
prop bulk 1e8 shear 7e7
fix x range x -0.1 0.1
def x_stress
x_stress = stress_inc * step
end
set stress_inc = -1e3
apply sxx 1.0 hist x_stress range x 5.9,6.1 y 0,6 z 0,2 hist zone sxx 6,0,0
hist x_stress
step 100
;网格显示
plot set rotation 30 0 30
plot grid red fap green
new
gen zone brick size 6 6 6
model elas
prop bulk 1e8 shear 7e7
fix x range x -0.1 0.1
table 1 0,0 100,-1e5
apply sxx 1.0 hist table 1 range x 5.9,6.1 y 0,6 z 0,2 hist zone sxx 6,0,0
step 100;网格显示
plot set rotation 30 0 30
plot grid red fap green
new
gen zone brick size 6,6,6 p1 6,0,-1
model elas
prop bulk 8e9 shear 5e9
apply sxx -2e6 range x -0.1 0.1
apply sxx -2e6 range x 5.9 6.1
step 500
plot grid disp
new
gen zone brick size 5,5,5
model elas
prop bulk 8e9 shear 5e9
set grav 0 0 -10
ini dens 1000
fix x range x -.1 .1
fix x range x 4.9 5.1
fix y range y -.1 .1
fix y range y 4.9 5.1
ini szz -5e4 grad 0 0 -1e4
app szz -5e4 range z -.1 .1
solve
model null range x 1,4 y 1 4 z 3 5
step 100
plot set plane dip 90 dd 180 origin 0,2.5,0
plot add surf plane behind yell
plot add vel plane behind black
plot set rotation 30 0 30
plot block group
plot show
4-12位移边界
new
gen zone brick size 4,4,4 p3 2,0,3.464
model elas
prop bulk 1e8 shear .3e8
macro left_boun ’plane dip 60 dd 270 origin 0.1,0,0 above’macro right_boun ’plane dip 60 dd 270 origin 3.9,0,0 below’apply nvel 0.1 plane dip 60 dd 270 range left_boun
apply nvel 0.1 plane dip 120 dd 90 range right_boun
step 1
;网格显示
plot set rotation 30 0 30
plot block group
plot grid vel
new
new
gen zone brick size 4 4 4
model elastic
prop bulk 1e8 shear .3e8
apply nstress -1e6 plane dip 0 dd 0 range y 3.9 4.1
apply nvel 0.0 plane dip 0 dd 0 range y -.1 .1
apply nvel 0.0 plane dip 0 dd 0 range x -.1 .1
step 10
;网格显示
plot set rotation 30 0 30
plot block group
plot grid vel
new
gen zone brick size 10 5 5
mod el
prop shear 1e8 bulk 2e8
fix x y z range x -.1 .1 y 0 5 z 0 5
fix x y z range x 0 10 y 0 5 z -.1 .1
fix x y z range x 9.9 10.1 y 0 5 z 0 5
table 1
def find_add
head = null
p_gp = gp_head
loop while p_gp # null
x_pos = gp_xpos(p_gp)
if x_pos = width then
new = get_mem(2)
mem(new) = head
mem(new+1) = p_gp
head = new
endif
p_gp = gp_next(p_gp)
endloop
end
set width=10.0
find_add
def apply_vel
while_stepping
ad = head
loop while ad # null
p_gp = mem(ad+1)
gp_xvel(p_gp) = vel_max * gp_zpos(p_gp) / height
gp_zvel(p_gp) = -vel_max * (gp_xpos(p_gp) - width) / height ad = mem(ad)
endloop
end
set large vel_max=1e-2 height=5.0
step 100
;网格显示
plot set rotation 30 0 30
plot block group
plot grid vel
4-13不考虑重力影响的均匀应力
new ;set grav 0 0 0
gen zone brick size 6 6 6
model elas
ini sxx=-5e6 syy=-1e7 szz=-2e7
apply sxx=-5e6 range x -0.1 0.1
apply sxx=-5e6 range x 5.9 6.1
apply syy=-1e7 range y -0.1 0.1
apply syy=-1e7 range y 5.9 6.1
apply szz=-2e7 range z -0.1 0.1
apply szz=-2e7 range z 5.9 6.1
prop bulk 8e9 shear 5e9
ini dens 2000
fix x range z -.1 .1
step 10
plot set rotation 30 0 30 ;网格三维显示
plot cont zdisp outl on
4-14考虑应力梯度的均匀材料
new
gen zone brick size 10 10 10 p1 20,0,0 p2 0,20,0 p3 0,0,20 model mohr
prop bulk 5e9 shear 3e9 fric 35
ini density 2500
set gravity 0,0,-10
fix x y z range z -0.1 0.1
ini szz -5.0e6 grad 0,0,2.5e4
ini syy -2.5e6 grad 0,0,1.25e4
ini sxx -2.5e6 grad 0,0,1.25e4
apply szz -4.5e6 range z 19.9 20.1
apply szz -5.0e6 range z -0.1 0.1
apply sxx -2.5e6 grad 0,0,1.25e4 range x -0.1 0.1
apply sxx -2.5e6 grad 0,0,1.25e4 range x 19.9 20.1
apply syy -2.5e6 grad 0,0,1.25e4 range y -0.1 0.1
apply syy -2.5e6 grad 0,0,1.25e4 range y 19.9 20.1
step 10
;网格三维显示
plot set rotation 30 0 30
plot cont zdisp outl on
4-15考虑应力梯度的非均匀材料
new
gen zone brick size 10 10 10 &
p0 0,-25,0 p1 20,-25,0 p2 0,0,0 p3 0,-25,20
model elas
prop bulk 5e9 shear 3e9
ini density 1600 range y -10,0
ini density 2000 range y -15,-10
ini density 2200 range y -25,-15
set gravity 0,-10,0
fix x range x -.1 .1
fix x range x 19.9 20.1
fix z range z -.1 .1
fix z range z 19.9 20.1
fix y range y -25.1 -24.9
ini syy 0.0 grad 0,1.6e4,0 range y -10,0
ini syy 4e4 grad 0,2.0e4,0 range y -15,-10
ini syy 7e4 grad 0,2.2e4,0 range y -15,-25
step 10
;网格三维显示
plot set rotation 30 0 30
plot cont zdisp outl on
4-16非均匀网格应力初始化
new
gen zone radcyl size 3 8 4 5 fill p1 10,0,0 p2 0,10,0 p3 0,0,10 mode elastic
prop shear 3e8 bulk 5e8
fix x range x -.1 .1
fix x range x 9.9 10.1
fix y range y -.1 .1
fix y range y 9.9 10.1
fix z range z -.1 .1
ini szz = -2.5e5 grad 0,0,2.5e4
ini density 2500
set grav 0,0,-10
step 10
;网格三维显示
plot set rotation 30 0 30
plot cont zdisp outl on
4-17不规则自由面应力初始化
gen zone brick size 15 15 10 p0 0,0,0 edge=100.0
model elastic
prop shear 3e8 bulk 5e8
def mountain
gp = gp_head
loop while gp # null
zz = sqrt(gp_xpos(gp)?2 + gp_ypos(gp)?2)
dz = 0.06 * sin(0.2 * zz + 100.0) ; Sum Fourier terms for
dz = dz + 0.06 * sin(0.22 * zz - 20.3) ; quasi-random surface
dz = dz - 0.04 * sin(0.33 * zz + 33.3) ; topology.
gp_zpos(gp) = 0.5 * gp_zpos(gp) * (1.0 + dz)
gp = gp_next(gp)
end_loop
end
mountain
fix x range x -.1 .1
fix x range x 99.9 100.1
fix y range y -.1 .1
fix y range y 99.9 100.1
fix z range z -.1 .1
set grav 0,0,-10
ini density=2000
ini szz=-2.0e6 (grad 0,0,2.0e4) sxx=-4.0e6 (grad 0,0,4.0e4) syy=-4.0e6 (grad 0,0,4.0e4) step 100
plot set rotation 30 0 30
plot cont zdisp outl on;网格三维显示
4-18非均网格内部压实
new
gen zone brick size 8 8 10 ratio 1.2 1 1
model mohr
ini dens 2000
prop bulk 2e8 shear 1e8
prop fric 30
fix x range x -.1 .1
fix x range x 7.9 8.1
fix y range y -.1 .1
fix y range y 7.9 8.1
fix z range z -.1 .1
八天学会Ansys命令流
八天学会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,,
FLAC命令流
FLAC命令流 ;**************************************** group命令 group soil range z 1 2 ;定义group prop bulk 7.8e6 shear 3.0e6 coh 10e3 fric 15 ran group soil ;给group赋值model null range group soil model elastic range group soil range命令 range name trench x 0 1 y 0 4 z 0 2 model null range trench model null range x=2,4 y=2,6 z=5,10 range name Big_Brick x -3 3 y -2 2 z -1 1 model elastic range Big_Brick prop bulk 1e8 shear 1e8 range Big_Brick range name Layer1 plane dip 0 dd 0 ori 0 0 0 above range name Layer2 plane dip 0 dd 0 ori 0 0 0 below macro命令 macro Sand 'bulk 1e8 shear 0.5e8 coh 0 tens 0 fric 35' macro Clay 'bulk 1e7 shear 0.3e7 coh 1e7 tens 0 fric 0' prop sand range Layer1 prop clay range layer2 macro Pt0 'p0 0 0 0' macro Pt1 'p1 add 10 0 0' macro Pt2 'p2 add 0 10 0' macro Pt3 'p3 add 0 0 10' macro Model_Size 'size 4 5 6' macro Big_Brick 'gen zone brick Pt0 Pt1 Pt2 Pt3 Model_Size' Big_Brick macro 'Pt0' 'p0 15 15 15' gen Big_Brick ; this will cause an error
ANSYS学习心得
一学习ANSYS需要认识到的几点 相对于其他应用型软件而言,ANSYS作为大型权威性的有限元分析软件,对提高解决问题的能力是一个全面的锻炼过程,是一门相当难学的软件,因而,要学好ANSYS,对学习者就提出了很高的要求,一方面,需要学习者有比较扎实的力学理论基础,对ANSYS分析结果能有个比较准确的预测和判断,可以说,理论水平的高低在很大程度上决定了ANSYS使用水平;另一方面,需要学习者不断摸索出软件的使用经验不断总结以提高解决问题的效率。在学习ANSYS的方法上,为了让初学者有一个比较好的把握,特提出以下五点建议:(1)将ANSYS的学习紧密与工程力学专业结合起来 毫无疑问,刚开始接触ANSYS时,如果对有限元,单元,节点,形函数等《有限元单元法及程序设计》中的基本概念没有清楚的了解话,那么学ANSYS很长一段时间都会感觉还没入门,只是在僵硬的模仿,即使已经了解了,在学ANSYS之前,也非常有必要先反复看几遍书,加深对有限元单元法及其基本概念的理解。 作为工程力学专业的学生,虽然力学理论知识学了很多,但对许多基本概念的理解许多人基本上是只停留于一个符号的认识上,理论认识不够,更没有太多的感性认识,比如一开始学ANSYS时可能很多人都不知道钢材应输入一个多大的弹性模量是合适的。而在进行有限元数值计算时,需要对相关参数的数值有很清楚的了解,比如材料常数,直接关系到结果的正确性,一定要准确。实际上在学ANSYS时,以前学的很多基本概念和力学理论知识都忘得差不多了,因而遇到有一
定理论难度的问题可能很难下手,特别是对结果的分析,需要用到《材料力学》,《弹性力学》和《塑性力学》里面的知识进行理论上的判断,所以在这种情况下,复习一下《材料力学》,《弹性力学》和《塑性力学》是非常有必要的,加深对基本概念的理解,实际上,适当的复习并不要花很多时间,效果却很明显,不仅能勾起遥远的回忆,加深理解,又能使遇到的问题得到顺利的解决。 在涉及到复杂的非线性问题时(比如接触问题),一方面,不同的问题对应着不同的数值计算方法,求解器的选择直接关系到程序的计算代价和问题是否能顺利解决;另一方面,需要对非线性的求解过程有比较清楚的了解,知道程序的求解是如何实现的。只有这样,才能在程序的求解过程中,对计算的情况做出正确的判断。因此,要能对具体的问题选择什么计算方法做出正确判断以及对计算过程进行适当控制,对《计算方法》里面的知识必须要相当熟悉,将其理解运用到ANSYS的计算过程中来,彼此相互加强理解。要知道ANSYS是基于有限元单元法与现代数值计算方法的发展而逐步发展起来的。因此,在解决非线性问题时,千万别忘了复习一下《计算方法》。此外,对《计算固体力学》也要有所了解(一门非常难学的课),ANSYS对非线性问题处理的理论基础就是基于《计算固体力学》里面所讲到的复杂理论。 作为学工程力学的学生,提高建模能力是非常急需加强的一个方面。在做偏向于理论的分析时,可能对建模能力要求不是很高,但对于实际的工程问题,有限元模型的建立可以说是一个最重要的问题,而后
ANSYS命令流解释大全
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 !选择所有单元
FLAC3D 实例命令流1
第1部分命令流按照顺序进行 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)
ANSYS中的APDL命令总结
在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: 单元类型号
ansys命令流解释
对ansys主要命令的解释 本文给出了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 ! 申明分析类型是静力分析(STATIC或者0)
Ansys命令流大全(整理)
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 在参数表达式中,为角度函数指定单位.
flac 命令流 参考 实用
;**************************************** group命令 group soil range z 1 2 ;定义group prop bulk 7.8e6 shear 3.0e6 coh 10e3 fric 15 ran group soil ;给group赋值 model null range group soil model elastic range group soil plot block group range group dam ;只显示dam range命令 range name trench x 0 1 y 0 4 z 0 2 model null range trench model null range x=2,4 y=2,6 z=5,10 range name Big_Brick x -3 3 y -2 2 z -1 1 model elastic range Big_Brick prop bulk 1e8 shear 1e8 range Big_Brick range name Layer1 plane dip 0 dd 0 ori 0 0 0 above range name Layer2 plane dip 0 dd 0 ori 0 0 0 below range cylinder end1 x1 y1 z1 end2 x2 y2 z2 radius r cylindrical range with one end of the cylinder axis (end1) at location (x1, y1, z1), the other end (end2) at location (x2, y2, z2), and with a cylinder radius of r '由(x1, y1, z1) 、(x2, y2, z2)两点确定旋转轴 macro命令 macro Sand 'bulk 1e8 shear 0.5e8 coh 0 tens 0 fric 35' macro Clay 'bulk 1e7 shear 0.3e7 coh 1e7 tens 0 fric 0' prop sand range Layer1 prop clay range layer2 macro Pt0 'p0 0 0 0' macro Pt1 'p1 add 10 0 0' macro Pt2 'p2 add 0 10 0' macro Pt3 'p3 add 0 0 10' macro Model_Size 'size 4 5 6'
个人总结ansys命令流
Q235 属性:弹性模量E=2.1e5 N/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.360docs.net/doc/8413103454.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 命令流解释大全
一、定义材料号及特性 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 !选择所有结点
FLAC3D基础命令流解释
;模型镜像 gen zone radcylinder size 25 1 25 25 gen zone reflect normal -1 0 0 origin x y z(面上一点);沿X轴镜像,通过对称平面法线向量确定对称面 gen zone reflect normal 0 0 -1 ;沿z轴镜像 ;绘图控制 pl
ANSYS命令流学习笔记workbench中命令流的一些应用
!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——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中生成等势切面云图动画)
flacd软件隧道支护与开挖命令流
f l a c d软件隧道支护与开 挖命令流 The latest revision on November 22, 2020
;----------------------------------------------------- ; ---- Excavation and Support for a Shallow Tunnel --- ;----------------------------------------------------- new ;新建项目 set fish autocreate off title 'Excavation and Support for a Shallow Tunnel' ;定义题目 ; generate primitive components of grid ; concrete liner - upper tunnel gen zon cshell p0 0 0 0 p1 7 0 0 p2 0 51 0 p3 0 0 5.5 & dim 5 5 5 5 size 2 51 10 group zone 'concrete liner' ; ; upper tunnel gen zon cylinder p0 0 0 0 p1 5 0 0 p2 0 51 0 p3 0 0 5 & size 5 51 10 group zone tunnel range group 'concrete liner' not ; ; lower tunnel & liner gen zone brick p0 0 0 -4.5 p1 add 7 0 0 p2 add 0 51 0 p3 add 0 0 4.5 & size 7 51 3 ; ; surrounding rock (8 primitives) gen zon radcyl p0 0 0 0 p1 27 0 0 p2 0 51 0 p3 0 0 25 & dim 7 5.5 7 5.5 size 5 51 10 8 rat 1 1 1 1.3 ; gen zone brick p0 7 0 -4.5 p1 27 0 -15 p2 add 0 51 0 p3 7 0 0 & p4 27 51 -15 p5 7 51 0 p6 27 0 0 p7 27 51 0 & size 8 51 3 ratio 1.3 1 1 ; gen zone brick p0 0 0 -15 p1 add 27 0 0 p2 add 0 51 0 p3 0 0 -4.5 & p4 27 51 -15 p5 0 51 -4.5 p6 7 0 -4.5 p7 7 51 -4.5 & size 7 51 8 rat 1 1 0.7692307692307692 ; gen zon brick p0 0 0 25 p1 add 27 0 0 p2 add 0 51 0 p3 add 0 0 10 & size 5 51 2 ; gen zon bric p0 27 0 25 p1 add 17 0 0 p2 add 0 51 0 p3 add 0 0 10 & size 2 51 2 rat 2 1 1 ;
ansys命令流
第一天目标: 熟悉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/TILE,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,,K,3,50,10,,K,4,10,10,,K,5,10,50,,K,6,0,50,,A,1,2,3,4,5,6,!由关键点生成面......!划分网格ESIZE,1,0,AMESH, 1......FINISH!前处理结束标识/SOLU!进入求解模块标识!施加约束和载荷DL,5,,ALLSFL,3,PRES,1000SFL,2,PRES, 1000......SOLVE!求解标识FINISH!求解模块结束标识/POST1!进入通用后处理器标识....../POST26!进入时间历程后处理器……/EXIT,SAVE!退出并存盘以下是日志文件中常出现的一些命令的标识说明,希望能给大家在整理LOG文件时有所帮助/ANGLE!指定绕轴旋转视图/DIST!说明对视图进行缩放/DEVICE!设置图例的显示,如: 风格,字体等/REPLOT!重新显示当前图例/RESET!恢复缺省的图形设置/VIEW!设置观察方向/ZOOM!对图形显示窗口的某一区域进行缩放第三天生成关键点和线部分 1.生成关键点K,关键点编号,X坐标,Y坐标,Z坐标例:
FLAC3D命令流(整理版)
实用标准文档 1、怎样查看模型? 答:plot grid 可以查看网格,plot grid num 可以查看节点号。 2、请问在圆柱体四周如何施加约束条件? 答:可以用fix ... ran cylinder end1 end2 radius r1 cylinder end1 end2 radius r2 not,其中r2