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

Go to the source code of this file.

Functions/Subroutines

subroutine sigeps110c (nel, ngl, nuparam, nuvar, npf, time, timestep, uparam, uvar, jthe, off, gs, rho, pla, dpla, epsp, soundsp, depsxx, depsyy, depsxy, depsyz, depszx, asrate, epspxx, epspyy, epspxy, epspyz, epspzx, sigoxx, sigoyy, sigoxy, sigoyz, sigozx, signxx, signyy, signxy, signyz, signzx, thkly, thk, sigy, et, tempel, temp, seq, tf, numtabl, itable, table, nvartmp, vartmp, siga, inloc, dplanl, loff, ioff_duct)

Function/Subroutine Documentation

◆ sigeps110c()

subroutine sigeps110c ( integer nel,
integer, dimension(nel), intent(in) ngl,
integer nuparam,
integer nuvar,
integer, dimension(*) npf,
time,
timestep,
intent(in) uparam,
intent(inout) uvar,
integer jthe,
intent(inout) off,
intent(in) gs,
intent(in) rho,
intent(inout) pla,
intent(inout) dpla,
intent(out) epsp,
intent(out) soundsp,
intent(in) depsxx,
intent(in) depsyy,
intent(in) depsxy,
intent(in) depsyz,
intent(in) depszx,
asrate,
intent(in) epspxx,
intent(in) epspyy,
intent(in) epspxy,
intent(in) epspyz,
intent(in) epspzx,
intent(in) sigoxx,
intent(in) sigoyy,
intent(in) sigoxy,
intent(in) sigoyz,
intent(in) sigozx,
intent(out) signxx,
intent(out) signyy,
intent(out) signxy,
intent(out) signyz,
intent(out) signzx,
intent(in) thkly,
intent(inout) thk,
intent(out) sigy,
intent(out) et,
intent(in) tempel,
intent(inout) temp,
intent(inout) seq,
tf,
integer numtabl,
integer, dimension(numtabl) itable,
type(ttable), dimension(ntable) table,
integer nvartmp,
integer, dimension(nel,nvartmp) vartmp,
intent(inout) siga,
integer inloc,
intent(in) dplanl,
intent(in) loff,
integer, dimension(nel), intent(inout) ioff_duct )

Definition at line 36 of file sigeps110c.F.

47 !=======================================================================
48 ! Modules
49 !=======================================================================
50 USE table_mod
51 USE interface_table_mod
52 !=======================================================================
53 ! Implicit types
54 !=======================================================================
55#include "implicit_f.inc"
56 !=======================================================================
57 ! Common
58 !=======================================================================
59#include "com04_c.inc"
60 !=======================================================================
61 ! Dummy arguments
62 !=======================================================================
63 INTEGER NEL,NUPARAM,NUVAR,JTHE,NUMTABL,ITABLE(NUMTABL),NVARTMP,NPF(*),INLOC
64 INTEGER,DIMENSION(NEL), INTENT(IN) :: NGL
65 INTEGER,DIMENSION(NEL), INTENT(INOUT) :: IOFF_DUCT
66 my_real
67 . time,timestep,asrate,tf(*)
68 INTEGER :: VARTMP(NEL,NVARTMP)
69 my_real,DIMENSION(NUPARAM), INTENT(IN) ::
70 . uparam
71 my_real,DIMENSION(NEL), INTENT(IN) ::
72 . rho,tempel,
73 . depsxx,depsyy,depsxy,depsyz,depszx,
74 . epspxx,epspyy,epspxy,epspyz,epspzx ,
75 . sigoxx,sigoyy,sigoxy,sigoyz,sigozx,
76 . gs,thkly,dplanl,loff
77c
78 my_real ,DIMENSION(NEL), INTENT(OUT) ::
79 . soundsp,sigy,et,epsp,
80 . signxx,signyy,signxy,signyz,signzx
81c
82 my_real ,DIMENSION(NEL), INTENT(INOUT) ::
83 . pla,dpla,off,thk,temp,seq
84 my_real ,DIMENSION(NEL,NUVAR), INTENT(INOUT) ::
85 . uvar
86 my_real ,DIMENSION(NEL,3) ,INTENT(INOUT) ::
87 . siga
88 TYPE(TTABLE), DIMENSION(NTABLE) :: TABLE
89C-----------------------------------------------
90C L o c a l V a r i a b l e s
91C-----------------------------------------------
92 INTEGER ICRIT,IRES
93C=======================================================================
94c UVAR(1)
95c UVAR(2)
96C DEPIJ = PLASTIC STRAIN TENSOR COMPONENT
97C DEPSIJ = TOTAL STRAIN TENSOR COMPONENT (EL+PL)
98C=======================================================================
99c
100 icrit = nint(uparam(27))! Vegter formulation flag
101 ! = 1 => Vegter classical
102 ! = 2 => Vegter standard
103 ! = 3 => Vegter 2017
104 ! = 4 => Vegter-Lite
105 ires = nint(uparam(29))! Plastic projection method
106 ! = 1 => Nice method
107 ! = 2 => Newton-iteration method
108c--------------------------
109 ! Choice of return mapping algorithm
110 SELECT CASE (ires)
111 ! Nice explicit return mapping
112 CASE(1)
113 ! Choice of Vegter formulation
114 SELECT CASE (icrit)
115 CASE(1,2,3)
116 CALL sigeps110c_nice(
117 1 nel ,ngl ,nuparam ,nuvar ,npf ,
118 2 time ,timestep,uparam ,uvar ,jthe ,off ,
119 3 gs ,rho ,pla ,dpla ,epsp ,soundsp ,
120 4 depsxx ,depsyy ,depsxy ,depsyz ,depszx ,asrate ,
121 5 epspxx ,epspyy ,epspxy ,epspyz ,epspzx ,
122 6 sigoxx ,sigoyy ,sigoxy ,sigoyz ,sigozx ,
123 7 signxx ,signyy ,signxy ,signyz ,signzx ,thkly ,
124 8 thk ,sigy ,et ,tempel ,temp ,seq ,
125 9 tf ,numtabl ,itable ,table ,nvartmp ,vartmp ,
126 a siga ,inloc ,dplanl ,loff )
127 CASE(4)
128 CALL sigeps110c_lite_nice(
129 1 nel ,ngl ,nuparam ,nuvar ,npf ,
130 2 time ,timestep,uparam ,uvar ,jthe ,off ,
131 3 gs ,rho ,pla ,dpla ,epsp ,soundsp ,
132 4 depsxx ,depsyy ,depsxy ,depsyz ,depszx ,asrate ,
133 5 epspxx ,epspyy ,epspxy ,epspyz ,epspzx ,
134 6 sigoxx ,sigoyy ,sigoxy ,sigoyz ,sigozx ,
135 7 signxx ,signyy ,signxy ,signyz ,signzx ,thkly ,
136 8 thk ,sigy ,et ,tempel ,temp ,seq ,
137 9 tf ,numtabl ,itable ,table ,nvartmp ,vartmp ,
138 a siga ,inloc ,dplanl ,loff )
139 END SELECT
140 ! Newton cutting plane return mapping
141 CASE(2)
142 ! Choice of Vegter formulation
143 SELECT CASE (icrit)
144 CASE(1,2,3)
145 CALL sigeps110c_newton(
146 1 nel ,ngl ,nuparam ,nuvar ,npf ,
147 2 time ,timestep,uparam ,uvar ,jthe ,off ,
148 3 gs ,rho ,pla ,dpla ,epsp ,soundsp ,
149 4 depsxx ,depsyy ,depsxy ,depsyz ,depszx ,asrate ,
150 5 epspxx ,epspyy ,epspxy ,epspyz ,epspzx ,
151 6 sigoxx ,sigoyy ,sigoxy ,sigoyz ,sigozx ,
152 7 signxx ,signyy ,signxy ,signyz ,signzx ,thkly ,
153 8 thk ,sigy ,et ,tempel ,temp ,seq ,
154 9 tf ,numtabl ,itable ,table ,nvartmp ,vartmp ,
155 a siga ,inloc ,dplanl ,loff )
156 CASE(4)
157 CALL sigeps110c_lite_newton(
158 1 nel ,ngl ,nuparam ,nuvar ,npf ,
159 2 time ,timestep,uparam ,uvar ,jthe ,off ,
160 3 gs ,rho ,pla ,dpla ,epsp ,soundsp ,
161 4 depsxx ,depsyy ,depsxy ,depsyz ,depszx ,asrate ,
162 5 epspxx ,epspyy ,epspxy ,epspyz ,epspzx ,
163 6 sigoxx ,sigoyy ,sigoxy ,sigoyz ,sigozx ,
164 7 signxx ,signyy ,signxy ,signyz ,signzx ,thkly ,
165 8 thk ,sigy ,et ,tempel ,temp ,seq ,
166 9 tf ,numtabl ,itable ,table ,nvartmp ,vartmp ,
167 a siga ,inloc ,dplanl ,loff )
168 END SELECT
169 END SELECT
170 ! Ductile failure activation
171 ioff_duct(1:nel) = 1
172c--------------------------
my_real
#define my_real
Definition cppsort.cpp:32
interface_table_mod
Definition table_mod.F:238
table_mod
Definition table_mod.F:112
sigeps110c_lite_newton
subroutine sigeps110c_lite_newton(nel, ngl, nuparam, nuvar, npf, time, timestep, uparam, uvar, jthe, off, gs, rho, pla, dpla, epsp, soundsp, depsxx, depsyy, depsxy, depsyz, depszx, asrate, epspxx, epspyy, epspxy, epspyz, epspzx, sigoxx, sigoyy, sigoxy, sigoyz, sigozx, signxx, signyy, signxy, signyz, signzx, thkly, thk, sigy, et, tempel, temp, seq, tf, numtabl, itable, table, nvartmp, vartmp, siga, inloc, dplanl, loff)
Definition sigeps110c_lite_newton.F:45
sigeps110c_lite_nice
subroutine sigeps110c_lite_nice(nel, ngl, nuparam, nuvar, npf, time, timestep, uparam, uvar, jthe, off, gs, rho, pla, dpla, epsp, soundsp, depsxx, depsyy, depsxy, depsyz, depszx, asrate, epspxx, epspyy, epspxy, epspyz, epspzx, sigoxx, sigoyy, sigoxy, sigoyz, sigozx, signxx, signyy, signxy, signyz, signzx, thkly, thk, sigy, et, tempel, temp, seq, tf, numtabl, itable, table, nvartmp, vartmp, siga, inloc, dplanl, loff)
Definition sigeps110c_lite_nice.F:45
sigeps110c_newton
subroutine sigeps110c_newton(nel, ngl, nuparam, nuvar, npf, time, timestep, uparam, uvar, jthe, off, gs, rho, pla, dpla, epsp, soundsp, depsxx, depsyy, depsxy, depsyz, depszx, asrate, epspxx, epspyy, epspxy, epspyz, epspzx, sigoxx, sigoyy, sigoxy, sigoyz, sigozx, signxx, signyy, signxy, signyz, signzx, thkly, thk, sigy, et, tempel, temp, seq, tf, numtabl, itable, table, nvartmp, vartmp, siga, inloc, dplanl, loff)
Definition sigeps110c_newton.F:45
sigeps110c_nice
subroutine sigeps110c_nice(nel, ngl, nuparam, nuvar, npf, time, timestep, uparam, uvar, jthe, off, gs, rho, pla, dpla, epsp, soundsp, depsxx, depsyy, depsxy, depsyz, depszx, asrate, epspxx, epspyy, epspxy, epspyz, epspzx, sigoxx, sigoyy, sigoxy, sigoyz, sigozx, signxx, signyy, signxy, signyz, signzx, thkly, thk, sigy, et, tempel, temp, seq, tf, numtabl, itable, table, nvartmp, vartmp, siga, inloc, dplanl, loff)
Definition sigeps110c_nice.F:45