本帖最后由 zsl303466906 于 2015-1-15 23:23 编辑
我做的是瞬态分析 用的加速度荷载:
finish
/cle
/filename,suidao
/title,anlysis of suidao in 3d
/PREP7
et,1,combin14
r,1,30e6
et,2,solid45
et,3,plane42
mp,ex,2,2.1E8
mp,prxy,2,0.3
mp,dens,2,7800
k,1,0,0
k,2,0,0.3
k,3,0,3.3
k,4,0,3.6
k,5,0.3,3.6
k,6,0.3,3.3
k,7,0.3,0.3
k,8,0.3,0
k,9,3.3,0
k,10,3.3,0.3
k,11,3.3,3.3
k,12,3.3,3.6
k,13,3.6,3.6
k,14,3.6,3.3
k,15,3.6,0.3
k,16,3.6,0
k,17,6.6,0
k,18,6.6,0.3
k,19,6.6,3.3
k,20,6.6,3.6
k,21,6.9,3.6
k,22,6.9,3.3
k,23,6.9,0.3
k,24,6.9,0
l,1,2 !先连短的
l,3,4
l,8,7
l,6,5
l,9,10
l,11,12
l,16,15
l,14,13
l,17,18
l,19,20
l,24,23
l,22,21
l,1,8
l,9,16
l,17,24
l,2,7
l,10,15
l,18,23
l,3,6
l,11,14
l,19,22
l,4,5
l,12,13
l,20,21
l,2,3 !在链接长的
l,7,6
l,10,11
l,15,14
l,18,19
l,23,22
l,8,9
l,16,17
l,7,10
l,15,18
l,6,11
l,14,19
l,5,12
l,13,20
al,1,13,3,16 !6个四方角
al,5,14,7,17
al,9,15,11,18
al,2,19,4,22
al,6,20,8,23
al,10,21,12,24
al,16,26,19,25 !长边墙
al,17,28,20,27
al,18,30,21,29
al,4,35,6,37
al,3,31,5,33
al,38,8,36,10
al,34,7,32,9
al,26,33,27,35 !四方地板
al,28,34,29,36
lsel,s,line,,1,24,1 !选线
lesize,all,,,2 !分段
lsel,s,line,,25,38,1
lesize,all,,,20
allsel
type,3
mat,2
amesh,all !划网格
allsel
k,100,0,0,0.15 !画拉伸线
l,1,100
lsel,s,line,,39,39,1
lesize,all,,,1
asel,s,loc,z,0 !选择Z平面所有
type,2
mat,2
vdrag,all,,,,,,39 !拉伸
allsel
lsel,s,line,,39 !删除线和点
ldele,all
allsel
kdele,100,100,1
allsel
k,100,0,0,2.55 !画拉伸线
l,25,100
lsel,s,line,,39,39,1
lesize,all,,,16
asel,s,area,,20
asel,a,area,,25
asel,a,area,,30
asel,a,area,,35
asel,a,area,,40
asel,a,area,,45
asel,a,area,,48
asel,a,area,,51
asel,a,area,,54
asel,a,area,,57
asel,a,area,,60
asel,a,area,,63
asel,a,area,,66
type,2
mat,2
vdrag,all,,,,,,39
allsel
lsel,s,line,,39 !删除线和点
ldele,all
allsel
kdele,100,100,1
allsel
al,150,158,153,157
al,152,162,155,161
asel,s,area,,120,121,1
type,3
mat,2
mshape,0,2d
mshkey,1
amesh,all
allsel
k,100,0,0,2.7
l,49,100
lesize,39,,,1
asel,s,loc,z,2.55
type,2
mat,2
vdrag,all,,,,,,39
lsel,s,line,,39
ldele,all
allsel
kdele,100,100,1
allsel !盒子建成
allsel
*get,npmax,node,,num,max !得到已经定义的最大结点数,存入npmax
asel,s,loc,z,0
nsla,s,1
*get,n_nodes,node,,count !得到已经定义的节点个数
Node_next=0
*do,i,1,n_nodes
Node_next=NDNEXT(Node_next)
x=nx(node_next)
y=ny(node_next)
z=nz(node_next)
npmax=npmax+1 !定义法向节点
n,npmax,x-0.1,y,z
type,1
real,1
e,node_next,npmax
d,npmax,all,0 !约束新生成的点
npmax=npmax+1
n,npmax,x,y-0.1,z !定义切向结点以便与边界点形成切向单元
type,1
real,1
e,node_next,npmax
d,npmax,all,0 !约束新生成的点
npmax=npmax+1
n,npmax,x,y,z-0.1 !定义底边界切向结点以便与边界点形成切向单元
type,1
real,1
e,node_next,npmax
d,npmax,all,0 !约束新生成的点
*enddo
Allsel
finish
!定义和读入时程曲线
/config,nres,10000000 !定义最大子步
NT=100
DT=0.05 !时间间隔
*dim,ff,array,NT,1
*vread,ff(1,1),eqf,txt,,IJK,NT,1 !地震波采用 EL-Centro波(东西方向)
(1f10.4)
finish
/solu
NSUBST,1,,,1 !1个子步
OUTRES,ALL,1 !输出每个子步的结果
ANTYPE,TRANS !时程分析
*do,i,1,NT
asel,s,loc,z,0
nsla,s,1
acel,,,ff(i,1)
TIME,i*DT
allsel,all
solve
*enddo
eqf.txt是一个文件里面是加速度值都在0.1左右。大神你可以自己建一个放在ansys里面。
|