OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
s4volnod_sm.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!|| s4volnod_sm ../engine/source/elements/solid/solide4_sfem/s4volnod_sm.F
25!||--- called by ------------------------------------------------------
26!|| s4lagsfem ../engine/source/elements/solid/solide4_sfem/s4lagsfem.F
27!||--- calls -----------------------------------------------------
28!|| foat_to_6_float ../engine/source/system/parit.F
29!||====================================================================
30 SUBROUTINE s4volnod_sm(
31 1 VOLNOD6, V, NC1, NC2,
32 2 NC3, NC4, VOL0, AMU,
33 3 OFFG, SAV, NEL, ISMSTR)
34C-----------------------------------------------
35C I m p l i c i t T y p e s
36C-----------------------------------------------
37#include "implicit_f.inc"
38C-----------------------------------------------
39C G l o b a l P a r a m e t e r s
40C-----------------------------------------------
41#include "mvsiz_p.inc"
42C-----------------------------------------------
43C D u m m y A r g u m e n t s
44C-----------------------------------------------
45 INTEGER, INTENT(IN) :: ISMSTR
46 INTEGER NC1(*),NC2(*),NC3(*),NC4(*),NEL
47 DOUBLE PRECISION
48 . volnod6(6,2*numnod),sav(nel,9)
49C REAL
50 my_real
51 . v(3,*),offg(*),vol0(*),amu(*)
52C-----------------------------------------------
53C C o m m o n B l o c k s
54C-----------------------------------------------
55#include "scr18_c.inc"
56#include "com04_c.inc"
57#include "com08_c.inc"
58C-----------------------------------------------
59C L o c a l V a r i a b l e s
60C-----------------------------------------------
61 INTEGER I, K,INDEX(MVSIZ),NE_SM,J
62C REAL
63
64 double precision
65 . x41, y41, z41, x42, y42, z42, x43, y43, z43,
66 . dvol,detdp,b1dp,c1dp,d1dp
67 double precision
68 . x1(mvsiz), x2(mvsiz), x3(mvsiz), x4(mvsiz),
69 . y1(mvsiz), y2(mvsiz), y3(mvsiz), y4(mvsiz),
70 . z1(mvsiz), z2(mvsiz), z3(mvsiz), z4(mvsiz)
71 my_real
72 . b1(mvsiz), b2(mvsiz), b3(mvsiz), b4(mvsiz),
73 . c1(mvsiz), c2(mvsiz), c3(mvsiz), c4(mvsiz),
74 . d1(mvsiz), d2(mvsiz), d3(mvsiz), d4(mvsiz),
75 . det(mvsiz),vol(mvsiz),volref(mvsiz),d,dxx,dyy,dzz,
76 . px1(mvsiz), px2(mvsiz), px3(mvsiz), px4(mvsiz),
77 . py1(mvsiz), py2(mvsiz), py3(mvsiz), py4(mvsiz),
78 . pz1(mvsiz), pz2(mvsiz), pz3(mvsiz), pz4(mvsiz),
79 . vx1(mvsiz), vx2(mvsiz), vx3(mvsiz), vx4(mvsiz),
80 . vy1(mvsiz), vy2(mvsiz), vy3(mvsiz), vy4(mvsiz),
81 . vz1(mvsiz), vz2(mvsiz), vz3(mvsiz), vz4(mvsiz)
82
83 double precision
84 . det6(6,mvsiz)
85C----------------------------
86C NODAL COORDINATES
87C----------------------------
88 ne_sm = 0
89C-----case ISMSTR=1,11,2,12
90 IF(ismstr==1.OR.ismstr==11)THEN
91 ne_sm = nel
92 DO i=1,nel
93 vx1(i)=v(1,nc1(i))
94 vy1(i)=v(2,nc1(i))
95 vz1(i)=v(3,nc1(i))
96 vx2(i)=v(1,nc2(i))
97 vy2(i)=v(2,nc2(i))
98 vz2(i)=v(3,nc2(i))
99 vx3(i)=v(1,nc3(i))
100 vy3(i)=v(2,nc3(i))
101 vz3(i)=v(3,nc3(i))
102 vx4(i)=v(1,nc4(i))
103 vy4(i)=v(2,nc4(i))
104 vz4(i)=v(3,nc4(i))
105 volref(i) = vol0(i)
106 index(i)=i
107 ENDDO
108 IF(ismstr==11)THEN
109 DO i=1,nel
110 x1(i)=sav(i,1)
111 y1(i)=sav(i,4)
112 z1(i)=sav(i,7)
113 x2(i)=sav(i,2)
114 y2(i)=sav(i,5)
115 z2(i)=sav(i,8)
116 x3(i)=sav(i,3)
117 y3(i)=sav(i,6)
118 z3(i)=sav(i,9)
119 x4(i)=zero
120 y4(i)=zero
121 z4(i)=zero
122 ENDDO
123 ELSE
124 DO i=1,nel
125 x1(i)=sav(i,1)
126 y1(i)=sav(i,2)
127 z1(i)=sav(i,3)
128 x2(i)=sav(i,4)
129 y2(i)=sav(i,5)
130 z2(i)=sav(i,6)
131 x3(i)=sav(i,7)
132 y3(i)=sav(i,8)
133 z3(i)=sav(i,9)
134 x4(i)=zero
135 y4(i)=zero
136 z4(i)=zero
137 ENDDO
138 END IF
139 ELSEIF((ismstr==2.OR.ismstr==12).AND.idtmin(1)==3) THEN
140 DO i=1,nel
141 IF(abs(offg(i))>one)THEN
142 ne_sm = ne_sm+1
143 index(ne_sm)=i
144 x1(ne_sm)=sav(i,1)
145 y1(ne_sm)=sav(i,2)
146 z1(ne_sm)=sav(i,3)
147 x2(ne_sm)=sav(i,4)
148 y2(ne_sm)=sav(i,5)
149 z2(ne_sm)=sav(i,6)
150 x3(ne_sm)=sav(i,7)
151 y3(ne_sm)=sav(i,8)
152 z3(ne_sm)=sav(i,9)
153 x4(ne_sm)=zero
154 y4(ne_sm)=zero
155 z4(ne_sm)=zero
156 volref(ne_sm) = vol0(i)
157 vx1(ne_sm)=v(1,nc1(i))
158 vy1(ne_sm)=v(2,nc1(i))
159 vz1(ne_sm)=v(3,nc1(i))
160 vx2(ne_sm)=v(1,nc2(i))
161 vy2(ne_sm)=v(2,nc2(i))
162 vz2(ne_sm)=v(3,nc2(i))
163 vx3(ne_sm)=v(1,nc3(i))
164 vy3(ne_sm)=v(2,nc3(i))
165 vz3(ne_sm)=v(3,nc3(i))
166 vx4(ne_sm)=v(1,nc4(i))
167 vy4(ne_sm)=v(2,nc4(i))
168 vz4(ne_sm)=v(3,nc4(i))
169 ENDIF
170 ENDDO
171 ENDIF
172 IF (ne_sm==0) RETURN
173C----------------------------
174C TETRAHEDRON equivalent volume in smalll strain options
175C----------------------------
176C-----dt1=0 at T=0, no use to do
177C mean value on rate of vol/vol_ref
178 IF (tt==zero) THEN
179 vol(1:ne_sm) = volref(1:ne_sm)
180 ELSE
181 DO i=1,ne_sm
182 x43 = x4(i) - x3(i)
183 y43 = y4(i) - y3(i)
184 z43 = z4(i) - z3(i)
185 x41 = x4(i) - x1(i)
186 y41 = y4(i) - y1(i)
187 z41 = z4(i) - z1(i)
188 x42 = x4(i) - x2(i)
189 y42 = y4(i) - y2(i)
190 z42 = z4(i) - z2(i)
191C
192 b1dp = y43*z42 - y42*z43
193 c1dp = z43*x42 - z42*x43
194 d1dp = x43*y42 - x42*y43
195 b1(i) = b1dp
196 c1(i) = c1dp
197 d1(i) = d1dp
198 b2(i) = y41*z43 - y43*z41
199 b3(i) = y42*z41 - y41*z42
200 b4(i) = -(b1(i) + b2(i) + b3(i))
201C
202 c2(i) = z41*x43 - z43*x41
203 c3(i) = z42*x41 - z41*x42
204 c4(i) = -(c1(i) + c2(i) + c3(i))
205C
206 d2(i) = x41*y43 - x43*y41
207 d3(i) = x42*y41 - x41*y42
208 d4(i) = -(d1(i) + d2(i) + d3(i))
209C
210 detdp= (x41*b1dp + y41*c1dp + z41*d1dp)*one_over_6
211 det(i)= detdp
212 ENDDO
213 DO i=1,ne_sm
214 d = one/det(i)/six
215 px1(i)=-b1(i)*d
216 py1(i)=-c1(i)*d
217 pz1(i)=-d1(i)*d
218 px2(i)=-b2(i)*d
219 py2(i)=-c2(i)*d
220 pz2(i)=-d2(i)*d
221 px3(i)=-b3(i)*d
222 py3(i)=-c3(i)*d
223 pz3(i)=-d3(i)*d
224 px4(i)=-b4(i)*d
225 py4(i)=-c4(i)*d
226 pz4(i)=-d4(i)*d
227 END DO
228C-----compute equivalent (V/V0)*V0=(rho0/rho)*V0
229C-----rho/rho0=1+AMU(t)-dvol
230 DO i=1,ne_sm
231 dxx=px1(i)*vx1(i)+px2(i)*vx2(i)+
232 . px3(i)*vx3(i)+px4(i)*vx4(i)
233 dyy=py1(i)*vy1(i)+py2(i)*vy2(i)+
234 . py3(i)*vy3(i)+py4(i)*vy4(i)
235 dzz=pz1(i)*vz1(i)+pz2(i)*vz2(i)+
236 . pz3(i)*vz3(i)+pz4(i)*vz4(i)
237 dvol = (dxx+dyy+dzz)*dt1
238 detdp = volref(i)/max(em20,(one+amu(i)-dvol))
239 vol(i) = detdp
240 ENDDO
241 END IF ! (TT==ZERO) THEN
242
243 !Parith-On treatment
244 CALL foat_to_6_float(1 ,ne_sm ,vol ,det6 )
245
246 IF(ismstr==1.OR.ismstr==11)THEN
247 DO i=1,nel
248 !Parith-On treatment
249 DO k=1,6
250 volnod6(k,nc1(i)) = volnod6(k,nc1(i)) + det6(k,i)
251 volnod6(k,nc2(i)) = volnod6(k,nc2(i)) + det6(k,i)
252 volnod6(k,nc3(i)) = volnod6(k,nc3(i)) + det6(k,i)
253 volnod6(k,nc4(i)) = volnod6(k,nc4(i)) + det6(k,i)
254 ENDDO
255 ENDDO
256 ELSE
257#include "vectorize.inc"
258 DO j=1,ne_sm
259 i = index(j)
260 !Parith-On treatment
261 DO k=1,6
262 volnod6(k,nc1(i)) = volnod6(k,nc1(i)) + det6(k,j)
263 volnod6(k,nc2(i)) = volnod6(k,nc2(i)) + det6(k,j)
264 volnod6(k,nc3(i)) = volnod6(k,nc3(i)) + det6(k,j)
265 volnod6(k,nc4(i)) = volnod6(k,nc4(i)) + det6(k,j)
266 ENDDO
267 ENDDO
268 ENDIF
269C
270 RETURN
271 END
my_real
#define my_real
Definition cppsort.cpp:32
max
#define max(a, b)
Definition macros.h:21
foat_to_6_float
subroutine foat_to_6_float(jft, jlt, f, f6)
Definition parit.F:225
s4volnod_sm
subroutine s4volnod_sm(volnod6, v, nc1, nc2, nc3, nc4, vol0, amu, offg, sav, nel, ismstr)
Definition s4volnod_sm.F:34