马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
出错后查看.err文件有提示:
*** ERROR *** CP = 5.922 TIME= 16:10:12
Present time 0 is less than or equal to the previous time in a
transient analysis.
因为是用命令流做得*do循环加载-求解过程,所以怀疑是循环逻辑有问题,尤其是time设置可能有错,但是自己没能想明白,想请大侠帮忙查查看错哪了。。先谢过!
附图是计划的每步加载示意图
下面是我的solve:
finish
/solu
!下面把已有的几个.dat文件数据,定义成数组load_3m--load_48m(100,2),.dat文件格式为999行*2列,第一列是1-999,第二列是小数
!之后计划找出模型3-48m高度各三个节点,定义载荷F
*dim,load_3m,array,100,2
*vread,load_3m(1,1),load_3m,dat,,ijk,100,2
(f3.0,f6.4)
*dim,load_6m,array,100,2
*vread,load_6m(1,1),load_6m,dat,,ijk,100,2
(f3.0,f6.4)
*dim,load_9m,array,100,2
*vread,load_9m(1,1),load_9m,dat,,ijk,100,2
(f3.0,f6.4)
*dim,load_12m,array,100,2
*vread,load_12m(1,1),load_12m,dat,,ijk,100,2
(f3.0,f6.4)
*dim,load_15m,array,100,2
*vread,load_15m(1,1),load_15m,dat,,ijk,100,2
(f3.0,f6.4)
*dim,load_18m,array,100,2
*vread,load_18m(1,1),load_18m,dat,,ijk,100,2
(f3.0,f6.4)
*dim,load_21m,array,100,2
*vread,load_21m(1,1),load_21m,dat,,ijk,100,2
(f3.0,f6.4)
*dim,load_24m,array,100,2
*vread,load_24m(1,1),load_24m,dat,,ijk,100,2
(f3.0,f6.4)
*dim,load_27m,array,100,2
*vread,load_27m(1,1),load_27m,dat,,ijk,100,2
(f3.0,f6.4)
*dim,load_30m,array,100,2
*vread,load_30m(1,1),load_30m,dat,,ijk,100,2
(f3.0,f6.4)
*dim,load_33m,array,100,2
*vread,load_33m(1,1),load_33m,dat,,ijk,100,2
(f3.0,f6.4)
*dim,load_36m,array,100,2
*vread,load_36m(1,1),load_36m,dat,,ijk,100,2
(f3.0,f6.4)
*dim,load_39m,array,100,2
*vread,load_39m(1,1),load_39m,dat,,ijk,100,2
(f3.0,f6.4)
*dim,load_42m,array,100,2
*vread,load_42m(1,1),load_42m,dat,,ijk,100,2
(f3.0,f6.4)
*dim,load_45m,array,100,2
*vread,load_45m(1,1),load_45m,dat,,ijk,100,2
(f3.0,f6.4)
*dim,load_48m,array,100,2
*vread,load_48m(1,1),load_48m,dat,,ijk,100,2
(f3.0,f6.4)
!分析设置
ANTYPE,4
TRNOPT,FULL
NLGEOM,1
OUTRES,ALL,ALL
AUTOTS,0
KBC,1 !!!!!!采用step载荷增加方式
!约束底部三个节点
d,1,all,0
d,6,all,0
d,11,all,0
!加载
!循环意图是想让time自己,去循环取数组矩阵的第一列(序号列)作为当前载荷步终点时间,然后循环取所有数组的第i行第2列数加载,可能错在这
*do,i,1,3
time,load_3m(i,1)
nsubst,2,0,0
f,2,FY,load_3m(i,2) !定义3m载荷
f,7,FY,load_3m(i,2)
f,12,FY,load_3m(i,2)
f,4,FY,load_6m(i,2) !定义6m载荷
f,9,FY,load_6m(i,2)
f,14,FY,load_6m(i,2)
f,46,FY,load_9m(i,2) !定义9m载荷
f,50,FY,load_9m(i,2)
f,54,FY,load_9m(i,2)
f,48,FY,load_12m(i,2) !定义12m载荷
f,52,FY,load_12m(i,2)
f,56,FY,load_12m(i,2)
f,88,FY,load_15m(i,2) !定义15m载荷
f,92,FY,load_15m(i,2)
f,96,FY,load_15m(i,2)
f,90,FY,load_18m(i,2) !定义18m载荷
f,94,FY,load_18m(i,2)
f,98,FY,load_18m(i,2)
f,132,FY,load_21m(i,2) !定义21m载荷
f,145,FY,load_21m(i,2)
f,152,FY,load_21m(i,2)
f,136,FY,load_24m(i,2) !定义24m载荷
f,137,FY,load_24m(i,2)
f,139,FY,load_24m(i,2)
f,219,FY,load_27m(i,2) !定义27m载荷
f,232,FY,load_27m(i,2)
f,239,FY,load_27m(i,2)
f,223,FY,load_30m(i,2) !定义30m载荷
f,224,FY,load_30m(i,2)
f,226,FY,load_30m(i,2)
f,306,FY,load_33m(i,2) !定义33m载荷
f,319,FY,load_33m(i,2)
f,326,FY,load_33m(i,2)
f,310,FY,load_36m(i,2) !定义36m载荷
f,311,FY,load_36m(i,2)
f,313,FY,load_36m(i,2)
f,395,FY,load_39m(i,2) !定义39m载荷
f,410,FY,load_39m(i,2)
f,419,FY,load_39m(i,2)
f,397,FY,load_42m(i,2) !定义42m载荷
f,398,FY,load_42m(i,2)
f,400,FY,load_42m(i,2)
f,503,FY,load_45m(i,2) !定义45m载荷
f,513,FY,load_45m(i,2)
f,523,FY,load_45m(i,2)
f,507,FY,load_48m(i,2) !定义48m载荷
f,517,FY,load_48m(i,2)
f,527,FY,load_48m(i,2)
solve
*enddo
finish |