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

wjx541297

自己比着别人的代码写了一个，加载后怎么都不行 请各位给检查检查看看啥问题

1.       subroutine vdload (
   
2. C Read only (unmodifiable)variables -
   
3.      1 nblock, ndim, stepTime, totalTime,
   
4.      2 amplitude, curCoords, velocity, dirCos, jltyp, sname,
   
5. C Write only (modifiable) variable -
   
6.      1 value )
   
7. C
   
8.       include 'vaba_param.inc'
   
9. C
   
10.       dimension curCoords(nblock,ndim), velocity(nblock,ndim),
    
11.      1  dirCos(nblock,ndim,ndim), value(nblock)
    
12.       character*80 sname
    
13. C      
    
14.       parameter(vel=16.67,dlen=0.2,pressure=0.5d6)
    
15.       integer k
    
16.         disp=vel*steptime
    
17.           zc=disp-0
    
18.           zmax=zc
    
19.           zmin=zmax-dlen
    
20.         IF (SNAME(1:15).EQ.'ASSEMBLY_SURF-5') THEN
    
21.           do  k=1, nblock
    
22.           if(curCoords(k,1).lt.zmax.and.curCoords(k,1).ge.zmin) then
    
23.                value(k)=pressure
    
24.             else
    
25.                value(k)=0.0
    
26.            end if
    
27.         enddo
    
28.       ENDIF
    
29.           return
    
30.       end
    

复制代码

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

wjx541297

我的模型如图所示

wjx541297

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

smokyeyes

value(k)=pressure

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

wjx541297

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

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

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

飒沓自峥嵘

先下载学习了，谢谢分享。