OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
sigeps85c_void.F
Go to the documentation of this file.
1Copyright> OpenRadioss
2Copyright> Copyright (C) 1986-2025 Altair Engineering Inc.
3Copyright>
4Copyright> This program is free software: you can redistribute it and/or modify
5Copyright> it under the terms of the GNU Affero General Public License as published by
6Copyright> the Free Software Foundation, either version 3 of the License, or
7Copyright> (at your option) any later version.
8Copyright>
9Copyright> This program is distributed in the hope that it will be useful,
10Copyright> but WITHOUT ANY WARRANTY; without even the implied warranty of
11Copyright> MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
12Copyright> GNU Affero General Public License for more details.
13Copyright>
14Copyright> You should have received a copy of the GNU Affero General Public License
15Copyright> along with this program. If not, see <https://www.gnu.org/licenses/>.
16Copyright>
17Copyright>
18Copyright> Commercial Alternative: Altair Radioss Software
19Copyright>
20Copyright> As an alternative to this open-source version, Altair also offers Altair Radioss
21Copyright> software under a commercial license. Contact Altair to discuss further if the
22Copyright> commercial version may interest you: https://www.altair.com/radioss/.
23!||====================================================================
24!|| sigeps85c_void ../engine/source/materials/mat/mat085/sigeps85c_void.F
25!||--- called by ------------------------------------------------------
26!|| mulawc ../engine/source/materials/mat_share/mulawc.F90
27!||--- calls -----------------------------------------------------
28!|| finter ../engine/source/tools/curve/finter.F
29!||====================================================================
30 SUBROUTINE sigeps85c_void(
31 1 NEL ,NUPARAM,NUVAR ,NFUNC ,IFUNC ,
32 2 NPF ,NPT ,IPT ,IFLAG ,
33 2 TF ,TIME ,TIMESTEP,UPARAM ,RHO0 ,
34 3 AREA ,EINT ,THKLY ,
35 4 EPSPXX ,EPSPYY ,EPSPXY ,EPSPYZ ,EPSPZX ,
36 5 DEPSXX ,DEPSYY ,DEPSXY ,DEPSYZ ,DEPSZX ,
37 6 EPSXX ,EPSYY ,EPSXY ,EPSYZ ,EPSZX ,
38 7 SIGOXX ,SIGOYY ,SIGOXY ,SIGOYZ ,SIGOZX ,
39 8 SIGNXX ,SIGNYY ,SIGNXY ,SIGNYZ ,SIGNZX ,
40 9 SIGVXX ,SIGVYY ,SIGVXY ,SIGVYZ ,SIGVZX ,
41 A SOUNDSP,VISCMAX,THK ,PLA ,UVAR ,
42 B OFF ,NGL ,SHF)
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
168 END
sigeps85c_void
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)
Definition sigeps85c_void.F:43