OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
hm_read_perturb_part_shell.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!|| hm_read_perturb_part_shell ../starter/source/general_controls/computation/hm_read_perturb_part_shell.F
25!||--- called by ------------------------------------------------------
26!|| hm_read_perturb ../starter/source/general_controls/computation/hm_read_perturb.F
27!||--- calls -----------------------------------------------------
28!|| ancmsg ../starter/source/output/message/message.F
29!|| hm_get_floatv ../starter/source/devtools/hm_reader/hm_get_floatv.F
30!|| hm_get_intv ../starter/source/devtools/hm_reader/hm_get_intv.F
31!|| hm_get_string ../starter/source/devtools/hm_reader/hm_get_string.F
32!|| hm_option_read_key ../starter/source/devtools/hm_reader/hm_option_read_key.F
33!|| hm_option_start ../starter/source/devtools/hm_reader/hm_option_start.F
34!|| plot_distrib ../starter/source/general_controls/computation/plot_distrib.F
35!||--- uses -----------------------------------------------------
36!|| format_mod ../starter/share/modules1/format_mod.F90
37!|| hm_option_read_mod ../starter/share/modules1/hm_option_read_mod.F
38!|| message_mod ../starter/share/message_module/message_mod.F
39!|| submodel_mod ../starter/share/modules1/submodel_mod.F
40!||====================================================================
41 SUBROUTINE hm_read_perturb_part_shell(
42 1 IPART ,RNOISE ,IPARTC ,IPARTG ,IGRPART ,
43 2 IPM ,PERTURB ,LSUBMODEL,UNITAB ,IDPERTURB,
44 3 INDEX ,INDEX_ITYP,NPART_SHELL,OFFS,QP_IPERTURB,
45 4 QP_RPERTURB)
46C-----------------------------------------------
47C M o d u l e s
48C-----------------------------------------------
49 USE message_mod
50 USE groupdef_mod
51 USE unitab_mod
52 USE submodel_mod
53 USE hm_option_read_mod
54 USE names_and_titles_mod , ONLY : nchartitle, ncharfield, ncharkey
55 USE format_mod , ONLY : lfield
56C-----------------------------------------------
57C I m p l i c i t T y p e s
58C-----------------------------------------------
59#include "implicit_f.inc"
60C-----------------------------------------------
61C C o m m o n B l o c k s
62C-----------------------------------------------
63#include "com04_c.inc"
64#include "scr17_c.inc"
65#include "units_c.inc"
66#include "param_c.inc"
67#include "sphcom.inc"
68C-----------------------------------------------
69C D u m m y A r g u m e n t s
70C-----------------------------------------------
71 TYPE (UNIT_TYPE_),INTENT(IN) ::UNITAB
72 my_real
73 . RNOISE(NPERTURB,NUMELC+NUMELTG+NUMELS+NUMSPH),
74 . QP_RPERTURB(NPERTURB,4)
75 INTEGER IPART(LIPART1,*),IPARTC(*),
76 . ipartg(*),ipm(npropmi,*),
77 . perturb(nperturb),offs,
78 . idperturb(nperturb),index(numelc+numeltg+numels+numsph),
79 . index_ityp(numelc+numeltg+numels+numsph),npart_shell,
80 . qp_iperturb(nperturb,6)
81 TYPE(submodel_data) LSUBMODEL(*)
82C-----------------------------------------------
83 TYPE (GROUP_) , DIMENSION(NGRPART) :: IGRPART
84C-----------------------------------------------
85C L o c a l V a r i a b l e s
86C-----------------------------------------------
87 INTEGER I,J,K,NUMA,I_METHOD,MAX_PART,
88 . CPT_PART,NB_RANDOM,I_SEED,DISTRIB(50),
89 . II,NB_INTERV,IGRPRT,N,IOK,SEED,SEED_RANDOM,
90 . ITYP,L,I_PERTURB_VAR,SIZEY,EMPTY
91 CHARACTER(LEN=NCHARTITLE) :: TITR
92 CHARACTER(LEN=NCHARKEY) :: KEY
93 CHARACTER MES*40
94 INTEGER, DIMENSION(:,:), ALLOCATABLE :: TAB_PART
95 INTEGER, DIMENSION(:), ALLOCATABLE :: A_SEED
96 INTEGER, DIMENSION(1:8) :: DT_SEED
97 my_real
98 . mean,sd,mean_input,sd_input,max_distrib,temp,min_value,
99 . max_value,interv,VALUE,max_value1,minval,maxval,bid
100 my_real, DIMENSION(:), ALLOCATABLE :: array
101 CHARACTER*100 CHAR(100)
102 CHARACTER*100 CHAR1(100)
103 CHARACTER*100 CHAR2
104 CHARACTER MESS*40
105 CHARACTER(LEN=NCHARFIELD)::CHVAR
106 LOGICAL IS_AVAILABLE
107C-----------------------------------------------
108C E x t e r n a l F u n c t i o n s
109C-----------------------------------------------
110 DATA mess/'PERTURBATION DEFINITION '/
111C=======================================================================
112C-----------------------------------------------
113C B e g i n n i n g o f s o u r c e
114C-----------------------------------------------
115 ! Initialization and allocation of tables
116 max_part = 0
117 ityp = 0
118 bid = 0
119 is_available = .false.
120 CALL hm_option_start('/PERTURB/PART/SHELL')
121 !----------------------------------------------------------------------
122 ! 1st Loop over /PERTURB/PART for computing table dimension
123 !----------------------------------------------------------------------
124 DO i=1+offs,npart_shell+offs
125C
126 ! Reading the option
127 titr = ''
128 CALL hm_option_read_key(lsubmodel,
129 . option_id = idperturb(i),
130 . option_titr = titr)
131C
132 ! Perturbation type
133 ityp = 1
134c
135 i_perturb_var = 0
136 cpt_part = 0
137 igrprt = 0
138 iok = 0
139 chvar = ''
140C
141 ! Reading the number of the Group of Part + the perturbed variable
142 CALL hm_get_intv('grpart_ID' ,igrprt ,is_available,lsubmodel)
143 CALL hm_get_string('chvar',chvar ,ncharfield,is_available)
144 IF (chvar(1:5) == 'thick' .OR. chvar(1:5) == 'THICK') i_perturb_var = 1
145c
146 empty = 1
147 DO k = 1,lfield
148 IF(chvar(k:k) /= ' ') empty = 0
149 ENDDO
150c
151 ! Checking the perturbed variable
152 IF (i_perturb_var /= 1 .AND. empty == 0) CALL ancmsg(msgid=1194,
153 . msgtype=msgerror,
154 . anmode=aninfo,
155 . i1=idperturb(i),
156 . c1=titr,
157 . c2=chvar)
158c
159 ! Checking Shell part group
160 IF (igrprt /= 0) THEN
161 DO n=1,ngrpart
162 IF (igrpart(n)%ID == igrprt) THEN
163 igrprt=n
164 iok = 1
165 ityp = 1
166 EXIT
167 END IF
168 END DO
169 ENDIF
170c
171 ! Saving the perturbation type
172 perturb(i) = ityp
173c
174 ! Error messages or counting
175 IF (iok == 0) THEN
176 CALL ancmsg(msgid=1137,
177 . msgtype=msgerror,
178 . anmode=aninfo,
179 . i1=idperturb(i),
180 . c1=titr,
181 . i2=igrprt,
182 . c2='GROUP OF PART')
183 ELSEIF(iok == 1)THEN
184 cpt_part = igrpart(igrprt)%NENTITY
185 ENDIF
186 max_part = max(max_part,cpt_part)
187 ENDDO
188 ! Allocation of tables
189 ALLOCATE(tab_part(nperturb,max_part))
190c
191 !----------------------------------------------------------------------
192 ! 2nd Loop over /PERTURB/PART for reading and computing perturbation
193 !----------------------------------------------------------------------
194 CALL hm_option_start('/PERTURB/PART/SHELL')
195 DO i = 1+offs,npart_shell+offs
196C
197 ! Resetting index tables
198 index(1:(numelc+numeltg+numels+numsph)) = 0
199 index_ityp(1:numelc+numeltg+numels+numsph) = 0
200C
201 ! Reading the option
202 titr = ''
203 CALL hm_option_read_key(lsubmodel,
204 . option_id = idperturb(i),
205 . option_titr = titr)
206C
207 ! Perturbation type
208 ityp = 1
209C
210 ! Reading the card
211 CALL hm_get_floatv('F_Mean' ,mean ,is_available, lsubmodel, unitab)
212 CALL hm_get_floatv('Deviation' ,sd ,is_available, lsubmodel, unitab)
213 CALL hm_get_floatv('Min_cut' ,minval ,is_available, lsubmodel, unitab)
214 CALL hm_get_floatv('Max_cut' ,maxval ,is_available, lsubmodel, unitab)
215 CALL hm_get_intv('Seed' ,seed ,is_available, lsubmodel)
216 CALL hm_get_intv('Idistri' ,i_method,is_available, lsubmodel)
217c
218 ! Default value
219 IF(i_method == 0) i_method = 2
220 IF(minval == zero .AND. maxval == zero) THEN
221 IF(i_method == 1) THEN
222 ELSEIF(i_method == 2)THEN
223 minval = -ep30
224 maxval = ep30
225 ENDIF
226 ENDIF
227 sd_input = sd
228 mean_input = mean
229c
230 ! QAPRINT table
231 qp_iperturb(i,1) = idperturb(i)
232 qp_iperturb(i,2) = ityp
233 qp_iperturb(i,3) = seed
234 qp_iperturb(i,4) = i_method
235 qp_rperturb(i,1) = mean
236 qp_rperturb(i,2) = sd
237 qp_rperturb(i,3) = minval
238 qp_rperturb(i,4) = maxval
239c
240 ! Initialization flag and counter
241 cpt_part = 0
242 igrprt = 0
243 iok = 0
244C
245 ! Reading the number of the Group of Part + the perturbed variable
246 CALL hm_get_intv('grpart_ID',igrprt,is_available,lsubmodel)
247 qp_iperturb(i,5) = igrprt
248 CALL hm_get_string('chvar',chvar,ncharfield,is_available)
249 IF (chvar(1:5) == 'thick' .OR. chvar(1:5) == 'THICK') qp_iperturb(i,6) = 1
250c
251 ! checking shell part group
252 IF (igrprt /= 0) THEN
253 DO n=1,ngrpart
254 IF (igrpart(n)%ID == igrprt) THEN
255 igrprt = n
256 iok = 1
257 ityp = 1
258 EXIT
259 END IF
260 END DO
261 ENDIF
262c
263 ! Tag the parts
264 IF (iok == 1) THEN
265 cpt_part = 0
266 DO j=1,igrpart(igrprt)%NENTITY
267 cpt_part = cpt_part + 1
268 numa = igrpart(igrprt)%ENTITY(j)
269 tab_part(i,cpt_part) = numa
270 END DO
271 ENDIF
272c
273 ! Printing out the information
274 IF(i_method == 2) THEN
275 WRITE (iout,1000)
276 . idperturb(i),'GAUSSIAN',mean_input,sd_input,seed
277 WRITE (iout,'(10I10)') ipart(4,tab_part(i,1:cpt_part))
278 WRITE(iout,*) ' '
279 WRITE(iout,*) ' '
280 ELSEIF(i_method == 1) THEN
281 WRITE (iout,1100)
282 . idperturb(i),'RANDOM',seed
283 WRITE (iout,'(10I10)') ipart(4,tab_part(i,1:cpt_part))
284 WRITE(iout,*) ' '
285 WRITE(iout,*) ' '
286 ENDIF
287c
288 ! Filling the index table
289 nb_random = 0
290 DO ii=1,numelc
291 DO k=1,cpt_part
292 IF (ipartc(ii) == tab_part(i,k)) THEN
293 nb_random = nb_random + 1
294 index(nb_random) = ii
295 index_ityp(nb_random) = 3
296 ENDIF
297 ENDDO
298 ENDDO
299 DO ii=1,numeltg
300 DO k=1,cpt_part
301 IF(ipartg(ii) == tab_part(i,k)) THEN
302 nb_random = nb_random + 1
303 index(nb_random) = ii
304 index_ityp(nb_random) = 7
305 ENDIF
306 ENDDO
307 ENDDO
308C
309 ! Set up random seed
310 IF( seed == 0 )THEN
311 CALL random_seed(size=i_seed)
312 ALLOCATE(a_seed(1:i_seed))
313 CALL random_seed(get=a_seed)
314 CALL date_and_time(values=dt_seed)
315 a_seed=dt_seed(8)*dt_seed(7)*dt_seed(6)
316 seed=dt_seed(8)*dt_seed(7)*dt_seed(6)
317 CALL random_seed(put=a_seed)
318 seed_random = 1
319 DEALLOCATE(a_seed)
320 ELSE
321 CALL random_seed(size=i_seed)
322 ALLOCATE(a_seed(1:i_seed))
323 a_seed=seed
324 CALL random_seed(put=a_seed)
325 seed_random = 0
326 DEALLOCATE(a_seed)
327 ENDIF
328C
329 ! Build uniform distribution
330 char=''
331 char1=''
332 char2=''
333 distrib(1:50) = 0
334 ALLOCATE(array(nb_random+2))
335 CALL random_number(array)
336C
337 ! Build normal distribution
338 max_value = -ep30
339 min_value = ep30
340 IF ( i_method == 2) THEN
341 DO ii = 1, nb_random+1, 2
342 temp = sd * sqrt(-2.0*log(array(ii))) * cos(2*pi*array(ii+1)) +
343 . mean
344 array(ii+1) =
345 . sd * sqrt(-2.0*log(array(ii))) * sin(2*pi*array(ii+1)) + mean
346 array(ii) = temp
347 END DO
348 DO ii = 1, nb_random
349 array(ii) = max(min(maxval,array(ii)),minval)
350 max_value = max(array(ii),max_value)
351 min_value = min(array(ii),min_value)
352 END DO
353 ELSEIF(i_method == 1)THEN
354 DO ii = 1, nb_random
355 array(ii) = array(ii)*(maxval-minval)+minval
356 max_value = max(array(ii),max_value)
357 min_value = min(array(ii),min_value)
358 END DO
359 ENDIF
360c
361 ! Filling RNOISE table
362 DO ii = 1, nb_random
363 IF (index_ityp(ii) == 3) THEN
364 rnoise(i,index(ii)) = array(ii)
365 ELSEIF (index_ityp(ii) == 7) THEN
366 rnoise(i,index(ii)+numelc) = array(ii)
367 ENDIF
368 ENDDO
369c
370 ! Check mean and standard deviation
371 mean = sum(array)/nb_random
372 sd = sqrt(sum((array - mean)**2)/nb_random)
373c
374 ! Plot the normal distribution
375 IF(i_method == 2) THEN
376 max_distrib = one /(sd*sqrt(two * pi))
377 ELSEIF(i_method == 1) THEN
378 max_distrib = one /(max_value-min_value)
379 ENDIF
380 WRITE (iout,3000)
381 WRITE(iout,*) ' '
382 nb_interv = 50
383 sizey = 20
384 IF (minval /= -ep30 .AND. maxval /= ep30)THEN
385 min_value = minval
386 max_value = maxval
387 ENDIF
388 CALL plot_distrib( array,nb_random, nb_interv,sizey,min_value,
389 . max_value,max_distrib,'#')
390 IF(i_method == 2) THEN
391 WRITE (iout,2000) mean,sd
392 ELSEIF(i_method == 1) THEN
393 WRITE (iout,2050) mean
394 ENDIF
395 IF(seed_random == 1) WRITE (iout,2100) seed
396 WRITE(iout,*) ' '
397 WRITE(iout,*) ' '
398 IF (ALLOCATED(array)) DEALLOCATE(array)
399 ENDDO
400 DEALLOCATE(tab_part)
401C-------------------------------------------------------------
402 1000 FORMAT(/' PERTURBATION ID',i10/
403 + ' ---------------'/
404 + ' TYPE . . . . . . . . . . . . . . .',a/
405 + ' INPUT MEAN VALUE . . . . . . . . .',1pg20.13/
406 + ' INPUT STANDARD DEVIATION . . . . .',1pg20.13/
407 + ' INPUT SEED VALUE . . . . . . . . .',i10/
408 + ' SHELL THICKNESSES, PARTS:')
409 1100 FORMAT(/' PERTURBATION ID',i10/
410 + ' ---------------'/
411 + ' TYPE . . . . . . . . . . . . . . .',a/
412 + ' INPUT SEED VALUE . . . . . . . . .',i10/
413 + ' SHELL THICKNESSES, PARTS:')
414C-------------------------------------------------------------
415 2000 FORMAT(/
416 + ' GENERATED MEAN VALUE . . . . . . .',1pg20.13/
417 + ' GENERATED STANDARD DEVIATION . . .',1pg20.13)
418 2050 FORMAT(/
419 + ' GENERATED MEAN VALUE . . . . . . .',1pg20.13)
420 2100 FORMAT(/
421 + ' GENERATED SEED VALUE . . . . . . .',i10/)
422C-------------------------------------------------------------
423 3000 FORMAT(/
424 + ' DISTRIBUTION OF SCALE FACTORS APPLIED TO THICKNESSES OF SHELLS')
425C------------------------------
426C------------------------------
427 END
hm_get_floatv
subroutine hm_get_floatv(name, rval, is_available, lsubmodel, unitab)
Definition hm_get_floatv.F:415
hm_get_intv
subroutine hm_get_intv(name, ival, is_available, lsubmodel)
Definition hm_get_intv.F:523
hm_get_string
subroutine hm_get_string(name, sval, size, is_available)
Definition hm_get_string.F:76
hm_option_start
subroutine hm_option_start(entity_type)
Definition hm_option_start.F:251
hm_read_perturb_part_shell
subroutine hm_read_perturb_part_shell(ipart, rnoise, ipartc, ipartg, igrpart, ipm, perturb, lsubmodel, unitab, idperturb, index, index_ityp, npart_shell, offs, qp_iperturb, qp_rperturb)
Definition hm_read_perturb_part_shell.F:46
hm_option_read_mod::hm_option_read_key
Definition hm_option_read_mod.F:410
min
#define min(a, b)
Definition macros.h:20
max
#define max(a, b)
Definition macros.h:21
groupdef_mod
Definition groupdef_mod.F:662
hm_option_read_mod
Definition hm_option_read_mod.F:408
message_mod
Definition message_mod.F:1249
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
names_and_titles_mod::ncharkey
integer, parameter ncharkey
Definition names_and_titles_mod.F:1000
names_and_titles_mod::ncharfield
integer, parameter ncharfield
Definition names_and_titles_mod.F:1005
submodel_mod
Definition submodel_mod.F:712
unitab_mod
Definition unitab_mod.F:513
plot_distrib
subroutine plot_distrib(array, s_array, nb_interv, sizey, x_minvalue, x_maxvalue, y_maxvalue, ecrit)
Definition plot_distrib.F:34
ancmsg
subroutine ancmsg(msgid, msgtype, anmode, i1, i2, i3, i4, i5, i6, i7, i8, i9, i10, i11, i12, i13, i14, i15, i16, i17, i18, i19, i20, r1, r2, r3, r4, r5, r6, r7, r8, r9, c1, c2, c3, c4, c5, c6, c7, c8, c9, prmode)
Definition message.F:889
submodel_mod::submodel_data
Definition submodel_mod.F:723