OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
roto.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!|| roto ../engine/source/airbag/roto.F
25!||--- called by ------------------------------------------------------
26!|| mulawglc ../engine/source/materials/mat_share/mulawglc.F
27!|| porfor5 ../engine/source/airbag/porfor5.F
28!|| porform5 ../engine/source/airbag/porfor5.F
29!||====================================================================
30 SUBROUTINE roto(JFT,JLT,SIG,DIR,NEL)
31C-----------------------------------------------
32C I m p l i c i t T y p e s
33C-----------------------------------------------
34#include "implicit_f.inc"
35C-----------------------------------------------
36C D u m m y A r g u m e n t s
37C-----------------------------------------------
38 INTEGER JFT,JLT,NEL
39 my_real
40 . sig(5,*), dir(nel,2)
41C-----------------------------------------------
42C L o c a l V a r i a b l e s
43C-----------------------------------------------
44 INTEGER I
45 my_real s1, s2, s3, s4, s5
46C-----------------------------------------------
47 DO i=jft,jlt
48 s1 = dir(i,1)*dir(i,1)*sig(1,i)
49 . +dir(i,2)*dir(i,2)*sig(2,i)
50 . + two*dir(i,1)*dir(i,2)*sig(3,i)
51 s2 = dir(i,2)*dir(i,2)*sig(1,i)
52 . +dir(i,1)*dir(i,1)*sig(2,i)
53 . -two*dir(i,2)*dir(i,1)*sig(3,i)
54 s3 =-dir(i,1)*dir(i,2)*sig(1,i)
55 . +dir(i,2)*dir(i,1)*sig(2,i)
56 . +(dir(i,1)*dir(i,1)-dir(i,2)*dir(i,2))*sig(3,i)
57 s4 =-dir(i,2)*sig(5,i)+dir(i,1)*sig(4,i)
58 s5 = dir(i,1)*sig(5,i)+dir(i,2)*sig(4,i)
59 sig(1,i)=s1
60 sig(2,i)=s2
61 sig(3,i)=s3
62 sig(4,i)=s4
63 sig(5,i)=s5
64 ENDDO
65C
66 RETURN
67 END
68!||====================================================================
69!|| rotov ../engine/source/airbag/roto.F
70!||--- called by ------------------------------------------------------
71!|| mat25_crasurv_c ../engine/source/materials/mat/mat025/mat25_crasurv_c.F90
72!|| mat25_tsaiwu_c ../engine/source/materials/mat/mat025/mat25_tsaiwu_c.f90
73!|| mulawc ../engine/source/materials/mat_share/mulawc.F90
74!|| usermat_shell ../engine/source/materials/mat_share/usermat_shell.F
75!||====================================================================
76 SUBROUTINE rotov(JFT,JLT,SIG,DIR,NEL)
77C-----------------------------------------------
78C I m p l i c i t T y p e s
79C-----------------------------------------------
80#include "implicit_f.inc"
81C-----------------------------------------------
82C G l o b a l P a r a m e t e r s
83C-----------------------------------------------
84#include "mvsiz_p.inc"
85C-----------------------------------------------
86C D u m m y A r g u m e n t s
87C-----------------------------------------------
88 INTEGER JFT,JLT,NEL
89 my_real
90 . sig(mvsiz,5), dir(nel,2)
91C-----------------------------------------------
92C L o c a l V a r i a b l e s
93C-----------------------------------------------
94 INTEGER I
95 my_real
96 . s1, s2, s3, s4, s5
97C-----------------------------------------------
98 DO i=jft,jlt
99 s1 = dir(i,1)*dir(i,1)*sig(i,1)
100 . +dir(i,2)*dir(i,2)*sig(i,2)
101 . + two*dir(i,1)*dir(i,2)*sig(i,3)
102 s2 = dir(i,2)*dir(i,2)*sig(i,1)
103 . +dir(i,1)*dir(i,1)*sig(i,2)
104 . -two*dir(i,2)*dir(i,1)*sig(i,3)
105 s3 =-dir(i,1)*dir(i,2)*sig(i,1)
106 . +dir(i,2)*dir(i,1)*sig(i,2)
107 . +(dir(i,1)*dir(i,1)-dir(i,2)*dir(i,2))*sig(i,3)
108 s4 =-dir(i,2)*sig(i,5)+dir(i,1)*sig(i,4)
109 s5 = dir(i,1)*sig(i,5)+dir(i,2)*sig(i,4)
110 sig(i,1)=s1
111 sig(i,2)=s2
112 sig(i,3)=s3
113 sig(i,4)=s4
114 sig(i,5)=s5
115 ENDDO
116C
117 RETURN
118 END
119!||====================================================================
120!|| roto_sig ../engine/source/airbag/roto.F
121!||--- called by ------------------------------------------------------
122!|| prony25c ../engine/source/materials/mat/mat025/prony25c.F
123!|| sigeps27c ../engine/source/materials/mat/mat027/sigeps27c.F
124!||====================================================================
125 SUBROUTINE roto_sig(JFT,JLT,SIG,DIR,NEL)
126C-----------------------------------------------
127C I m p l i c i t T y p e s
128C-----------------------------------------------
129#include "implicit_f.inc"
130C-----------------------------------------------
131C D u m m y A r g u m e n t s
132C-----------------------------------------------
133 INTEGER JFT,JLT,NEL
134 my_real
135 . sig(nel,5), dir(nel,2)
136C-----------------------------------------------
137C L o c a l V a r i a b l e s
138C-----------------------------------------------
139 INTEGER I
140 my_real
141 . s1, s2, s3, s4, s5
142C-----------------------------------------------
143 DO i=jft,jlt
144 s1 = dir(i,1)*dir(i,1)*sig(i,1)
145 . +dir(i,2)*dir(i,2)*sig(i,2)
146 .+ two*dir(i,1)*dir(i,2)*sig(i,3)
147 s2 = dir(i,2)*dir(i,2)*sig(i,1)
148 . +dir(i,1)*dir(i,1)*sig(i,2)
149 .- two*dir(i,2)*dir(i,1)*sig(i,3)
150 s3 =-dir(i,1)*dir(i,2)*sig(i,1)
151 . +dir(i,2)*dir(i,1)*sig(i,2)
152 . +(dir(i,1)*dir(i,1)-dir(i,2)*dir(i,2))*sig(i,3)
153 s4 =-dir(i,2)*sig(i,5)+dir(i,1)*sig(i,4)
154 s5 = dir(i,1)*sig(i,5)+dir(i,2)*sig(i,4)
155 sig(i,1)=s1
156 sig(i,2)=s2
157 sig(i,3)=s3
158 sig(i,4)=s4
159 sig(i,5)=s5
160 ENDDO
161C
162 RETURN
163 END
my_real
#define my_real
Definition cppsort.cpp:32
roto
subroutine roto(jft, jlt, sig, dir, nel)
Definition roto.F:31
roto_sig
subroutine roto_sig(jft, jlt, sig, dir, nel)
Definition roto.F:126
rotov
subroutine rotov(jft, jlt, sig, dir, nel)
Definition roto.F:77