新手,起步阶段。 c,O;B_}M]
准备做路基在列车动荷载作用下的动力响应。 \&~YFj B
先试着建了一个简单的10米*10米*10米的地基上作用一个正弦荷载的例子 uUV"86B_
问题不少。。看了陈育民的书,收获不少,但问题更多 #eX<=H]
1.首先动力计算时的dyntime每一步的计算时间步怎么放大或缩小?陈的书建议不好修改,但太小了,计算时间太长。 1Wiz0X/
2.自己试算的模型最大不平衡力太大,不收敛。找不出问题。 \98N8p;,I
3.动力边界条件不知设置对否? /DP0K
@%
4.动力荷载的施加,面荷载可以仿照书上的,但点荷载不会施加,apply szz -1 hist wave range 去掉range 直接写上坐标位置,程序显示不对,加了后,也不知道荷载具体在哪个位置。不想ABAQUS那样明显易得。 HA7%8R*.2i
"Q:h[) a
new h-RhmQA=Iz
;第一步:静力计算阶段 ,fqM>Q
config dyn 29CzG0?B
set dyn off aGfp"NtL
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米一个网格 luyu7`
TD%L`Gk
model mohr ,LP^v'[V7
prop bulk 5.33e7 shear 3.077e7 coh 5.41e4 fri 36.7 M< .1U?_#
ini dens 2000 &,=FPlTC=
set grav 0 0 -9.81 O&&_)
cx ("F/Jm
fix x y z range z -0.01 0.01 ;地基底部固定 hgzNEx%^q
fix x range x -0.1 0.1 xRdx`
YY u
fix x range x 9.99 10.01 ;地基左右两侧x方向固定 X;2LK!x;y
fix y range y -0.01 0.01 m
OE!`fd
fix y range y 9.99 10.01 ;地基前后纵向y方向固定 U=7nz|
hist unbal gFrNk
Uqp
solve =9W\;xE S
ini xdis 0 ydis 0 zdis 0 ;将节点位移清零 La3rX
ini xvel 0 yvel 0 zvel 0 ;将节点速度清零 4Gh\T`=
;第二步:动力计算阶段 mlCg&fnDB
set dyn on $4~Z]-38#A
def wave F_m'
9KX4E
wave=sin(2*pi*15*dytime) Jz3,vVfQ:
end 0 q1x+
;定义一个正弦波,频率15HZ bM5o-U#^ C
free x y z range z -0.1 0.1 ;去掉模型底部的静力条件 rZcSG(d`53
apply nquiet squiet dquiet ran z -0.1 0.1 ;静态边界条件 VumM`SH
apply szz -1 hist wave range x 5 6 y 5 6 z 9.99 10.01 ;地基面施加一个动力荷载 mwyB~,[d+W
apply ff ;施加自由场边界条件 Jp%5qBS^
set dyn time = 0 ;设置动力计算从0s开始 hvQXYo>TZx
hist reset ;清空已有的历史信息 6/?onEL9_
hist unbal v;U5[
hist dytime gu%i|-}
hist wave V0#Ocq,
solve age 5 f8aY6o"i
save 11.sav