新手,起步阶段。 #b)`as?!1
准备做路基在列车动荷载作用下的动力响应。 ~@Bw(!
先试着建了一个简单的10米*10米*10米的地基上作用一个正弦荷载的例子 Nc[N 11?O
问题不少。。看了陈育民的书,收获不少,但问题更多 ;<T,W[3J
1.首先动力计算时的dyntime每一步的计算时间步怎么放大或缩小?陈的书建议不好修改,但太小了,计算时间太长。 `S4*~Xx
2.自己试算的模型最大不平衡力太大,不收敛。找不出问题。 ]gnEo.R
3.动力边界条件不知设置对否? 3rHn?
4.动力荷载的施加,面荷载可以仿照书上的,但点荷载不会施加,apply szz -1 hist wave range 去掉range 直接写上坐标位置,程序显示不对,加了后,也不知道荷载具体在哪个位置。不想ABAQUS那样明显易得。 icIWv
' e!WZvr
new C .B=E"e
;第一步:静力计算阶段 Z
Q*hrgQ
config dyn N8kNi4$mp=
set dyn off e, 2/3jO
gen zone brick p0 0 0 0 p1 10 0 0 p2 0 10 0 p3 0 0 10 size 20 20 20;建立一个10*10*10的地基模型,0.5米一个网格 V'dw=W17V
60ciI,_`
model mohr m##!sF^k~J
prop bulk 5.33e7 shear 3.077e7 coh 5.41e4 fri 36.7 A\9LJ#E
ini dens 2000 NLx TiyQy
set grav 0 0 -9.81 0uM&F[.x@g
fyT|xI`iD
fix x y z range z -0.01 0.01 ;地基底部固定 uJ0'`Q?6R9
fix x range x -0.1 0.1 JJg;X :p
fix x range x 9.99 10.01 ;地基左右两侧x方向固定 nvwf!iU6
fix y range y -0.01 0.01 { Dm@_&
fix y range y 9.99 10.01 ;地基前后纵向y方向固定 >L "+8N6
hist unbal b?,%M^9\`
solve Z 1wtOL
ini xdis 0 ydis 0 zdis 0 ;将节点位移清零 xj8z*fC;
ini xvel 0 yvel 0 zvel 0 ;将节点速度清零 B:UM2Jl
;第二步:动力计算阶段 qgfP6W$
set dyn on KlS#f
def wave [r Nd7-j <
wave=sin(2*pi*15*dytime) GB}=
end t~4Cf])
;定义一个正弦波,频率15HZ Fkpaou
free x y z range z -0.1 0.1 ;去掉模型底部的静力条件 4bw4!z9G
apply nquiet squiet dquiet ran z -0.1 0.1 ;静态边界条件 Yd~J(
apply szz -1 hist wave range x 5 6 y 5 6 z 9.99 10.01 ;地基面施加一个动力荷载 nJYIkfdA
apply ff ;施加自由场边界条件 Q1yXdw
set dyn time = 0 ;设置动力计算从0s开始 ! N!pvK;
hist reset ;清空已有的历史信息 (xTGt",_Jo
hist unbal 3'i(wI~<[
hist dytime {fV$\^c
hist wave $'eY-U8q
solve age 5 ,H.5TQ#
save 11.sav