*** ERROR *** CP = 124.891 TIME= 21:02:01
Present time 1.E-06 is less than or equal to the previous time in a
transient analysis.
*** ERROR *** CP = 124.922 TIME= 21:02:11
The step data was checked and there were errors found.
Please check output or errors file ( E:\UIDL\file.err ) for messages.
(以下是命令流)
finish
/clear
/filname,kangzhenzuoye
/title,shuangtaxielaqiao
!加速度谱计算
tg=0.65
gra=9.8
kes=0.05
ci=0.34
cs=1.2
cd=1
a=0.15*gra
smax=2.25*ci*cs*cd*a
*dim,ttt,,22
*dim,fre,,22
*dim,ace,,22
ttt(1)=0,0.1,0.15,0.2,0.25,0.3,0.4,0.5,0.6,0.65
ttt(11)=0.7,0.8,0.9,1.0,1.5,2.0,2.5,3.0,4.0,6.0,8.0,10.0
*do,i,1,22
*if,ttt(i),le,0.1,then
ace(i)=smax*(5.5*ttt(i)+0.45)
*elseif,ttt(i),gt,0.1,and,ttt(i),le,0.3,then
ace(i)=smax*(1.0-2.5*(ttt(i)-0.1))
*elseif,ttt(i),gt,0.3,and,ttt(i),le,tg,then
ace(i)=smax*0.5
*elseif,ttt(i),gt,tg,then
ace(i)=smax*(tg/ttt(i))
*endif
*enddo
!求频率,从小到大排列
fre(22)=1e5
*do,i,2,22
fre(23-i)=1/ttt(i)
*enddo
*do,i,1,22
ttt(23-i)=ace(i)
*enddo
*do,i,1,22
ace(i)=ttt(i)
*enddo
!建模
/prep7
et,1,beam3
et,2,link1
mp,ex,1,2.1e11 !1号材料为梁的材料
mp,prxy,1,0.2
mp,dens,1,27343.75
mp,ex,2,2.1e11 !2号材料为塔的材料
mp,prxy,2,0.2
mp,dens,2,24000
mp,ex,3,2.1e11 !3号材料为索的材料
mp,prxy,3,0.3
mp,dens,3,7673.81
r,1,0.32,1.13,1 !主梁实常数
r,2,0.25,0.35,1 !主塔实常数
r,3,0.042,7000/0.042/2.1e8 !外索实常数
r,4,0.042,6000/0.042/2.1e8 !中索实常数
r,5,0.042,5000/0.042/2.1e8 !内索实常数
k,1
k,2,15
k,3,45
k,4,75
k,5,90
k,6,105
k,7,135
k,8,165
k,9,90,-25
k,10,90,30
k,11,90,45
k,12,90,60
*do,i,1,7
l,i,i+1
*enddo
*do,i,9,11
l,i,i+1
*enddo
l,2,12
l,3,11
l,4,10
l,6,10
l,7,11
l,8,12
lsymm,x,all
numcmp,kp
nummrg,kp
!主梁网分控制
lsel,s,loc,y,0
latt,1,1,1
lesize,all,1
!主塔网分控制
lsel,s,loc,x,90
lsel,a,loc,x,-90
latt,2,2,1
lesize,all,1
!外索网分控制
lsel,s,line,,11,16,5
lsel,a,line,,27,32,5
latt,3,3,2
lesize,all,,,1
!中索网分控制
lsel,s,line,,12,15,3
lsel,a,line,,28,31,3
latt,3,4,2
lesize,all,,,1
!内索网分控制
lsel,s,line,,13,14,1
lsel,a,line,,29,30,1
latt,3,5,2
lesize,all,,,1
!网分
allsel,all
lmesh,all
!施加约束
ksel,s,loc,y,-25
dk,all,all
ksel,s,loc,x,165
ksel,a,loc,x,-165
dk,all,uy
allsel,all
finish
!获得模态解
/solu
antype,modal
modopt,lanb,30
solve
finish
!时程分析
/CONFIG,NRES,20000
/CONFIG,STATU
/SOLU
allsel
ANTYPE,TRANS
DAMPRATIO=0.05
PI=3.1416
FREQ1=0.191132*2*PI
FREQ2=1.007*2*PI
ALPHAD,2*DAMPRATIO*FREQ1*FREQ2/(FREQ1+FREQ2) !定义质量阻尼系数
BETAD,2*DAMPRATIO/(FREQ1+FREQ2) !定义刚度阻尼系数
TRNOPT,FULL !完全法
TIMINT,OFF !关闭时间积分效应
NLGEOM,ON
TIME,1E-8
NSUBST,2
KBC,1
ACEL,0,9.8
SOLVE
TIMINT,ON
nt=1000
*dim,aca,,nt
*dim,act,,nt
*creat,dizhenbo,mac
*vread,ACa(1),tianjin,txt,,IJK,1000
(f16.6)
*vread,ACt(1),time,txt,,IJK,1000
(f16.6)
*end
/input,dizhenbo,mac
nsubst,10
outpr,all,none
outres,all,none
outres,nsol,all
*do,i,1,nt
acel,0,aca(i)
time,act(i)
solve
*enddo
/POST1
*dim,Tmax,,330
*dim,Tmin,,330
*DIM,T1,,1000,1
*DO,I,1,330
*DO,J,1,1000,1
SET,j
*GET,T1(j),ELEM,i,SMISC,6 !第J号单元的i端弯矩
*ENDDO
*VSCFUN,Tmax(I),MAX,T1
*VSCFUN,Tmin(I),Min,T1
*ENDDO |