OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
write_matparam.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!|| write_matparam ../engine/source/output/restart/write_matparam.f
25!||--- called by ------------------------------------------------------
26!|| wrrestp ../engine/source/output/restart/wrrestp.F
27!||--- calls -----------------------------------------------------
28!|| write_ale_rezoning_param ../engine/source/output/restart/write_ale_rezoning_param.F90
29!|| write_c_c ../common_source/tools/input_output/write_routtines.c
30!|| write_db ../common_source/tools/input_output/write_db.F
31!|| write_eosparam ../engine/source/output/restart/write_eosparam.F90
32!|| write_failparam ../engine/source/output/restart/write_failparam.F
33!|| write_i_c ../common_source/tools/input_output/write_routtines.c
34!|| write_mat_table ../engine/source/materials/tools/write_mat_table.F
35!|| write_thermparam ../engine/source/output/restart/write_thermparam.F90
36!|| write_viscparam ../engine/source/output/restart/write_viscparam.F
37!||--- uses -----------------------------------------------------
38!|| mat_elem_mod ../common_source/modules/mat_elem/mat_elem_mod.F90
39!|| names_and_titles_mod ../common_source/modules/names_and_titles_mod.F
40!|| write_ale_rezoning_param_mod ../engine/source/output/restart/write_ale_rezoning_param.F90
41!|| write_therpmaram_mod ../engine/source/output/restart/write_thermparam.F90
42!||====================================================================
43 SUBROUTINE write_matparam(MAT_ELEM)
44C-----------------------------------------------
45C M o d u l e s
46C-----------------------------------------------
47 USE mat_elem_mod
48 USE names_and_titles_mod
49 use write_therpmaram_mod
50 use write_ale_rezoning_param_mod , only : write_ale_rezoning_param
51C-----------------------------------------------
52C I m p l i c i t T y p e s
53C-----------------------------------------------
54#include "implicit_f.inc"
55C-----------------------------------------------
56C D u m m y A r g u m e n t s
57C-----------------------------------------------
58 TYPE(mat_elem_) ,INTENT(IN) :: MAT_ELEM
59C-----------------------------------------------
60C L o c a l V a r i a b l e s
61C-----------------------------------------------
62 INTEGER :: I,NUMMAT,IMAT,NUPARAM,NIPARAM,NUMTABL,NFAIL,ITHERM
63 INTEGER :: IAD,NFIX,NFIXR,LEN,LENR,NMOD,MOD
64 INTEGER ,DIMENSION(NCHARTITLE) :: NAME
65 INTEGER ,DIMENSION(:) ,ALLOCATABLE :: IBUF
66 INTEGER :: NBSUBMAT
67 my_real ,DIMENSION(:) ,ALLOCATABLE :: rbuf
68C=======================================================================
69 nfix = 24
70 nummat = mat_elem%NUMMAT
71 len = nfix*nummat
72 ALLOCATE (ibuf(len) )
73
74 ! write integer parameters and flags
75 iad = 0
76 DO imat = 1,nummat
77 ibuf(iad+ 1) = mat_elem%MAT_PARAM(imat)%ILAW
78 ibuf(iad+ 2) = mat_elem%MAT_PARAM(imat)%MAT_ID
79 ibuf(iad+ 3) = mat_elem%MAT_PARAM(imat)%NUPARAM
80 ibuf(iad+ 4) = mat_elem%MAT_PARAM(imat)%NIPARAM
81 ibuf(iad+ 5) = mat_elem%MAT_PARAM(imat)%NFUNC
82 ibuf(iad+ 6) = mat_elem%MAT_PARAM(imat)%NTABLE
83 ibuf(iad+ 7) = mat_elem%MAT_PARAM(imat)%NSUBMAT
84 ibuf(iad+ 8) = mat_elem%MAT_PARAM(imat)%NFAIL
85 ibuf(iad+ 9) = mat_elem%MAT_PARAM(imat)%IVISC
86 ibuf(iad+10) = mat_elem%MAT_PARAM(imat)%IEOS
87 ibuf(iad+11) = mat_elem%MAT_PARAM(imat)%ITHERM
88 ibuf(iad+12) = mat_elem%MAT_PARAM(imat)%IEXPAN
89 ibuf(iad+13) = mat_elem%MAT_PARAM(imat)%IALE
90 ibuf(iad+14) = mat_elem%MAT_PARAM(imat)%ITURB
91 ibuf(iad+15) = mat_elem%MAT_PARAM(imat)%HEAT_FLAG
92 ibuf(iad+16) = mat_elem%MAT_PARAM(imat)%COMPRESSIBILITY
93 ibuf(iad+17) = mat_elem%MAT_PARAM(imat)%SMSTR
94 ibuf(iad+18) = mat_elem%MAT_PARAM(imat)%STRAIN_FORMULATION
95 ibuf(iad+19) = mat_elem%MAT_PARAM(imat)%IPRES
96 ibuf(iad+20) = mat_elem%MAT_PARAM(imat)%ORTHOTROPY
97 ibuf(iad+21) = mat_elem%MAT_PARAM(imat)%NLOC
98 ibuf(iad+22) = mat_elem%MAT_PARAM(imat)%IFAILWAVE
99 ibuf(iad+23) = mat_elem%MAT_PARAM(imat)%IXFEM
100 ibuf(iad+24) = mat_elem%MAT_PARAM(imat)%NMOD
101c
102 iad = iad + nfix
103 END DO
104
105 CALL write_i_c(ibuf,len)
106 DEALLOCATE(ibuf)
107
108 ! write real parameters
109 nfixr = 9
110 lenr = nfixr*nummat
111 iad = 0
112 ALLOCATE (rbuf(lenr) )
113!
114 DO imat = 1,nummat
115 rbuf(iad + 1) = mat_elem%MAT_PARAM(imat)%RHO
116 rbuf(iad + 2) = mat_elem%MAT_PARAM(imat)%RHO0
117 rbuf(iad + 3) = mat_elem%MAT_PARAM(imat)%YOUNG
118 rbuf(iad + 4) = mat_elem%MAT_PARAM(imat)%BULK
119 rbuf(iad + 5) = mat_elem%MAT_PARAM(imat)%SHEAR
120 rbuf(iad + 6) = mat_elem%MAT_PARAM(imat)%NU
121 rbuf(iad + 7) = mat_elem%MAT_PARAM(imat)%STIFF_CONTACT
122 rbuf(iad + 8) = mat_elem%MAT_PARAM(imat)%STIFF_HGLASS
123 rbuf(iad + 9) = mat_elem%MAT_PARAM(imat)%STIFF_TSTEP
124 iad = iad + nfixr
125 END DO
126 CALL write_db(rbuf ,lenr)
127 DEALLOCATE(rbuf)
128
129 ! write material title
130 DO imat = 1,nummat
131 DO i=1,nchartitle
132 name(i) = ichar(mat_elem%MAT_PARAM(imat)%TITLE(i:i))
133 END DO
134 CALL write_c_c(name,nchartitle)
135 END DO
136
137 ! write material parameter array
138 DO imat = 1,nummat
139 nuparam = mat_elem%MAT_PARAM(imat)%NUPARAM
140 niparam = mat_elem%MAT_PARAM(imat)%NIPARAM
141 IF (nuparam > 0) THEN
142 CALL write_db(mat_elem%MAT_PARAM(imat)%UPARAM ,nuparam)
143 END IF
144 IF (niparam > 0) THEN
145 CALL write_i_c(mat_elem%MAT_PARAM(imat)%IPARAM ,niparam)
146 END IF
147 END DO
148
149 ! write local 4D function tables
150 DO imat = 1,nummat
151 numtabl = mat_elem%MAT_PARAM(imat)%NTABLE
152 IF (numtabl > 0) THEN
153 CALL write_mat_table(mat_elem%MAT_PARAM(imat)%TABLE, numtabl)
154 END IF
155 END DO
156
157 ! write viscosity model parameters
158 DO imat = 1,nummat
159 IF (mat_elem%MAT_PARAM(imat)%IVISC > 0) THEN
160 CALL write_viscparam(mat_elem%MAT_PARAM(imat)%VISC)
161 END IF
162 END DO
163
164c
165c write thermal parameters
166c
167 DO imat = 1,nummat
168 itherm = mat_elem%mat_param(imat)%itherm
169 if (itherm > 0) then
170 CALL write_thermparam(mat_elem%MAT_PARAM(imat)%THERM)
171 END IF
172 END DO
173c
174c write parameters of failure models per material
175c
176 DO imat = 1,nummat
177 nfail = mat_elem%MAT_PARAM(imat)%NFAIL
178 IF (nfail > 0) THEN
179 DO i = 1,nfail
180 CALL write_failparam(mat_elem%MAT_PARAM(imat)%FAIL(i))
181 END DO
182 END IF
183 END DO
184
185 ! write damage modes
186 DO imat = 1,nummat
187 nmod = mat_elem%MAT_PARAM(imat)%NMOD
188 IF (nmod > 0) THEN
189 DO mod = 1,nmod
190 DO i=1,nchartitle
191 name(i) = ichar(mat_elem%MAT_PARAM(imat)%MODE(mod)(i:i))
192 END DO
193 CALL write_c_c(name,nchartitle)
194 ENDDO
195 ENDIF
196 ENDDO
197
198 ! Write multimaterial buffer
199 ALLOCATE (ibuf(21) ) !< mex is 21 submaterials
200 ALLOCATE (rbuf(21) )
201 DO imat = 1,nummat
202 nbsubmat = mat_elem%MAT_PARAM(imat)%MULTIMAT%NB
203 ibuf(1) = nbsubmat
204 CALL write_i_c(ibuf, 1)
205 IF(nbsubmat > 0)THEN
206 ibuf(1:nbsubmat) = mat_elem%MAT_PARAM(imat)%MULTIMAT%mid
207 rbuf(1:nbsubmat) = mat_elem%MAT_PARAM(imat)%MULTIMAT%vfrac
208 CALL write_i_c(ibuf, nbsubmat)
209 CALL write_db(rbuf, nbsubmat)
210 ENDIF
211 ENDDO
212 DEALLOCATE(ibuf)
213 DEALLOCATE(rbuf)
214
215 ! write eos model parameters
216 DO imat = 1,nummat
217 IF (mat_elem%MAT_PARAM(imat)%IEOS > 0) THEN
218 CALL write_eosparam(mat_elem%MAT_PARAM(imat)%EOS)
219 END IF
220 END DO
221
222 ! write ALE rezoning parameters
223 DO imat = 1,nummat
224 CALL write_ale_rezoning_param(mat_elem%MAT_PARAM(imat)%REZON)
225 END DO
226!-----------
227 RETURN
228 END
my_real
#define my_real
Definition cppsort.cpp:32
write_matparam
subroutine write_matparam(mat_elem)
Definition write_matparam.F:44
names_and_titles_mod
Definition names_and_titles_mod.F:997
names_and_titles_mod::nchartitle
integer, parameter nchartitle
Definition names_and_titles_mod.F:1003
write_failparam
subroutine write_failparam(fail, len)
Definition write_failparam.F:31
write_viscparam
subroutine write_viscparam(visc, len)
Definition write_viscparam.F:32
write_mat_table
subroutine write_mat_table(table, numtabl)
Definition write_mat_table.F:33
write_db
subroutine write_db(a, n)
Definition write_db.F:140
write_i_c
void write_i_c(int *w, int *len)
Definition write_routtines.c:645
write_c_c
void write_c_c(int *w, int *len)
Definition write_routtines.c:591