OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
sigeps85c_void.F File Reference
#include "implicit_f.inc"

Go to the source code of this file.

Functions/Subroutines

subroutine sigeps85c_void (nel, nuparam, nuvar, nfunc, ifunc, npf, npt, ipt, iflag, tf, time, timestep, uparam, rho0, area, eint, thkly, epspxx, epspyy, epspxy, epspyz, epspzx, depsxx, depsyy, depsxy, depsyz, depszx, epsxx, epsyy, epsxy, epsyz, epszx, sigoxx, sigoyy, sigoxy, sigoyz, sigozx, signxx, signyy, signxy, signyz, signzx, sigvxx, sigvyy, sigvxy, sigvyz, sigvzx, soundsp, viscmax, thk, pla, uvar, off, ngl, shf)

Function/Subroutine Documentation

◆ sigeps85c_void()

subroutine sigeps85c_void ( integer nel,
integer nuparam,
integer nuvar,
integer nfunc,
integer, dimension(nfunc) ifunc,
integer, dimension(*) npf,
integer npt,
integer ipt,
integer, dimension(*) iflag,
tf,
time,
timestep,
uparam,
rho0,
area,
eint,
thkly,
epspxx,
epspyy,
epspxy,
epspyz,
epspzx,
depsxx,
depsyy,
depsxy,
depsyz,
depszx,
epsxx,
epsyy,
epsxy,
epsyz,
epszx,
sigoxx,
sigoyy,
sigoxy,
sigoyz,
sigozx,
signxx,
signyy,
signxy,
signyz,
signzx,
sigvxx,
sigvyy,
sigvxy,
sigvyz,
sigvzx,
soundsp,
viscmax,
thk,
pla,
uvar,
off,
integer, dimension(nel) ngl,
shf )

Definition at line 30 of file sigeps85c_void.F.

43C-----------------------------------------------
44C I m p l i c i t T y p e s
45C-----------------------------------------------
46#include "implicit_f.inc"
47C-----------------------------------------------
48C G l o b a l P a r a m e t e r s
49C-----------------------------------------------
50C---------+---------+---+---+--------------------------------------------
51C VAR | SIZE |TYP| RW| DEFINITION
52C---------+---------+---+---+--------------------------------------------
53C NEL | 1 | I | R | SIZE OF THE ELEMENT GROUP NEL
54C NUPARAM | 1 | I | R | SIZE OF THE USER PARAMETER ARRAY
55C NUVAR | 1 | I | R | NUMBER OF USER ELEMENT VARIABLES
56C---------+---------+---+---+--------------------------------------------
57C NFUNC | 1 | I | R | NUMBER FUNCTION USED FOR THIS USER LAW
58C IFUNC | NFUNC | I | R | FUNCTION INDEX
59C NPF | * | I | R | FUNCTION ARRAY
60C NPT | 1 | I | R | NUMBER OF LAYERS OR INTEGRATION POINTS
61C IPT | 1 | I | R | LAYER OR INTEGRATION POINT NUMBER
62C IFLAG | * | I | R | GEOMETRICAL FLAGS
63C TF | * | F | R | FUNCTION ARRAY
64C NGL | NEL | I | R | ELEMENT NUMBER
65C SHF | NEL | F | R | SHEAR FACTOR
66C---------+---------+---+---+--------------------------------------------
67C TIME | 1 | F | R | CURRENT TIME
68C TIMESTEP| 1 | F | R | CURRENT TIME STEP
69C UPARAM | NUPARAM | F | R | USER MATERIAL PARAMETER ARRAY
70C RHO0 | NEL | F | R | INITIAL DENSITY
71C AREA | NEL | F | R | AREA
72C EINT | 2*NEL | F | R | INTERNAL ENERGY(MEMBRANE,BENDING)
73C THKLY | NEL | F | R | LAYER THICKNESS
74C EPSPXX | NEL | F | R | STRAIN RATE XX
75C EPSPYY | NEL | F | R | STRAIN RATE YY
76C ... | | | |
77C DEPSXX | NEL | F | R | STRAIN INCREMENT XX
78C DEPSYY | NEL | F | R | STRAIN INCREMENT YY
79C ... | | | |
80C EPSXX | NEL | F | R | STRAIN XX
81C EPSYY | NEL | F | R | STRAIN YY
82C ... | | | |
83C SIGOXX | NEL | F | R | OLD ELASTO PLASTIC STRESS XX
84C SIGOYY | NEL | F | R | OLD ELASTO PLASTIC STRESS YY
85C ... | | | |
86C---------+---------+---+---+--------------------------------------------
87C SIGNXX | NEL | F | W | NEW ELASTO PLASTIC STRESS XX
88C SIGNYY | NEL | F | W | NEW ELASTO PLASTIC STRESS YY
89C ... | | | |
90C SIGVXX | NEL | F | W | VISCOUS STRESS XX
91C SIGVYY | NEL | F | W | VISCOUS STRESS YY
92C ... | | | |
93C SOUNDSP | NEL | F | W | SOUND SPEED (NEEDED FOR TIME STEP)
94C VISCMAX | NEL | F | W | MAXIMUM DAMPING MODULUS(NEEDED FOR TIME STEP)
95C---------+---------+---+---+--------------------------------------------
96C THK | NEL | F |R/W| THICKNESS
97C PLA | NEL | F |R/W| PLASTIC STRAIN
98C UVAR |NEL*NUVAR| F |R/W| USER ELEMENT VARIABLE ARRAY
99C OFF | NEL | F |R/W| DELETED ELEMENT FLAG (=1. ON, =0. OFF)
100C---------+---------+---+---+--------------------------------------------
101C-----------------------------------------------
102C I N P U T A r g u m e n t s
103C-----------------------------------------------
104C
105 INTEGER NEL, NUPARAM, NUVAR, NPT, IPT,IFLAG(*),
106 . NGL(NEL)
107 my_real time,timestep,uparam(nuparam),
108 . area(nel),rho0(nel),eint(nel,2),
109 . thkly(nel),pla(nel),shf(nel),
110 . epspxx(nel),epspyy(nel),
111 . epspxy(nel),epspyz(nel),epspzx(nel),
112 . depsxx(nel),depsyy(nel),
113 . depsxy(nel),depsyz(nel),depszx(nel),
114 . epsxx(nel) ,epsyy(nel) ,
115 . epsxy(nel) ,epsyz(nel) ,epszx(nel) ,
116 . sigoxx(nel),sigoyy(nel),
117 . sigoxy(nel),sigoyz(nel),sigozx(nel)
118C-----------------------------------------------
119C O U T P U T A r g u m e n t s
120C-----------------------------------------------
121 my_real
122 . signxx(nel),signyy(nel),
123 . signxy(nel),signyz(nel),signzx(nel),
124 . sigvxx(nel),sigvyy(nel),
125 . sigvxy(nel),sigvyz(nel),sigvzx(nel),
126 . soundsp(nel),viscmax(nel)
127C-----------------------------------------------
128C I N P U T O U T P U T A r g u m e n t s
129C-----------------------------------------------
130 my_real uvar(nel,nuvar), off(nel),thk(nel)
131C-----------------------------------------------
132C VARIABLES FOR FUNCTION INTERPOLATION
133C-----------------------------------------------
134 INTEGER NPF(*), NFUNC, IFUNC(NFUNC)
135 my_real finter ,tf(*)
136 EXTERNAL finter
137C Y = FINTER(IFUNC(J),X,NPF,TF,DYDX)
138C Y : y = f(x)
139C X : x
140C DYDX : f'(x) = dy/dx
141C IFUNC(J): FUNCTION INDEX
142C J : FIRST(J=1), SECOND(J=2) .. FUNCTION USED FOR THIS LAW
143C NPF,TF : FUNCTION PARAMETER
144C-----------------------------------------------
145C L o c a l V a r i a b l e s
146C-----------------------------------------------
147 INTEGER I
148 my_real a1,a2,g
149C=======================================================================
150 a1 = uparam(4)
151 a2 = uparam(5)
152 g = uparam(6)
153C
154 DO i=1,nel
155C
156 signxx(i)=sigoxx(i)+a1*depsxx(i)+a2*depsyy(i)
157 signyy(i)=sigoyy(i)+a2*depsxx(i)+a1*depsyy(i)
158 signxy(i)=sigoxy(i)+g *depsxy(i)
159 signyz(i)=sigoyz(i)+g *depsyz(i)
160 signzx(i)=sigozx(i)+g *depszx(i)
161C
162 soundsp(i) = sqrt(a1/rho0(i))
163 viscmax(i) = zero
164C
165 ENDDO
166C-----------
167 RETURN
#define my_real
Definition cppsort.cpp:32
subroutine area(d1, x, x2, y, y2, eint, stif0)