OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
szsvm_or.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!|| szsvm_or ../engine/source/elements/solid/solidez/szsvm_or.F
25!||--- called by ------------------------------------------------------
26!|| szhour3_or ../engine/source/elements/solid/solidez/szhour3_or.F
27!||====================================================================
28 SUBROUTINE szsvm_or(
29 1 JR0, JS0, JT0, CC,
30 2 CG, G33, FHOUR, SIGY,
31 3 SIG0, NU, SVM1, SVM2,
32 4 NEL, IINT)
33C-----------------------------------------------
34C I m p l i c i t T y p e s
35C-----------------------------------------------
36#include "implicit_f.inc"
37C-----------------------------------------------
38C G l o b a l P a r a m e t e r s
39C-----------------------------------------------
40#include "mvsiz_p.inc"
41C-----------------------------------------------
42C C o m m o n B l o c k s
43C-----------------------------------------------
44C-----------------------------------------------
45C D u m m y A r g u m e n t s
46C-----------------------------------------------
47 INTEGER, INTENT(IN) :: IINT
48 INTEGER NEL
49 my_real
50 . FHOUR(NEL,3,4),JR0(*),JS0(*),JT0(*) ,
51 . sigy(*) ,sig0(nel,6),svm1(*),svm2(*),nu,
52 . cc(mvsiz,3,3),cg(mvsiz,3,3),g33(mvsiz,3,3)
53C-----------------------------------------------
54C L o c a l V a r i a b l e s
55C-----------------------------------------------
56 INTEGER I, J,K,IKT
57 my_real
58 . S1,S2,S3,SVM0,SR1,SR2,SR3,SR4,
59 . SS1,SS2,SS3,SS4,ST1,ST2,ST3,ST4,SVMR,SVMS,SVMT,
60 . svm,svmrst,coef,coef1, jr_1,js_1,jt_1,nu1,nu2,t1,t2,
61 . rs,st,rt,vt1(4),vt2(4),max1,min1,max2,min2,max0,min0,
62 . sr5,sr6,ss5,ss6, st5,st6
63C-----------------------------------------------
64C coef= (3*8*0.75)^2 coef1 =1/0.75
65 coef= threehundred24
66 coef1= onep333
67 ikt = iint
68 IF (ikt ==2) THEN
69 coef1= coef/(ten*ikt)
70 coef= max(coef1,em01)
71 ELSEIF (ikt > 2) THEN
72 coef1= coef/(ten*ikt)
73 coef= max(coef1,em01)
74 ENDIF
75 nu1 =one/(one - nu)
76 nu2 =nu*nu1
77 DO i=1,nel
78 svm1(i) = zero
79 svm2(i) = zero
80 IF (sigy(i)<zep9ep30) THEN
81 s1 =sig0(i,1)-sig0(i,2)
82 s2 =sig0(i,2)-sig0(i,3)
83 s3 =sig0(i,1)-sig0(i,3)
84 svm0 = (s1*s1+s2*s2+s3*s3)*half + three*(sig0(i,4)*sig0(i,4)
85 . +sig0(i,5)*sig0(i,5)+sig0(i,6)*sig0(i,6))
86 jr_1 = one/max(em20,jr0(i))
87 js_1 = one/max(em20,js0(i))
88 jt_1 = one/max(em20,jt0(i))
89 j=2
90 k=1
91 t1 = nu2*(cc(i,j,2)-cc(i,k,2))
92 t2 = jr0(i)*jt_1*fhour(i,1,1)+jt0(i)*jr_1*fhour(i,3,3)
93 sr1 =(cc(i,j,1)-cc(i,k,1)-t1)*fhour(i,1,3)+
94 . (cc(i,j,3)-cc(i,k,3)-t1)*fhour(i,3,1)+
95 . (cg(i,j,2)-cg(i,k,2))*t2
96 j=3
97 k=2
98 t1 = nu2*(cc(i,j,2)-cc(i,k,2))
99 sr2 =(cc(i,j,1)-cc(i,k,1)-t1)*fhour(i,1,3)+
100 . (cc(i,j,3)-cc(i,k,3)-t1)*fhour(i,3,1)+
101 . (cg(i,j,2)-cg(i,k,2))*t2
102 j=3
103 k=1
104 t1 = nu2*(cc(i,j,2)-cc(i,k,2))
105 sr3 =(cc(i,j,1)-cc(i,k,1)-t1)*fhour(i,1,3)+
106 . (cc(i,j,3)-cc(i,k,3)-t1)*fhour(i,3,1)+
107 . (cg(i,j,2)-cg(i,k,2))*t2
108 j=1
109 t1 = nu2*cg(i,2,j)
110 sr4 =(cg(i,1,j)-t1)*fhour(i,1,3)+(cg(i,3,j)-t1)*fhour(i,3,1)+
111 . g33(i,j,2)*t2
112 j=2
113 t1 = nu2*cg(i,2,j)
114 sr5 =(cg(i,1,j)-t1)*fhour(i,1,3)+(cg(i,3,j)-t1)*fhour(i,3,1)+
115 . g33(i,j,2)*t2
116 j=3
117 t1 = nu2*cg(i,2,j)
118 sr6 =(cg(i,1,j)-t1)*fhour(i,1,3)+(cg(i,3,j)-t1)*fhour(i,3,1)+
119 . g33(i,j,2)*t2
120 svmr =half*(sr1*sr1+sr2*sr2+sr3*sr3)+three*(sr4*sr4+sr5*sr5+sr6*sr6)
121 j=2
122 k=1
123 t1 = nu2*(cc(i,j,3)-cc(i,k,3))
124 t2 = jr0(i)*js_1*fhour(i,1,1)+js0(i)*jr_1*fhour(i,2,2)
125 ss1 =(cc(i,j,1)-cc(i,k,1)-t1)*fhour(i,1,2)+
126 . (cc(i,j,2)-cc(i,k,2)-t1)*fhour(i,2,1)+
127 . (cg(i,j,1)-cg(i,k,1))*t2
128 j=3
129 k=2
130 t1 = nu2*(cc(i,j,3)-cc(i,k,3))
131 ss2 =(cc(i,j,1)-cc(i,k,1)-t1)*fhour(i,1,2)+
132 . (cc(i,j,2)-cc(i,k,2)-t1)*fhour(i,2,1)+
133 . (cg(i,j,1)-cg(i,k,1))*t2
134 j=3
135 k=1
136 t1 = nu2*(cc(i,j,3)-cc(i,k,3))
137 ss3 =(cc(i,j,1)-cc(i,k,1)-t1)*fhour(i,1,2)+
138 . (cc(i,j,2)-cc(i,k,2)-t1)*fhour(i,2,1)+
139 . (cg(i,j,1)-cg(i,k,1))*t2
140 j=1
141 t1 = nu2*cg(i,3,j)
142 ss4 =(cg(i,1,j)-t1)*fhour(i,1,2)+(cg(i,2,j)-t1)*fhour(i,2,1)+
143 . g33(i,j,1)*t2
144 j=2
145 t1 = nu2*cg(i,3,j)
146 ss5 =(cg(i,1,j)-t1)*fhour(i,1,2)+(cg(i,2,j)-t1)*fhour(i,2,1)+
147 . g33(i,j,1)*t2
148 j=3
149 t1 = nu2*cg(i,3,j)
150 ss6 =(cg(i,1,j)-t1)*fhour(i,1,2)+(cg(i,2,j)-t1)*fhour(i,2,1)+
151 . g33(i,j,1)*t2
152 svms =half*(ss1*ss1+ss2*ss2+ss3*ss3)+three*(ss4*ss4+ss5*ss5+ss6*ss6)
153 j=2
154 k=1
155 t1 = nu2*(cc(i,j,1)-cc(i,k,1))
156 t2 = js0(i)*jt_1*fhour(i,2,2)+jt0(i)*js_1*fhour(i,3,3)
157 st1 =(cc(i,j,2)-cc(i,k,2)-t1)*fhour(i,2,3)+
158 . (cc(i,j,3)-cc(i,k,3)-t1)*fhour(i,3,2)+
159 . (cg(i,j,3)-cg(i,k,3))*t2
160 j=3
161 k=2
162 t1 = nu2*(cc(i,j,1)-cc(i,k,1))
163 st2 =(cc(i,j,2)-cc(i,k,2)-t1)*fhour(i,2,3)+
164 . (cc(i,j,3)-cc(i,k,3)-t1)*fhour(i,3,2)+
165 . (cg(i,j,3)-cg(i,k,3))*t2
166 j=3
167 k=1
168 t1 = nu2*(cc(i,j,1)-cc(i,k,1))
169 st3 =(cc(i,j,2)-cc(i,k,2)-t1)*fhour(i,2,3)+
170 . (cc(i,j,3)-cc(i,k,3)-t1)*fhour(i,3,2)+
171 . (cg(i,j,3)-cg(i,k,3))*t2
172 j=1
173 t1 = nu2*cg(i,1,j)
174 st4 =(cg(i,2,j)-t1)*fhour(i,2,3)+(cg(i,3,j)-t1)*fhour(i,3,2)+
175 . g33(i,j,3)*t2
176 j=2
177 t1 = nu2*cg(i,1,j)
178 st5 =(cg(i,2,j)-t1)*fhour(i,2,3)+(cg(i,3,j)-t1)*fhour(i,3,2)+
179 . g33(i,j,3)*t2
180 j=3
181 t1 = nu2*cg(i,1,j)
182 st6 =(cg(i,2,j)-t1)*fhour(i,2,3)+(cg(i,3,j)-t1)*fhour(i,3,2)+
183 . g33(i,j,3)*t2
184 svmt =half*(st1*st1+st2*st2+st3*st3)+three*(st4*st4+st5*st5+st6*st6)
185C
186 rs =half*sr1*ss1
187 rt =half*sr1*st1
188 st =half*ss1*st1
189 vt2(1) = rs+rt+st
190 vt2(2) = rs-rt-st
191 vt2(3) = -rs+rt-st
192 vt2(4) = -rs-rt+st
193 max0 =max(vt2(1),vt2(2),vt2(3),vt2(4))
194 min0 =min(vt2(1),vt2(2),vt2(3),vt2(4))
195 rs =half*sr2*ss2
196 rt =half*sr2*st2
197 st =half*ss2*st2
198 vt2(1) = rs+rt+st
199 vt2(2) = rs-rt-st
200 vt2(3) = -rs+rt-st
201 vt2(4) = -rs-rt+st
202 max0 =max0+max(vt2(1),vt2(2),vt2(3),vt2(4))
203 min0 =min0+min(vt2(1),vt2(2),vt2(3),vt2(4))
204 rs =half*sr3*ss3
205 rt =half*sr3*st3
206 st =half*ss3*st3
207 vt2(1) = rs+rt+st
208 vt2(2) = rs-rt-st
209 vt2(3) = -rs+rt-st
210 vt2(4) = -rs-rt+st
211 max0 =max0+max(vt2(1),vt2(2),vt2(3),vt2(4))
212 min0 =min0+min(vt2(1),vt2(2),vt2(3),vt2(4))
213 svmrst = svm0+coef*(svmr+svms+svmt)
214 svm1(i) = sqrt(abs(svmrst+coef*max0))
215 svm2(i) = sqrt(max(svmrst+coef*min0,zero))
216 svm2(i) = min(svm2(i),svm0)
217 ENDIF
218 ENDDO
219C
220 RETURN
221 END
cg
subroutine cg(dim, mat, rhs, sol, max_iter, tol)
Definition conjugate_gradient.F:29
min
#define min(a, b)
Definition macros.h:20
max
#define max(a, b)
Definition macros.h:21
szsvm_or
subroutine szsvm_or(jr0, js0, jt0, cc, cg, g33, fhour, sigy, sig0, nu, svm1, svm2, nel, iint)
Definition szsvm_or.F:33