OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
ccoef3.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!|| ccoef3 ../engine/source/elements/shell/coque/ccoef3.F
25!||--- called by ------------------------------------------------------
26!|| cforc3 ../engine/source/elements/shell/coque/cforc3.F
27!|| cforc3_crk ../engine/source/elements/xfem/cforc3_crk.F
28!||====================================================================
29 SUBROUTINE ccoef3(
30 1 JFT ,JLT ,PM ,MAT ,GEO ,
31 2 PID ,OFF ,AREA ,STI ,STIR ,
32 3 SHF ,THK0 ,THK02 ,NU ,
33 4 G ,YM ,A11 ,A12 ,THK ,
34 5 SSP ,RHO ,H1 ,H2 ,H3 ,
35 6 VOL0 ,VOL00 ,ALPE ,GS ,MTN ,
36 7 ITHK ,ISMSTR ,NPT ,KFTS ,
37 8 SRH1 ,SRH2 ,SRH3 ,IGEO ,
38 9 A11R ,ISUBSTACK, PM_STACK )
39C-----------------------------------------------
40C I m p l i c i t T y p e s
41C-----------------------------------------------
42#include "implicit_f.inc"
43C-----------------------------------------------
44C G l o b a l P a r a m e t e r s
45C-----------------------------------------------
46#include "mvsiz_p.inc"
47C-----------------------------------------------
48C C o m m o n B l o c k s
49C-----------------------------------------------
50#include "param_c.inc"
51C-----------------------------------------------
52C D u m m y A r g u m e n t s
53C-----------------------------------------------
54 INTEGER JFT, JLT, MTN,ITHK,
55 . ISMSTR,NPT ,KFTS ,ISUBSTACK
56 INTEGER MAT(MVSIZ), PID(MVSIZ), IGEO(NPROPGI,*)
57C REAL
58 my_real
59 . GEO(NPROPG,*), PM(NPROPM,*), OFF(*), AREA(*),
60 . STI(*),STIR(*),SHF(*),THK0(*),THK02(*),THK(*),
61 . NU(*),G(*),YM(*),A11(*),A12(*),GS(*),
62 . H1(*),H2(*),H3(*),VOL0(*),SSP(*),RHO(*),ALPE(*),VOL00(*),
63 . srh1(*) ,srh2(*) ,srh3(*),a11r(*),pm_stack(20,*)
64C-----------------------------------------------
65C L o c a l V a r i a b l e s
66C-----------------------------------------------
67 INTEGER I,ISH,MX,IPID,J,IGTYP,IPGMAT,IGMAT
68C REAL
69 my_real
70 . FAC1(MVSIZ),FSH, FAC1TMP,
71 . RHO_1, YM_1, NU_1, G_1, A11_1, A12_1, SSP_1
72C=======================================================================
73 IF(ITHK>0.AND.ISMSTR/=3)then
74 DO i=jft,jlt
75 alpe(i)=one
76C --- A CORRIGER : VOL00(I) = THK0(I)*AREA0(I)
77 vol00(i) = thk0(i)*area(i)
78 thk0(i)=thk(i)
79 vol0(i) = thk0(i)*area(i)
80 thk02(i) = thk0(i)*thk0(i)
81 ENDDO
82 ELSE
83 DO i=jft,jlt
84 alpe(i)=one
85C --- A CORRIGER : VOL00(I) = THK0(I)*AREA0(I)
86 vol00(i) = thk0(i)*area(i)
87 vol0(i) = vol00(i)
88 thk02(i) = thk0(i)*thk0(i)
89 ENDDO
90 ENDIF
91C
92 igtyp = igeo(11,pid(1))
93 igmat = igeo(98,pid(1))
94 ipgmat = 700
95
96 IF(igtyp == 11 .AND. igmat > 0) THEN
97 mx = mat(jft)
98 DO i=jft,jlt
99 mx = pid(i)
100 rho(i) = geo(ipgmat +1 ,mx)
101 ym(i) = geo(ipgmat +2 ,mx)
102 nu(i) = geo(ipgmat +3 ,mx)
103 g(i) = geo(ipgmat +4 ,mx)
104 a11(i) = geo(ipgmat +5 ,mx)
105 a12(i) = geo(ipgmat +6 ,mx)
106 a11r(i)= geo(ipgmat +7 ,mx)
107 ssp(i) = geo(ipgmat +9 ,mx)
108 ENDDO
109 ELSEIF(igtyp == 52 .OR.
110 . ((igtyp == 17 .OR. igtyp == 51) .AND. igmat > 0))THEN
111 DO i=jft,jlt
112 rho(i) = pm_stack(1 ,isubstack)
113 ym(i) = pm_stack(2 ,isubstack)
114 nu(i) = pm_stack(3 ,isubstack)
115 g(i) = pm_stack(4 ,isubstack)
116 a11(i) = pm_stack(5 ,isubstack)
117 a12(i) = pm_stack(6 ,isubstack)
118 a11r(i)= pm_stack(7 ,isubstack)
119 ssp(i) = pm_stack(9 ,isubstack)
120 ENDDO
121
122 ELSE
123 mx = mat(jft)
124 rho_1=pm(1,mx)
125 ym_1 =pm(20,mx)
126 nu_1 =pm(21,mx)
127 g_1 =pm(22,mx)
128 a11_1 =pm(24,mx)
129 a12_1 =pm(25,mx)
130 ssp_1 =pm(27,mx)
131 DO i=jft,jlt
132 rho(i)=rho_1
133 ym(i) =ym_1
134 nu(i) =nu_1
135 g(i) =g_1
136 a11(i) =a11_1
137 a12(i) =a12_1
138 ssp(i) =ssp_1
139 ENDDO
140 ENDIF !igtyp
141C
142 DO i=jft,jlt
143 ipid=pid(i)
144 h1(i) =geo(13,ipid)
145 h2(i) =geo(14,ipid)
146 h3(i) =geo(15,ipid)
147 srh1(i)=geo(18,ipid)
148 srh2(i)=geo(19,ipid)
149 srh3(i)=geo(20,ipid)
150 ENDDO
151C
152 IF(npt==1) THEN
153 DO i=jft,jlt
154 shf(i)=zero
155 ENDDO
156 ELSE
157 DO i=jft,jlt
158C ds021 41i 26/9/00 +++
159 fac1tmp =2.*(1.+nu(i))*thk02(i)
160 ish =nint(geo(37,pid(i)))
161 fsh =geo(38,pid(i))
162 shf(i)=fsh*(1. - ish +
163 . ish*fac1tmp / (fsh*area(i)+fac1tmp) )
164 fac1(i) = fac1tmp
165C ds021 41i 26/9/00 ---
166 ENDDO
167 ENDIF
168 DO i=jft,jlt
169 gs(i)=g(i)*shf(i)
170 ENDDO
171C
172 RETURN
173 END
ccoef3
subroutine ccoef3(jft, jlt, pm, mat, geo, pid, off, area, sti, stir, shf, thk0, thk02, nu, g, ym, a11, a12, thk, ssp, rho, h1, h2, h3, vol0, vol00, alpe, gs, mtn, ithk, ismstr, npt, kfts, srh1, srh2, srh3, igeo, a11r, isubstack, pm_stack)
Definition ccoef3.F:39