[求助]恳请牛人进来看一下我的动力分析哪里出错了。谢谢 [复制链接]

新手，起步阶段。 #U14-
^7  
准备做路基在列车动荷载作用下的动力响应。 &jmRA';sK  
先试着建了一个简单的10米*10米*10米的地基上作用一个正弦荷载的例子 O%bEB g  
问题不少。。看了陈育民的书，收获不少，但问题更多
:OuA)f  
1.首先动力计算时的dyntime每一步的计算时间步怎么放大或缩小？陈的书建议不好修改，但太小了，计算时间太长。 d3xmtG {i  
2.自己试算的模型最大不平衡力太大，不收敛。找不出问题。 =?!wXOg_  
3.动力边界条件不知设置对否？ ~@}Bi@*  
4.动力荷载的施加，面荷载可以仿照书上的，但点荷载不会施加，apply szz -1 hist wave range  去掉range 直接写上坐标位置，程序显示不对，加了后，也不知道荷载具体在哪个位置。不想ABAQUS那样明显易得。 % >=!p  
!r<pmr3f@7  
new Dh=9Gns9  
;第一步：静力计算阶段 YPxM<Gfa8  
config dyn 9 AJ(&qY(  
set dyn off <7~'; K  
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米一个网格 z4N*b"QF  
kGl~GOB a  
model mohr .[_L=_.  
prop bulk 5.33e7 shear 3.077e7 coh 5.41e4 fri 36.7 CB^U6ZS  
ini dens 2000
@{25xTt  
set grav 0 0 -9.81 B6={&7U2  
'dn]rV0(C  
fix x y z  range z -0.01 0.01 ;地基底部固定 Hl,W=2N  
fix x range x -0.1 0.1 m;,N)<~  
fix x range x 9.99 10.01      ;地基左右两侧x方向固定 ueUuJxq)  
fix y range y -0.01 0.01 }~L.qG  
fix y range y 9.99 10.01      ;地基前后纵向y方向固定 E 7{U|\  
hist unbal DA\2rLs  
solve j:v@p
zTD  
ini xdis 0 ydis 0 zdis 0                          ;将节点位移清零 +L;e^#>d  
ini xvel 0 yvel 0 zvel 0                          ;将节点速度清零 J\b^)  
;第二步：动力计算阶段 A6thXs2  
set dyn on .6Pw|xu`Pw  
def wave ln6d<; M5  
    wave=sin(2*pi*15*dytime) M!o##* *`  
end a^I\ /&aw'  
;定义一个正弦波，频率15HZ XuFYYx~ ^3  
free x y z range z -0.1 0.1                 ;去掉模型底部的静力条件             )P sY($ &  
apply nquiet squiet dquiet ran z -0.1 0.1   ;静态边界条件 jWA(C;W  
apply szz -1 hist wave range x 5 6 y 5 6 z 9.99 10.01  ;地基面施加一个动力荷载 .:F%_dS D  
apply ff                                    ;施加自由场边界条件 X9V*UXTc  
set dyn time = 0                            ;设置动力计算从0s开始 $|@ (  
hist  reset                                 ;清空已有的历史信息 HMNLa*CL'  
hist  unbal "]} bFO7C  
hist  dytime dl.p\t(1  
hist  wave , K~}\CR  
solve age 5 ~UP[A'9jJ  
save 11.sav

1#< '&Lr  
Re:恳请牛人进来看一下我的动力分析哪里出错了。谢谢