wjx541297 发表于 2019-4-28 14:13:44

使用vdload对路面结构施加车辆荷载，运算错误，请指教。

自己比着别人的代码写了一个，加载后怎么都不行 请各位给检查检查看看啥问题      subroutine vdload (
C Read only (unmodifiable)variables -
   1 nblock, ndim, stepTime, totalTime,
   2 amplitude, curCoords, velocity, dirCos, jltyp, sname,
C Write only (modifiable) variable -
   1 value )
C
      include 'vaba_param.inc'
C
      dimension curCoords(nblock,ndim), velocity(nblock,ndim),
   1dirCos(nblock,ndim,ndim), value(nblock)
      character*80 sname
C      
      parameter(vel=16.67,dlen=0.2,pressure=0.5d6)
      integer k
      disp=vel*steptime
          zc=disp-0
          zmax=zc
          zmin=zmax-dlen
        IF (SNAME(1:15).EQ.'ASSEMBLY_SURF-5') THEN
          dok=1, nblock
          if(curCoords(k,1).lt.zmax.and.curCoords(k,1).ge.zmin) then
             value(k)=pressure
          else
             value(k)=0.0
           end if
      enddo
      ENDIF
        return
      end

在请教一个基础性的问题，代码是不是直到parameter那行之前都是vdload的格式，无论模型是啥都是一样的？

wjx541297 发表于 2019-4-28 14:20:38

我的模型如图所示

wjx541297 发表于 2019-4-28 14:34:43

我的CAE文件和vdload文件如下，懂得劳烦一下了

smokyeyes 发表于 2019-4-28 22:12:52

value(k)=pressure

这里，pressure的值是瞬间加上去的，导致网格扭曲，因此计算停止了。
将pressure减小一点就能完成计算了。或者将pressure与时间相关一点点加上去也可以。

wjx541297 发表于 2019-4-29 08:45:22

smokyeyes 发表于 2019-4-28 22:12
value(k)=pressure

这里，pressure的值是瞬间加上去的，导致网格扭曲，因此计算停止了。


我把荷载调成1.0还是不行
IF (SNAME(1:15).EQ.'ASSEMBLY_SURF-5') THEN
          dok=1, nblock
这各是啥意思，能解答一下么？

飒沓自峥嵘 发表于 2022-3-12 13:53:07


先下载学习了，谢谢分享。
页: [1]
查看完整版本: 使用vdload对路面结构施加车辆荷载，运算错误，请指教。