OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
s8cforc3.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!|| s8cforc3 ../engine/source/elements/thickshell/solide8c/s8cforc3.F
25!||--- called by ------------------------------------------------------
26!|| forint ../engine/source/elements/forint.F
27!||--- calls -----------------------------------------------------
28!|| basis8 ../engine/source/elements/solid/solide8/basis8.F
29!|| boltst ../engine/source/elements/solid/solide/boltst.F
30!|| csmall3 ../engine/source/elements/solid/solide/csmall3.F
31!|| degenes8 ../engine/source/elements/solid/solide/degenes8.F
32!|| mmain ../engine/source/materials/mat_share/mmain.F90
33!|| s8cdefo3 ../engine/source/elements/thickshell/solide8c/s8cdefo3.F
34!|| s8cfint_reg ../engine/source/elements/thickshell/solide8c/s8cfint_reg.F
35!|| s8csigp3 ../engine/source/elements/thickshell/solide8c/s8csigp3.F
36!|| s8etherm ../engine/source/elements/solid/solide8e/s8etherm.F
37!|| s8for_distor ../engine/source/elements/solid/solide/s8for_distor.F
38!|| s8get_x3 ../engine/source/elements/solid/solide/s8get_x3.F
39!|| s8sav3 ../engine/source/elements/solid/solide/s8sav3.F
40!|| s8zderi3 ../engine/source/elements/solid/solide8z/s8zderi3.F
41!|| s8zderic3 ../engine/source/elements/solid/solide8z/s8zderic3.F
42!|| s8zfint3 ../engine/source/elements/solid/solide8z/s8zfint3.F
43!|| s8zfintp3 ../engine/source/elements/solid/solide8z/s8zfintp3.F
44!|| s8ztempel ../engine/source/elements/solid/solide8z/s8ztempel.F
45!|| s8zzero3 ../engine/source/elements/solid/solide8z/s8zzero3.F
46!|| sbilan ../engine/source/elements/solid/solide/sbilan.F
47!|| scoor3 ../engine/source/elements/solid/solide/scoor3.F
48!|| sctorth3 ../engine/source/elements/thickshell/solide8c/sctorth3.F
49!|| scumu3 ../engine/source/elements/solid/solide/scumu3.F
50!|| scumu3p ../engine/source/elements/solid/solide/scumu3p.F
51!|| scumualpha ../engine/source/elements/thickshell/solidec/scumualpha.F
52!|| sdistor_ini ../engine/source/elements/solid/solide/sdistror_ini.F90
53!|| sdlen_dege ../engine/source/elements/solid/solide/sdlen_dege.F
54!|| sdlensh14 ../engine/source/elements/thickshell/solide8c/sdlensh14.F
55!|| sdlensh8 ../engine/source/elements/thickshell/solide8c/sdlensh8.F
56!|| sfillopt ../engine/source/elements/solid/solide/sfillopt.F
57!|| sgparav3 ../engine/source/elements/solid/solide/sgparav3.F
58!|| smalla3 ../engine/source/elements/solid/solide/smalla3.F
59!|| smallb3 ../engine/source/elements/solid/solide/smallb3.F
60!|| srbilan ../engine/source/elements/solid/solide/srbilan.F
61!|| srcoor3 ../engine/source/elements/solid/solide/srcoor3.F
62!|| srho3 ../engine/source/elements/solid/solide/srho3.F
63!|| srota3 ../engine/source/elements/solid/solide/srota3.F
64!|| sroto3 ../engine/source/elements/solid/solidez/sroto3.F
65!|| srrota3 ../engine/source/elements/solid/solide/srrota3.f
66!|| sstra3 ../engine/source/elements/solid/solide/sstra3.F
67!|| startime ../engine/source/system/timer_mod.F90
68!|| stoptime ../engine/source/system/timer_mod.F90
69!|| storth3 ../engine/source/elements/solid/solidez/storth3.F
70!|| szordef3 ../engine/source/elements/solid/solidez/szordef3.F
71!|| tshgeodel3 ../engine/source/elements/thickshell/solidec/tshgeodel3.F
72!||--- uses -----------------------------------------------------
73!|| ale_connectivity_mod ../common_source/modules/ale/ale_connectivity_mod.f
74!|| dt_mod ../engine/source/modules/dt_mod.F
75!|| elbufdef_mod ../common_source/modules/mat_elem/elbufdef_mod.F90
76!|| glob_therm_mod ../common_source/modules/mat_elem/glob_therm_mod.F90
77!|| mat_elem_mod ../common_source/modules/mat_elem/mat_elem_mod.F90
78!|| mmain_mod ../engine/source/materials/mat_share/mmain.F90
79!|| nlocal_reg_mod ../common_source/modules/nlocal_reg_mod.F
80!|| output_mod ../common_source/modules/output/output_mod.F90
81!|| sdistor_ini_mod ../engine/source/elements/solid/solide/sdistror_ini.F90
82!|| sensor_mod ../common_source/modules/sensor_mod.F90
83!|| table_mod ../engine/share/modules/table_mod.F
84!|| timer_mod ../engine/source/system/timer_mod.F90
85!||====================================================================
86 SUBROUTINE s8cforc3(TIMERS, OUTPUT, ELBUF_TAB,NG ,
87 1 PM ,GEO ,IXS ,X ,
88 2 A ,V ,MS ,W ,FLUX ,
89 3 FLU1 ,VEUL ,FV ,ALE_CONNECT ,IPARG ,
90 4 TF ,NPF ,BUFMAT ,PARTSAV,NLOC_DMG,
91 5 DT2T ,NELTST ,ITYPTST,STIFN ,FSKY ,
92 6 IADS ,OFFSET ,EANI ,IPARTS ,
93 7 F11 ,F21 ,F31 ,F12 ,F22 ,
94 8 F32 ,F13 ,F23 ,F33 ,F14 ,
95 9 F24 ,F34 ,F15 ,F25 ,F35 ,
96 A F16 ,F26 ,F36 ,F17 ,F27 ,
97 B F37 ,F18 ,F28 ,F38 ,NEL ,
98 C ICP ,
99 F ICSIG ,SMR ,SMS ,SMT ,MFXX ,
100 G MFXY ,MFXZ ,MFYX ,MFYY ,MFYZ ,
101 H MFZX ,MFZY ,MFZZ ,NVC ,IPM ,
102 I ITASK ,ISTRAIN,TEMP ,FTHE ,
103 J FTHESKY ,IEXPAN ,IGEO ,NNPT ,GRESAV ,
104 K GRTH ,IGRTH ,MSSA ,DMELS ,TABLE ,
105 L XDP ,VOLN ,CONDN ,CONDNSKY,SENSORS,
106 M IOUTPRT ,MAT_ELEM,H3D_STRAIN,DT ,NODADT ,
107 N DTFAC1 ,DTMIN1, IDTMIN,SNPC,STF,SBUFMAT,SVIS,
108 O NSVOIS ,IDTMINS,IRESP,MAXFUNC,USERL_AVAIL,
109 V GLOB_THERM,impl_s,idyna)
110C-----------------------------------------------
111C M o d u l e s
112C-----------------------------------------------
113 USE timer_mod
114 USE output_mod, only : output_
115 USE mmain_mod
116 USE table_mod
117 USE mat_elem_mod
118 USE nlocal_reg_mod
119 USE ale_connectivity_mod
120 USE sensor_mod
121 USE dt_mod
122 USE elbufdef_mod
123 USE sdistor_ini_mod, ONLY : sdistor_ini
124 use glob_therm_mod
125C-----------------------------------------------
126C I m p l i c i t T y p e s
127C-----------------------------------------------
128#include "implicit_f.inc"
129C-----------------------------------------------
130C G l o b a l P a r a m e t e r s
131C-----------------------------------------------
132#include "mvsiz_p.inc"
133C-----------------------------------------------
134C C o m m o n B l o c k s
135C-----------------------------------------------
136#include "vect01_c.inc"
137#include "com01_c.inc"
138#include "com04_c.inc"
139#include "com08_c.inc"
140#include "scr17_c.inc"
141#include "parit_c.inc"
142#include "param_c.inc"
143#include "timeri_c.inc"
144C-----------------------------------------------
145C D u m m y A r g u m e n t s
146C-----------------------------------------------
147 TYPE(timer_), INTENT(inout) :: TIMERS
148 TYPE(OUTPUT_), INTENT(inout) :: OUTPUT
149 INTEGER, INTENT(IN) :: SNPC
150 INTEGER, INTENT(IN) :: STF
151 INTEGER, INTENT(IN) :: SBUFMAT
152 INTEGER, INTENT(IN) :: NSVOIS
153 INTEGER, INTENT(IN) :: IDTMINS
154 INTEGER ,INTENT(IN) :: IRESP
155 INTEGER ,INTENT(IN) :: MAXFUNC
156 INTEGER, INTENT(IN) :: IMPL_S
157 INTEGER, INTENT(IN) :: IDYNA
158 INTEGER, INTENT(IN) :: USERL_AVAIL
159 INTEGER IXS(NIXS,*), IPARG(NPARG,NGROUP),NPF(*),
160 . IPARTS(*), IPM(NPROPMI,*),IGEO(NPROPGI,*),GRTH(*),
161 . IGRTH(*),IADS(8,*),IOUTPRT,H3D_STRAIN
162 INTEGER, INTENT(IN) :: NODADT,IDTMIN(102)
163 my_real, INTENT(IN) :: DTFAC1(102),DTMIN1(102)
164C
165 INTEGER NELTST,ITYPTST,OFFSET,NEL,ICP,
166 . ICSIG, NVC,ITASK, ISTRAIN, IEXPAN ,NNPT,NG
167
168 DOUBLE PRECISION
169 . XDP(3,*)
170
171 my_real
172 . dt2t
173C
174 my_real
175 . pm(npropm,*),geo(npropg,*),x(*),a(*),v(3,*),ms(*),w(*),
176 . flux(6,*),flu1(*),veul(*),fv(*),tf(*),bufmat(*),
177 . partsav(*),stifn(*), fsky(*),eani(*),
178 . f11(mvsiz),f21(mvsiz),f31(mvsiz),
179 . f12(mvsiz),f22(mvsiz),f32(mvsiz),
180 . f13(mvsiz),f23(mvsiz),f33(mvsiz),
181 . f14(mvsiz),f24(mvsiz),f34(mvsiz),
182 . f15(mvsiz),f25(mvsiz),f35(mvsiz),
183 . f16(mvsiz),f26(mvsiz),f36(mvsiz),
184 . f17(mvsiz),f27(mvsiz),f37(mvsiz),
185 . f18(mvsiz),f28(mvsiz),f38(mvsiz)
186 my_real
187 . smr(mvsiz,729),sms(mvsiz,729),smt(mvsiz,729),
188 . mfxx(mvsiz,729),mfxy(mvsiz,729),mfyx(mvsiz,729),
189 . mfyy(mvsiz,729),mfyz(mvsiz,729),mfzy(mvsiz,729),
190 . mfzz(mvsiz,729),mfzx(mvsiz,729),mfxz(mvsiz,729),
191 . temp(*), fthe(*), fthesky(*),gresav(*), mssa(*),
192 . dmels(*), voln(mvsiz),condn(*),condnsky(*)
193 my_real, DIMENSION(MVSIZ,6), INTENT(INOUT) :: svis
194 TYPE(ttable) TABLE(*)
195 TYPE (ELBUF_STRUCT_), DIMENSION(NGROUP), TARGET :: ELBUF_TAB
196 TYPE (NLOCAL_STR_) , TARGET :: NLOC_DMG
197 TYPE(t_ale_connectivity), INTENT(IN) :: ALE_CONNECT
198 TYPE (MAT_ELEM_) ,INTENT(INOUT) :: MAT_ELEM
199 TYPE (SENSORS_) ,INTENT(INOUT) :: SENSORS
200 TYPE(DT_), INTENT(INOUT) :: DT
201 type (glob_therm_) ,intent(inout) :: glob_therm
202C-----------------------------------------------
203C L o c a l V a r i a b l e s
204C-----------------------------------------------
205 INTEGER I,J,N,LCO,NF1,IFLAG,IAD0,NLAY,L_PLA,L_EPSD,
206 . IP,IR,IS,IT,IL,NPTR,NPTS,NPTT,ICR,ICS,ICT,PID,MTN0,IPTHK,
207 . IPPOS,IPMAT,NLYMAX,MID,IPANG,IOR_TSH,NPTR0,NPTS0,NPTT0,IBID,MX,
208 . II(6),CURRENT_LAYER
209 INTEGER MXT(MVSIZ),NGL(MVSIZ),NGEO(MVSIZ),MXT0(MVSIZ),IBIDON(1),
210 . IDEG(MVSIZ)
211
212 DOUBLE PRECISION
213 . X0(MVSIZ,8),Y0(MVSIZ,8),Z0(MVSIZ,8),
214 . XD1(MVSIZ), XD2(MVSIZ), XD3(MVSIZ), XD4(MVSIZ),
215 . XD5(MVSIZ), XD6(MVSIZ), XD7(MVSIZ), XD8(MVSIZ),
216 . yd1(mvsiz), yd2(mvsiz), yd3(mvsiz), yd4(mvsiz),
217 . yd5(mvsiz), yd6(mvsiz), yd7(mvsiz), yd8(mvsiz),
218 . zd1(mvsiz), zd2(mvsiz), zd3(mvsiz), zd4(mvsiz),
219 . zd5(mvsiz), zd6(mvsiz), zd7(mvsiz), zd8(mvsiz),voldp(mvsiz)
220
221 my_real
222 . vd2(mvsiz) , dvol(mvsiz),deltax(mvsiz),
223 . vis(mvsiz) , qvis(mvsiz), cxx(mvsiz) ,
224 . s1(mvsiz) , s2(mvsiz) , s3(mvsiz) ,
225 . s4(mvsiz) , s5(mvsiz) , s6(mvsiz) ,
226 . dxx(mvsiz) , dyy(mvsiz) , dzz(mvsiz) ,
227 . d4(mvsiz) , d5(mvsiz) , d6(mvsiz) ,
228 . ajc1(mvsiz) , ajc2(mvsiz) , ajc3(mvsiz) ,
229 . ajc4(mvsiz) , ajc5(mvsiz) , ajc6(mvsiz) ,
230 . ajc7(mvsiz) , ajc8(mvsiz) , ajc9(mvsiz) ,
231 . aj1(mvsiz) , aj2(mvsiz) , aj3(mvsiz) ,
232 . aj4(mvsiz) , aj5(mvsiz) , aj6(mvsiz) ,
233 . vdx(mvsiz) , vdy(mvsiz) , vdz(mvsiz),ssp_eq(mvsiz),aire(mvsiz),
234 . e0(mvsiz),c1,fac(mvsiz),them(mvsiz,8),tempel(mvsiz),ni(8),
235 . die(mvsiz),conde(mvsiz),conden(mvsiz),amu(mvsiz)
236C-----
237C Variables utilisees en argument par les materiaux.
238 my_real
239 . sti(mvsiz),stin(mvsiz),gama(mvsiz,6),
240 . wxx(mvsiz) , wyy(mvsiz) , wzz(mvsiz)
241C Variables utilisees en argument par les materiaux si SPH uniquement.
242 my_real
243 . muvoid(mvsiz)
244 INTEGER IOFFS,IFVM22,NN_DEL,IPRES,ISCTL,ISTAB(MVSIZ)
245 my_real
246 . OFFS(MVSIZ),DSV(MVSIZ)
247C-----
248C Variables utilisees dans les routines solides uniquement (en arguments).
249 INTEGER NC1(MVSIZ), NC2(MVSIZ), NC3(MVSIZ), NC4(MVSIZ),
250 . NC5(MVSIZ), NC6(MVSIZ), NC7(MVSIZ), NC8(MVSIZ)
251 my_real
252 . OFF(MVSIZ) ,OFFL(MVSIZ) , RHOO(MVSIZ),
253 . X1(MVSIZ), X2(MVSIZ), X3(MVSIZ), X4(MVSIZ),
254 . X5(MVSIZ), X6(MVSIZ), X7(MVSIZ), X8(MVSIZ),
255 . Y1(MVSIZ), Y2(MVSIZ), Y3(MVSIZ), Y4(MVSIZ),
256 . Y5(MVSIZ), Y6(MVSIZ), Y7(MVSIZ), Y8(MVSIZ),
257 . Z1(MVSIZ), Z2(MVSIZ), Z3(MVSIZ), Z4(MVSIZ),
258 . Z5(MVSIZ), Z6(MVSIZ), Z7(MVSIZ), Z8(MVSIZ),
259 . VX1(MVSIZ),VX2(MVSIZ),VX3(MVSIZ),VX4(MVSIZ),
260 . VX5(MVSIZ),VX6(MVSIZ),VX7(MVSIZ),VX8(MVSIZ),
261 . vy1(mvsiz),vy2(mvsiz),vy3(mvsiz),vy4(mvsiz),
262 . vy5(mvsiz),vy6(mvsiz),vy7(mvsiz),vy8(mvsiz),
263 . vz1(mvsiz),vz2(mvsiz),vz3(mvsiz),vz4(mvsiz),
264 . vz5(mvsiz),vz6(mvsiz),vz7(mvsiz),vz8(mvsiz),
265 . hx1(mvsiz),hx2(mvsiz),hx3(mvsiz),hx4(mvsiz),
266 . hy1(mvsiz),hy2(mvsiz),hy3(mvsiz),hy4(mvsiz),
267 . hz1(mvsiz),hz2(mvsiz),hz3(mvsiz),hz4(mvsiz),
268 . px1(mvsiz),px2(mvsiz),px3(mvsiz),px4(mvsiz),
269 . px5(mvsiz),px6(mvsiz),px7(mvsiz),px8(mvsiz),
270 . py1(mvsiz),py2(mvsiz),py3(mvsiz),py4(mvsiz),
271 . py5(mvsiz),py6(mvsiz),py7(mvsiz),py8(mvsiz),
272 . pz1(mvsiz),pz2(mvsiz),pz3(mvsiz),pz4(mvsiz),
273 . pz5(mvsiz),pz6(mvsiz),pz7(mvsiz),pz8(mvsiz),
274 . pxy1(mvsiz),pxy2(mvsiz),pxy3(mvsiz),pxy4(mvsiz),
275 . pxy5(mvsiz),pxy6(mvsiz),pxy7(mvsiz),pxy8(mvsiz),
276 . pyx1(mvsiz),pyx2(mvsiz),pyx3(mvsiz),pyx4(mvsiz),
277 . pyx5(mvsiz),pyx6(mvsiz),pyx7(mvsiz),pyx8(mvsiz),
278 . pxz1(mvsiz),pxz2(mvsiz),pxz3(mvsiz),pxz4(mvsiz),
279 . pxz5(mvsiz),pxz6(mvsiz),pxz7(mvsiz),pxz8(mvsiz),
280 . pzx1(mvsiz),pzx2(mvsiz),pzx3(mvsiz),pzx4(mvsiz),
281 . pzx5(mvsiz),pzx6(mvsiz),pzx7(mvsiz),pzx8(mvsiz),
282 . pyz1(mvsiz),pyz2(mvsiz),pyz3(mvsiz),pyz4(mvsiz),
283 . pyz5(mvsiz),pyz6(mvsiz),pyz7(mvsiz),pyz8(mvsiz),
284 . pzy1(mvsiz),pzy2(mvsiz),pzy3(mvsiz),pzy4(mvsiz),
285 . pzy5(mvsiz),pzy6(mvsiz),pzy7(mvsiz),pzy8(mvsiz),
286 . pxc1(mvsiz),pxc2(mvsiz),pxc3(mvsiz),pxc4(mvsiz),
287 . pyc1(mvsiz),pyc2(mvsiz),pyc3(mvsiz),pyc4(mvsiz),
288 . pzc1(mvsiz),pzc2(mvsiz),pzc3(mvsiz),pzc4(mvsiz),
289 . bxy1(mvsiz),bxy2(mvsiz),bxy3(mvsiz),bxy4(mvsiz),
290 . bxy5(mvsiz),bxy6(mvsiz),bxy7(mvsiz),bxy8(mvsiz),
291 . byx1(mvsiz),byx2(mvsiz),byx3(mvsiz),byx4(mvsiz),
292 . byx5(mvsiz),byx6(mvsiz),byx7(mvsiz),byx8(mvsiz),
293 . bxz1(mvsiz),bxz2(mvsiz),bxz3(mvsiz),bxz4(mvsiz),
294 . bxz5(mvsiz),bxz6(mvsiz),bxz7(mvsiz),bxz8(mvsiz),
295 . bzx1(mvsiz),bzx2(mvsiz),bzx3(mvsiz),bzx4(mvsiz),
296 . bzx5(mvsiz),bzx6(mvsiz),bzx7(mvsiz),bzx8(mvsiz),
297 . byz1(mvsiz),byz2(mvsiz),byz3(mvsiz),byz4(mvsiz),
298 . byz5(mvsiz),byz6(mvsiz),byz7(mvsiz),byz8(mvsiz),
299 . bzy1(mvsiz),bzy2(mvsiz),bzy3(mvsiz),bzy4(mvsiz),
300 . bzy5(mvsiz),bzy6(mvsiz),bzy7(mvsiz),bzy8(mvsiz),
301 . px1h1(mvsiz),px2h1(mvsiz),px3h1(mvsiz),px4h1(mvsiz),
302 . px1h2(mvsiz),px2h2(mvsiz),px3h2(mvsiz),px4h2(mvsiz),
303 . px1h3(mvsiz),px2h3(mvsiz),px3h3(mvsiz),px4h3(mvsiz),
304 . px1h4(mvsiz),px2h4(mvsiz),px3h4(mvsiz),px4h4(mvsiz),
305 . vdx1(mvsiz),vdx2(mvsiz),vdx3(mvsiz),vdx4(mvsiz),
306 . vdx5(mvsiz),vdx6(mvsiz),vdx7(mvsiz),vdx8(mvsiz),
307 . vdy1(mvsiz),vdy2(mvsiz),vdy3(mvsiz),vdy4(mvsiz),
308 . vdy5(mvsiz),vdy6(mvsiz),vdy7(mvsiz),vdy8(mvsiz),
309 . vdz1(mvsiz),vdz2(mvsiz),vdz3(mvsiz),vdz4(mvsiz),
310 . vdz5(mvsiz),vdz6(mvsiz),vdz7(mvsiz),vdz8(mvsiz),
311 . vgxa(mvsiz),vgya(mvsiz),vgza(mvsiz), vga2(mvsiz),
312 . hx(mvsiz,4), hy(mvsiz,4), hz(mvsiz,4),
313 . xgxa(mvsiz),xgya(mvsiz),xgza(mvsiz),
314 . xgxya(mvsiz),xgyza(mvsiz),xgzxa(mvsiz),
315 . xgxa2(mvsiz),xgya2(mvsiz),xgza2(mvsiz)
316 my_real
317 . dxy(mvsiz),dyx(mvsiz),
318 . dyz(mvsiz),dzy(mvsiz),
319 . dzx(mvsiz),dxz(mvsiz),divde(mvsiz)
320 my_real
321 . r11(mvsiz),r12(mvsiz),r13(mvsiz),
322 . r21(mvsiz),r22(mvsiz),r23(mvsiz),
323 . r31(mvsiz),r32(mvsiz),r33(mvsiz),
324 . g1x(mvsiz),g2x(mvsiz),g3x(mvsiz),
325 . g1y(mvsiz),g2y(mvsiz),g3y(mvsiz),
326 . g1z(mvsiz),g2z(mvsiz),g3z(mvsiz)
327 my_real
328 . wi,smax(mvsiz),volg(mvsiz),nu(mvsiz),pp(mvsiz),usb(mvsiz),
329 . volm(mvsiz),sigm(mvsiz),dti,bid(mvsiz)
330 my_real
331 . sigy(mvsiz), sign(nel,6),et(mvsiz), nu1(mvsiz),
332 . r1_free(mvsiz),r3_free(mvsiz),r4_free(mvsiz)
333 my_real
334 . vx0(mvsiz,8),vy0(mvsiz,8),vz0(mvsiz,8),
335 . alpha_e(mvsiz),llsh(mvsiz),llsmin(mvsiz),area(mvsiz)
336 DOUBLE PRECISION
337 . FACDP
338 my_real
339 . shf(mvsiz),zr,zs,zt,wt,zz,
340 . rx(mvsiz), ry(mvsiz), rz(mvsiz),
341 . sx(mvsiz), sy(mvsiz), sz(mvsiz),
342 . tx(mvsiz), ty(mvsiz), tz(mvsiz)
343 my_real, DIMENSION(MVSIZ,NNPT) :: smtt
344 my_real,
345 . DIMENSION(:), POINTER :: eint
346 my_real, dimension(mvsiz) :: fheat
347C-----
348 my_real,
349 . DIMENSION(:,:), ALLOCATABLE ::
350 . var_reg,pxx1,pxx2,pxx3,pxx4,pxx5,pxx6,pxx7,pxx8,
351 . pyy1,pyy2,pyy3,pyy4,pyy5,pyy6,pyy7,pyy8,
352 . pzz1,pzz2,pzz3,pzz4,pzz5,pzz6,pzz7,pzz8
353 INTEGER :: INLOC, L_NLOC, INOD(8), IPOS(8)
354 my_real,
355 . DIMENSION(:), POINTER :: DNL
356 my_real
357 . h(8),pr(8),ps(8),pt(8)
358 my_real
359 . ev(mvsiz),evd(mvsiz) !!!! BIDON !!!!!
360 TYPE(g_bufel_) ,POINTER :: GBUF
361 TYPE(l_bufel_) ,POINTER :: LBUF
362 INTEGER SZ_IX
363 my_real :: STI_C(MVSIZ),LL(MVSIZ),FLD(MVSIZ),
364 . CNS2,FQMAX
365c-----------------------------------------------------
366 parameter(nlymax = 200,ipmat = 100,ipang = 200)
367c-----------------------------------------------------
368 my_real
369 . w_gauss(9,9),a_gauss(9,9)
370 DATA w_gauss /
371 1 2. ,0. ,0. ,
372 1 0. ,0. ,0. ,
373 1 0. ,0. ,0. ,
374 2 1. ,1. ,0. ,
375 2 0. ,0. ,0. ,
376 2 0. ,0. ,0. ,
377 3 0.555555555555556,0.888888888888889,0.555555555555556,
378 3 0. ,0. ,0. ,
379 3 0. ,0. ,0. ,
380 4 0.347854845137454,0.652145154862546,0.652145154862546,
381 4 0.347854845137454,0. ,0. ,
382 4 0. ,0. ,0. ,
383 5 0.236926885056189,0.478628670499366,0.568888888888889,
384 5 0.478628670499366,0.236926885056189,0. ,
385 5 0. ,0. ,0. ,
386 6 0.171324492379170,0.360761573048139,0.467913934572691,
387 6 0.467913934572691,0.360761573048139,0.171324492379170,
388 6 0. ,0. ,0. ,
389 7 0.129484966168870,0.279705391489277,0.381830050505119,
390 7 0.417959183673469,0.381830050505119,0.279705391489277,
391 7 0.129484966168870,0. ,0. ,
392 8 0.101228536290376,0.222381034453374,0.313706645877887,
393 8 0.362683783378362,0.362683783378362,0.313706645877887,
394 8 0.222381034453374,0.101228536290376,0. ,
395 9 0.081274388361574,0.180648160694857,0.260610696402935,
396 9 0.312347077040003,0.330239355001260,0.312347077040003,
397 9 0.260610696402935,0.180648160694857,0.081274388361574/
398 DATA a_gauss /
399 1 0. ,0. ,0. ,
400 1 0. ,0. ,0. ,
401 1 0. ,0. ,0. ,
402 2 -.577350269189626,0.577350269189626,0. ,
403 2 0. ,0. ,0. ,
404 2 0. ,0. ,0. ,
405 3 -.774596669241483,0. ,0.774596669241483,
406 3 0. ,0. ,0. ,
407 3 0. ,0. ,0. ,
408 4 -.861136311594053,-.339981043584856,0.339981043584856,
409 4 0.861136311594053,0. ,0. ,
410 4 0. ,0. ,0. ,
411 5 -.906179845938664,-.538469310105683,0. ,
412 5 0.538469310105683,0.906179845938664,0. ,
413 5 0. ,0. ,0. ,
414 6 -.932469514203152,-.661209386466265,-.238619186083197,
415 6 0.238619186083197,0.661209386466265,0.932469514203152,
416 6 0. ,0. ,0. ,
417 7 -.949107912342759,-.741531185599394,-.405845151377397,
418 7 0. ,0.405845151377397,0.741531185599394,
419 7 0.949107912342759,0. ,0. ,
420 8 -.960289856497536,-.796666477413627,-.525532409916329,
421 8 -.183434642495650,0.183434642495650,0.525532409916329,
422 8 0.796666477413627,0.960289856497536,0. ,
423 9 -.968160239507626,-.836031107326636,-.613371432700590,
424 9 -.324253423403809,0. ,0.324253423403809,
425 9 0.613371432700590,0.836031107326636,0.968160239507626/
426c-----------------------------------------------------
427c Flag Bolt Preloading
428 INTEGER IBOLTP,NBPRELD
429 INTEGER SZ_R1_FREE
430 my_real,
431 . DIMENSION(:), POINTER :: BPRELD
432C-----------------------------------------------
433C S o u r c e L i n e s
434C=======================================================================
435 SZ_IX=numelq+numels+nsvois
436 sz_r1_free=mvsiz
437 gbuf => elbuf_tab(ng)%GBUF
438 tempel(:) = zero
439 fheat(:) = zero
440!
441 DO i=1,6
442 ii(i) = nel*(i-1)
443 ENDDO
444!
445 iboltp = iparg(72,ng)
446 nbpreld = gbuf%G_BPRELD
447 bpreld =>gbuf%BPRELD(1:nbpreld*nel)
448C-----
449 ibid = 0
450 ibidon(1) = 0
451 bid = zero
452 IF (jcvt==1 .AND. isorth==1) jcvt=2
453 iad0 = 1
454 IF (igtyp /= 22) THEN
455 IF (isorth > 0) iad0 = 1 + 6*nel
456 isorthg = 0
457 END IF
458 ior_tsh = 0
459 IF (igtyp == 21.OR.igtyp == 22) ior_tsh = 1
460C
461C-----------
462 nf1=nft+1
463C--------------------------
464 IF (ior_tsh > 0) THEN
465 CALL sgparav3(
466 1 8, x, ixs(1,nf1),rx,
467 2 ry, rz, sx, sy,
468 3 sz, tx, ty, tz,
469 4 nel)
470 ENDIF
471C-----------
472C GATHER NODAL VARIABLES AND COMPUTE INTRINSIC ROTATION.
473C-----------
474 IF (jcvt == 0) THEN
475C GATHER NODAL VARIABLES
476 CALL scoor3(x,ixs(1,nf1),v,w,gbuf%GAMA,gama,
477 2 x1, x2, x3, x4, x5, x6, x7, x8,
478 3 y1, y2, y3, y4, y5, y6, y7, y8,
479 4 z1, z2, z3, z4, z5, z6, z7, z8,
480 5 vx1, vx2, vx3, vx4, vx5, vx6, vx7, vx8,
481 6 vy1, vy2, vy3, vy4, vy5, vy6, vy7, vy8,
482 7 vz1, vz2, vz3, vz4, vz5, vz6, vz7, vz8,
483 8 vdx1, vdx2, vdx3, vdx4, vdx5, vdx6, vdx7, vdx8,
484 9 vdy1, vdy2, vdy3, vdy4, vdy5, vdy6, vdy7, vdy8,
485 a vdz1, vdz2, vdz3, vdz4, vdz5, vdz6, vdz7, vdz8,
486 b vdx,vdy,vdz,vd2,vis,gbuf%OFF,offl,gbuf%SMSTR,
487 c gbuf%RHO,rhoo,nc1,nc2,nc3,nc4,nc5,nc6,
488 c nc7,nc8,ngl,mxt,ngeo,
489 e bid,bid,bid,bid,bid,bid,
490 f xd1, xd2 , xd3, xd4, xd5, xd6, xd7, xd8,
491 g yd1, yd2 , yd3, yd4, yd5, yd6, yd7, yd8,
492 h zd1, zd2 , zd3, zd4, zd5, zd6, zd7, zd8,
493 i xdp, iparg, ng , nel)
494 ELSE
495C GATHER NODAL VARIABLES AND COMPUTE INTRINSIC ROTATION.
496 CALL srcoor3(x,ixs(1,nf1),v,w,gbuf%GAMA,gama,
497 . x1, x2, x3, x4, x5, x6, x7, x8,
498 . y1, y2, y3, y4, y5, y6, y7, y8,
499 . z1, z2, z3, z4, z5, z6, z7, z8,
500 . vx1, vx2, vx3, vx4, vx5, vx6, vx7, vx8,
501 . vy1, vy2, vy3, vy4, vy5, vy6, vy7, vy8,
502 . vz1, vz2, vz3, vz4, vz5, vz6, vz7, vz8,
503 . vd2,vis,gbuf%OFF,offl,gbuf%SMSTR,gbuf%RHO,rhoo,
504 . r11, r12, r13, r21, r22, r23, r31, r32, r33,
505 . nc1,nc2,nc3,nc4,nc5,nc6,nc7,nc8,
506 . ngl,mxt,ngeo,ioutprt, vgxa, vgya, vgza, vga2,
507 . xd1, xd2, xd3, xd4, xd5, xd6, xd7, xd8,
508 . yd1, yd2, yd3, yd4, yd5, yd6, yd7, yd8,
509 . zd1, zd2, zd3, zd4, zd5, zd6, zd7, zd8,
510 . xdp, x0 , y0 , z0 , nel, xgxa, xgya, xgza,
511 . xgxa2,xgya2,xgza2,xgxya,xgyza,xgzxa,iparg(1,ng),
512 . gbuf%GAMA_R)
513 ENDIF
514 nn_del = 0
515 pid = ngeo(1)
516 IF (geo(190,pid)+geo(191,pid)+geo(192,pid)+geo(192,pid)>zero)
517 . nn_del=8
518 IF (nn_del ==0 .AND. dt%IDEL_BRICK>0) nn_del=8
519 isctl = igeo(97,pid)
520c
521 IF (jcvt == 2 .AND. ior_tsh == 0) THEN
522 CALL storth3(
523 1 lft, llt, nel, g1x,
524 2 g1y, g1z, g2x, g2y,
525 3 g2z, g3x, g3y, g3z,
526 4 gbuf%GAMA,isorth)
527 ENDIF
528c
529 mx = mxt(1)
530 ipres = mat_elem%MAT_PARAM(mx)%IPRES
531 DO i=lft,llt
532 nu(i)=min(half,pm(21,mxt(i)))
533 deltax(i)=ep30
534 c1 =pm(32,mxt(i))
535 e0(i) =three*(one-two*nu(i))*c1
536 usb(i)=em01/c1
537 ideg(i)=0
538 ENDDO
539c
540C---------R,S,T->S,T,R-------------
541 ics=icsig/100
542 ict=mod(icsig/10,10)
543 icr=mod(icsig,10)
544
545 nlay=elbuf_tab(ng)%NLAY
546 nptr=elbuf_tab(ng)%NPTR
547 npts=elbuf_tab(ng)%NPTS
548 nptt=elbuf_tab(ng)%NPTT
549 it = 1
550 inloc = iparg(78,ng)
551 ALLOCATE(var_reg(nel,nptr*npts*nlay))
552 IF (inloc > 0) ALLOCATE(
553 . pxx1(nel,nptr*npts*nlay),pxx2(nel,nptr*npts*nlay),pxx3(nel,nptr*npts*nlay),
554 . pxx4(nel,nptr*npts*nlay),pxx5(nel,nptr*npts*nlay),pxx6(nel,nptr*npts*nlay),
555 . pxx7(nel,nptr*npts*nlay),pxx8(nel,nptr*npts*nlay),pyy1(nel,nptr*npts*nlay),
556 . pyy2(nel,nptr*npts*nlay),pyy3(nel,nptr*npts*nlay),pyy4(nel,nptr*npts*nlay),
557 . pyy5(nel,nptr*npts*nlay),pyy6(nel,nptr*npts*nlay),pyy7(nel,nptr*npts*nlay),
558 . pyy8(nel,nptr*npts*nlay),pzz1(nel,nptr*npts*nlay),pzz2(nel,nptr*npts*nlay),
559 . pzz3(nel,nptr*npts*nlay),pzz4(nel,nptr*npts*nlay),pzz5(nel,nptr*npts*nlay),
560 . pzz6(nel,nptr*npts*nlay),pzz7(nel,nptr*npts*nlay),pzz8(nel,nptr*npts*nlay))
561c---------------------------
562 CALL s8zderic3(
563 1 offl, volg, ngl, xd1,
564 2 xd2, xd3, xd4, xd5,
565 3 xd6, xd7, xd8, yd1,
566 4 yd2, yd3, yd4, yd5,
567 5 yd6, yd7, yd8, zd1,
568 6 zd2, zd3, zd4, zd5,
569 7 zd6, zd7, zd8, pxc1,
570 8 pxc2, pxc3, pxc4, pyc1,
571 9 pyc2, pyc3, pyc4, pzc1,
572 a pzc2, pzc3, pzc4, px1h1,
573 b px1h2, px1h3, px1h4, px2h1,
574 c px2h2, px2h3, px2h4, px3h1,
575 d px3h2, px3h3, px3h4, px4h1,
576 e px4h2, px4h3, px4h4, hx,
577 f hy, hz, ajc1, ajc2,
578 g ajc3, ajc4, ajc5, ajc6,
579 h ajc7, ajc8, ajc9, smax,
580 i gbuf%SMSTR,gbuf%OFF, nel, ismstr,
581 j jlag)
582 IF (idts6>0) THEN
583 CALL sdlen_dege(
584 1 volg, deltax, x1, x2,
585 2 x3, x4, x5, x6,
586 3 x7, x8, y1, y2,
587 4 y3, y4, y5, y6,
588 5 y7, y8, z1, z2,
589 6 z3, z4, z5, z6,
590 7 z7, z8, ixs(1,nf1),ideg,
591 8 nel)
592 END IF
593 IF (ntsheg > 0 .AND. isctl == 0) THEN
594 CALL sdlensh8(voln,llsh,area ,
595 . x1, x2, x3, x4, x5, x6, x7, x8,
596 . y1, y2, y3, y4, y5, y6, y7, y8,
597 . z1, z2, z3, z4, z5, z6, z7, z8,icsig,nel)
598 llsmin(lft:llt) = deltax(lft:llt)
599 END IF
600C --------------------------
601C --- UPDATE REF CONFIGURATION (possible future change to small strain option)
602C --- ! Total strain option doesn't change the Ref CONF.
603C --------------------------
604 IF (ismstr <= 3.OR.(ismstr==4.AND.jlag>0)) THEN
605 CALL s8sav3(
606 1 gbuf%OFF, gbuf%SMSTR,xd1, xd2,
607 2 xd3, xd4, xd5, xd6,
608 3 xd7, xd8, yd1, yd2,
609 4 yd3, yd4, yd5, yd6,
610 5 yd7, yd8, zd1, zd2,
611 6 zd3, zd4, zd5, zd6,
612 7 zd7, zd8, nel)
613 END IF !(ISMSTR == 2) THEN
614C
615 IF (icp == 2) THEN
616 CALL s8csigp3(gbuf%SIG,e0 ,gbuf%PLA,fac,gbuf%G_PLA,nel)
617 DO i=lft,llt
618 nu1(i)=nu(i)+(half-nu(i))*fac(i)
619 ENDDO
620 ELSE
621 DO i=lft,llt
622 nu1(i)=nu(i)
623 ENDDO
624 ENDIF
625C-----don't do it w/ degenerated elm
626 IF (icp >0 .AND. idts6==0) CALL degenes8(
627 1 ixs(1,nf1),ideg, nel)
628C-----------
629C INITIALIZATION---Before increment-----
630C-----------
631 CALL s8zzero3(
632 1 f11, f21, f31, f12,
633 2 f22, f32, f13, f23,
634 3 f33, f14, f24, f34,
635 4 f15, f25, f35, f16,
636 5 f26, f36, f17, f27,
637 6 f37, f18, f28, f38,
638 7 gbuf%SIG, gbuf%EINT, gbuf%RHO, gbuf%QVIS,
639 8 gbuf%PLA, gbuf%EPSD, stin, pp,
640 9 gbuf%G_PLA, gbuf%G_EPSD,iexpan, gbuf%EINTTH,
641 a nel, conden)
642C-------------------------------------------
643C COMPUTE AVERAGE TEMPERATURE IN ELEMENT
644C-------------------------------------------
645 IF(jthe < 0) THEN
646 DO i=lft,llt
647 IF(gbuf%OFF(i) == zero) cycle
648 tempel(i) = one_over_8 *( temp(nc1(i)) + temp(nc2(i))
649 . + temp(nc3(i)) + temp(nc4(i))
650 . + temp(nc5(i)) + temp(nc6(i))
651 . + temp(nc7(i)) + temp(nc8(i)))
652 gbuf%TEMP(i) = tempel(i)
653 ENDDO
654 ENDIF
655C
656C-----------Constant SIG traitement-----
657C-----------dans ce cas unique constant direction- EV(NB2) est dans ortho. sys------
658 pid=ngeo(lft)
659 DO ir=1,nptr
660 DO is=1,npts
661 DO i=lft,llt
662 sigm(i) = zero
663 volm(i) = zero
664 ENDDO
665 DO il=1,nlay
666 lbuf => elbuf_tab(ng)%BUFLY(il)%LBUF(ir,is,it)
667 DO i=lft,llt
668 sigm(i) = sigm(i)+lbuf%SIG(ii(3)+i)*lbuf%VOL(i)
669 volm(i) = volm(i)+lbuf%VOL(i)
670 ENDDO
671 ENDDO
672 DO il=1,nlay
673 lbuf => elbuf_tab(ng)%BUFLY(il)%LBUF(ir,is,it)
674 ip = ir + ( (is-1) + (il-1)*npts )*nptr
675 DO i=lft,llt
676 smtt(i,ip)=(lbuf%SIG(ii(3)+i)-sigm(i)/volm(i))*usb(i)
677 ENDDO
678 ENDDO
679 ENDDO
680 ENDDO
681c------------------
682C
683 IF (igtyp == 21) THEN
684 CALL sctorth3(
685 1 lft, llt, icsig, nel,
686 2 rx, ry, rz, sx,
687 3 sy, sz, tx, ty,
688 4 tz, r11, r21, r31,
689 5 r12, r22, r32, r13,
690 6 r23, r33, g1x, g1y,
691 7 g1z, g2x, g2y, g2z,
692 8 g3x, g3y, g3z, gbuf%GAMA,
693 9 irep)
694c . G1X,G1Y,G1Z,G2X,G2Y,G2Z,G3X,G3Y,G3Z,EV(NBGAMA))
695 ELSEIF (igtyp == 22) THEN
696 ipthk = ipang+nlymax
697 ippos = ipthk+nlymax
698 mtn0=mtn
699 DO i=lft,llt
700 mxt0(i)=mxt(i)
701 shf(i)=geo(38,ngeo(i))
702 ENDDO
703 ENDIF
704C---------------------------------
705 IF ((ict+icr+ics) >= 1) THEN
706 IF (dt1 == zero) THEN
707 dti =zero
708 ELSE
709 dti = one/dt1
710 ENDIF
711 ENDIF
712C
713 ioffs=0
714 DO i=lft,llt
715 offs(i)=ep20
716C
717 fac(i) = one
718 END DO
719 IF(jthe < 0) them(lft:llt,1:8) = zero
720C-----------Matrix B----------
721C-----------Begin integrating points-----
722c
723 DO ir=1,nptr
724 DO is=1,npts
725 DO il=1,nlay
726 ip = ir + ( (is-1) + (il-1)*npts )*nptr
727 lbuf => elbuf_tab(ng)%BUFLY(il)%LBUF(ir,is,it)
728C-----avoid multi-print
729 IF (ioffs == 1)THEN
730 DO i=lft,llt
731 IF (offs(i)<=two) lbuf%OFF(i)=offs(i)
732 ENDDO
733 END IF
734C-----------
735c IF (IFAILURE > 0) ILAY = IP
736C----------------
737 IF (igtyp == 22) THEN
738 wt = geo(ipthk+il,pid)
739 zz = geo(ippos+il,pid)
740 mid=igeo(ipmat+il,pid)
741 mtn=nint(pm(19,mid))
742 DO i=lft,llt
743 mxt(i)=mid
744 ENDDO
745 ELSE
746 zz = a_gauss(il,nlay)
747 wt = w_gauss(il,nlay)
748 ENDIF
749 IF (ict == 1) THEN
750 zr = a_gauss(ir,nptr)
751 zs = a_gauss(is,npts)
752 zt = zz
753 ELSEIF (ics == 1) THEN
754 zr = a_gauss(ir,nptr)
755 zs = zz
756 zt = a_gauss(is,npts)
757 ELSEIF (icr == 1) THEN
758 zr = zz
759 zs = a_gauss(ir,nptr)
760 zt = a_gauss(is,npts)
761 ENDIF
762 wi = w_gauss(ir,nptr)*w_gauss(is,npts)*wt
763c
764c-------------------------------------------
765c COMPUTE Regularized non local variable in Gauss point
766c-------------------------------------------
767 IF (inloc > 0) THEN
768 l_nloc = nloc_dmg%L_NLOC
769 dnl => nloc_dmg%DNL(1:l_nloc) ! DNL = non local variable increment
770 CALL basis8 (zr,zs,zt,h,pr,ps,pt)
771 DO i=lft,llt
772 inod(1) = nloc_dmg%IDXI(nc1(i))
773 inod(2) = nloc_dmg%IDXI(nc2(i))
774 inod(3) = nloc_dmg%IDXI(nc3(i))
775 inod(4) = nloc_dmg%IDXI(nc4(i))
776 inod(5) = nloc_dmg%IDXI(nc5(i))
777 inod(6) = nloc_dmg%IDXI(nc6(i))
778 inod(7) = nloc_dmg%IDXI(nc7(i))
779 inod(8) = nloc_dmg%IDXI(nc8(i))
780 DO j = 1, 8
781 ipos(j) = nloc_dmg%POSI(inod(j))+il-1
782 ENDDO
783 var_reg(i,ip) = h(1)*dnl(ipos(1)) + h(2)*dnl(ipos(2)) + h(3)*dnl(ipos(3))
784 . + h(4)*dnl(ipos(4)) + h(5)*dnl(ipos(5)) + h(6)*dnl(ipos(6))
785 . + h(7)*dnl(ipos(7)) + h(8)*dnl(ipos(8))
786 ENDDO
787 ENDIF
788c
789 CALL s8zderi3(
790 1 offl, off, voln, ngl,
791 2 zr, zs, zt, wi,
792 3 pxc1, pxc2, pxc3, pxc4,
793 4 pyc1, pyc2, pyc3, pyc4,
794 5 pzc1, pzc2, pzc3, pzc4,
795 6 px1h1, px1h2, px1h3, px1h4,
796 7 px2h1, px2h2, px2h3, px2h4,
797 8 px3h1, px3h2, px3h3, px3h4,
798 9 px4h1, px4h2, px4h3, px4h4,
799 a hx, hy, hz, px1,
800 b px2, px3, px4, px5,
801 c px6, px7, px8, py1,
802 d py2, py3, py4, py5,
803 e py6, py7, py8, pz1,
804 f pz2, pz3, pz4, pz5,
805 g pz6, pz7, pz8, pxy1,
806 h pxy2, pxy3, pxy4, pxy5,
807 i pxy6, pxy7, pxy8, pyx1,
808 j pyx2, pyx3, pyx4, pyx5,
809 k pyx6, pyx7, pyx8, pxz1,
810 l pxz2, pxz3, pxz4, pxz5,
811 m pxz6, pxz7, pxz8, pzx1,
812 n pzx2, pzx3, pzx4, pzx5,
813 o pzx6, pzx7, pzx8, pyz1,
814 p pyz2, pyz3, pyz4, pyz5,
815 q pyz6, pyz7, pyz8, pzy1,
816 r pzy2, pzy3, pzy4, pzy5,
817 s pzy6, pzy7, pzy8, bxy1,
818 t bxy2, bxy3, bxy4, bxy5,
819 u bxy6, bxy7, bxy8, byx1,
820 v byx2, byx3, byx4, byx5,
821 w byx6, byx7, byx8, bxz1,
822 x bxz2, bxz3, bxz4, bxz5,
823 y bxz6, bxz7, bxz8, bzx1,
824 z bzx2, bzx3, bzx4, bzx5,
825 1 bzx6, bzx7, bzx8, byz1,
826 2 byz2, byz3, byz4, byz5,
827 3 byz6, byz7, byz8, bzy1,
828 4 bzy2, bzy3, bzy4, bzy5,
829 5 bzy6, bzy7, bzy8, ajc1,
830 6 ajc2, ajc3, ajc4, ajc5,
831 7 ajc6, ajc7, ajc8, ajc9,
832 8 aj1, aj2, aj3, aj4,
833 9 aj5, aj6, smax, deltax,
834 a nu1, icp, ideg, voldp,
835 b nel, mtn, ismstr, jhbe)
836C
837 IF (ict == 1) THEN
838 smt(lft:llt,ip) = smtt(lft:llt,ip)
839 ELSEIF (ics == 1) THEN
840 sms(lft:llt,ip) = smtt(lft:llt,ip)
841 ELSEIF (icr == 1) THEN
842 smr(lft:llt,ip) = smtt(lft:llt,ip)
843 ENDIF
844 CALL s8cdefo3(
845 1 pxc1, pxc2, pxc3, pxc4,
846 2 pyc1, pyc2, pyc3, pyc4,
847 3 pzc1, pzc2, pzc3, pzc4,
848 4 px1, px2, px3, px4,
849 5 px5, px6, px7, px8,
850 6 py1, py2, py3, py4,
851 7 py5, py6, py7, py8,
852 8 pz1, pz2, pz3, pz4,
853 9 pz5, pz6, pz7, pz8,
854 a pxy1, pxy2, pxy3, pxy4,
855 b pxy5, pxy6, pxy7, pxy8,
856 c pyx1, pyx2, pyx3, pyx4,
857 d pyx5, pyx6, pyx7, pyx8,
858 e pxz1, pxz2, pxz3, pxz4,
859 f pxz5, pxz6, pxz7, pxz8,
860 g pzx1, pzx2, pzx3, pzx4,
861 h pzx5, pzx6, pzx7, pzx8,
862 i pyz1, pyz2, pyz3, pyz4,
863 j pyz5, pyz6, pyz7, pyz8,
864 k pzy1, pzy2, pzy3, pzy4,
865 l pzy5, pzy6, pzy7, pzy8,
866 m bxy1, bxy2, bxy3, bxy4,
867 n bxy5, bxy6, bxy7, bxy8,
868 o byx1, byx2, byx3, byx4,
869 p byx5, byx6, byx7, byx8,
870 q bxz1, bxz2, bxz3, bxz4,
871 r bxz5, bxz6, bxz7, bxz8,
872 s bzx1, bzx2, bzx3, bzx4,
873 t bzx5, bzx6, bzx7, bzx8,
874 u byz1, byz2, byz3, byz4,
875 v byz5, byz6, byz7, byz8,
876 w bzy1, bzy2, bzy3, bzy4,
877 x bzy5, bzy6, bzy7, bzy8,
878 y vx1, vx2, vx3, vx4,
879 z vx5, vx6, vx7, vx8,
880 1 vy1, vy2, vy3, vy4,
881 2 vy5, vy6, vy7, vy8,
882 3 vz1, vz2, vz3, vz4,
883 4 vz5, vz6, vz7, vz8,
884 5 dxx, dxy, dxz, dyx,
885 6 dyy, dyz, dzx, dzy,
886 7 dzz, d4, d5, d6,
887 8 wxx, wyy, wzz, lbuf%VOL,
888 9 off, lbuf%EINT, gbuf%OFF, dsv,
889 a icp, icr, ics, ict,
890 b dti, smr(1,ip), sms(1,ip), smt(1,ip),
891 c ngl, ideg, lbuf%VOL0DP,nel,
892 d ismstr, ipres)
893C
894 mx = mxt(lft)
895 DO i=lft,llt
896 rhoo(i) = lbuf%RHO(i)
897 tempel(i) = pm(79,mx)
898 ENDDO
899 IF (igtyp == 22)
900 . CALL sctorth3(
901 1 lft, llt, icsig, nel,
902 2 rx, ry, rz, sx,
903 3 sy, sz, tx, ty,
904 4 tz, r11, r21, r31,
905 5 r12, r22, r32, r13,
906 6 r23, r33, g1x, g1y,
907 7 g1z, g2x, g2y, g2z,
908 8 g3x, g3y, g3z, lbuf%GAMA,
909 9 irep)
910 IF (jcvt == 2) THEN
911 CALL szordef3(lft,llt,dxx,dyy,dzz,d4,d5,d6,
912 . g1x, g1y, g1z, g2x, g2y, g2z, g3x, g3y, g3z)
913 ENDIF
914 IF (igtyp == 22) THEN
915 IF (ics == 1) THEN
916 DO i=lft,llt
917 d4(i)=shf(i)*d4(i)
918 d5(i)=shf(i)*d5(i)
919 ENDDO
920 ELSEIF (icr == 1) THEN
921 DO i=lft,llt
922 d4(i)=shf(i)*d4(i)
923 d6(i)=shf(i)*d6(i)
924 ENDDO
925 ELSEIF (ict == 1) THEN
926 DO i=lft,llt
927 d5(i)=shf(i)*d5(i)
928 d6(i)=shf(i)*d6(i)
929 ENDDO
930 ENDIF
931 ENDIF
932C
933 IF (ntsheg > 0 .AND. isctl == 0) THEN
934 DO i=lft,llt
935 IF (gbuf%IDT_TSH(i)<=0) cycle
936 llsmin(i) = min(llsmin(i),deltax(i))
937 deltax(i) = max(llsh(i),deltax(i))
938 ENDDO
939 END IF
940 divde(1:nel) = dt1*(dxx(1:nel)+ dyy(1:nel)+ dzz(1:nel))+dsv(1:nel)
941 CALL srho3(
942 1 pm, lbuf%VOL, lbuf%RHO, lbuf%EINT,
943 2 divde, flux(1,nf1),flu1(nf1), voln,
944 3 dvol, ngl, mxt, off,
945 4 0, gbuf%TAG22, voldp, lbuf%VOL0DP,
946 5 amu, gbuf%OFF, nel, mtn,
947 6 jale, ismstr, jeul, jlag)
948C
949 IF (jcvt == 0) THEN
950 CALL srota3(
951 1 lbuf%SIG,s1, s2, s3,
952 2 s4, s5, s6, wxx,
953 3 wyy, wzz, nel, mtn,
954 4 ismstr)
955C-----------------------------
956C SMALL STRAIN
957C-----------------------------
958 CALL smalla3(
959 1 gbuf%SMSTR,gbuf%OFF, off, wxx,
960 2 wyy, wzz, nel, ismstr,
961 3 jlag)
962C
963 ELSE
964C-----------------------------
965C EXTRACT STRESSES + SMALL STRAIN
966C-----------------------------
967 CALL csmall3(lbuf%SIG,s1,s2,s3,s4,s5,s6,
968 . gbuf%OFF,off,nel)
969 ENDIF
970C
971C for heat transfert
972C
973 IF (jthe < 0 ) THEN
974 CALL s8ztempel(lft,llt,ni,
975 1 zr,zs,zt,
976 2 nc1,nc2,nc3,nc4,nc5,nc6,nc7,nc8,
977 3 temp,tempel)
978 ENDIF
979C------------------------------------------------------
980C CALCUL DES CONTRAINTES SUIVANT LOIS CONSTITUTIVES
981C TIME STEP EST CALCULE ICI DT2T
982C------------------------------------------------------
983 IF ((itask==0).AND.(imon_mat==1)) CALL startime(timers,35)
984c
985 IF(iboltp /= 0) CALL boltst(
986 . ip, bpreld, lbuf%SIG,tt ,
987 . nel ,npt ,sensors%NSENSOR,sensors%SENSOR_TAB)
988
989 current_layer=il
990 CALL mmain(timers, output,
991 1 elbuf_tab, ng, pm, geo,
992 2 ale_connect, ixs, iparg,
993 3 v, tf, npf, bufmat,
994 4 sti, x, dt2t, neltst,
995 5 ityptst, offset, nel, w,
996 6 off, ngeo, mxt, ngl,
997 7 voln, vd2, dvol, deltax,
998 8 vis, qvis, cxx, s1,
999 9 s2, s3, s4, s5,
1000 a s6, dxx, dyy, dzz,
1001 b d4, d5, d6, wxx,
1002 c wyy, wzz, aj1, aj2,
1003 d aj3, aj4, aj5, aj6,
1004 e vdx, vdy, vdz, muvoid,
1005 f ssp_eq, aire, sigy, et,
1006 g r1_free, lbuf%PLA, r3_free, amu,
1007 h mfxx(1,ip), mfxy(1,ip), mfxz(1,ip), mfyx(1,ip),
1008 i mfyy(1,ip), mfyz(1,ip), mfzx(1,ip), mfzy(1,ip),
1009 j mfzz(1,ip), ipm, gama, bid,
1010 k bid, bid, bid, bid,
1011 l bid, bid, istrain, tempel,
1012 m die, iexpan, current_layer,mssa,
1013 n dmels, ir, is, it,
1014 o table, bid, bid, bid,
1015 p bid, iparg(1,ng), igeo, conde,
1016 q itask, nloc_dmg, var_reg(1,ip),mat_elem,
1017 r h3d_strain, jplasol, jsph, sz_r1_free,
1018 * snpc, stf, sbufmat, glob_therm,
1019 * svis, sz_ix, iresp,
1020 * n2d, th_strain, ngroup, tt,
1021 . dt1, ntable, numelq, nummat,
1022 . numgeo, numnod, numels,
1023 . idel7nok, idtmin, maxfunc,
1024 . imon_mat, userl_avail, impl_s,
1025 . idyna, dt, fheat ,sensors, opt_mtn=mtn, opt_jcvt=jcvt,
1026 * opt_isorth=isorth,opt_isorthg=isorthg)
1027C
1028 IF (istrain == 1) THEN
1029 CALL sstra3(
1030 1 dxx, dyy, dzz, d4,
1031 2 d5, d6, lbuf%STRA,wxx,
1032 3 wyy, wzz, off, nel,
1033 4 jcvt)
1034 ENDIF
1035C
1036 IF ((itask==0).AND.(imon_mat==1)) CALL stoptime(timers,35)
1037 l_pla = elbuf_tab(ng)%BUFLY(il)%L_PLA
1038 l_epsd = elbuf_tab(ng)%BUFLY(il)%L_EPSD
1039C
1040 IF (jcvt == 2) THEN
1041 CALL sroto3(lft,llt,lbuf%SIG,sign,
1042 . g1x, g2x, g3x, g1y, g2y, g3y, g1z, g2z, g3z,nel)
1043 CALL s8zfint3(
1044 1 sign, sign, px1, px2,
1045 2 px3, px4, py1, py2,
1046 3 py3, py4, pz1, pz2,
1047 4 pz3, pz4, px5, px6,
1048 5 px7, px8, py5, py6,
1049 6 py7, py8, pz5, pz6,
1050 7 pz7, pz8, pxy1, pxy2,
1051 8 pxy3, pxy4, pxy5, pxy6,
1052 9 pxy7, pxy8, pyx1, pyx2,
1053 a pyx3, pyx4, pyx5, pyx6,
1054 b pyx7, pyx8, pxz1, pxz2,
1055 c pxz3, pxz4, pxz5, pxz6,
1056 d pxz7, pxz8, pzx1, pzx2,
1057 e pzx3, pzx4, pzx5, pzx6,
1058 f pzx7, pzx8, pyz1, pyz2,
1059 g pyz3, pyz4, pyz5, pyz6,
1060 h pyz7, pyz8, pzy1, pzy2,
1061 i pzy3, pzy4, pzy5, pzy6,
1062 j pzy7, pzy8, bxy1, bxy2,
1063 k bxy3, bxy4, bxy5, bxy6,
1064 l bxy7, bxy8, byx1, byx2,
1065 m byx3, byx4, byx5, byx6,
1066 n byx7, byx8, bxz1, bxz2,
1067 o bxz3, bxz4, bxz5, bxz6,
1068 p bxz7, bxz8, bzx1, bzx2,
1069 q bzx3, bzx4, bzx5, bzx6,
1070 r bzx7, bzx8, byz1, byz2,
1071 s byz3, byz4, byz5, byz6,
1072 t byz7, byz8, bzy1, bzy2,
1073 u bzy3, bzy4, bzy5, bzy6,
1074 v bzy7, bzy8, f11, f21,
1075 w f31, f12, f22, f32,
1076 x f13, f23, f33, f14,
1077 y f24, f34, f15, f25,
1078 z f35, f16, f26, f36,
1079 1 f17, f27, f37, f18,
1080 2 f28, f38, voln, qvis,
1081 3 pp, lbuf%EINT, lbuf%RHO, lbuf%QVIS,
1082 4 lbuf%PLA, lbuf%EPSD, gbuf%EPSD, gbuf%SIG,
1083 5 gbuf%EINT, gbuf%RHO, gbuf%QVIS, gbuf%PLA,
1084 6 volg, sti, stin, icp,
1085 7 off, lbuf%VOL, gbuf%VOL, gbuf%G_PLA*l_pla,
1086 8 l_epsd, fac, lbuf%EINTTH, gbuf%EINTTH,
1087 9 iexpan, nel, ideg, conde,
1088 a conden, mtn, ismstr, svis ,glob_therm%NODADT_THERM,
1089 b gbuf%WPLA, lbuf%WPLA, gbuf%G_WPLA )
1090 ELSE
1091C----------------------------
1092C INTERNAL FORCES
1093C----------------------------
1094 CALL s8zfint3(
1095 1 lbuf%SIG, lbuf%SIG, px1, px2,
1096 2 px3, px4, py1, py2,
1097 3 py3, py4, pz1, pz2,
1098 4 pz3, pz4, px5, px6,
1099 5 px7, px8, py5, py6,
1100 6 py7, py8, pz5, pz6,
1101 7 pz7, pz8, pxy1, pxy2,
1102 8 pxy3, pxy4, pxy5, pxy6,
1103 9 pxy7, pxy8, pyx1, pyx2,
1104 a pyx3, pyx4, pyx5, pyx6,
1105 b pyx7, pyx8, pxz1, pxz2,
1106 c pxz3, pxz4, pxz5, pxz6,
1107 d pxz7, pxz8, pzx1, pzx2,
1108 e pzx3, pzx4, pzx5, pzx6,
1109 f pzx7, pzx8, pyz1, pyz2,
1110 g pyz3, pyz4, pyz5, pyz6,
1111 h pyz7, pyz8, pzy1, pzy2,
1112 i pzy3, pzy4, pzy5, pzy6,
1113 j pzy7, pzy8, bxy1, bxy2,
1114 k bxy3, bxy4, bxy5, bxy6,
1115 l bxy7, bxy8, byx1, byx2,
1116 m byx3, byx4, byx5, byx6,
1117 n byx7, byx8, bxz1, bxz2,
1118 o bxz3, bxz4, bxz5, bxz6,
1119 p bxz7, bxz8, bzx1, bzx2,
1120 q bzx3, bzx4, bzx5, bzx6,
1121 r bzx7, bzx8, byz1, byz2,
1122 s byz3, byz4, byz5, byz6,
1123 t byz7, byz8, bzy1, bzy2,
1124 u bzy3, bzy4, bzy5, bzy6,
1125 v bzy7, bzy8, f11, f21,
1126 w f31, f12, f22, f32,
1127 x f13, f23, f33, f14,
1128 y f24, f34, f15, f25,
1129 z f35, f16, f26, f36,
1130 1 f17, f27, f37, f18,
1131 2 f28, f38, voln, qvis,
1132 3 pp, lbuf%EINT, lbuf%RHO, lbuf%QVIS,
1133 4 lbuf%PLA, lbuf%EPSD, gbuf%EPSD, gbuf%SIG,
1134 5 gbuf%EINT, gbuf%RHO, gbuf%QVIS, gbuf%PLA,
1135 6 volg, sti, stin, icp,
1136 7 off, lbuf%VOL, gbuf%VOL, l_pla,
1137 8 l_epsd, fac, lbuf%EINTTH,gbuf%EINTTH,
1138 9 iexpan, nel, ideg, conde,
1139 a conden, mtn, ismstr, svis,glob_therm%NODADT_THERM,
1140 b gbuf%WPLA, lbuf%WPLA, gbuf%G_WPLA )
1141 ENDIF
1142 DO i=lft,llt
1143 offl(i)=min(offl(i),off(i))
1144 IF (lbuf%OFF(i) > one .AND. gbuf%OFF(i) == one) THEN
1145 offs(i)=min(lbuf%OFF(i),offs(i))
1146 ioffs =1
1147 END IF
1148 ENDDO
1149CC-------------------------
1150c finite element heat transfert
1151C--------------------------
1152 IF (jthe < 0) THEN
1153 IF (mat_elem%MAT_PARAM(mx)%HEAT_FLAG == 1) THEN
1154 CALL s8etherm(
1155 1 pm, mx , voln, ni,
1156 2 nc1, nc2, nc3, nc4,
1157 3 nc5, nc6, nc7, nc8,
1158 4 px1, px2, px3, px4,
1159 5 py1, py2, py3, py4,
1160 6 pz1, pz2, pz3, pz4,
1161 7 px5, px6, px7, px8,
1162 8 py5, py6, py7, py8,
1163 9 pz5, pz6, pz7, pz8,
1164 a dt1, temp, tempel, fheat,
1165 b them, gbuf%OFF,lbuf%OFF,partsav,
1166 c iparts, lbuf%VOL,nel,glob_therm%THEACCFACT)
1167 ELSE
1168 CALL s8etherm(
1169 1 pm, mx , voln, ni,
1170 2 nc1, nc2, nc3, nc4,
1171 3 nc5, nc6, nc7, nc8,
1172 4 px1, px2, px3, px4,
1173 5 py1, py2, py3, py4,
1174 6 pz1, pz2, pz3, pz4,
1175 7 px5, px6, px7, px8,
1176 8 py5, py6, py7, py8,
1177 9 pz5, pz6, pz7, pz8,
1178 a dt1, temp, tempel, die,
1179 b them, gbuf%OFF,lbuf%OFF,partsav,
1180 c iparts, lbuf%VOL,nel,glob_therm%THEACCFACT)
1181 END IF
1182 ENDIF
1183c
1184C-------------------------
1185c Save shape function derivatives value for non-local
1186C--------------------------
1187 IF (inloc > 0) THEN
1188 pxx1(1:nel,ip) = px1(1:nel)
1189 pxx2(1:nel,ip) = px2(1:nel)
1190 pxx3(1:nel,ip) = px3(1:nel)
1191 pxx4(1:nel,ip) = px4(1:nel)
1192 pxx5(1:nel,ip) = px5(1:nel)
1193 pxx6(1:nel,ip) = px6(1:nel)
1194 pxx7(1:nel,ip) = px7(1:nel)
1195 pxx8(1:nel,ip) = px8(1:nel)
1196 pyy1(1:nel,ip) = py1(1:nel)
1197 pyy2(1:nel,ip) = py2(1:nel)
1198 pyy3(1:nel,ip) = py3(1:nel)
1199 pyy4(1:nel,ip) = py4(1:nel)
1200 pyy5(1:nel,ip) = py5(1:nel)
1201 pyy6(1:nel,ip) = py6(1:nel)
1202 pyy7(1:nel,ip) = py7(1:nel)
1203 pyy8(1:nel,ip) = py8(1:nel)
1204 pzz1(1:nel,ip) = pz1(1:nel)
1205 pzz2(1:nel,ip) = pz2(1:nel)
1206 pzz3(1:nel,ip) = pz3(1:nel)
1207 pzz4(1:nel,ip) = pz4(1:nel)
1208 pzz5(1:nel,ip) = pz5(1:nel)
1209 pzz6(1:nel,ip) = pz6(1:nel)
1210 pzz7(1:nel,ip) = pz7(1:nel)
1211 pzz8(1:nel,ip) = pz8(1:nel)
1212 ENDIF
1213 ENDDO ! IL=1,NLAY
1214 ENDDO ! IS=1,NPTS
1215 ENDDO ! IR=1,NPTR
1216C-----------End integrating points-----
1217c
1218c-------------------------
1219c Virtual internal forces of regularized non local ddl
1220c--------------------------
1221 IF (inloc > 0) THEN
1222 CALL sdlensh8(volg,llsh,area ,
1223 . x1, x2, x3, x4, x5, x6, x7, x8,
1224 . y1, y2, y3, y4, y5, y6, y7, y8,
1225 . z1, z2, z3, z4, z5, z6, z7, z8,icsig,nel)
1226 DO ir = 1,nptr
1227 DO is = 1,npts
1228 CALL s8cfint_reg(
1229 1 nloc_dmg, var_reg, nel, gbuf%OFF,
1230 2 volg, nc1, nc2, nc3,
1231 3 nc4, nc5, nc6, nc7,
1232 4 nc8, pxx1, pxx2, pxx3,
1233 5 pxx4, pxx5, pxx6, pxx7,
1234 6 pxx8, pyy1, pyy2, pyy3,
1235 7 pyy4, pyy5, pyy6, pyy7,
1236 8 pyy8, pzz1, pzz2, pzz3,
1237 9 pzz4, pzz5, pzz6, pzz7,
1238 a pzz8, mxt(lft), itask, dt2t,
1239 b gbuf%VOL, nft, nlay, nptr,
1240 c npts, ir, is, w_gauss,
1241 d a_gauss, icr, ics, ict,
1242 e elbuf_tab(ng)%NLOCTS(ir,is),area, nodadt,
1243 f dt1, dt12, dtfac1, dtmin1,
1244 g idtmin)
1245 ENDDO
1246 ENDDO
1247 ENDIF
1248c
1249 IF (ioffs == 1) THEN
1250 DO i=lft,llt
1251 IF (offs(i) <= two) gbuf%OFF(i)=offs(i)
1252 END DO
1253c
1254 DO ir=1,nptr
1255 DO is=1,npts
1256 DO il=1,nlay
1257 lbuf => elbuf_tab(ng)%BUFLY(il)%LBUF(ir,is,it)
1258 DO i=lft,llt
1259 IF (gbuf%OFF(i) > one) lbuf%OFF(i)=gbuf%OFF(i)
1260 END DO
1261 END DO
1262 END DO
1263 END DO
1264 END IF
1265c----------------------
1266 IF (icp == 1) THEN
1267 CALL s8zfintp3(
1268 1 pxc1, pxc2, pxc3, pxc4,
1269 2 pyc1, pyc2, pyc3, pyc4,
1270 3 pzc1, pzc2, pzc3, pzc4,
1271 4 f11, f21, f31, f12,
1272 5 f22, f32, f13, f23,
1273 6 f33, f14, f24, f34,
1274 7 f15, f25, f35, f16,
1275 8 f26, f36, f17, f27,
1276 9 f37, f18, f28, f38,
1277 a volg, pp, ideg, nel)
1278 ENDIF
1279 IF (igtyp == 22) THEN
1280 mtn = mtn0
1281 DO i=lft,llt
1282 mxt(i)=mxt0(i)
1283 ENDDO
1284 ENDIF
1285 IF ( nn_del> 0) THEN
1286 CALL sdlensh14(volg,llsh,area ,
1287 . x1, x2, x3, x4, x5, x6, x7, x8,
1288 . y1, y2, y3, y4, y5, y6, y7, y8,
1289 . z1, z2, z3, z4, z5, z6, z7, z8,icsig,nel)
1290 CALL tshgeodel3(ngl,gbuf%OFF,volg,area,gbuf%VOL,
1291 . llsh,geo(1,pid),nn_del,dt,nel )
1292 ENDIF
1293C-----------------------------
1294C SMALL STRAIN
1295C-----------------------------
1296 CALL smallb3(
1297 1 gbuf%OFF,offl, nel, ismstr)
1298C--------------------------
1299C BILANS PAR MATERIAU
1300C--------------------------
1301 iflag=mod(ncycle,ncpri)
1302 IF(ioutprt>0)THEN
1303 IF (jcvt == 0) THEN
1304 ifvm22 = 0
1305 CALL sbilan(partsav,gbuf%EINT,gbuf%RHO,gbuf%RK,gbuf%VOL,
1306 . vx1, vx2, vx3, vx4, vx5, vx6, vx7, vx8,
1307 . vy1, vy2, vy3, vy4, vy5, vy6, vy7, vy8,
1308 . vz1, vz2, vz3, vz4, vz5, vz6, vz7, vz8,
1309 . volg,iparts , gresav,grth,igrth,gbuf%OFF,
1310 . iexpan,gbuf%EINTTH,gbuf%FILL,gbuf%MOM,nel,ifvm22,
1311 . x1, x2, x3, x4, x5, x6, x7, x8,
1312 . y1, y2, y3, y4, y5, y6, y7, y8,
1313 . z1, z2, z3, z4, z5, z6, z7, z8,
1314 . itask,iparg(1,ng),sensors,gbuf%G_WPLA,gbuf%WPLA)
1315 ELSE
1316 CALL srbilan(partsav,gbuf%EINT,gbuf%RHO,gbuf%RK,gbuf%VOL,
1317 . vgxa, vgya, vgza, vga2, volg,iparts,
1318 . gresav,grth,igrth,gbuf%OFF,iexpan,gbuf%EINTTH,
1319 . gbuf%FILL,xgxa, xgya, xgza,
1320 . xgxa2,xgya2,xgza2,xgxya,xgyza,xgzxa,itask,iparg(1,ng),sensors,
1321 . nel,gbuf%G_WPLA,gbuf%WPLA)
1322 ENDIF
1323 ENDIF
1324C----------------------------
1325C CONVECTE --> GLOBAL.
1326C----------------------------
1327 IF (jcvt/=0) THEN
1328 CALL srrota3(
1329 1 r11, r21, r31, r12,
1330 2 r22, r32, r13, r23,
1331 3 r33, f11, f12, f13,
1332 4 f14, f15, f16, f17,
1333 5 f18, f21, f22, f23,
1334 6 f24, f25, f26, f27,
1335 7 f28, f31, f32, f33,
1336 8 f34, f35, f36, f37,
1337 9 f38, nel)
1338 ENDIF
1339C----------------------------
1340C distortion control
1341C----------------------------
1342 IF (isctl > 0) THEN
1343 alpha_e(1:nel) = one ! incompatibale
1344 CALL sdistor_ini(
1345 1 nel ,sti_c ,npropm ,nummat ,
1346 2 ismstr ,mxt ,istab ,pm ,
1347 3 gbuf%SIG ,gbuf%RHO ,cxx ,offl ,
1348 4 gbuf%OFF ,ll ,volg ,fld ,
1349 5 cns2 ,fqmax )
1350! all in global system
1351 CALL s8get_x3(
1352 . nc1, nc2, nc3, nc4,
1353 . nc5, nc6, nc7, nc8,
1354 . x1, x2, x3, x4,
1355 . x5, x6, x7, x8,
1356 . y1, y2, y3, y4,
1357 . y5, y6, y7, y8,
1358 . z1, z2, z3, z4,
1359 . z5, z6, z7, z8,
1360 . x, xdp, numnod, nel )
1361 DO i=1,nel
1362 vx1(i)=v(1,nc1(i))
1363 vy1(i)=v(2,nc1(i))
1364 vz1(i)=v(3,nc1(i))
1365 vx2(i)=v(1,nc2(i))
1366 vy2(i)=v(2,nc2(i))
1367 vz2(i)=v(3,nc2(i))
1368 vx3(i)=v(1,nc3(i))
1369 vy3(i)=v(2,nc3(i))
1370 vz3(i)=v(3,nc3(i))
1371 vx4(i)=v(1,nc4(i))
1372 vy4(i)=v(2,nc4(i))
1373 vz4(i)=v(3,nc4(i))
1374 vx5(i)=v(1,nc5(i))
1375 vy5(i)=v(2,nc5(i))
1376 vz5(i)=v(3,nc5(i))
1377 vx6(i)=v(1,nc6(i))
1378 vy6(i)=v(2,nc6(i))
1379 vz6(i)=v(3,nc6(i))
1380 vx7(i)=v(1,nc7(i))
1381 vy7(i)=v(2,nc7(i))
1382 vz7(i)=v(3,nc7(i))
1383 vx8(i)=v(1,nc8(i))
1384 vy8(i)=v(2,nc8(i))
1385 vz8(i)=v(3,nc8(i))
1386 ENDDO
1387 CALL s8for_distor(
1388 . x1, x2, x3, x4,
1389 . x5, x6, x7, x8,
1390 . y1, y2, y3, y4,
1391 . y5, y6, y7, y8,
1392 . z1, z2, z3, z4,
1393 . z5, z6, z7, z8,
1394 . vx1, vx2, vx3, vx4,
1395 . vx5, vx6, vx7, vx8,
1396 . vy1, vy2, vy3, vy4,
1397 . vy5, vy6, vy7, vy8,
1398 . vz1, vz2, vz3, vz4,
1399 . vz5, vz6, vz7, vz8,
1400 . f11, f12, f13, f14,
1401 . f15, f16, f17, f18,
1402 . f21, f22, f23, f24,
1403 . f25, f26, f27, f28,
1404 . f31, f32, f33, f34,
1405 . f35, f36, f37, f38,
1406 . stin, sti_c, fld , cns2,
1407 . ll , istab, fqmax, nel,
1408 . gbuf%EINT_DISTOR, dt1)
1409 ENDIF
1410C----------------------------
1411 IF(nfilsol/=0) CALL sfillopt(
1412 1 gbuf%FILL,stin, f11, f21,
1413 2 f31, f12, f22, f32,
1414 3 f13, f23, f33, f14,
1415 4 f24, f34, f15, f25,
1416 5 f35, f16, f26, f36,
1417 6 f17, f27, f37, f18,
1418 7 f28, f38, nel)
1419C----------------------------
1420 IF(iparit == 0)THEN
1421 CALL scumu3(
1422 1 gbuf%OFF,a, nc1, nc2,
1423 2 nc3, nc4, nc5, nc6,
1424 3 nc7, nc8, stifn, stin,
1425 4 f11, f21, f31, f12,
1426 5 f22, f32, f13, f23,
1427 6 f33, f14, f24, f34,
1428 7 f15, f25, f35, f16,
1429 8 f26, f36, f17, f27,
1430 9 f37, f18, f28, f38,
1431 a nvc, bid, bid, bid,
1432 b bid, bid, bid, bid,
1433 c bid, bid, bid, bid,
1434 d bid, bid, bid, bid,
1435 e bid, bid, bid, bid,
1436 f bid, bid, bid, bid,
1437 g bid, bid, bid, bid,
1438 h them, fthe, condn, conden,
1439 i nel, jthe, isrot, ipartsph,glob_therm%NODADT_THERM)
1440 ELSE
1441 CALL scumu3p(
1442 1 gbuf%OFF,stin, fsky, fsky,
1443 2 iads, f11, f21, f31,
1444 3 f12, f22, f32, f13,
1445 4 f23, f33, f14, f24,
1446 5 f34, f15, f25, f35,
1447 6 f16, f26, f36, f17,
1448 7 f27, f37, f18, f28,
1449 8 f38, nc1, nc2, nc3,
1450 9 nc4, nc5, nc6, nc7,
1451 a nc8, bid, bid, bid,
1452 b bid, bid, bid, bid,
1453 c bid, bid, bid, bid,
1454 d bid, bid, bid, bid,
1455 e bid, bid, bid, bid,
1456 f bid, bid, bid, bid,
1457 g bid, bid, bid, bid,
1458 h them, fthesky, condnsky,conden,
1459 i nel, nft, jthe, isrot,
1460 j ipartsph,glob_therm%NODADT_THERM)
1461 ENDIF
1462 IF (ntsheg > 0) THEN
1463 alpha_e(lft:llt) = one
1464 DO i=lft,llt
1465 IF (gbuf%IDT_TSH(i)<=0) cycle
1466 facdp = 1.343*llsh(i)/llsmin(i)
1467 alpha_e(i) = facdp*facdp
1468 ENDDO
1469 CALL scumualpha(
1470 1 gbuf%OFF,alpha_e, nc1, nc2,
1471 2 nc3, nc4, nc5, nc6,
1472 3 nc7, nc8, nel)
1473 END IF
1474C-----------
1475
1476 IF (ALLOCATED(var_reg)) DEALLOCATE(var_reg)
1477 IF (ALLOCATED(pxx1)) DEALLOCATE(pxx1)
1478 IF (ALLOCATED(pxx2)) DEALLOCATE(pxx2)
1479 IF (ALLOCATED(pxx3)) DEALLOCATE(pxx3)
1480 IF (ALLOCATED(pxx4)) DEALLOCATE(pxx4)
1481 IF (ALLOCATED(pxx5)) DEALLOCATE(pxx5)
1482 IF (ALLOCATED(pxx6)) DEALLOCATE(pxx6)
1483 IF (ALLOCATED(pxx7)) DEALLOCATE(pxx7)
1484 IF (ALLOCATED(pxx8)) DEALLOCATE(pxx8)
1485 IF (ALLOCATED(pyy1)) DEALLOCATE(pyy1)
1486 IF (ALLOCATED(pyy2)) DEALLOCATE(pyy2)
1487 IF (ALLOCATED(pyy3)) DEALLOCATE(pyy3)
1488 IF (ALLOCATED(pyy4)) DEALLOCATE(pyy4)
1489 IF (ALLOCATED(pyy5)) DEALLOCATE(pyy5)
1490 IF (ALLOCATED(pyy6)) DEALLOCATE(pyy6)
1491 IF (ALLOCATED(pyy7)) DEALLOCATE(pyy7)
1492 IF (ALLOCATED(pyy8)) DEALLOCATE(pyy8)
1493 IF (ALLOCATED(pzz1)) DEALLOCATE(pzz1)
1494 IF (ALLOCATED(pzz2)) DEALLOCATE(pzz2)
1495 IF (ALLOCATED(pzz3)) DEALLOCATE(pzz3)
1496 IF (ALLOCATED(pzz4)) DEALLOCATE(pzz4)
1497 IF (ALLOCATED(pzz5)) DEALLOCATE(pzz5)
1498 IF (ALLOCATED(pzz6)) DEALLOCATE(pzz6)
1499 IF (ALLOCATED(pzz7)) DEALLOCATE(pzz7)
1500 IF (ALLOCATED(pzz8)) DEALLOCATE(pzz8)
1501 RETURN
1502 END
boltst
subroutine boltst(ip, bpreld, sig, tt, nel, npt, nsensor, sensor_tab)
Definition boltst.F:38
csmall3
subroutine csmall3(sig, s1, s2, s3, s4, s5, s6, offg, off, nel)
Definition csmall3.F:35
degenes8
subroutine degenes8(ixs, idege, nel)
Definition degenes8.F:37
area
subroutine area(d1, x, x2, y, y2, eint, stif0)
Definition hm_read_prop32.F:567
min
#define min(a, b)
Definition macros.h:20
max
#define max(a, b)
Definition macros.h:21
mmain
subroutine mmain(pm, elbuf_str, ix, nix, x, geo, iparg, nel, skew, bufmat, ipart, ipartel, nummat, matparam, imat, ipm, ngl, pid, npf, tf, mfxx, mfxy, mfxz, mfyx, mfyy, mfyz, mfzx, mfzy, mfzz, rx, ry, rz, sx, sy, sz, gama, voln, dvol, s1, s2, s3, s4, s5, s6, dxx, dyy, dzz, d4, d5, d6, wxx, wyy, wzz)
Definition mmain.F:43
ale_connectivity_mod
Definition ale_connectivity_mod.F:223
dt_mod
Definition dt_mod.F:79
nlocal_reg_mod
Definition nlocal_reg_mod.F:117
table_mod
Definition table_mod.F:112
s8cdefo3
subroutine s8cdefo3(pxc1, pxc2, pxc3, pxc4, pyc1, pyc2, pyc3, pyc4, pzc1, pzc2, pzc3, pzc4, px1, px2, px3, px4, px5, px6, px7, px8, py1, py2, py3, py4, py5, py6, py7, py8, pz1, pz2, pz3, pz4, pz5, pz6, pz7, pz8, pxy1, pxy2, pxy3, pxy4, pxy5, pxy6, pxy7, pxy8, pyx1, pyx2, pyx3, pyx4, pyx5, pyx6, pyx7, pyx8, pxz1, pxz2, pxz3, pxz4, pxz5, pxz6, pxz7, pxz8, pzx1, pzx2, pzx3, pzx4, pzx5, pzx6, pzx7, pzx8, pyz1, pyz2, pyz3, pyz4, pyz5, pyz6, pyz7, pyz8, pzy1, pzy2, pzy3, pzy4, pzy5, pzy6, pzy7, pzy8, bxy1, bxy2, bxy3, bxy4, bxy5, bxy6, bxy7, bxy8, byx1, byx2, byx3, byx4, byx5, byx6, byx7, byx8, bxz1, bxz2, bxz3, bxz4, bxz5, bxz6, bxz7, bxz8, bzx1, bzx2, bzx3, bzx4, bzx5, bzx6, bzx7, bzx8, byz1, byz2, byz3, byz4, byz5, byz6, byz7, byz8, bzy1, bzy2, bzy3, bzy4, bzy5, bzy6, bzy7, bzy8, vx1, vx2, vx3, vx4, vx5, vx6, vx7, vx8, vy1, vy2, vy3, vy4, vy5, vy6, vy7, vy8, vz1, vz2, vz3, vz4, vz5, vz6, vz7, vz8, dxx, dxy, dxz, dyx, dyy, dyz, dzx, dzy, dzz, d4, d5, d6, wxx, wyy, wzz, volo, off, eint, offs, dsv, icp, icr, ics, ict, dti, dsr, dss, dst, ngl, ideg, vol0dp, nel, ismstr, ipres)
Definition s8cdefo3.F:77
s8cfint_reg
subroutine s8cfint_reg(nloc_dmg, var_reg, nel, offg, vol, nc1, nc2, nc3, nc4, nc5, nc6, nc7, nc8, px1, px2, px3, px4, px5, px6, px7, px8, py1, py2, py3, py4, py5, py6, py7, py8, pz1, pz2, pz3, pz4, pz5, pz6, pz7, pz8, imat, itask, dt2t, vol0, nft, nlay, nptr, npts, ir, is, ws, as, icr, ics, ict, bufnlts, area, nodadt, dt1, dt12, dtfac1, dtmin1, idtmin)
Definition s8cfint_reg.F:49
s8cforc3
subroutine s8cforc3(timers, output, elbuf_tab, ng, pm, geo, ixs, x, a, v, ms, w, flux, flu1, veul, fv, ale_connect, iparg, tf, npf, bufmat, partsav, nloc_dmg, dt2t, neltst, ityptst, stifn, fsky, iads, offset, eani, iparts, f11, f21, f31, f12, f22, f32, f13, f23, f33, f14, f24, f34, f15, f25, f35, f16, f26, f36, f17, f27, f37, f18, f28, f38, nel, icp, icsig, smr, sms, smt, mfxx, mfxy, mfxz, mfyx, mfyy, mfyz, mfzx, mfzy, mfzz, nvc, ipm, itask, istrain, temp, fthe, fthesky, iexpan, igeo, nnpt, gresav, grth, igrth, mssa, dmels, table, xdp, voln, condn, condnsky, sensors, ioutprt, mat_elem, h3d_strain, dt, nodadt, dtfac1, dtmin1, idtmin, snpc, stf, sbufmat, svis, nsvois, idtmins, iresp, maxfunc, userl_avail, glob_therm, impl_s, idyna)
Definition s8cforc3.F:110
s8csigp3
subroutine s8csigp3(sig, e0, defp, fac, g_pla, nel)
Definition s8csigp3.F:32
s8etherm
subroutine s8etherm(pm, imat, vol, ni, nc1, nc2, nc3, nc4, nc5, nc6, nc7, nc8, px1, px2, px3, px4, py1, py2, py3, py4, pz1, pz2, pz3, pz4, px5, px6, px7, px8, py5, py6, py7, py8, pz5, pz6, pz7, pz8, dt1, tempnc, tel, heat, fphi, offg, off, partsav, iparts, vol0, nel, theaccfact)
Definition s8etherm.F:43
s8for_distor
subroutine s8for_distor(x1, x2, x3, x4, x5, x6, x7, x8, y1, y2, y3, y4, y5, y6, y7, y8, z1, z2, z3, z4, z5, z6, z7, z8, vx1, vx2, vx3, vx4, vx5, vx6, vx7, vx8, vy1, vy2, vy3, vy4, vy5, vy6, vy7, vy8, vz1, vz2, vz3, vz4, vz5, vz6, vz7, vz8, f11, f12, f13, f14, f15, f16, f17, f18, f21, f22, f23, f24, f25, f26, f27, f28, f31, f32, f33, f34, f35, f36, f37, f38, sti, sti_c, fld, mu, ll, istab, fqmax, nel, e_distor, dt1)
Definition s8for_distor.F:60
s8get_x3
subroutine s8get_x3(nc1, nc2, nc3, nc4, nc5, nc6, nc7, nc8, x1, x2, x3, x4, x5, x6, x7, x8, y1, y2, y3, y4, y5, y6, y7, y8, z1, z2, z3, z4, z5, z6, z7, z8, x, xdp, numnod, nel)
Definition s8get_x3.F:42
s8sav3
subroutine s8sav3(offg, sav, xd1, xd2, xd3, xd4, xd5, xd6, xd7, xd8, yd1, yd2, yd3, yd4, yd5, yd6, yd7, yd8, zd1, zd2, zd3, zd4, zd5, zd6, zd7, zd8, nel)
Definition s8sav3.F:41
s8zfint3
subroutine s8zfint3(sig, sigor, px1, px2, px3, px4, py1, py2, py3, py4, pz1, pz2, pz3, pz4, px5, px6, px7, px8, py5, py6, py7, py8, pz5, pz6, pz7, pz8, pxy1, pxy2, pxy3, pxy4, pxy5, pxy6, pxy7, pxy8, pyx1, pyx2, pyx3, pyx4, pyx5, pyx6, pyx7, pyx8, pxz1, pxz2, pxz3, pxz4, pxz5, pxz6, pxz7, pxz8, pzx1, pzx2, pzx3, pzx4, pzx5, pzx6, pzx7, pzx8, pyz1, pyz2, pyz3, pyz4, pyz5, pyz6, pyz7, pyz8, pzy1, pzy2, pzy3, pzy4, pzy5, pzy6, pzy7, pzy8, bxy1, bxy2, bxy3, bxy4, bxy5, bxy6, bxy7, bxy8, byx1, byx2, byx3, byx4, byx5, byx6, byx7, byx8, bxz1, bxz2, bxz3, bxz4, bxz5, bxz6, bxz7, bxz8, bzx1, bzx2, bzx3, bzx4, bzx5, bzx6, bzx7, bzx8, byz1, byz2, byz3, byz4, byz5, byz6, byz7, byz8, bzy1, bzy2, bzy3, bzy4, bzy5, bzy6, bzy7, bzy8, f11, f21, f31, f12, f22, f32, f13, f23, f33, f14, f24, f34, f15, f25, f35, f16, f26, f36, f17, f27, f37, f18, f28, f38, vol, qvis, pp, eint, rho, q, defp, epsd, epsdm, sigm, eintm, rhom, qm, defpm, volg, sti, stin, icp, off, vol0, vol0g, g_pla, g_epsd, jfac, eintth, eintthm, iexpan, nel, ideg, conde, conden, mtn, ismstr, svis, nodadt_therm, g_wpla, l_wpla, g_wpla_flag)
Definition s8zfint3.F:76
s8zfintp3
subroutine s8zfintp3(px1, px2, px3, px4, py1, py2, py3, py4, pz1, pz2, pz3, pz4, f11, f21, f31, f12, f22, f32, f13, f23, f33, f14, f24, f34, f15, f25, f35, f16, f26, f36, f17, f27, f37, f18, f28, f38, vol, pp, idege, nel)
Definition s8zfintp3.F:42
s8ztempel
subroutine s8ztempel(lft, llt, ni, ksi, eta, zeta, nc1, nc2, nc3, nc4, nc5, nc6, nc7, nc8, temp, tempel)
Definition s8ztempel.F:34
s8zzero3
subroutine s8zzero3(fx1, fy1, fz1, fx2, fy2, fz2, fx3, fy3, fz3, fx4, fy4, fz4, fx5, fy5, fz5, fx6, fy6, fz6, fx7, fy7, fz7, fx8, fy8, fz8, sigm, eintm, rhom, qm, defpm, epsdm, stin, pp, g_pla, g_epsd, iexpan, eintthm, nel, conden)
Definition s8zzero3.F:42
sbilan
subroutine sbilan(partsav, eint, rho, rk, vol, vx1, vx2, vx3, vx4, vx5, vx6, vx7, vx8, vy1, vy2, vy3, vy4, vy5, vy6, vy7, vy8, vz1, vz2, vz3, vz4, vz5, vz6, vz7, vz8, vnew, iparts, gresav, grth, igrth, off, iexpan, eintth, fill, mom, nel, ifvm22, x1, x2, x3, x4, x5, x6, x7, x8, y1, y2, y3, y4, y5, y6, y7, y8, z1, z2, z3, z4, z5, z6, z7, z8, itask, iparg, sensors, g_wpla, wpla)
Definition sbilan.F:47
sctorth3
subroutine sctorth3(jft, jlt, icstr, nel, rx, ry, rz, sx, sy, sz, tx, ty, tz, e1x, e1y, e1z, e2x, e2y, e2z, e3x, e3y, e3z, g1x, g1y, g1z, g2x, g2y, g2z, g3x, g3y, g3z, gama, irep)
Definition sctorth3.F:39
scumu3
subroutine scumu3(offg, e, nc1, nc2, nc3, nc4, nc5, nc6, nc7, nc8, stifn, sti, f11, f21, f31, f12, f22, f32, f13, f23, f33, f14, f24, f34, f15, f25, f35, f16, f26, f36, f17, f27, f37, f18, f28, f38, nvc, ar, fr_wave, fr_wav, mx1, my1, mz1, mx2, my2, mz2, mx3, my3, mz3, mx4, my4, mz4, mx5, my5, mz5, mx6, my6, mz6, mx7, my7, mz7, mx8, my8, mz8, them, fthe, condn, conde, nel, jthe, isrot, ipartsph, nodadt_therm)
Definition scumu3.F:55
scumu3p
subroutine scumu3p(offg, sti, fsky, fskyv, iads, f11, f21, f31, f12, f22, f32, f13, f23, f33, f14, f24, f34, f15, f25, f35, f16, f26, f36, f17, f27, f37, f18, f28, f38, nc1, nc2, nc3, nc4, nc5, nc6, nc7, nc8, ar, fr_wave, fr_wav, mx1, my1, mz1, mx2, my2, mz2, mx3, my3, mz3, mx4, my4, mz4, mx5, my5, mz5, mx6, my6, mz6, mx7, my7, mz7, mx8, my8, mz8, them, fthesky, condnsky, conde, nel, nft, jthe, isrot, ipartsph, nodadt_therm)
Definition scumu3p.F:56
scumualpha
subroutine scumualpha(offg, alpha_e, nc1, nc2, nc3, nc4, nc5, nc6, nc7, nc8, nel)
Definition scumualpha.F:35
sdlen_dege
subroutine sdlen_dege(volg, lat, x1, x2, x3, x4, x5, x6, x7, x8, y1, y2, y3, y4, y5, y6, y7, y8, z1, z2, z3, z4, z5, z6, z7, z8, ixs, idege, nel)
Definition sdlen_dege.F:45
sdlensh8
subroutine sdlensh8(voln, llsh, area, x1, x2, x3, x4, x5, x6, x7, x8, y1, y2, y3, y4, y5, y6, y7, y8, z1, z2, z3, z4, z5, z6, z7, z8, ics, nel)
Definition sdlensh8.F:38
sfillopt
subroutine sfillopt(fill, sti, f11, f21, f31, f12, f22, f32, f13, f23, f33, f14, f24, f34, f15, f25, f35, f16, f26, f36, f17, f27, f37, f18, f28, f38, nel)
Definition sfillopt.F:43
sgparav3
subroutine sgparav3(npe, x, ixs, rx, ry, rz, sx, sy, sz, tx, ty, tz, nel)
Definition sgparav3.F:35
smalla3
subroutine smalla3(sav, offg, off, wxx, wyy, wzz, nel, ismstr, jlag)
Definition smalla3.F:35
smallb3
subroutine smallb3(offg, off, nel, ismstr)
Definition smallb3.F:44
srbilan
subroutine srbilan(partsav, eint, rho, rk, vol, vxa, vya, vza, va2, vnew, iparts, gresav, grth, igrth, off, iexpan, eintth, fill, xx, yy, zz, xx2, yy2, zz2, xy, yz, zx, itask, iparg, sensors, nel, g_wpla, wpla)
Definition srbilan.F:45
srota3
subroutine srota3(sig, s1, s2, s3, s4, s5, s6, wxx, wyy, wzz, nel, mtn, ismstr)
Definition srota3.F:43
sstra3
subroutine sstra3(dxx, dyy, dzz, d4, d5, d6, strain, wxx, wyy, wzz, off, nel, jcvt)
Definition sstra3.F:46
basis8
subroutine basis8(r, s, t, h, pr, ps, pt)
Definition basis8.F:30
s8zderi3
subroutine s8zderi3(vol, veul, geo, ksi, eta, zeta, wi, hx, hy, hz, cj1, cj2, cj3, cj4, cj5, cj6, cj7, cj8, cj9, smax, deltax, ngl, voldp)
Definition s8zderi3.F:40
s8zderic3
subroutine s8zderic3(vol, hx, hy, hz, jac1, jac2, jac3, jac4, jac5, jac6, jac7, jac8, jac9, smax, det, ngl, xd1, xd2, xd3, xd4, xd5, xd6, xd7, xd8, yd1, yd2, yd3, yd4, yd5, yd6, yd7, yd8, zd1, zd2, zd3, zd4, zd5, zd6, zd7, zd8)
Definition s8zderi3.F:142
scoor3
subroutine scoor3(x, xrefs, ixs, geo, mxt, ngeo, ngl, ix1, ix2, ix3, ix4, ix5, ix6, ix7, ix8, x1, x2, x3, x4, x5, x6, x7, x8, y1, y2, y3, y4, y5, y6, y7, y8, z1, z2, z3, z4, z5, z6, z7, z8, rx, ry, rz, sx, sy, sz, tx, ty, tz, e1x, e1y, e1z, e2x, e2y, e2z, e3x, e3y, e3z, f1x, f1y, f1z, f2x, f2y, f2z, temp0, temp, nintemp, xd1, xd2, xd3, xd4, xd5, xd6, xd7, xd8, yd1, yd2, yd3, yd4, yd5, yd6, yd7, yd8, zd1, zd2, zd3, zd4, zd5, zd6, zd7, zd8)
Definition scoor3.F:52
srcoor3
subroutine srcoor3(x, xrefs, ixs, geo, mxt, ngeo, ngl, jhbe, ix1, ix2, ix3, ix4, ix5, ix6, ix7, ix8, x1, x2, x3, x4, x5, x6, x7, x8, y1, y2, y3, y4, y5, y6, y7, y8, z1, z2, z3, z4, z5, z6, z7, z8, rx, ry, rz, sx, sy, sz, tx, ty, tz, e1x, e1y, e1z, e2x, e2y, e2z, e3x, e3y, e3z, f1x, f1y, f1z, f2x, f2y, f2z, temp0, temp, nintemp, xd1, xd2, xd3, xd4, xd5, xd6, xd7, xd8, yd1, yd2, yd3, yd4, yd5, yd6, yd7, yd8, zd1, zd2, zd3, zd4, zd5, zd6, zd7, zd8)
Definition srcoor3.F:52
srho3
subroutine srho3(pm, volo, rhon, eint, dxx, dyy, dzz, voln, dvol, mat)
Definition srho3.F:31
srrota3
subroutine srrota3(r11, r12, r13, r21, r22, r23, r31, r32, r33, x1, x2, x3, x4, x5, x6, x7, x8, y1, y2, y3, y4, y5, y6, y7, y8, z1, z2, z3, z4, z5, z6, z7, z8)
Definition srrota3.F:33
sdlensh14
subroutine sdlensh14(nel, llsh, x1, x2, x3, x4, x5, x6, x7, x8, y1, y2, y3, y4, y5, y6, y7, y8, z1, z2, z3, z4, z5, z6, z7, z8, ics, idt1sol)
Definition sdlensh14.F:38
startime
subroutine startime(event, itask)
Definition timer.F:93
stoptime
subroutine stoptime(event, itask)
Definition timer.F:135
ale_connectivity_mod::t_ale_connectivity
Definition ale_connectivity_mod.F:254
table_mod::ttable
Definition table_mod.F:125
storth3
subroutine storth3(isorth, nel, g1x, g1y, g1z, g2x, g2y, g2z, g3x, g3y, g3z, gama)
Definition szorth3.F:30
sroto3
subroutine sroto3(nel, sig, sign, g1x, g1y, g1z, g2x, g2y, g2z, g3x, g3y, g3z)
Definition szorth3.F:136
szordef3
subroutine szordef3(nel, dxx, dyy, dzz, d4, d5, d6, g1x, g1y, g1z, g2x, g2y, g2z, g3x, g3y, g3z)
Definition szorth3.F:84
tshgeodel3
subroutine tshgeodel3(ngl, offg, volg, area, volg0, l_sh, geo, nnod, dt, nel)
Definition tshgeodel3.F:36