发一个拱桥地震响应的命令流,附带请教一问题

西风独自凉... 发表于 2010-5-31 08:25:50 | 显示全部楼层 | 阅读模式
本帖最后由 西风独自凉 于 2010-5-31 08:32 编辑

        小弟对ansys不是太熟练,做完硕士课题后,感觉学好ansys还是必要的,再说离开学还有一段时间,小弟就较为系统地学了下ansys,所用参考教材为:《AlgorAnsys在桥梁工程中的应用方法与实例》,《ANSYS.结构有限元高级分析方法与范例应用(第一版)》,《ANSYS.结构有限元高级分析方法与范例应用(第二版)》,《ANSYS土木工程应用实例》,《ansys高级工程有限元分析范例精选——祝效华》。这些教材,基本上学习ansys最好的教材了,但遗憾的是这些书中都有错误。如能将这些错误改正,对自己和别人都大有益处。
这些书中的错误有些地方小弟可以自己看出来,有些地方就无能为力了,在此,特向各位师长,学长,同仁请教。
首先,将自己发现并改正的例子发出来,供大家讨论。
ANSYS.结构有限元高级分析方法与范例应用(第一版)》第8章悬索拱桥地震荷载响应分析有错误。本人作了改动,模型和地震响应分析的完整命令流:
FINISH
/CLEAR
/FILENAME,bridge harmonic analysis
/TITLE,bridge harmonic analysis
/PREP7
et,1,beam4
ET,2,LINK10
et,3,shell63
r,1,0.6667,0.1667,1,2
r,2,0.0314,0.1
r,3,0.5,0.5,0.5,0.5,,,
!定义材料模型及参数
MP,EX,1,2.0e11
!
弹性模量
MP,PRXY,1,0.3

!
泊松比
MP,DENS,1,7800
!
密度
MP,EX,2,3.0e10
!
弹性模量
MP,PRXY,2,0.2
!
泊松比
MP,DENS,2,3000
!
密度
k,1,-60,-10,0
k,2,-50,0,0
k,3,-30,15,0
k,4,-10,20,0
k,5,10,20,0
k,6,30,15,0
k,7,50,0,0
k,8,60,-10,0
k,9,-70,0,
k,10,-60
k,11,-30
k,12,-10
k,13,10
k,14,30
k,15,60
k,16,70
kgen,2,all,,,,,-20,,0
*do,ii,1,7
l,ii,ii+1
l,ii+16,ii+17
*enddo
*do,ii,2,7
l,ii,ii+16
*enddo
l,1,9
l,1,10
l,17,25
l,17,26
l,8,15
l,8,16
l,24,31
l,24,32
l,1,18
l,2,17
l,8,23
l,7,24
*do,ii,3,6
l,ii,ii+8
l,ii+16,ii+24
*enddo
a,9,10,26,25
a,10,2,18,26
a,2,11,27,18
a,11,12,28,27
a,12,13,29,28
a,13,14,30,29
a,14,7,23,30
a,7,15,31,23
a,15,16,32,31
lsel,s,,,1,32
latt,1,1,1,,,,
lesize,all,2
lmesh,all
lsel,s,,,33,40
latt,1,2,2
lesize,all,,,1
lmesh,all
aatt,2,3,3
aesize,all,2
amesh,all
nsel,s,loc,y,-10
d,all,all
nsel,s,loc,x,-70
nsel,a,loc,x,70
d,all,ux
allsel,all
fini
!计算模态解
/sol
antype,modal
modopt,lanb,10
mxpand,10
solve
fini
!谱分析求解
/sol
antype,spectrum
spopt,sprs,10,0
svtyp,3,1,
sed,1,0,0
freq,0.5,1.0,1.5,2.0,2.5,0,0,0,0
sv,0.8,0.5,0.6,0.9,1.2,
solve
fini
!合并模态
/sol
antype,spectrum
srss,0.01,disp
solve
fini
!后处理
/post1
set,list
set,last
prnsol,dof
plnsol,u,x,0,1 !显示X方向位移云图
plnsol,u,y,0,1 !显示y方向位移云图
plnsol,u,z,0,1 !显示z方向位移云图
fini

                               
登录/注册后可看大图


                               
登录/注册后可看大图

                               
登录/注册后可看大图


ANSYS.结构有限元高级分析方法与范例应用(第一版)》第9章钢筋混凝土梁的弹塑性分析的错误本人始终找不出来,希望大家人尽其才,帮忙看看。
FINISH
/CLEAR
!****************************************************************
!钢筋混凝土梁的弹塑性分析命令流
!****************************************************************
/FILENAME,rc-beam
/TITLE,analysis OF rc-beam
/PREP7
et,1,link8
et,2,solid65
r,1,28.3e-6
r,2,113.1e-6
r,3,314.1e-6
r,4,
!4)定义材料模型及参数
MP,EX,1,2.1e11
!
弹性模量
MP,PRXY,1,0.3

!
泊松比
tb,bkin,1,1,2,1
tbdata,,2.1e8,0
MP,EX,2,2.55e10
!
弹性模量
MP,PRXY,2,0.3

!
泊松比
tb,kinh,2,1,8,
tbtemp,0
tbpt,,0.0001,0.936e6
tbpt,,0.0003,2.664e6
tbpt,,0.0006,4.896e6
tbpt,,0.0009,6.696e6
tbpt,,0.0012,8.064e6
tbpt,,0.0016,9.216e6
tbpt,,0.002,9.6e6
tbpt,,0.0033,9.6e6
tb,conc,2,1,9,
tbdata,,0.3,0.55,1.5e6,-1
block,0.0,3,0,0.4,0,-0.2
/view,1,1,2,3
/rep,fast
lplot
/pnum,lie,1
/replot
vatt,2,4,2
lsel,s,line,,4,5,1
lsel,a,line,,2,7,5
lesize,all,0.05
lsel,inve
lesize,all,0.025
lsel,all
vmesh,all
nsel,s,loc,x,0.125,0.175
nplot
/view,1,1
/rep,fast
/pnum,node,1
/replot
type,1
mat,1
real,1
*do,i,3349,3354,1
e,i,i+1
*enddo
*do,i,3355,3446,7
e,i,i+7
*enddo
*do,i,3452,3447,-1
e,i+1,i
*enddo
*do,i,3447,3356,-7
e,i,i-7
*enddo
esel,s,type,,1
egen,19,315,all,,,,,,,,0.15,0,0,
allsel,all
type,1
mat,1
real,2
*do,i,3145,9130,105
e,i,i+105
*enddo
esel,r,real,,2
egen,3,-3,all,,,,,,,,0,0,-0.075,
allsel,all
type,1
mat,1
real,3
*do,i,3237,9222,105
e,i,i+105
*enddo
esel,r,real,,3
egen,3,-3,all,,,,,,,,0,0,0.15,
allsel,all
esel,s,type,,1
eplot
/view,1,1,2,3
/rep
nsel,s,loc,y,0.19,0.21
nsel,r,loc,x,-0.01,0.01
d,all,all
allsel,all
nsel,s,loc,y,0.19,0.21
nsel,r,loc,x,2.99,3.01
d,all,all
allsel,all
fini
/solu
antype,0
nlgeom,1
nropt,full
eqslv,spar,,0,
time,1.0
autots,1
nsubst,100,200,50,1
kbc,0
cnvtol,u,,0.03,0,,
neqit,50,
pred,on,,on
outres,all,all
solve
fini
/post1
set,last
/dscale,1,100
/replot
pldisp,2
plnsol,u,y,0,1
plcrack,0,0
fini
精彩评论倒序浏览

4120查看2评论

西风独自凉... 发表于 2010-5-31 08:49:28
本帖最后由 西风独自凉 于 2010-5-31 08:51 编辑

有些地方小弟可以自己看出来,有些地方就无能为力了,在此,特向各位师长,学长,同仁请教。
举报 回复
幵开 发表于 2018-11-27 23:25:07
来学习一下
举报 回复
您需要登录后才可以回帖 登录 | 注册

本版积分规则

 
 
  • QQ:56984982
  • 点击这里给我发消息
    电话:13527553862
    站务咨询群桥头堡站务咨询桥梁专业交流群:
    中国桥梁专业领袖群
    工作时间
    8:00-18:00