Hello,
I am trying to run a simulation using the dload subroutine. However, the simulation exited with an error and the error file does not exists.
I checked whether abaqus is properly linked with intel fortran, and it was properly linked. However, when i tried running verification, both standard and explicit using subroutines displayed error. Can someone help me solve this.
Details;
abaqus2022
Microsoft visual studio 2022
Intel oneAPI base tool kit and HPC tool kit 2024.1
Any help is highly appreciated.
Images and also DLOAD code is attached below.
DLOAD code
SUBROUTINE DLOAD(F,KSTEP,KINC,TIME,NOEL,NPT,LAYER,KSPT,
1 COORDS,JLTYP,SNAME)
C
INCLUDE 'ABA_PARAM.INC'
C
DIMENSION TIME(2), COORDS (3)
CHARACTER*80 SNAME
double precision x01, x02, x03, x04, y01, y02, y03, y04
double precision r0, td, x, y, r1, r2, r3, r4, P_max
double precision t, pt, pr1, pr2, pr3, pr4, pr_tot
r0=1.125
x01=3.0*r0
x02=5.0*r0
x03=3.0*r0
x04=5.0*r0
y01=5.0*r0
y02=5.0*r0
y03=3.0*r0
y04=3.0*r0
td=10.0**(-8.0)
x=coords(1)
y=coords(2)
t=time(1)
P_max=5170.0
r1=((x-x01)**2.0+(y-y01)**2.0)**0.5
r2=((x-x02)**2.0+(y-y02)**2.0)**0.5
r3=((x-x03)**2.0+(y-y03)**2.0)**0.5
r4=((x-x04)**2.0+(y-y04)**2.0)**0.5
pt=0.0
if (0.0<=t .and. t<=td)then
pt=(t/td)*P_max
end if
if (td<=t .and. t<=2.0*td)then
pt=(2.0-t/td)*P_max
end if
if (t>2.0*td)then
pt=0.0
end if
pr1=exp((-r1**2.0)/(2.0*r0**2.0))
pr2=exp((-r2**2.0)/(2.0*r0**2.0))
pr3=exp((-r3**2.0)/(2.0*r0**2.0))
pr4=exp((-r4**2.0)/(2.0*r0**2.0))
pr_tot=pr1+pr2+pr3+pr4
F=pt*pr_tot
RETURN
END
Images
