wangguobo4220 发表于 2008-6-2 21:48:36

移动荷载作用于梁上的分析

在陈文化老师写的那本ansys参考书上看到关于车辆过桥的问题，我用FLAC3D写了简单程序，也想试试，可问题多多，现把命令流贴上，敬请大家讨论，指点迷津！
sel beam id=1 begin=( 1, 0, 0) end=( 9, 0, 0) nseg=9
sel beam id=1 prop emod=3e10 nu=0.20 density=2500 xcarea=1 xcj=0 xciy=0.083 xciz=0.083
sel nodefix x y z xr yr zr range id=1
sel nodefix x y z xr yr zr range id=2   ;为简单见，就直接全部约束了两端
config dynamic
set dyn on
;---------定义列车荷载-----------
def load
if dytime=0.02 then
    command
      sel node apply force=(0,0,0)
      sel node apply force=(0,0,-1e6) range id=1
   endcommand
end_if
if dytime=0.04 then
    command
      sel node apply force=(0,0,0)
      sel node apply force=(0,0,-1e6) range id=3
   endcommand
end_if
if dytime=0.06 then
    command
      sel node apply force=(0,0,0)
      sel node apply force=(0,0,-1e6) range id=4
   endcommand
end_if
if dytime=0.08 then
    command
      sel node apply force=(0,0,0)
      sel node apply force=(0,0,-1e6) range id=5
   endcommand
end_if
if dytime=0.1 then
    command
   sel node apply force=(0,0,0)
      sel node apply force=(0,0,-1e6) range id=6
   endcommand
end_if
if dytime=0.12 then
    command
      sel node apply force=(0,0,0)
      sel node apply force=(0,0,-1e6) range id=7
   endcommand
end_if
if dytime=0.14 then
    command
      sel node apply force=(0,0,0)
      sel node apply force=(0,0,-1e6) range id=8
   endcommand
end_if
if dytime=0.16 then
    command
      sel node apply force=(0,0,0)
      sel node apply force=(0,0,-1e6) range id=9
   endcommand
end_if
if dytime=0.18 then
    command
      sel node apply force=(0,0,0)
      sel node apply force=(0,0,-1e6) range id=10
   endcommand
end_if
if dytime=0.22 then
    command
   sel node apply force=(0,0,0)
      sel node apply force=(0,0,-1e6) range id=2
   endcommand
end_if
end
;load

;－－－选取的几个比较典型的节点进行记录－－－－－
hist reset
hist unbal               
hist dytime
hist sel node zvel id=5;好像不能监测加速度！
solve age 0.3
save dynamic.sav

命令流记事本也附上


忘记了把出现的问题附上： 最后监测的节点5的速度始终为零，于是就以为荷载没加上去，但在call这个记事本时，将load前面的分号去掉，窗口显示load的过程，本人又觉得荷载施加上去了。具体的我也不敢肯定了。

[ 本帖最后由 wangguobo4220 于 2008-6-3 07:58 编辑 ]

hoholiuwenyi 发表于 2008-6-2 22:43:29

觉得好复杂的问题！

wangguobo4220 发表于 2008-6-3 08:41:16

本人又尝试了下，可能是动力计算前没有进行静力计算的缘故，于是我先给了重力场进行静力计算，然后进行动力计算，这是node 5有速度，但很小，与ansys比起来简直不是一个数量级。期待继续讨论....修改后的命令流：sel beam id=1 begin=( 1, 0, 0) end=( 9, 0, 0) nseg=9sel beam id=1 prop emod=3e10 nu=0.20 density=2500 xcarea=0.1 xcj=0 xciy=0.000083 xciz=0.0083sel nodefix x y z xr yr zr range id=1sel nodefix x y z xr yr zr range id=2   ;为简单见，就直接全部约束了两端set grav 0 0 -10hist unbalsolvesave static.sav res static.savconfig dynamicsel node init xdisp 0sel node init ydisp 0sel node init zdisp 0sel node init xrdisp 0sel node init yrdisp 0sel node init zrdisp 0sel node init xvel 0sel node init yvel 0sel node init zvel 0sel node init xrvel 0sel node init yrvel 0sel node init zrvel 0set dyn on;---------定义列车荷载-----------def load if dytime=0.02 then    command      sel node apply force=(0,0,0)      sel node apply force=(0,0,-1e6) range id=1   endcommandend_if if dytime=0.04 then    command      sel node apply force=(0,0,0)      sel node apply force=(0,0,-1e6) range id=3   endcommandend_if if dytime=0.06 then    command      sel node apply force=(0,0,0)      sel node apply force=(0,0,-1e6) range id=4   endcommandend_if if dytime=0.08 then    command      sel node apply force=(0,0,0)      sel node apply force=(0,0,-1e6) range id=5   endcommandend_if if dytime=0.1 then    command   sel node apply force=(0,0,0)      sel node apply force=(0,0,-1e6) range id=6   endcommandend_if if dytime=0.12 then    command      sel node apply force=(0,0,0)      sel node apply force=(0,0,-1e6) range id=7   endcommandend_if if dytime=0.14 then    command      sel node apply force=(0,0,0)      sel node apply force=(0,0,-1e6) range id=8   endcommandend_if if dytime=0.16 then    command      sel node apply force=(0,0,0)      sel node apply force=(0,0,-1e6) range id=9   endcommandend_if if dytime=0.18 then    command      sel node apply force=(0,0,0)      sel node apply force=(0,0,-1e6) range id=10   endcommandend_if if dytime=0.22 then    command   sel node apply force=(0,0,0)      sel node apply force=(0,0,-1e6) range id=2   endcommandend_ifendload;－－－选取的几个比较典型的节点进行记录－－－－－hist resethist unbal                hist dytimehist sel node zvel id=5;好像不能监测加速度！solve age 0.3save dynamic.sav另外，用ansys计算的部分如下：/prep7ET,1,BEAM3R,1,0.1,0.0001/12,0.1, , , ,MPDATA,EX,1,,3e8MPDATA,PRXY,1,,0.2MPDATA,DENS,1,,2500 n,1,1   n,2,2   n,3,3   n,4,4   n,5,5   n,6,6   n,7,7   n,8,8   n,9,9   n,10,10 e,1,2   e,2,3   e,3,4   e,4,5   e,5,6   e,6,7   e,7,8   e,8,9   e,9,10 d,1,all,0d,10,all,0ANTYPE,4TRNOPT,FULL LUMPM,0 *do,i,1,10delt,0.02time,i*0.02fdele,i-1,allf, i, fy,-1e6allselsolveoutres,all,all*enddo附图分别为两种情况下的速度图

[ 本帖最后由 wangguobo4220 于 2008-6-3 08:42 编辑 ]

wangguobo4220 发表于 2008-6-3 08:45:50

怎么这么乱啊！还是用附件吧

wangguobo4220 发表于 2008-6-3 09:54:06

中间定义和在可以简化，
;---------定义列车荷载-----------
def load
   loop i (1,8)
    t1=0.02*i
    t2=0.02*(i-1)
    if dytime<=t1 then
      if dytime>t2 then
      command
          sel node apply force=(0,0,-1e6) range id=i+2
          load=force
      endcommand
      endif
    endif
   endloop
end

但最后print fish时显示如图：怎么荷载为0了？看来得到的节点速度曲线是重力作用引起的自振，我去掉这个荷载计算得到的曲线一致，验证了这种猜想，那就是说荷载还是没有加上去，这才是问题的关键。


[ 本帖最后由 wangguobo4220 于 2008-6-3 10:18 编辑 ]

wzn 发表于 2009-3-1 19:45:35

顶起来，请大家都讨论下，我也想做路面移动荷载的模拟，正不知怎么弄呢，请楼主和大家指教呀！谢谢！

gcchfut 发表于 2009-3-18 15:07:30

在学习中

ldhltw 发表于 2009-4-8 11:11:42

请问当速度增大时，动力响应为什么会减小呢

xiaozhou350321 发表于 2009-4-12 22:29:02

依wangguobo4220 兄的推论，那荷载没加上去的原因是什么？怎么解决啊？请高人指点！

s20070882 发表于 2009-4-13 09:15:52

我的同学做出来也是这个结果，高速时要小，他用的是移动质量模型，可能是高速时车与梁接触不完全造成的，车的质量来不及完全加到梁上

xiaozhou350321 发表于 2009-4-30 09:21:24

类似问题一直困扰我啊

oldfish2006 发表于 2009-5-4 22:23:34

好问题，顶起来，望高人相助

victory85 发表于 2009-7-16 09:53:32

你也是做路面的移动荷载模拟啊？相互讨论学习下吧，qq：286307176，（flac3d模拟）
6# wzn
页: [1]
查看完整版本: 移动荷载作用于梁上的分析