Hallo an alle. Ich bin ein kompletter Anfänger was die ABAQUS subroutinen angeht. Ich habs versucht mir die einfachste Subroutine zu überlegen, die so aussieht:
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
Wenn ich den Job laufen lasse, kommt eine Fehlermeldung: "Problem during compilation". ABAQUS ist 100% richtig mit Fortran bzw. mit Visual Studio verknüpft da ich einen DLOAD subroutine geschrieben habe und hat problemlos gelaufen.
So sieht meine Input Datei:
*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
Danke euch im Voraus für ihre Unterstützung!
Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP