Search
Close this search box.
Notifications
Clear all

DFLUX Subroutine

1 Posts
1 Users
0 Reactions
714 Views
Posts: 1
Topic starter
(@kill25874)
New Member
Joined: 10 years ago

Hello guys. I am a complete newbie regarding the Abaqus subroutines. Ive tried to write the simplest DFLUX subroutine that I could have think of:

SUBROUTINE DFLUX(FLUX,SOL,KSTEP,KINC,TIME,NOEL,NPT,COORDS,

1 JLTYP,TEMP,PRESS,SNAME)

C

INCLUDE ABA_PARAM.INC

C

DIMENSION FLUX(2), TIME(2), COORDS(3)

CHARACTER*80 SNAME

C

x=COORDS(1)

y=COORDS(2)

z=COORDS(3)

FLUX(2)=0.0

FLUX(1)=0.0

T=TIME(2)

Q=200.0

C

IF ((JLTYP .EQ. 1.0) .AND. (T .LE. 1)) THEN

IF (x .LE. 5.0 .AND. y .LE. 5.0) THEN

IF (x .GE. -5.0 .AND. y .GE. -5.0) THEN

FLUX(2)=Q

ELSE FLUX(2)=0.0

ENDIF

ELSE FLUX(2)=0.0

ENDIF

ELSE FLUX(2)=0.0

ENDIF

RETURN

END

The problem is that when I want to submit the job an error occurs: Problem during compilation. The ABAQUS is fully and correctly coupled with Fortran and VS cause Ive tested a DLOAD subroutine and it worked just fine.

This is how my my INPUT file looks like:

*Heading

** Job name: dflux_V2 Model name: Model-1

** Generated by: Abaqus/CAE 6.13-2

*Preprint, echo=NO, model=NO, history=NO, contact=NO

**

** PARTS

**

*Part, name=Part-1

*End Part

**

**

** ASSEMBLY

**

*Assembly, name=Assembly

**

*Instance, name=Part-1-1, part=Part-1

*Node

1, 50., 50.

2, 25., 50.

3, 0., 50.

4, -25., 50.

5, -50., 50.

6, 50., 25.

7, 25., 25.

8, 0., 25.

9, -25., 25.

10, -50., 25.

11, 50., 0.

12, 25., 0.

13, 0., 0.

14, -25., 0.

15, -50., 0.

16, 50., -25.

17, 25., -25.

18, 0., -25.

19, -25., -25.

20, -50., -25.

21, 50., -50.

22, 25., -50.

23, 0., -50.

24, -25., -50.

25, -50., -50.

*Element, type=CPE4T

1, 1, 2, 7, 6

2, 2, 3, 8, 7

3, 3, 4, 9, 8

4, 4, 5, 10, 9

5, 6, 7, 12, 11

6, 7, 8, 13, 12

7, 8, 9, 14, 13

8, 9, 10, 15, 14

9, 11, 12, 17, 16

10, 12, 13, 18, 17

11, 13, 14, 19, 18

12, 14, 15, 20, 19

13, 16, 17, 22, 21

14, 17, 18, 23, 22

15, 18, 19, 24, 23

16, 19, 20, 25, 24

*Nset, nset=_PickedSet2, internal, generate

1, 25, 1

*Elset, elset=_PickedSet2, internal, generate

1, 16, 1

** Section: Section-1

*Solid Section, elset=_PickedSet2, material=Material-Stahl

,

*End Instance

**

*Nset, nset=_PickedSet4, internal, instance=Part-1-1

1, 2, 3, 4, 5, 21, 22, 23, 24, 25

*Elset, elset=_PickedSet4, internal, instance=Part-1-1

1, 2, 3, 4, 13, 14, 15, 16

*Nset, nset=_PickedSet5, internal, instance=Part-1-1, generate

1, 25, 1

*Elset, elset=_PickedSet5, internal, instance=Part-1-1, generate

1, 16, 1

*Nset, nset=_PickedSet10, internal, instance=Part-1-1, generate

1, 25, 1

*Elset, elset=_PickedSet10, internal, instance=Part-1-1, generate

1, 16, 1

*Nset, nset=_PickedSet12, internal, instance=Part-1-1, generate

1, 25, 1

*Elset, elset=_PickedSet12, internal, instance=Part-1-1, generate

1, 16, 1

*End Assembly

**

** MATERIALS

**

*Material, name=Material-Stahl

*Conductivity

25.3,

*Density

7.78e-06,

*Elastic

210000., 0.287

*Expansion

1.1e-05,

*Specific Heat

4.5e+08,

**

** BOUNDARY CONDITIONS

**

** Name: BC-1 Type: Symmetry/Antisymmetry/Encastre

*Boundary

_PickedSet4, ENCASTRE

**

** PREDEFINED FIELDS

**

** Name: Predefined Field-1 Type: Temperature

*Initial Conditions, type=TEMPERATURE

_PickedSet10, 20.

** ----------------------------------------------------------------

**

** STEP: Step-1

**

*Step, name=Step-1, nlgeom=YES

*Coupled Temperature-displacement, creep=none

1., 1.,

**

** LOADS

**

** Name: Load-1 Type: Body heat flux

*Dflux

_PickedSet5, BFNU, 1.

**

** OUTPUT REQUESTS

**

*Restart, write, frequency=0

**

** FIELD OUTPUT: F-Output-1

**

*Output, field, variable=PRESELECT

**

** HISTORY OUTPUT: H-Output-1

**

*Output, history, variable=PRESELECT

*End Step

Thanks for your suggestions and your support!

Topic Tags
Share: