dload子程序

abaqus ⽤户dload⼦程序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),1 dirCos(nblock,ndim,ndim), value(nblock)character*80 snameC

parameter(vel=5.0d0,dlen=0.1d0,pressure=500000.0d0)integer k

disp=vel*steptimezc=dispzmax=zc

zmin=zmax-dlenzmax1=zmax+0.5zmin1=zmin+0.5zmax2=zmax+1.5zmin2=zmin+1.5zmax3=zmax+2zmin3=zmin+2

IF (SNAME(1:15).EQ.'ASSEMBLY_SURF-1') THENdo k=1, nblock

if(curCoords(k,3).lt.zmax.and.curCoords(k,3).ge.zmin) thenvalue(k)=amplitude * pressure

else if(curCoords(k,3).lt.zmax1.and.curCoords(k,3).ge.zmin1) thenvalue(k)=amplitude * pressure

else if(curCoords(k,3).lt.zmax2.and.curCoords(k,3).ge.zmin2) thenvalue(k)=amplitude * pressure

else if(curCoords(k,3).lt.zmax3.and.curCoords(k,3).ge.zmin3) thenvalue(k)=amplitude * pressureelse

value(k)=0.0end ifenddoENDIFreturnend