OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
h3d_fld_strain.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/.
23C
24!||====================================================================
25!|| h3d_fld_strain ../engine/source/output/h3d/h3d_results/h3d_fld_strain.F
26!||--- called by ------------------------------------------------------
27!|| h3d_skin_scalar ../engine/source/output/h3d/h3d_results/h3d_skin_scalar.F
28!||--- calls -----------------------------------------------------
29!|| roto_sig2d ../engine/source/output/h3d/h3d_results/h3d_skin_tensor.F
30!|| tsh_dir2 ../engine/source/output/h3d/h3d_results/h3d_skin_tensor.F
31!||--- uses -----------------------------------------------------
32!|| elbufdef_mod ../common_source/modules/mat_elem/elbufdef_mod.F90
33!|| element_mod ../common_source/modules/elements/element_mod.F90
34!||====================================================================
35 SUBROUTINE h3d_fld_strain(ELBUF_TAB,X ,IXS ,
36 . JHBE,MLWI,ILAY,KCVT,IOR_TSH,
37 . ICSTR,NPTR,NPTS,NEL,F_EXP,EVAR )
38C-----------------------------------------------
39C M o d u l e s
40C-----------------------------------------------
41 USE elbufdef_mod
42 use element_mod , only : nixs
43C-----------------------------------------------
44C I m p l i c i t T y p e s
45C-----------------------------------------------
46#include "implicit_f.inc"
47C-----------------------------------------------
48C C o m m o n B l o c k s
49C-----------------------------------------------
50#include "mvsiz_p.inc"
51#include "com04_c.inc"
52C-----------------------------------------------
53C D u m m y A r g u m e n t s
54C-----------------------------------------------
55 INTEGER , INTENT(IN) :: JHBE,ILAY,MLWI,KCVT,IOR_TSH,
56 . NPTR,NPTS,ICSTR,NEL
57 INTEGER ,DIMENSION(NIXS,NUMELS), INTENT(IN) :: IXS
58 my_real , INTENT(IN) :: F_EXP
59 my_real ,DIMENSION(3,MVSIZ), INTENT(OUT) :: evar
60 my_real ,DIMENSION(3,NUMNOD), INTENT(IN) :: x
61 TYPE (ELBUF_STRUCT_), TARGET :: ELBUF_TAB
62C-----------------------------------------------
63C L o c a l V a r i a b l e s
64C-----------------------------------------------
65 my_real
66 . dir(mvsiz,2),dirb(mvsiz,2)
67 INTEGER I,IR,IS,IT,JJ(4)
68
69 TYPE(g_bufel_) ,POINTER :: GBUF
70 TYPE(L_BUFEL_) ,POINTER :: LBUF
71C-----------------------------------------------
72 evar(1:3,1:nel) = zero
73 it = 1
74 ir = 0
75 is = 0
76 DO i=1,4
77 jj(i) = nel*(i-1)
78 ENDDO
79 gbuf => elbuf_tab%GBUF
80 IF (jhbe==15) THEN
81 ir = 1
82 is = 1
83 lbuf => elbuf_tab%BUFLY(ilay)%LBUF(ir,is,it)
84 IF (mlwi == 12 .OR. mlwi == 14) THEN
85 DO i=1,nel
86 evar(1:2,i) = lbuf%EPE(jj(1:2) + i)
87 evar(3,i) = half*lbuf%EPE(jj(4) + i)
88 ENDDO
89 ELSEIF (mlwi /= 49 ) THEN
90 DO i=1,nel
91 evar(1:2,i) = lbuf%STRA(jj(1:2) + i)
92 evar(3,i) = half*lbuf%STRA(jj(4) + i)
93 ENDDO
94 END IF
95C------
96 ELSE ! 14,16
97 IF (mlwi == 12 .OR. mlwi == 14) THEN
98 DO ir=1,nptr
99 DO is=1,npts
100 lbuf => elbuf_tab%BUFLY(ilay)%LBUF(ir,is,it)
101 DO i=1,nel
102 evar(1:2,i) = evar(1:2,i)+lbuf%EPE(jj(1:2) + i)
103 evar(3,i) = evar(3,i)+half*lbuf%EPE(jj(4) + i)
104 ENDDO
105 END DO
106 END DO
107 ELSEIF (mlwi /= 49 ) THEN
108 DO ir=1,nptr
109 DO is=1,npts
110 lbuf => elbuf_tab%BUFLY(ilay)%LBUF(ir,is,it)
111 DO i=1,nel
112 evar(1:2,i) = evar(1:2,i)+lbuf%STRA(jj(1:2) + i)
113 evar(3,i) = evar(3,i)+ half*lbuf%STRA(jj(4) + i)
114 ENDDO
115 ENDDO
116 END DO
117 END IF
118 END IF
119 evar(1:3,1:nel) = f_exp*evar(1:3,1:nel)
120C------
121 IF (kcvt==2) THEN
122 IF(ior_tsh==1)THEN
123 DO i=1,nel
124 dir(i,1:2)= gbuf%GAMA(jj(1:2) + i)
125 ENDDO
126 ELSEIF(ior_tsh==2)THEN
127 IF(jhbe==14)THEN
128 ir = 1
129 is = 1
130 END IF
131 lbuf => elbuf_tab%BUFLY(ilay)%LBUF(ir,is,it)
132 DO i=1,nel
133 dir(i,1:2)= lbuf%GAMA(jj(1:2) + i)
134 ENDDO
135 END IF
136 CALL tsh_dir2(x,ixs,dir,dirb,icstr,nel)
137 CALL roto_sig2d(1,nel,evar,dirb)
138 END IF !(KCVT==2) THEN
139C
140 RETURN
141 END SUBROUTINE h3d_fld_strain
h3d_fld_strain
subroutine h3d_fld_strain(elbuf_tab, x, ixs, jhbe, mlwi, ilay, kcvt, ior_tsh, icstr, nptr, npts, nel, f_exp, evar)
Definition h3d_fld_strain.F:38
tsh_dir2
subroutine tsh_dir2(x, ixs, dir, dirb, icstr, nel)
Definition h3d_skin_tensor.F:610
roto_sig2d
subroutine roto_sig2d(jft, jlt, sig, dir)
Definition h3d_skin_tensor.F:779