- 积分
- 0
- 注册时间
- 2007-5-13
- 仿真币
-
- 最后登录
- 1970-1-1
|
本帖最后由 cloud328 于 2010-5-11 12:40 编辑
运行结果出现:Timestep rejected by module
自我感觉与网格划分和荷载施加无多大关系
请大伙帮忙分析一下原因,多谢了:)
;-------------------------------
问题已经解决,谢谢!因为在动力模式下,结构单元必须要考虑密度!
以下命令流正确!
;---------------------------------
config dyn
;
def iniwave
per=0.01
end
def wave
wave=0.5*(1.0-cos(2*pi*dytime/per))
end
iniwave
;
gen zone brick p0 0 0 0 p1 10 0 0 p2 0 1 0 p3 0 0 1 size 20 8 8
;
model elastic
prop bulk 66667 shear 40000
ini dens 0.0025
;
sel pile id=1 begin 4,0.5,0.5 end 6,0.5,0.5 nseg=7
sel pile id=1 begin 6,0.5,0.5 end 10,0.5,0.5 nseg=1
sel pile id=1 prop dens=1100 Emod=7.17e10 Nu=0.33 pmoment 3.04e3 XCArea=28.26e-6 &
XCJ=61.27e-10 XCIy=6.36e-11 XCIz=6.36e-11 Per=1.88e-2 &
CS_sK=1.48e6 CS_sCoh=8.0e8 CS_sFric=25.5 &
CS_nK=1.82e6 CS_nCoh=5.70e5 CS_nFric= 0.0 CS_nGap=on slide=on range cid 1,7
sel pile id=1 prop dens=1100 Emod=7.17e10 Nu=0.33 pmoment 3.04e3 XCArea=28.26e-6 &
XCJ=61.27e-10 XCIy=6.36e-11 XCIz=6.36e-11Per=1.88e-2 &
CS_sK=1 CS_sCoh=1 CS_sFric=0.0 &
CS_nK=1 CS_nCoh=1 CS_nFric=0.0 CS_nGap=on slide=on range cid 8
sel delete link id=9
sel link id=10 9 target zone
sel link attach xdir=rigid xrdir=rigid ydir=rigid yrdir=rigid &
zdir=rigid zrdir=rigid
sel link constit nydeform 6 area=1.0 k=2.84e12 &
ycomp=4.17e3 ytens=4.17e3 range id=10
sel node apply force (-100,0,0) range id=9
set grav 0 0 -10
fix x range x -0.01 0.01
fix x range x 0.99 1.01
fix y range y -0.01 0.01
fix y range y 0.99 1.01
fix z range z -0.01 0.01
;
set dyn off
;
step 100
;solve ratio 1e-5
;save ini.sav
;
set dyn on
;
free x y z range z -0.01 0.01
apply nquiet squiet dquiet range z -0.01 0.01
apply dstress 1.0 hist wave range z -0.01 0.01
;
set dyn time=0
;
solve age=0.01
save pile.sav
[ 本帖最后由 whutzll 于 2008-5-5 11:35 编辑 ] |
|