OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
m14gtf.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!|| m14gtf ../engine/source/materials/mat/mat014/m14gtf.F
25!||--- called by ------------------------------------------------------
26!|| m12law ../engine/source/materials/mat/mat012/m12law.F
27!|| m14law ../engine/source/materials/mat/mat014/m14law.F
28!||====================================================================
29 SUBROUTINE m14gtf(SIG,AX ,AY ,AZ ,BX ,BY,
30 2 BZ ,CX ,CY ,CZ ,D1 ,D2,
31 3 D3 ,D4 ,D5 ,D6 ,T1 ,T2,
32 4 T3 ,T4 ,T5 ,T6 ,E1 ,E2,
33 5 E3 ,E4 ,E5 ,E6 ,NEL)
34C-----------------------------------------------
35C I m p l i c i t T y p e s
36C-----------------------------------------------
37#include "implicit_f.inc"
38C-----------------------------------------------
39C C o m m o n B l o c k s
40C-----------------------------------------------
41C-----------------------------------------------
42C D u m m y A r g u m e n t s
43C-----------------------------------------------
44 INTEGER NEL
45C REAL
46 my_real
47 . SIG(NEL,6),
48 . ax(*),ay(*),az(*),bx(*),by(*),bz(*),cx(*),cy(*),cz(*)
49C REAL
50 my_real
51 . t1(*), t2(*), t3(*), t4(*), t5(*), t6(*), d1(*), d2(*),
52 . d3(*), d4(*), d5(*), d6(*), e1(*), e2(*), e3(*), e4(*), e5(*),
53 . e6(*)
54C-----------------------------------------------
55C L o c a l V a r i a b l e s
56C-----------------------------------------------
57 INTEGER I
58C-----------------------------------------------
59 DO 200 I=1,nel
60 e1(i)=ax(i)*ax(i)*d1(i)
61 . +ay(i)*ay(i)*d2(i)
62 . +az(i)*az(i)*d3(i)
63 . +ax(i)*ay(i)*d4(i)
64 . +ay(i)*az(i)*d5(i)
65 . +az(i)*ax(i)*d6(i)
66 e2(i)=bx(i)*bx(i)*d1(i)
67 . +by(i)*by(i)*d2(i)
68 . +bz(i)*bz(i)*d3(i)
69 . +bx(i)*by(i)*d4(i)
70 . +by(i)*bz(i)*d5(i)
71 . +bz(i)*bx(i)*d6(i)
72 e3(i)=cx(i)*cx(i)*d1(i)
73 . +cy(i)*cy(i)*d2(i)
74 . +cz(i)*cz(i)*d3(i)
75 . +cx(i)*cy(i)*d4(i)
76 . +cy(i)*cz(i)*d5(i)
77 . +cz(i)*cx(i)*d6(i)
78 e4(i)=ax(i)*bx(i)*d1(i)
79 . +ay(i)*by(i)*d2(i)
80 . +az(i)*bz(i)*d3(i)
81 . +(ax(i)*by(i)+ay(i)*bx(i))*d4(i)*half
82 . +(ay(i)*bz(i)+az(i)*by(i))*d5(i)*half
83 . +(az(i)*bx(i)+ax(i)*bz(i))*d6(i)*half
84 e5(i)=bx(i)*cx(i)*d1(i)
85 . +by(i)*cy(i)*d2(i)
86 . +bz(i)*cz(i)*d3(i)
87 . +(bx(i)*cy(i)+by(i)*cx(i))*d4(i)*half
88 . +(by(i)*cz(i)+bz(i)*cy(i))*d5(i)*half
89 . +(bz(i)*cx(i)+bx(i)*cz(i))*d6(i)*half
90 e6(i)=cx(i)*ax(i)*d1(i)
91 . +cy(i)*ay(i)*d2(i)
92 . +cz(i)*az(i)*d3(i)
93 . +(cx(i)*ay(i)+cy(i)*ax(i))*d4(i)*half
94 . +(cy(i)*az(i)+cz(i)*ay(i))*d5(i)*half
95 . +(cz(i)*ax(i)+cx(i)*az(i))*d6(i)*half
96 e4(i)=e4(i)*two
97 e5(i)=e5(i)*two
98 e6(i)=e6(i)*two
99 200 CONTINUE
100C
101 DO 400 i=1,nel
102 t1(i)=ax(i)*ax(i)*sig(i,1)
103 . +ay(i)*ay(i)*sig(i,2)
104 . +az(i)*az(i)*sig(i,3)
105 . +two*ax(i)*ay(i)*sig(i,4)
106 . +two*ay(i)*az(i)*sig(i,5)
107 . +two*az(i)*ax(i)*sig(i,6)
108 t2(i)=bx(i)*bx(i)*sig(i,1)
109 . +by(i)*by(i)*sig(i,2)
110 . +bz(i)*bz(i)*sig(i,3)
111 . +two*bx(i)*by(i)*sig(i,4)
112 . +two*by(i)*bz(i)*sig(i,5)
113 . +two*bz(i)*bx(i)*sig(i,6)
114 t3(i)=cx(i)*cx(i)*sig(i,1)
115 . +cy(i)*cy(i)*sig(i,2)
116 . +cz(i)*cz(i)*sig(i,3)
117 . +two*cx(i)*cy(i)*sig(i,4)
118 . +two*cy(i)*cz(i)*sig(i,5)
119 . +two*cz(i)*cx(i)*sig(i,6)
120 t4(i)=ax(i)*bx(i)*sig(i,1)
121 . +ay(i)*by(i)*sig(i,2)
122 . +az(i)*bz(i)*sig(i,3)
123 . +(ax(i)*by(i)+ay(i)*bx(i))*sig(i,4)
124 . +(ay(i)*bz(i)+az(i)*by(i))*sig(i,5)
125 . +(az(i)*bx(i)+ax(i)*bz(i))*sig(i,6)
126 t5(i)=bx(i)*cx(i)*sig(i,1)
127 . +by(i)*cy(i)*sig(i,2)
128 . +bz(i)*cz(i)*sig(i,3)
129 . +(bx(i)*cy(i)+by(i)*cx(i))*sig(i,4)
130 . +(by(i)*cz(i)+bz(i)*cy(i))*sig(i,5)
131 . +(bz(i)*cx(i)+bx(i)*cz(i))*sig(i,6)
132 t6(i)=cx(i)*ax(i)*sig(i,1)
133 . +cy(i)*ay(i)*sig(i,2)
134 . +cz(i)*az(i)*sig(i,3)
135 . +(cx(i)*ay(i)+cy(i)*ax(i))*sig(i,4)
136 . +(cy(i)*az(i)+cz(i)*ay(i))*sig(i,5)
137 . +(cz(i)*ax(i)+cx(i)*az(i))*sig(i,6)
138 400 CONTINUE
139C
140 RETURN
141 END
m14gtf
subroutine m14gtf(sig, ax, ay, az, bx, by, bz, cx, cy, cz, d1, d2, d3, d4, d5, d6, t1, t2, t3, t4, t5, t6, e1, e2, e3, e4, e5, e6, nel)
Definition m14gtf.F:34