Dave1231
Civil/Environmental
- Mar 2, 2011
- 9
I'm relatively new to ABAQUS and currently trying to model a wheel's movement over a beam element. I have succeeded in building a model using the AMPLITUDE function, now however, I want to build one using a FORTRAN subroutine. I have designed a DLOAD subroutine - enclosed below. I have attached the subroutine to ABAQUS and the job returns with ABAQUS error - Problem during compilation - C:\DLOAD.for. Does anyone have any suggestions?
I use ABAQUS 6.10-1 and Microsoft Visual Studio 2010 - Intel Visual Fortran.
I've modelled a simple 50m long beam element with a rectangular profile - height=0.153m and width=0.078m. The subroutine is attempting to model a 0N load over the beam element except from 2 uniform moving loads spaced at 1.8m apart. The moving loads are each the circumference of a wheel=2.388m. I have tried numerous alterations to my subroutine, all of which return the same error. Any help would be much appreciated. A step-by-step guide would be a treat. Regards.
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
C
C----------------------ROAD LENGTH
RDL=50.
C
C----------------------WHEEL LOADING Vehicle : 1 or 2. F - front R - Rear
WLD1F=52345.
WLD1R=52345.
C
C----------------------FRONT-REAR WHEEL SPACING
WLSP1=1.8
C
IROAD=1
C
IF(IROAD.EQ.1) THEN
ST1R=0.+TIME(1)*(RDL-2.388)
FI1R=2.388.+TIME(1)*(RDL-2.388)
C
ST1F=ST1R+WLSP1
FI1F=FI1R+WLSP1
ENDIF
C
F=0.
C
C--------------- VEHICLE 1 REAR wheel load
IF(COORDS(1).GE.ST1R.AND.COORDS(1).LE.FI1R)THEN
F=WLD1R
ST=ST1R
FI=FI1R
C
C---------------- VEHICLE 1 FRONT wheel load
ELSE IF(COORDS(1).GE.ST1F.AND.COORDS(1).LE.FI1F) THEN
C
C-------------FOR UNIFORM LOADING over the loaded region
F=WLD1F
ST=ST1F
FI=FI1F
C
ENDIF
C
RETURN
END
I use ABAQUS 6.10-1 and Microsoft Visual Studio 2010 - Intel Visual Fortran.
I've modelled a simple 50m long beam element with a rectangular profile - height=0.153m and width=0.078m. The subroutine is attempting to model a 0N load over the beam element except from 2 uniform moving loads spaced at 1.8m apart. The moving loads are each the circumference of a wheel=2.388m. I have tried numerous alterations to my subroutine, all of which return the same error. Any help would be much appreciated. A step-by-step guide would be a treat. Regards.
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
C
C----------------------ROAD LENGTH
RDL=50.
C
C----------------------WHEEL LOADING Vehicle : 1 or 2. F - front R - Rear
WLD1F=52345.
WLD1R=52345.
C
C----------------------FRONT-REAR WHEEL SPACING
WLSP1=1.8
C
IROAD=1
C
IF(IROAD.EQ.1) THEN
ST1R=0.+TIME(1)*(RDL-2.388)
FI1R=2.388.+TIME(1)*(RDL-2.388)
C
ST1F=ST1R+WLSP1
FI1F=FI1R+WLSP1
ENDIF
C
F=0.
C
C--------------- VEHICLE 1 REAR wheel load
IF(COORDS(1).GE.ST1R.AND.COORDS(1).LE.FI1R)THEN
F=WLD1R
ST=ST1R
FI=FI1R
C
C---------------- VEHICLE 1 FRONT wheel load
ELSE IF(COORDS(1).GE.ST1F.AND.COORDS(1).LE.FI1F) THEN
C
C-------------FOR UNIFORM LOADING over the loaded region
F=WLD1F
ST=ST1F
FI=FI1F
C
ENDIF
C
RETURN
END