OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
cfint_reg_ini.F File Reference
#include "implicit_f.inc"
#include "param_c.inc"
#include "com01_c.inc"
#include "com04_c.inc"
#include "scr03_c.inc"
#include "vectorize.inc"

Go to the source code of this file.

Functions/Subroutines

subroutine cfint_reg_ini (elbuf_tab, nloc_dmg, area, ixc, dt_nl, x, xrefc, nft, nel, ng, ipm, bufmat, time, failure)

Function/Subroutine Documentation

◆ cfint_reg_ini()

subroutine cfint_reg_ini ( type(elbuf_struct_), dimension(ngroup), target elbuf_tab,
type (nlocal_str_), target nloc_dmg,
intent(in) area,
integer, dimension(nixc,*) ixc,
dt_nl,
x,
xrefc,
integer nft,
integer nel,
integer ng,
integer, dimension(npropmi,*) ipm,
bufmat,
time,
logical failure )

Definition at line 33 of file cfint_reg_ini.F.

37C-----------------------------------------------
38C M o d u l e s
39C-----------------------------------------------
40 USE elbufdef_mod
41 USE message_mod
42 USE nlocal_reg_mod
43C-----------------------------------------------
44C I m p l i c i t T y p e s
45C-----------------------------------------------
46#include "implicit_f.inc"
47C-----------------------------------------------
48C C o m m o n B l o c k s
49C-----------------------------------------------
50#include "param_c.inc"
51#include "com01_c.inc"
52#include "com04_c.inc"
53#include "scr03_c.inc"
54C-----------------------------------------------
55C D u m m y A r g u m e n t s
56C-----------------------------------------------
57 TYPE(ELBUF_STRUCT_), TARGET, DIMENSION(NGROUP) :: ELBUF_TAB
58 TYPE (NLOCAL_STR_) , TARGET :: NLOC_DMG
59 INTEGER IXC(NIXC,*),NFT,NEL,NG,IPM(NPROPMI,*)
60 my_real ,DIMENSION(NUMELC+NUMELTG),INTENT(IN) ::
61 . area
62 my_real
63 . x(3,*),xrefc(4,3,*),dt_nl,bufmat(*),time
64 LOGICAL :: FAILURE
65C-----------------------------------------------
66C L o c a l V a r i a b l e s
67C-----------------------------------------------
68 INTEGER :: IMAT,NDOF,L_NLOC,N1,N2,N3,N4,K,I,NDNOD
69 INTEGER, DIMENSION(NEL) :: POS1,POS2,POS3,POS4
70 my_real
71 . le_min,len,damp,dens,ntn_unl,ntn_vnl,
72 . ntvar,z01(11,11),wf1(11,11),zn1(12,11),b1,b2,
73 . b3,b4,nth1,nth2,bth1,bth2,k1,k12,k2,sspnl,le_max
74 my_real, DIMENSION(:,:), ALLOCATABLE :: var_reg,vpred
75 my_real,
76 . DIMENSION(:), POINTER :: fnl,unl,vnl,dnl,mnl,thck
77 my_real, DIMENSION(NEL) :: x1,x2,x3,x4,
78 . y1,y2,y3,y4,px1,px2,py1,py2,e1x,e2x,e3x,
79 . e1y,e2y,e3y,e1z,e2z,e3z,x2l,y2l,x3l,y3l,
80 . x4l,y4l,z1,z2,z3,z4,surf,offg,vols,btb11,
81 . btb12,btb22
82 TYPE(BUF_NLOC_),POINTER :: BUFNL
83 my_real, DIMENSION(:,:), POINTER ::
84 . massth,fnlth,vnlth,unlth
85 ! Position of integration points in the thickness
86 DATA z01/
87 1 0. ,0. ,0. ,0. ,0. ,
88 1 0. ,0. ,0. ,0. ,0. ,0. ,
89 2 -.5 ,0.5 ,0. ,0. ,0. ,
90 2 0. ,0. ,0. ,0. ,0. ,0. ,
91 3 -.5 ,0. ,0.5 ,0. ,0. ,
92 3 0. ,0. ,0. ,0. ,0. ,0. ,
93 4 -.5 ,-.1666667,0.1666667,0.5 ,0. ,
94 4 0. ,0. ,0. ,0. ,0. ,0. ,
95 5 -.5 ,-.25 ,0. ,0.25 ,0.5 ,
96 5 0. ,0. ,0. ,0. ,0. ,0. ,
97 6 -.5 ,-.3 ,-.1 ,0.1 ,0.3 ,
98 6 0.5 ,0. ,0. ,0. ,0. ,0. ,
99 7 -.5 ,-.3333333,-.1666667,0.0 ,0.1666667,
100 7 0.3333333,0.5 ,0. ,0. ,0. ,0. ,
101 8 -.5 ,-.3571429,-.2142857,-.0714286,0.0714286,
102 8 0.2142857,0.3571429,0.5 ,0. ,0. ,0. ,
103 9 -.5 ,-.375 ,-.25 ,-.125 ,0.0 ,
104 9 0.125 ,0.25 ,0.375 ,0.5 ,0. ,0. ,
105 a -.5 ,-.3888889,-.2777778,-.1666667,-.0555555,
106 a 0.0555555,0.1666667,0.2777778,0.3888889,0.5 ,0. ,
107 b -.5 ,-.4 ,-.3 ,-.2 ,-.1 ,
108 b 0. ,0.1 ,0.2 ,0.3 ,0.4 ,0.5 /
109 ! Weight of integration in the thickness
110 DATA wf1/
111 1 1. ,0. ,0. ,0. ,0. ,
112 1 0. ,0. ,0. ,0. ,0. ,0. ,
113 2 0.5 ,0.5 ,0. ,0. ,0. ,
114 2 0. ,0. ,0. ,0. ,0. ,0. ,
115 3 0.25 ,0.5 ,0.25 ,0. ,0. ,
116 3 0. ,0. ,0. ,0. ,0. ,0. ,
117 4 0.1666667,0.3333333,0.3333333,0.1666667,0. ,
118 4 0. ,0. ,0. ,0. ,0. ,0. ,
119 5 0.125 ,0.25 ,0.25 ,0.25 ,0.125 ,
120 5 0. ,0. ,0. ,0. ,0. ,0. ,
121 6 0.1 ,0.2 ,0.2 ,0.2 ,0.2 ,
122 6 0.1 ,0. ,0. ,0. ,0. ,0. ,
123 7 0.0833333,0.1666667,0.1666667,0.1666667,0.1666667,
124 7 0.1666667,0.0833333,0. ,0. ,0. ,0. ,
125 8 0.0714286,0.1428571,0.1428571,0.1428571,0.1428571,
126 8 0.1428571,0.1428571,0.0714286,0. ,0. ,0. ,
127 9 0.0625 ,0.125 ,0.125 ,0.125 ,0.125 ,
128 9 0.125 ,0.125 ,0.125 ,0.0625 ,0. ,0. ,
129 a 0.0555556,0.1111111,0.1111111,0.1111111,0.1111111,
130 a 0.1111111,0.1111111,0.1111111,0.1111111,0.0555556,0. ,
131 b 0.05 ,0.1 ,0.1 ,0.1 ,0.1 ,
132 b 0.1 ,0.1 ,0.1 ,0.1 ,0.1 ,0.05 /
133 ! Position of nodes in the shell thickness
134 DATA zn1/
135 1 0. ,0. ,0. ,0. ,0. ,0. ,
136 1 0. ,0. ,0. ,0. ,0. ,0. ,
137 2 -.5 ,0.5 ,0. ,0. ,0. ,0. ,
138 2 0. ,0. ,0. ,0. ,0. ,0. ,
139 3 -.5 ,-.25 ,0.25 ,0.5 ,0. ,0. ,
140 3 0. ,0. ,0. ,0. ,0. ,0. ,
141 4 -.5 ,-.3333333,0. ,0.3333333,0.5 ,0. ,
142 4 0. ,0. ,0. ,0. ,0. ,0. ,
143 5 -.5 ,-.375 ,-0.125 ,0.125 ,0.375 ,0.5 ,
144 5 0. ,0. ,0. ,0. ,0. ,0. ,
145 6 -.5 ,-.4 ,-.2 ,0.0 ,0.2 ,0.4 ,
146 6 0.5 ,0. ,0. ,0. ,0. ,0. ,
147 7 -.5 ,-.4166667,-.25 ,-.0833333,0.0833333,0.25 ,
148 7 0.4166667,0.5 ,0. ,0. ,0. ,0. ,
149 8 -.5 ,-.4285715,-.2857143,-.1428572,0.0 ,0.1428572,
150 8 0.2857143,0.4285715,0.5 ,0. ,0. ,0. ,
151 9 -.5 ,-.4375 ,-.3125 ,-.1875 ,-.0625 ,0.0625 ,
152 9 0.1875 ,0.3125 ,0.4375 ,0.5 ,0. ,0. ,
153 a -.5 ,-.4444444,-.3333333,-.2222222,-.1111111,0. ,
154 a 0.1111111,0.2222222,0.3333333,0.4444444,0.5 ,0. ,
155 b -.5 ,-.45 ,-.35 ,-.25 ,-.15 ,-.05 ,
156 b 0.05 ,0.15 ,0.25 ,0.35 ,0.45 ,0.5 /
157C=======================================================================
158 ! Size of the non-local vectors
159 l_nloc = nloc_dmg%L_NLOC
160 ! Pointing the non-local forces vector
161 fnl => nloc_dmg%FNL(1:l_nloc,1)
162 vnl => nloc_dmg%VNL(1:l_nloc)
163 dnl => nloc_dmg%DNL(1:l_nloc)
164 unl => nloc_dmg%UNL(1:l_nloc)
165 mnl => nloc_dmg%MASS(1:l_nloc)
166 ! Number of the material law
167 imat = ixc(1,1+nft)
168 ! Minimal length
169 le_min = sqrt(minval(area(nft+1:nft+nel)))
170 ! Number of integration points
171 ndof = elbuf_tab(ng)%BUFLY(1)%NPTT
172 ! Thickness of the shell
173 thck => elbuf_tab(ng)%GBUF%THK(1:nel)
174 ! Global minimal length
175 IF (ndof>1) THEN
176 le_min = min(le_min,minval(thck(1:nel))/ndof)
177 ENDIF
178 ! Non-local internal length
179 len = nloc_dmg%LEN(imat)
180 ! Maximal length of convergence
181 le_max = nloc_dmg%LE_MAX(imat)
182 ! Non-local damping
183 damp = nloc_dmg%DAMP(imat)
184 ! Non-local density
185 dens = nloc_dmg%DENS(imat)
186 ! Non-local sound speed
187 sspnl = nloc_dmg%SSPNL(imat)
188 ! Non-local timestep
189 dt_nl = min(dt_nl,0.5d0*((two*min(le_min,le_max)*sqrt(three*dens))/
190 . (sqrt(twelve*(len**2)+(min(le_min,le_max)**2)))))
191 ! Allocation of the velocities predictor
192 IF (ndof>1) THEN
193 IF (ndof > 2) THEN
194 ALLOCATE(vpred(nel,ndof+1))
195 ndnod = ndof + 1
196 ELSE
197 ALLOCATE(vpred(nel,ndof))
198 ndnod = ndof
199 ENDIF
200 ENDIF
201 ! Variable to regularize
202 IF (.NOT.ALLOCATED(var_reg)) ALLOCATE(var_reg(nel,ndof))
203c
204 ! Computation of kinematical data
205# include "vectorize.inc"
206 DO i = 1,nel
207 ! Coordinates of the nodes of the element
208 IF (nxref == 0) THEN
209 x1(i)=x(1,ixc(2,nft+i))
210 y1(i)=x(2,ixc(2,nft+i))
211 z1(i)=x(3,ixc(2,nft+i))
212 x2(i)=x(1,ixc(3,nft+i))
213 y2(i)=x(2,ixc(3,nft+i))
214 z2(i)=x(3,ixc(3,nft+i))
215 x3(i)=x(1,ixc(4,nft+i))
216 y3(i)=x(2,ixc(4,nft+i))
217 z3(i)=x(3,ixc(4,nft+i))
218 x4(i)=x(1,ixc(5,nft+i))
219 y4(i)=x(2,ixc(5,nft+i))
220 z4(i)=x(3,ixc(5,nft+i))
221 ELSE
222 x1(i)=xrefc(1,1,nft+i)
223 y1(i)=xrefc(1,2,nft+i)
224 z1(i)=xrefc(1,3,nft+i)
225 x2(i)=xrefc(2,1,nft+i)
226 y2(i)=xrefc(2,2,nft+i)
227 z2(i)=xrefc(2,3,nft+i)
228 x3(i)=xrefc(3,1,nft+i)
229 y3(i)=xrefc(3,2,nft+i)
230 z3(i)=xrefc(3,3,nft+i)
231 x4(i)=xrefc(4,1,nft+i)
232 y4(i)=xrefc(4,2,nft+i)
233 z4(i)=xrefc(4,3,nft+i)
234 ENDIF
235C
236 ! Recovering the nodes of the non-local element
237 n1 = nloc_dmg%IDXI(ixc(2,nft+i))
238 n2 = nloc_dmg%IDXI(ixc(3,nft+i))
239 n3 = nloc_dmg%IDXI(ixc(4,nft+i))
240 n4 = nloc_dmg%IDXI(ixc(5,nft+i))
241 ! Recovering the positions of the first d.o.fs of each nodes
242 pos1(i) = nloc_dmg%POSI(n1)
243 pos2(i) = nloc_dmg%POSI(n2)
244 pos3(i) = nloc_dmg%POSI(n3)
245 pos4(i) = nloc_dmg%POSI(n4)
246 ENDDO
247c
248 ! Non-local variable transfer at Gauss point
249 ! Loop over Gauss points in the thickness
250 DO k = 1,ndof
251 ! Loop over element
252 DO i = 1,nel
253 var_reg(i,k) = fourth*(dnl(pos1(i)+k-1) + dnl(pos2(i)+k-1)
254 . + dnl(pos3(i)+k-1) + dnl(pos4(i)+k-1))
255 ENDDO
256 ENDDO
257c
258 CALL ceveci(1,nel ,surf,
259 . x1 ,x2 ,x3 ,x4 ,y1 ,y2 ,
260 . y3 ,y4 ,z1 ,z2 ,z3 ,z4 ,
261 . e1x ,e2x ,e3x ,e1y ,e2y ,e3y ,
262 . e1z ,e2z ,e3z )
263C
264 ! Filling internal variable data of the non-local material
265 CALL cnloc_matini(elbuf_tab(ng),nel ,ipm ,
266 . bufmat ,time ,var_reg ,
267 . failure )
268C
269 !-----------------------------------------------------------------------
270 ! Computation of the element volume and the BtB matrix product
271 !-----------------------------------------------------------------------
272 ! Loop over elements
273# include "vectorize.inc"
274 DO i=1,nel
275c
276 ! Computation of shape functions derivatives
277 x2l(i) = e1x(i)*(x2(i)-x1(i))+e1y(i)*(y2(i)-y1(i))+e1z(i)*(z2(i)-z1(i))
278 y2l(i) = e2x(i)*(x2(i)-x1(i))+e2y(i)*(y2(i)-y1(i))+e2z(i)*(z2(i)-z1(i))
279 x3l(i) = e1x(i)*(x3(i)-x1(i))+e1y(i)*(y3(i)-y1(i))+e1z(i)*(z3(i)-z1(i))
280 y3l(i) = e2x(i)*(x3(i)-x1(i))+e2y(i)*(y3(i)-y1(i))+e2z(i)*(z3(i)-z1(i))
281 x4l(i) = e1x(i)*(x4(i)-x1(i))+e1y(i)*(y4(i)-y1(i))+e1z(i)*(z4(i)-z1(i))
282 y4l(i) = e2x(i)*(x4(i)-x1(i))+e2y(i)*(y4(i)-y1(i))+e2z(i)*(z4(i)-z1(i))
283c
284 px1(i) = half *(y2l(i)-y4l(i))
285 px2(i) = half * y3l(i)
286 py1(i) = -half *(x2l(i)-x4l(i))
287 py2(i) = -half * x3l(i)
288c
289 ! Computation of the product BtxB
290 btb11(i) = px1(i)**2 + py1(i)**2
291 btb12(i) = px1(i)*px2(i) + py1(i)*py2(i)
292 btb22(i) = px2(i)**2 + py2(i)**2
293c
294 ! Computation of the element volume
295 vols(i) = area(nft+i)*thck(i)
296c
297 ! To check if element is not broken
298 offg(i) = elbuf_tab(ng)%GBUF%OFF(i)
299c
300 ENDDO
301C
302 !-----------------------------------------------------------------------
303 ! Pre-treatment non-local regularization in the thickness
304 !-----------------------------------------------------------------------
305 ! Only if NDOF > 1
306 IF ((ndof > 1).AND.(len>zero)) THEN
307c
308 ! Pointing the non-local values in the thickness of the corresponding element
309 bufnl => elbuf_tab(ng)%NLOC(1,1)
310c
311 ! Pointing the non-local values in the thickness of the corresponding element
312 massth => bufnl%MASSTH(1:nel,1:ndnod)
313 unlth => bufnl%UNLTH(1:nel ,1:ndnod)
314 vnlth => bufnl%VNLTH(1:nel ,1:ndnod)
315 fnlth => bufnl%FNLTH(1:nel ,1:ndnod)
316c
317 DO k = 1,ndnod
318 DO i = 1,nel
319 ! Prediction of the velocities
320 vpred(i,k) = vnlth(i,k) - (fnlth(i,k)/massth(i,k))*(dt_nl/two)
321 ENDDO
322 ENDDO
323 DO k = 1,ndnod
324 DO i = 1,nel
325 ! Resetting non-local forces
326 fnlth(i,k) = zero
327 ENDDO
328 ENDDO
329c
330 ! Computation of non-local forces in the shell thickness
331 DO k = 1, ndof
332c
333 ! Computation of shape functions value
334 IF ((ndof==2).AND.(k==2)) THEN
335 nth1 = (z01(k,ndof) - zn1(k,ndof))/(zn1(k-1,ndof) - zn1(k,ndof))
336 nth2 = (z01(k,ndof) - zn1(k-1,ndof)) /(zn1(k,ndof) - zn1(k-1,ndof))
337 ELSE
338 nth1 = (z01(k,ndof) - zn1(k+1,ndof))/(zn1(k,ndof) - zn1(k+1,ndof))
339 nth2 = (z01(k,ndof) - zn1(k,ndof)) /(zn1(k+1,ndof) - zn1(k,ndof))
340 ENDIF
341c
342 ! Loop over elements
343 DO i = 1,nel
344 ! Computation of B-matrix values
345 IF ((ndof==2).AND.(k==2)) THEN
346 bth1 = (one/(zn1(k-1,ndof) - zn1(k,ndof)))*(one/thck(i))
347 bth2 = (one/(zn1(k,ndof) - zn1(k-1,ndof)))*(one/thck(i))
348 ELSE
349 bth1 = (one/(zn1(k,ndof) - zn1(k+1,ndof)))*(one/thck(i))
350 bth2 = (one/(zn1(k+1,ndof) - zn1(k,ndof)))*(one/thck(i))
351 ENDIF
352c
353 ! Computation of the non-local K matrix
354 k1 = (len**2)*(bth1**2) + nth1**2
355 k12 = (len**2)*(bth1*bth2)+ (nth1*nth2)
356 k2 = (len**2)*(bth2**2) + nth2**2
357c
358 ! Computation of the non-local forces
359 IF ((ndof==2).AND.(k==2)) THEN
360 fnlth(i,k-1) = fnlth(i,k-1) + (k1*unlth(i,k-1) + k12*unlth(i,k)
361 . + damp*((nth1**2)*vpred(i,k-1)
362 . + (nth1*nth2)*vpred(i,k))
363 . - (nth1*var_reg(i,k)))*vols(i)*wf1(k,ndof)
364 fnlth(i,k) = fnlth(i,k) + (k12*unlth(i,k-1) + k2*unlth(i,k)
365 . + damp*(nth1*nth2*vpred(i,k-1)
366 . + (nth2**2)*vpred(i,k))
367 . - nth2*var_reg(i,k))*vols(i)*wf1(k,ndof)
368 ELSE
369 fnlth(i,k) = fnlth(i,k) + (k1*unlth(i,k) + k12*unlth(i,k+1)
370 . + damp*((nth1**2)*vpred(i,k)
371 . + (nth1*nth2)*vpred(i,k+1))
372 . - (nth1*var_reg(i,k)))*vols(i)*wf1(k,ndof)
373 fnlth(i,k+1) = fnlth(i,k+1) + (k12*unlth(i,k) + k2*unlth(i,k+1)
374 . + damp*(nth1*nth2*vpred(i,k)
375 . + (nth2**2)*vpred(i,k+1))
376 . - nth2*var_reg(i,k))*vols(i)*wf1(k,ndof)
377 ENDIF
378 ENDDO
379 ENDDO
380c
381 DO k = 1,ndnod
382 DO i = 1,nel
383 ! Updating the non-local in-thickness velocities
384 vnlth(i,k) = vnlth(i,k) - (fnlth(i,k)/massth(i,k))*dt_nl
385 ENDDO
386 ENDDO
387c
388 DO k = 1,ndnod
389 DO i = 1,nel
390 ! Computing the non-local in-thickness cumulated values
391 unlth(i,k) = unlth(i,k) + vnlth(i,k)*dt_nl
392 ENDDO
393 ENDDO
394c
395 ! Transfer at the integration point
396 DO k = 1, ndof
397 !Computation of shape functions value
398 IF ((ndof==2).AND.(k==2)) THEN
399 nth1 = (z01(k,ndof) - zn1(k,ndof))/(zn1(k-1,ndof) - zn1(k,ndof))
400 nth2 = (z01(k,ndof) - zn1(k-1,ndof)) /(zn1(k,ndof) - zn1(k-1,ndof))
401 ELSE
402 nth1 = (z01(k,ndof) - zn1(k+1,ndof))/(zn1(k,ndof) - zn1(k+1,ndof))
403 nth2 = (z01(k,ndof) - zn1(k,ndof)) /(zn1(k+1,ndof) - zn1(k,ndof))
404 ENDIF
405 ! Loop over elements
406 DO i = 1,nel
407 !Integration points non-local variables
408 IF ((ndof==2).AND.(k==2)) THEN
409 var_reg(i,k) = nth1*unlth(i,k-1) + nth2*unlth(i,k)
410 ELSE
411 var_reg(i,k) = nth1*unlth(i,k) + nth2*unlth(i,k+1)
412 ENDIF
413 ENDDO
414 ENDDO
415 ENDIF
416c
417 !-----------------------------------------------------------------------
418 ! Computation of the elementary non-local forces
419 !-----------------------------------------------------------------------
420 ! Loop over integration points in the thickness
421 DO k = 1,ndof
422c
423 ! Computation of non-local forces
424# include "vectorize.inc"
425 DO i = 1,nel
426c
427 ! Computing the elementary non-local forces
428 IF (offg(i) > zero) THEN
429 ! Computing the product BtB*UNL
430 b1 = ((len**2)/vols(i))*wf1(k,ndof)*(btb11(i)*unl(pos1(i)+k-1) + btb12(i)*unl(pos2(i)+k-1)
431 . - btb11(i)*unl(pos3(i)+k-1) - btb12(i)*unl(pos4(i)+k-1))
432 b2 = ((len**2)/vols(i))*wf1(k,ndof)*(btb12(i)*unl(pos1(i)+k-1) + btb22(i)*unl(pos2(i)+k-1)
433 . - btb12(i)*unl(pos3(i)+k-1) - btb22(i)*unl(pos4(i)+k-1))
434 b3 = ((len**2)/vols(i))*wf1(k,ndof)*(-btb11(i)*unl(pos1(i)+k-1) - btb12(i)*unl(pos2(i)+k-1)
435 . + btb11(i)*unl(pos3(i)+k-1) + btb12(i)*unl(pos4(i)+k-1))
436 b4 = ((len**2)/vols(i))*wf1(k,ndof)*(-btb12(i)*unl(pos1(i)+k-1) - btb22(i)*unl(pos2(i)+k-1)
437 . + btb12(i)*unl(pos3(i)+k-1) + btb22(i)*unl(pos4(i)+k-1))
438 ! Computing the product NtN*UNL
439 ntn_unl = ((unl(pos1(i)+k-1) + unl(pos2(i)+k-1) + unl(pos3(i)+k-1) +
440 . unl(pos4(i)+k-1))*fourth*fourth)*vols(i)*wf1(k,ndof)
441 ! Computing the product DAMP*NtN*VNL!
442 ntn_vnl = ((vnl(pos1(i)+k-1) + vnl(pos2(i)+k-1) + vnl(pos3(i)+k-1) +
443 . vnl(pos4(i)+k-1))*fourth*fourth)*damp*vols(i)*wf1(k,ndof)
444 ! Introducing the internal variable to be regularized
445 ntvar = var_reg(i,k)*fourth*vols(i)*wf1(k,ndof)
446 ! Computing the elementary non-local forces
447 fnl(pos1(i)+k-1) = fnl(pos1(i)+k-1) - (ntn_unl + ntn_vnl - ntvar + b1)
448 fnl(pos2(i)+k-1) = fnl(pos2(i)+k-1) - (ntn_unl + ntn_vnl - ntvar + b2)
449 fnl(pos3(i)+k-1) = fnl(pos3(i)+k-1) - (ntn_unl + ntn_vnl - ntvar + b3)
450 fnl(pos4(i)+k-1) = fnl(pos4(i)+k-1) - (ntn_unl + ntn_vnl - ntvar + b4)
451 ! If the element is broken, the non-local wave is absorbed
452 ELSE
453 ! Non-local absorbing forces
454 fnl(pos1(i)+k-1) = fnl(pos1(i)+k-1) - wf1(k,ndof)*dens*sspnl*vnl(pos1(i)+k-1)*le_max*thck(i)
455 fnl(pos2(i)+k-1) = fnl(pos2(i)+k-1) - wf1(k,ndof)*dens*sspnl*vnl(pos2(i)+k-1)*le_max*thck(i)
456 fnl(pos3(i)+k-1) = fnl(pos3(i)+k-1) - wf1(k,ndof)*dens*sspnl*vnl(pos3(i)+k-1)*le_max*thck(i)
457 fnl(pos4(i)+k-1) = fnl(pos4(i)+k-1) - wf1(k,ndof)*dens*sspnl*vnl(pos4(i)+k-1)*le_max*thck(i)
458 ENDIF
459 ENDDO
460 ENDDO
461c -------------------
462 IF (ALLOCATED(var_reg)) DEALLOCATE(var_reg)
463 IF (ALLOCATED(vpred)) DEALLOCATE(vpred)
ceveci
subroutine ceveci(jft, jlt, area, x1, x2, x3, x4, y1, y2, y3, y4, z1, z2, z3, z4, e1x, e2x, e3x, e1y, e2y, e3y, e1z, e2z, e3z)
Definition ceveci.F:37
cnloc_matini
subroutine cnloc_matini(elbuf_str, nel, ipm, bufmat, time, varnl, failure)
Definition cnloc_matini.F:38
my_real
#define my_real
Definition cppsort.cpp:32
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
message_mod
Definition message_mod.F:1249
nlocal_reg_mod
Definition nlocal_reg_mod.F:117