OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
write_buf_law51.F File Reference
#include "implicit_f.inc"
#include "comlock.inc"
#include "com01_c.inc"
#include "mvsiz_p.inc"
#include "inter22.inc"

Go to the source code of this file.

Functions/Subroutines

subroutine write_buf_law51 (ixs, nft, nuvar, nel, uvar, i, unephase, dd, dbvold, dbvold_f, volume, vold, epspxx, epspyy, epspzz, tag22, bfrac, rho01, rho02, rho03, rho04)

Function/Subroutine Documentation

◆ write_buf_law51()

subroutine write_buf_law51 ( integer, dimension(nixs,*) ixs,
integer nft,
integer nuvar,
integer nel,
uvar,
integer i,
integer unephase,
dd,
dbvold,
dbvold_f,
volume,
vold,
epspxx,
epspyy,
epspzz,
tag22,
bfrac,
rho01,
rho02,
rho03,
rho04 )

Definition at line 30 of file write_buf_law51.F.

36C-----------------------------------------------
37C M o d u l e s
38C-----------------------------------------------
39 USE elbufdef_mod
40 USE i22bufbric_mod
41 USE multimat_param_mod , ONLY : m51_n0phas, m51_nvphas
42C-----------------------------------------------
43C I m p l i c i t T y p e s
44C-----------------------------------------------
45#include "implicit_f.inc"
46#include "comlock.inc"
47C-----------------------------------------------
48C C o m m o n B l o c k s
49C-----------------------------------------------
50#include "com01_c.inc"
51#include "mvsiz_p.inc"
52#include "inter22.inc"
53C-----------------------------------------------
54C D u m m y A r g u m e n t s
55C-----------------------------------------------
56 INTEGER :: NEL, NUPARAM, NUVAR,IXS(NIXS,*)
57 INTEGER :: NFT, UNEPHASE
58 my_real :: uvar(nel,nuvar), dd, dbvold(4),dbvold_f(4),tag22(mvsiz),volume(mvsiz),vold, var
59 my_real :: epspxx(mvsiz), epspyy(mvsiz), epspzz(mvsiz),bfrac
60 my_real :: rho01, rho02, rho03, rho04, rho0(4)
61
62C-----------------------------------------------
63C L o c a l V a r i a b l e s
64C-----------------------------------------------
65 INTEGER :: I,II, ITRIMAT, IPRESENT(4), Q, KK
66 my_real :: p
67C-----------------------------------------------
68C S o u r c e L i n e s
69C-----------------------------------------------
70
71 rho0 = ( / rho01, rho02, rho03, rho04 / )
72
73 ii = i+ nft
74 write (*,fmt='(A16,I10)' ) "LAW51- brick ID=",ixs(nixs,i+nft)
75
76 var=zero
77 DO itrimat = 1,trimat
78 kk = m51_n0phas + (itrimat-1)*m51_nvphas
79 var = var + uvar(i,11+kk)
80 ENDDO
81
82 IF(int22>0)THEN
83 IF(tag22(i)>zero)THEN
84 write (*,fmt='(A)')
85 ." +---------SuperCell-----------+---------ELEMENT-----------+-------SUM(submat)-----------+-------SUM(submat
86 .)+flux-----+"
87 write (*,fmt='(A16,4F30.16)') " Vold=",brick_list(1,nint(tag22(i)))%Vold_SCell,vold ,sum(dbvold(1:4)),
88 . sum(dbvold_f(1:4))
89 write (*,fmt='(A16,3F30.16)') " Vnew=",brick_list(1,nint(tag22(i)))%Vnew_SCell,volume(i) ,var
90 ELSE
91 write (*,fmt='(A)')
92 ." +-----------ELEMENT-----------+"
93 write (*,fmt='(A16,1F30.16)') " Vold=",vold
94 write (*,fmt='(A16,1F30.16)') " Vnew=",volume(i)
95 ENDIF
96 ELSE
97 write (*,fmt='(A)')
98 ." +-----------ELEMENT-----------+"
99 write (*,fmt='(A16,1F30.16)') " Vold=",vold
100 write (*,fmt='(A16,1F30.16)') " Vnew=",volume(i)
101 ENDIF
102 write (*,fmt='(A16,3F30.16)') " EPSii=",epspxx(i),epspyy(i), epspzz(i)
103 write (*,fmt='(A16,1F30.16)') " DD=", dd
104 write (*,fmt='(A16,1F30.16)') " BFRAC=", bfrac
105
106 !submat4
107 q = unephase / 8
108 ipresent(4) = q
109 unephase = unephase - q*8
110 !submat3
111 q = unephase / 4
112 ipresent(3) = q
113 unephase = unephase - q*4
114 !submat2
115 q = unephase / 2
116 ipresent(2) = q
117 unephase = unephase - q*2
118 !submat1
119 ipresent(1) = unephase
120
121 var = zero
122
123 DO itrimat = 1,trimat
124 !=====================!
125 ! submat_id = ITRIMAT !
126 !=====================!
127 kk = m51_n0phas + (itrimat-1)*m51_nvphas
128 write (*,fmt='(A16,I1)') " +-----submat=",itrimat
129 IF(ipresent(itrimat)==0)THEN
130 write (*,fmt='(A16 )') " empty"
131 cycle
132 ENDIF
133 write (*,fmt='(A16,E30.16)') " 1. Vfrac=", uvar(i,1+kk)
134 write (*,fmt='(A16,E30.16)') " 2. Sd_xx=", uvar(i,2+kk)
135 write (*,fmt='(A16,E30.16)') " 3. Sd_yy=", uvar(i,3+kk)
136 write (*,fmt='(A16,E30.16)') " 4. Sd_zz=", uvar(i,4+kk)
137 write (*,fmt='(A16,E30.16)') " 5. Sd_xy=", uvar(i,5+kk)
138 write (*,fmt='(A16,E30.16)') " 6. Sd_yz=", uvar(i,6+kk)
139 write (*,fmt='(A16,E30.16)') " 7. Sd_zx=", uvar(i,7+kk)
140 write (*,fmt='(A16,E30.16)') " 8. E=", uvar(i,8+kk)
141 write (*,fmt='(A16,E30.16)') " 9. RHO=", uvar(i,9+kk)
142 write (*,fmt='(A16,E30.16)') " 8*11 EINT=", uvar(i,8+kk)*uvar(i,11+kk)
143 write (*,fmt='(A16,E30.16)') " 9*11 MASS=", uvar(i,9+kk)*uvar(i,11+kk)
144 write (*,fmt='(A16,E30.16)') " 10. Q=", uvar(i,10+kk)
145 write (*,fmt='(A16,E30.16)') " 11. Vnew=", uvar(i,11+kk)
146 write (*,fmt='(A16,E30.16)') " Vold=", dbvold(itrimat)
147 write (*,fmt='(A16,E30.16)') " 12. RHO=", uvar(i,12+kk)
148 write (*,fmt='(A16,E30.16)') " 13. ddvol=", UVAR(I,13+KK)
149 write (*,FMT='(A16,E30.16)') " 14. ssp=", UVAR(I,14+KK)
150 write (*,FMT='(A16,E30.16)') " 15. plas=", UVAR(I,15+KK)
151 write (*,FMT='(A16,E30.16)') " 16. t=", UVAR(I,16+KK)
152 write (*,FMT='(A16,E30.16)') " 17. edif/v=", UVAR(I,17+KK)
153 write (*,FMT='(A16,E30.16)') " 18. p=", UVAR(I,18+KK)
154 write (*,FMT='(A16,E30.16)') " 19. epx=", UVAR(I,19+KK)
155 write (*,FMT='(A16,E30.16)') " . mu=", UVAR(I,9+KK)/RHO0(ITRIMAT) - ONE
156 VAR = VAR + UVAR(I,11+KK)
157 ENDDO
158
159 RETURN
my_real
#define my_real
Definition cppsort.cpp:32
i22bufbric_mod
Definition cut-cell-search_mod.F:95
i22bufbric_mod::brick_list
type(brick_entity), dimension(:,:), allocatable, target brick_list
Definition cut-cell-search_mod.F:105