动力分析的例子

huangdaxian

;-----dynamic-----
set dyn on
set fluid off
set large
set dyn multi on
ini xvel 0 yvel 0 zvel 0
ini xdisp 0 ydisp 0 zdisp 0
ini state 0
model null range group pile
;定义钢护筒
model elas range group hutong
prop bulk 1.75e11 shear 0.808e11 range group hutong
ini dens 2900 range group hutong

interface 1 face range cylin end1(0,0,0) end2(0,6.1,0) rad 0.81 &
                      cylin end1(0,0,0) end2(0,6.1,0) rad 0.79

interface 1 prop kn 1e8 ks 1e8 fric 20 coh 30000

;定义受力面
range name mian group hutong y -0.1 0.1
free x y z range mian   
;定义波形            
def setup
  omega=2.0*pi*freq
  pulse=1.0/freq
end
def wave
   if dytime>pulse
      wave=0.0
   else
      wave=5*(1.0-cos(omega*dytime))
   endif
end
set freq=2.0
apply syy=2e6 hist=waverange mian ；请问我这动荷载加上了吗？我自己不敢肯定，因为我怎么改变syy值其加速度都不边且很小？
apply nquiet dquiet squiet range y 14.9 15.1
apply nquiet dquiet squiet range bianjie
fix pp 0 range group hutong
solve age 5.0

huangdaxian

:) 顶

huangdaxian

尊敬的陈育民博士帮忙看一下吧 :handshake

mymeiy

没有阻尼

yanqing63

定义波形             6 D$ Z# s/ h) ]: t; k- I
def setup- ~/ A  ]( I1 i1 f
  omega=2.0*pi*freq
* w* c" v9 {, X  [7 L  pulse=1.0/freq! p- Z9 L  W# t! n- x
end' T! y" s" l6 j3 O
def wave7 Q) h$ n3 a9 F1 \% R% T
   if dytime>pulse7 g$ X' v1 s- V! Q5 S
      wave=0.0- T; }+ b; K# F; B, B/ _, o7 s0 Y: }
   else
. Y( B4 V/ L# M1 V+ ?      wave=5*(1.0-cos(omega*dytime))

这段是输入的荷载嘛？  有些是好像是  0   0   0
                                                           1    -   -  类似的   这两种方法是怎么回事呢？ 请教

happysky_24

貌似又在静力边界条件上施加了加速度时程，错啦

只想风流不想累

好深奥没有看懂唉 怎么学啊 这么难

zhangxdo

看看似乎不懂，加紧学习

furnace110

好深奥没有看懂唉 怎么学啊 这么难

半葫酒

最近也在做这方面的，我QQ442434339。要请教楼主

yanqisha

静态边界上不能用加速度时程吧，需要用应力时程