OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
hm_read_prop26.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_prop26 ../starter/source/properties/spring/hm_read_prop26.F
25!||--- called by ------------------------------------------------------
26!|| hm_read_properties ../starter/source/properties/hm_read_properties.F
27!||--- calls -----------------------------------------------------
28!|| ancmsg ../starter/source/output/message/message.F
29!|| fretitl2 ../starter/source/starter/freform.F
30!|| hm_get_float_array_index ../starter/source/devtools/hm_reader/hm_get_float_array_index.F
31!|| hm_get_float_array_index_dim ../starter/source/devtools/hm_reader/hm_get_float_array_index_dim.F
32!|| hm_get_floatv ../starter/source/devtools/hm_reader/hm_get_floatv.F
33!|| hm_get_floatv_dim ../starter/source/devtools/hm_reader/hm_get_floatv_dim.F
34!|| hm_get_int_array_index ../starter/source/devtools/hm_reader/hm_get_int_array_index.F
35!|| hm_get_intv ../starter/source/devtools/hm_reader/hm_get_intv.F
36!|| hm_option_is_encrypted ../starter/source/devtools/hm_reader/hm_option_is_encrypted.f
37!||--- uses -----------------------------------------------------
38!|| elbuftag_mod ../starter/share/modules1/elbuftag_mod.F
39!|| message_mod ../starter/share/message_module/message_mod.F
40!|| submodel_mod ../starter/share/modules1/submodel_mod.F
41!||====================================================================
42 SUBROUTINE hm_read_prop26(GEO , IGEO , UNITAB , IG , IGTYP,
43 . PROP_TAG, IDTITL, LSUBMODEL)
44C-----------------------------------------------
45 USE unitab_mod
46 USE elbuftag_mod
47 USE message_mod
48 USE submodel_mod
49 USE names_and_titles_mod , ONLY : nchartitle
50C-----------------------------------------------
51C I m p l i c i t T y p e s
52C-----------------------------------------------
53#include "implicit_f.inc"
54C-----------------------------------------------
55C C o m m o n B l o c k s
56C-----------------------------------------------
57#include "scr17_c.inc"
58#include "units_c.inc"
59#include "param_c.inc"
60#include "tablen_c.inc"
61C-----------------------------------------------
62C D u m m y A r g u m e n t s
63C-----------------------------------------------
64 TYPE (UNIT_TYPE_),INTENT(IN) ::UNITAB
65 INTEGER IGEO(NPROPGI),IG,IGTYP,FLAG_FMT
66C REAL
67 my_real
68 . geo(npropg)
69 CHARACTER(LEN=NCHARTITLE)::IDTITL
70 TYPE(prop_tag_) , DIMENSION(0:MAXPROP) :: PROP_TAG
71 TYPE(submodel_data),INTENT(IN)::LSUBMODEL(*)
72C-----------------------------------------------
73C L o c a l V a r i a b l e s
74C-----------------------------------------------
75 INTEGER I,J,NFUNC,NFUND,IFUN,IAD,ISENS,IFL,ILENG,IRTYP
76C REAL
77 my_real
78 . mass,kmax,dmax,xfac,yfac,rate,alpha,dmin,
79 . pun,yfac_dim,xfac_dim
80 CHARACTER(LEN=NCHARTITLE) :: TITR
81 LOGICAL IS_AVAILABLE, IS_ENCRYPTED
82C=======================================================================
83 pun = em01
84C
85 is_encrypted = .false.
86 is_available = .false.
87C
88 igeo( 1)=ig
89 igeo(11)=igtyp
90 geo(12) =igtyp+pun
91C
92C--------------------------------------------------
93C EXTRACT DATA (IS OPTION CRYPTED)
94C--------------------------------------------------
95 CALL hm_option_is_encrypted(is_encrypted)
96C--------------------------------------------------
97C EXTRACT DATAS (INTEGER VALUES)
98C--------------------------------------------------
99 CALL hm_get_intv('ISFLAG',ifl,is_available,lsubmodel)
100 CALL hm_get_intv('ISENSOR',isens,is_available,lsubmodel)
101 CALL hm_get_intv('Ileng',ileng,is_available,lsubmodel)
102
103 CALL hm_get_intv('nfunc',NFUNC,IS_AVAILABLE,LSUBMODEL)
104 CALL HM_GET_INTV('nraten',NFUND,IS_AVAILABLE,LSUBMODEL)
105 CALL HM_GET_FLOATV('dmin',DMIN,IS_AVAILABLE,LSUBMODEL,UNITAB)
106C--------------------------------------------------
107C EXTRACT DATAS (REAL VALUES)
108C--------------------------------------------------
109 CALL HM_GET_FLOATV('m_coeff',MASS,IS_AVAILABLE,LSUBMODEL,UNITAB)
110 CALL HM_GET_FLOATV('scale',XFAC,IS_AVAILABLE,LSUBMODEL,UNITAB)
111 CALL HM_GET_FLOATV_DIM('scale',XFAC_DIM,IS_AVAILABLE,LSUBMODEL,UNITAB)
112 CALL HM_GET_FLOATV('stiff0',KMAX,IS_AVAILABLE,LSUBMODEL,UNITAB)
113 CALL HM_GET_FLOATV('dmax',DMAX,IS_AVAILABLE,LSUBMODEL,UNITAB)
114 CALL HM_GET_FLOATV('alpha1',ALPHA,IS_AVAILABLE,LSUBMODEL,UNITAB)
115C----
116 IRTYP = 7
117 CALL FRETITL2(TITR,IGEO(NPROPGI-LTITR+1),LTITR)
118c
119 ! Check if at least one loading curve is defined
120 IF (NFUNC <= 0) THEN
121 CALL ANCMSG(MSGID=2078,
122 . MSGTYPE=MSGERROR,
123 . ANMODE=ANINFO_BLIND_1,
124 . I1=IG,
125 . C1=TITR)
126 ENDIF
127c
128 IF (ALPHA == ZERO) ALPHA = ONE
129 IF (XFAC == ZERO) XFAC = ONE * XFAC_DIM
130 DMIN = -ABS(DMIN)
131 DMAX = ABS(DMAX)
132 IF (DMIN == ZERO) DMIN = -INFINITY
133 IF (DMAX == ZERO) DMAX = INFINITY
134 IF (ILENG == 1) XFAC = ONE
135c
136C--- Loading curves
137 IAD = 100
138 DO I = 1, NFUNC
139 CALL HM_GET_INT_ARRAY_INDEX('fun_load',IFUN,I,IS_AVAILABLE,LSUBMODEL)
140 CALL HM_GET_FLOAT_ARRAY_INDEX('scale_load',yfac,i,is_available,lsubmodel,unitab)
141 CALL hm_get_float_array_index('STRAINRATE_LOAD',rate,i,is_available,lsubmodel,unitab)
142 CALL hm_get_float_array_index_dim('SCALE_LOAD',yfac_dim,i,is_available,lsubmodel,unitab)
143C
144 IF (ifun <= 0) THEN
145 CALL ancmsg(msgid=862,
146 . msgtype=msgerror,
147 . anmode=aninfo_blind_1,
148 . i1=ig,
149 . c1=titr)
150 EXIT
151 ENDIF
152 IF(i > 1 .AND. rate < geo(iad+100+i-1)) THEN
153 CALL ancmsg(msgid=861,
154 . msgtype=msgerror,
155 . anmode=aninfo_blind_1,
156 . i1=ig,
157 . c1=titr)
158 EXIT
159 ENDIF
160 IF (yfac == zero) yfac = one * yfac_dim
161C
162 igeo(iad+i) = ifun
163 geo(iad+100+i) = rate
164 geo(iad+200+i) = yfac
165 ENDDO
166C--- Unloading curves
167 iad = 100+nfunc
168 ! -> defined by user
169 IF (nfund > 0) THEN
170 DO i = 1, nfund
171 CALL hm_get_int_array_index('FUN_UNLOAD',ifun,i,is_available,lsubmodel)
172 CALL hm_get_float_array_index('SCALE_UNLOAD',yfac,i,is_available,lsubmodel,unitab)
173 CALL hm_get_float_array_index('STRAINRATE_UNLOAD',rate,i,is_available,lsubmodel,unitab)
174 CALL hm_get_float_array_index_dim('SCALE_UNLOAD',yfac_dim,i,is_available,lsubmodel,unitab)
175C
176 IF (ifun <= 0) THEN
177 CALL ancmsg(msgid=862,
178 . msgtype=msgerror,
179 . anmode=aninfo_blind_1,
180 . i1=ig,
181 . c1=titr)
182 EXIT
183 ENDIF
184 IF(i > 1 .AND. rate < geo(iad+100+i-1)) THEN
185 CALL ancmsg(msgid=861,
186 . msgtype=msgerror,
187 . anmode=aninfo_blind_1,
188 . i1=ig,
189 . c1=titr)
190 EXIT
191 ENDIF
192 IF (yfac == zero) yfac = one * yfac_dim
193C
194 igeo(iad+i) = ifun
195 geo(iad+100+i) = rate
196 geo(iad+200+i) = yfac
197 ENDDO
198 ! -> defined by default
199 ELSE
200 CALL ancmsg(msgid=2079,
201 . msgtype=msgwarning,
202 . anmode=aninfo_blind_1,
203 . i1=ig,
204 . c1=titr)
205 nfund = nfunc
206 DO i = 1,nfund
207 igeo(iad+i) = igeo(100+i)
208 geo(iad+100+i) = geo(200+i)
209 geo(iad+200+i) = geo(300+i)
210 ENDDO
211 ENDIF
212C
213 igeo(20) = nfunc
214 igeo(21) = nfund
215 geo(1) = mass
216 geo(2) = kmax
217 geo(4) = alpha
218 geo(5) = xfac
219 geo(8) = irtyp + em20
220 geo(15) = dmin
221 geo(16) = dmax
222C
223 IF (mass < em15)THEN
224 CALL ancmsg(msgid=229,
225 . msgtype=msgerror,
226 . anmode=aninfo_blind_1,
227 . i1=ig,
228 . c1=titr)
229 ENDIF
230 IF (ifl == 0)THEN
231 igeo(3)=isens
232 ELSEIF (ifl == 1)THEN
233 igeo(3)=-isens
234 ELSEIF (ifl == 2)THEN
235 igeo(3)=isens
236 ENDIF
237 geo(80)=ifl
238 geo(93)=ileng
239C
240 IF(is_encrypted)THEN
241 WRITE(iout,1000)ig
242 ELSE
243 WRITE(iout,1500)ig,mass,kmax,nfunc,nfund,dmin,dmax,alpha,xfac,ileng
244 iad = 100
245 DO i=1,nfunc
246 WRITE(iout,1700) igeo(iad+i),geo(iad+200+i),geo(iad+100+i)
247 ENDDO
248 iad = 100+nfunc
249 DO i=1,nfund
250 WRITE(iout,1800) igeo(iad+i),geo(iad+200+i),geo(iad+100+i)
251 ENDDO
252 ENDIF
253C
254 prop_tag(igtyp)%G_EINT = 1
255 prop_tag(igtyp)%G_FOR = 1
256 prop_tag(igtyp)%G_LENGTH = 1
257 prop_tag(igtyp)%G_TOTDEPL = 1
258 prop_tag(igtyp)%G_FOREP = 1
259 prop_tag(igtyp)%G_DEP_IN_COMP = 1
260 prop_tag(igtyp)%G_POSX = 5 ! just temp - not really used -
261 prop_tag(igtyp)%G_LENGTH_ERR = 1
262 prop_tag(igtyp)%G_DV = 1
263 prop_tag(igtyp)%G_RUPTCRIT = 1
264C
265C-----------
266 RETURN
267C-----------
268 1000 FORMAT(
269 & 5x,'TABULATED ELASTO-PLASTIC SPRING PROPERTY SET'/,
270 & 5x,'-------------------'/,
271 & 5x,'PROPERTY SET NUMBER . . . . . . . . . .=',i10/,
272 & 5x,'CONFIDENTIAL DATA'//)
273 1500 FORMAT(
274 & 5x,'TABULATED ELASTIC SPRING PROPERTY SET'/,
275 & 5x,'-------------------------------------'/,
276 & 5x,'PROPERTY SET NUMBER . . . . . . . . . .=',i10/,
277 & 5x,'SPRING MASS . . . . . . . . . . . . . .=',1pg20.13/,
278 & 5x,'MAXIMUM STIFFNESS . . . . . . . . . . .=',1pg20.13/,
279 & 5x,'NUMBER OF LOADING CURVES . . . . . . .=',i10/,
280 & 5x,'NUMBER OF UNLOADING CURVES. . . . . . .=',i10/,
281 & 5x,'FAILURE DISPLACEMENT IN COMPRESSION . .=',1pg20.13/,
282 & 5x,'FAILURE DISPLACEMENT IN TENSION . . . .=',1pg20.13/,
283 & 5x,'STRAIN RATE FILTERING FACTOR . . . . .=',1pg20.13/,
284 & 5x,'ABSCISSA SCALE FACTOR . . . . .=',1pg20.13/,
285 & 5x,'UNIT LENGTH FLAG. . . . . . . . . . . .=',i10/,
286 & 5x,'IF=1 UNIT LENGTH MASS,STIFFNESS AND INPUT',/,
287 & 5x,' CURVE ARE STRAIN DEPENDING',/)
288 1700 FORMAT(
289 & 5x,'YIELD STRESS FUNCTION NUMBER . . . . . =',i10/
290 & 7x,'SCALE FACTOR. . . . . . . . . . . . . . =',1pg20.13/
291 & 7x,'STRAIN RATE . . . . . . . . . . . . . . =',1pg20.13)
292 1800 FORMAT(
293 & 5x,'UNLOADING FUNCTION NUMBER . . . . . . . =',i10/
294 & 7x,'SCALE FACTOR. . . . . . . . . . . . . . =',1pg20.13/
295 & 7x,'STRAIN RATE . . . . . . . . . . . . . . =',1pg20.13)
296c sensor and sensor flag not used
297C-----------
298 RETURN
299 END SUBROUTINE hm_read_prop26
my_real
#define my_real
Definition cppsort.cpp:32
alpha
#define alpha
Definition eval.h:35
hm_get_float_array_index
subroutine hm_get_float_array_index(name, rval, index, is_available, lsubmodel, unitab)
Definition hm_get_float_array_index.F:105
hm_get_float_array_index_dim
subroutine hm_get_float_array_index_dim(name, dim_fac, index, is_available, lsubmodel, unitab)
Definition hm_get_float_array_index_dim.F:47
hm_get_int_array_index
subroutine hm_get_int_array_index(name, ival, index, is_available, lsubmodel)
Definition hm_get_int_array_index.F:164
hm_get_intv
subroutine hm_get_intv(name, ival, is_available, lsubmodel)
Definition hm_get_intv.F:523
hm_option_is_encrypted
subroutine hm_option_is_encrypted(is_encrypted)
Definition hm_option_is_encrypted.F:284
hm_read_prop26
subroutine hm_read_prop26(geo, igeo, unitab, ig, igtyp, prop_tag, idtitl, lsubmodel)
Definition hm_read_prop26.F:44
elbuftag_mod
Definition elbuftag_mod.F:227
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
submodel_mod
Definition submodel_mod.F:712
unitab_mod
Definition unitab_mod.F:513
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
starter
program starter
Definition starter.F:39
elbuftag_mod::prop_tag_
Definition elbuftag_mod.F:374
submodel_mod::submodel_data
Definition submodel_mod.F:723