OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
sigeps01gpinch.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!|| sigeps01gpinch ../engine/source/materials/mat/mat001/sigeps01gpinch.F
25!||--- called by ------------------------------------------------------
26!|| mulawglcpinch ../engine/source/elements/shell/coqueba/mulawglcpinch.F
27!||====================================================================
28 SUBROUTINE sigeps01gpinch(
29 1 JFT ,JLT ,G ,THK ,OFF ,
30 2 GS ,A1 ,A2 ,NU ,THK0 ,
31 3 NEL ,SSP ,RHO ,
32 4 DEPSXX ,DEPSYY ,DEPSZZ ,
33 5 DEPSXY ,DEPSYZ ,DEPSZX ,
34 6 DEPBXX ,DEPBYY ,DEPBXY ,
35 7 DEPPXZ ,DEPPYZ ,
36 8 SIGOXX ,SIGOYY ,SIGOZZ ,
37 9 SIGOXY ,SIGOYZ ,SIGOZX ,
38 A MOMOXX ,MOMOYY ,MOMOXY ,
39 B MOMOPXZ ,MOMOPYZ ,
40 C SIGNXX ,SIGNYY ,SIGNZZ ,
41 D SIGNXY ,SIGNYZ ,SIGNZX ,
42 E MOMNXX ,MOMNYY ,MOMNXY ,
43 F MOMNPXZ ,MOMNPYZ )
44C-----------------------------------------------
45C I m p l i c i t T y p e s
46C-----------------------------------------------
47#include "implicit_f.inc"
48C-----------------------------------------------
49C G l o b a l P a r a m e t e r s
50C-----------------------------------------------
51#include "mvsiz_p.inc"
52C-----------------------------------------------
53C I N P U T A R G U M E N T S
54C-----------------------------------------------
55 INTEGER JFT,JLT,NEL,ISMSTR
56C REAL
57 my_real
58 . OFF(*),GS(*),THK(*)
59 my_real
60 . NU(MVSIZ),G(MVSIZ),E(MVSIZ),A1(MVSIZ),A2(MVSIZ),
61 . DEPSXX(NEL),DEPSYY(NEL),DEPSZZ(NEL),
62 . DEPSXY(NEL),DEPSYZ(NEL),DEPSZX(NEL),
63 . DEPBXX(NEL),DEPBYY(NEL),DEPBXY(NEL),
64 . DEPPXZ(NEL),DEPPYZ(NEL),
65 . SIGOXX(NEL),SIGOYY(NEL),SIGOZZ(NEL),
66 . SIGOXY(NEL),SIGOYZ(NEL),SIGOZX(NEL),
67 . MOMOXX(NEL),MOMOYY(NEL),MOMOXY(NEL),
68 . MOMOPXZ(NEL),MOMOPYZ(NEL)
69C-----------------------------------------------
70C O U T P U T A R G U M E N T S
71C-----------------------------------------------
72 my_real
73 . signxx(nel),signyy(nel),signzz(nel),
74 . signxy(nel),signyz(nel),signzx(nel),
75 . momnxx(nel),momnyy(nel),momnxy(nel),
76 . momnpxz(nel),momnpyz(nel),
77 . ssp(nel),rho(nel)
78C-----------------------------------------------
79C L o c a l V a r i a b l e s
80C-----------------------------------------------
81 INTEGER I,MX,J
82C REAL
83 my_real
84 . PA1(MVSIZ),PA2(MVSIZ),PA3(MVSIZ),PA4(MVSIZ),PA5(MVSIZ),
85 . B1(MVSIZ),B2(MVSIZ),B3(MVSIZ),B4(MVSIZ),B5(MVSIZ),
86 . THK0(MVSIZ),THK08(MVSIZ)
87C-----------------------------------------------
88C
89 DO i=jft,jlt
90C [PM] reverse engineer E and ANU from A11 and A12
91C [PM] can't use A11/A12 since they correspond to shell (plane stress)
92 e(i) = a1(i)*(one-nu(i)**2)
93C
94C PA1 = C1111, C2222, C3333
95C PA2 = C1122, C1133, C2233
96C PA3 = C1212, C1313, C2323
97C PA4 = C1111*, C2222*
98C PA5 = C1122*
99 pa1(i) = e(i)*(one-nu(i))/(one+nu(i))/(one-two*nu(i))
100 pa2(i) = e(i)*nu(i)/(one+nu(i))/(one-two*nu(i))
101 pa3(i) = e(i)/two/(one+nu(i))
102 pa4(i) = (one+nu(i))*(one-two*nu(i))/(one-nu(i)**2)/(one-nu(i))*pa1(i)
103 pa5(i) = (one+nu(i))*(one-two*nu(i))/(one-nu(i)**2)*pa2(i)
104C
105 thk08(i) = thk0(i)*one_over_12
106 b1(i) = pa1(i)*thk08(i)
107 b2(i) = pa2(i)*thk08(i)
108 b3(i) = pa3(i)*thk08(i)
109 b4(i) = pa4(i)*thk08(i)
110 b5(i) = pa5(i)*thk08(i)
111 ENDDO
112C
113 DO i=jft,jlt
114C
115 signxx(i)=sigoxx(i)+pa1(i)*depsxx(i)+pa2(i)*depsyy(i)+pa2(i)*depszz(i)
116 signyy(i)=sigoyy(i)+pa2(i)*depsxx(i)+pa1(i)*depsyy(i)+pa2(i)*depszz(i)
117 signzz(i)=sigozz(i)+pa2(i)*depsxx(i)+pa2(i)*depsyy(i)+pa1(i)*depszz(i)
118C
119 signxy(i)=sigoxy(i)+g(i)*depsxy(i)
120 signyz(i)=sigoyz(i)+gs(i)*(depsyz(i)+zero*deppyz(i))
121 signzx(i)=sigozx(i)+gs(i)*(depszx(i)+zero*deppxz(i))
122C
123 momnxx(i)=momoxx(i)+b4(i)*depbxx(i)+b5(i)*depbyy(i)
124 momnyy(i)=momoyy(i)+b5(i)*depbxx(i)+b4(i)*depbyy(i)
125 momnxy(i)=momoxy(i)+b3(i)*depbxy(i)
126C
127 momnpxz(i)=momopxz(i)+b3(i)*deppxz(i)
128 momnpyz(i)=momopyz(i)+b3(i)*deppyz(i)
129C
130C Sound speed
131C SSP(I) = SQRT(E(I)/RHO(I)) ! historic
132C SSP(I) = SQRT(E(I)/RHO(I)/(ONE-ANU(I)*ANU(I))) ! for true shells
133 ssp(i) = sqrt(pa1(i)/rho(i)) ! like solids
134C
135 ENDDO
136 RETURN
137 END
sigeps01gpinch
subroutine sigeps01gpinch(jft, jlt, g, thk, off, gs, a1, a2, nu, thk0, nel, ssp, rho, depsxx, depsyy, depszz, depsxy, depsyz, depszx, depbxx, depbyy, depbxy, deppxz, deppyz, sigoxx, sigoyy, sigozz, sigoxy, sigoyz, sigozx, momoxx, momoyy, momoxy, momopxz, momopyz, signxx, signyy, signzz, signxy, signyz, signzx, momnxx, momnyy, momnxy, momnpxz, momnpyz)
Definition sigeps01gpinch.F:44