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 31 of file write_buf_law51.F.

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