找回密码
 注册
Simdroid-非首页
查看: 115|回复: 2

[其他] 在地震时程求解时出错,求解答!

[复制链接]
发表于 2011-6-26 21:10:42 | 显示全部楼层 |阅读模式 来自 湖南长沙
*** 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
发表于 2011-9-7 13:20:10 | 显示全部楼层 来自 广东东莞
Simdroid开发平台
没见到 tianjin.txt 和 time.txt 两个文件
自己写了随便写了一个 初始时间>1e-6
运行完全没有问题
请楼主注意时间的这个序列必须是升序就可以了
回复 不支持

使用道具 举报

发表于 2013-3-20 14:31:22 | 显示全部楼层 来自 重庆沙坪坝区
TIME,1E-8
有问题,从0开始试一下。
回复 不支持

使用道具 举报

您需要登录后才可以回帖 登录 | 注册

本版积分规则

Simapps系列直播

Archiver|小黑屋|联系我们|仿真互动网 ( 京ICP备15048925号-7 )

GMT+8, 2024-9-29 13:21 , Processed in 0.033506 second(s), 14 queries , Gzip On, MemCache On.

Powered by Discuz! X3.5 Licensed

© 2001-2024 Discuz! Team.

快速回复 返回顶部 返回列表