OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
r13ke3.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!|| r13ke3 ../engine/source/elements/spring/r13ke3.f
25!||--- called by ------------------------------------------------------
26!|| imp_glob_k ../engine/source/implicit/imp_glob_k.F
27!|| imp_glob_k0 ../engine/source/implicit/imp_glob_k.F
28!||--- calls -----------------------------------------------------
29!|| assem_p ../engine/source/implicit/assem_p.F
30!|| peoff ../engine/source/elements/beam/peoff.F
31!|| pkeg3 ../engine/source/elements/beam/pke3.F
32!|| r13mat3 ../engine/source/elements/spring/r13mat3.F
33!|| r13sumg3 ../engine/source/elements/spring/r13sumg3.F
34!|| r4coork3 ../engine/source/elements/spring/r4coork3.F
35!||--- uses -----------------------------------------------------
36!|| elbufdef_mod ../common_source/modules/mat_elem/elbufdef_mod.F90
37!||====================================================================
38 SUBROUTINE r13ke3(JFT ,JLT ,NEL ,MTN ,PM ,
39 1 GEO ,IXR ,X ,ELBUF_TAB,NPF ,
40 2 TF ,SKEW ,OFFSET,FR_WAVE,IKGEO ,IGEO,
41 3 ETAG , IDDL ,NDOF ,K_DIAG ,K_LT ,
42 4 IADK ,JDIK )
43
44C-----------------------------------------------
45C M o d u l e s
46C-----------------------------------------------
47 USE elbufdef_mod
48C-----------------------------------------------
49C I m p l i c i t T y p e s
50C-----------------------------------------------
51#include "implicit_f.inc"
52C-----------------------------------------------
53C G l o b a l P a r a m e t e r s
54C-----------------------------------------------
55#include "mvsiz_p.inc"
56C-----------------------------------------------
57C C o m m o n B l o c k s
58C-----------------------------------------------
59#include "com04_c.inc"
60#include "param_c.inc"
61C-----------------------------------------------
62C D u m m y A r g u m e n t s
63C-----------------------------------------------
64 INTEGER IXR(NIXR,*), NPF(*),JFT,JLT,OFFSET,NEL,MTN,IKGEO
65 INTEGER IGEO(*),
66 . ETAG(*),IDDL(*) ,NDOF(*) ,IADK(*) ,JDIK(*)
67
68C REAL
69 my_real
70 . geo(npropg,*),x(*),tf(*),skew(lskew,*),
71 . fr_wave(*),pm(*),ke11(36,mvsiz),ke12(36,mvsiz),ke22(36,mvsiz),
72 . off(mvsiz),k_diag(*) ,k_lt(*)
73C
74 TYPE (elbuf_struct_), TARGET:: elbuf_tab
75C-----------------------------------------------
76C L o c a l V a r i a b l e s
77C-----------------------------------------------
78C REAL
79 my_real
80 . al(mvsiz),kx(mvsiz),ky(mvsiz) ,kz(mvsiz),
81 . mx(mvsiz) ,my(mvsiz) ,mz(mvsiz),
82 . r11(mvsiz),r12(mvsiz),r13(mvsiz),
83 . r21(mvsiz),r22(mvsiz),r23(mvsiz),
84 . r31(mvsiz),r32(mvsiz),r33(mvsiz)
85 INTEGER PID(MVSIZ),NGL(MVSIZ)
86 INTEGER NB1, NB2, NB3, NB4, NB5, NB6, NB7, NB8, NB9, NB10,
87 . NB11, NB2A, NB2B, NB4A, NB4B, NB5A, NB5B, NB6A, NB6B, NB7A,
88 . NB7B, NB8A, NB8B, NB9A, NB9B, NB10A, NB10B, NB11A, NB11B,
89 . nb12, nb12a, nb12b, nb13, nb13a, nb13b, nb14, nel3, nel4,
90 . nb15, nb15a, nb15b, nb15c, nb15d, nb15e, nb15f, nb16,
91 . nb8c, nb8d,offset3,offset4,nbfi,igtyp,i0,i,ii(3)
92C
93 TYPE(g_bufel_) ,POINTER :: GBUF
94C-----------------------------------------------
95 GBUF => elbuf_tab%GBUF
96C
97C
98 DO i=1,3
99 ii(i) = (i-1)*nel + 1
100 ENDDO
101C
102C I0 = IXR(1,1)
103C IGTYP = NINT(GEO(12,I0))
104c NB1 =1
105c NB2 =NB1 + NEL
106c NB1 =NB1 + OFFSET
107c NEL3 = 3 * NEL
108c NEL4 = 4 * NEL
109c OFFSET3 = 3 * OFFSET
110c OFFSET4 = 4 * OFFSET
111C
112C=======================================================================
113C IF(IGTYP==13)THEN
114C=======================================================================
115c NB3 =NB2 + NEL3
116c NB4 =NB3 + NEL
117c NB5 =NB4 + NEL3
118c IF (SCODVER<101) THEN
119c NB6 =NB5 + NEL3
120c ELSE
121c NB6 =NB5 + 2*NEL3
122c END IF
123c NB7 =NB6 + NEL3
124c NB8 =NB7 + NEL3
125c NB9 =NB8 + NEL3
126c NB10=NB9 + NEL3
127c NB11=NB10 + NEL3
128c NB12=NB11 + NEL3
129c NB13=NB12 + NEL3
130c NB14=NB13 + NEL3
131c IF (SCODVER<101) THEN
132c NB15 = NB14 + NEL3
133c ELSE
134c NB15 = NB14 + 2*NEL3
135c END IF
136c NBFI = NB15 + NEL4*6
137c NB2A=NB2 + NEL
138c NB2B=NB2A + NEL
139c NB4A=NB4 + NEL
140c NB4B=NB4A + NEL
141c NB5A=NB5 + NEL
142c NB5B=NB5A + NEL
143c NB6A=NB6 + NEL
144c NB6B=NB6A + NEL
145c NB7A=NB7 + NEL
146c NB7B=NB7A + NEL
147c NB8A=NB8 + NEL
148c NB8B=NB8A + NEL
149c NB9A=NB9 + NEL
150c NB9B=NB9A + NEL
151c NB10A=NB10 + NEL
152c NB10B=NB10A + NEL
153c NB11A=NB11 + NEL
154c NB11B=NB11A+ NEL
155c NB12A=NB12 + NEL
156c NB12B=NB12A+ NEL
157c NB13A=NB13 + NEL
158c NB13B=NB13A+ NEL
159c NB15A=NB15 + NEL4
160c NB15B=NB15A+ NEL4
161c NB15C=NB15B+ NEL4
162c NB15D=NB15C+ NEL4
163c NB15E=NB15D+ NEL4
164c NB16 =NB15E+ NEL4
165c NBFI =NB16 + 6*NEL
166C
167c NB2 =NB2 + OFFSET
168c NB3 =NB3 + OFFSET
169c NB4 =NB4 + OFFSET
170c NB5 =NB5 + OFFSET
171c NB6 =NB6 + OFFSET
172c NB7 =NB7 + OFFSET
173c NB8 =NB8 + OFFSET
174c NB9 =NB9 + OFFSET
175c NB10=NB10 + OFFSET
176c NB11=NB11 + OFFSET
177c NB12=NB12 + OFFSET
178c NB13=NB13 + OFFSET
179c NB14=NB14 + OFFSET3
180c NB2A=NB2A + OFFSET
181c NB2B=NB2B + OFFSET
182c NB4A=NB4A + OFFSET
183c NB4B=NB4B + OFFSET
184c NB5A=NB5A + OFFSET
185c NB5B=NB5B + OFFSET
186c NB6A=NB6A + OFFSET
187c NB6B=NB6B + OFFSET
188c NB7A=NB7A + OFFSET
189c NB7B=NB7B + OFFSET
190c NB8A=NB8A + OFFSET
191c NB8B=NB8B + OFFSET
192c NB9A=NB9A + OFFSET
193c NB9B=NB9B + OFFSET
194c NB10A=NB10A + OFFSET
195c NB10B=NB10B + OFFSET
196c NB11A=NB11A+ OFFSET
197c NB11B=NB11B+ OFFSET
198c NB12A=NB12A+ OFFSET
199c NB12B=NB12B+ OFFSET
200c NB13A=NB13A+ OFFSET
201c NB13B=NB13B+ OFFSET
202c NB15 =NB15 + OFFSET4
203c NB15A=NB15A + OFFSET4
204c NB15B=NB15B + OFFSET4
205c NB15C=NB15C + OFFSET4
206c NB15D=NB15D + OFFSET4
207c NB15E=NB15E + OFFSET4
208C FX,FY,FZ -> NB2,NB2A,NB2B E -> NB3
209C DX,DY,DZ -> NB4,NB4A,NB4B OFF -> NB1
210C DPX,DPY,DPZ -> NB6,NB6A,NB6B DPX2,DPY2,DPZ2 -> NB7,NB7A,NB7B
211C FXEP,FYEP,FZEP -> NB12,NB12A,NB12B X0,Y0,Z0 -> NB5,NB5A,NB5B
212C XMOM,YMOM,ZMOM -> NB8,NB8A,NB8B RX,RY,RZ -> NB9,NB9A,NB9B
213C RPX,RPY,RPZ -> NB10,NB10A,NB10B XMEP,YMEP,ZMEP -> NB11,NB11A,NB11B
214C RPX2,RPY2,RPZ2 -> NB13,NB13A,NB13B
215C RLOC -> NB14
216 CALL r4coork3(jft ,jlt ,x ,ixr ,pid ,
217 2 ngl ,al ,r11 ,r12 ,r13 ,
218 3 r21 ,r22 ,r23 ,r31 ,r32 ,
219 4 r33 ,gbuf%SKEW,gbuf%OFF,off )
220 CALL r13mat3(jft ,jlt ,geo ,kx ,ky ,
221 2 kz ,mx ,my ,mz ,pid ,
222 3 gbuf%LENGTH(ii(1)),gbuf%FOR(ii(1)),gbuf%FOR(ii(2)),gbuf%FOR(ii(3)),gbuf%TOTDEPL(ii(1)),
223 4 gbuf%TOTDEPL(ii(2)),gbuf%TOTDEPL(ii(3)),gbuf%MOM(ii(1)),gbuf%MOM(ii(2)),gbuf%MOM(ii(3)),
224 5 gbuf%TOTROT(ii(1)),gbuf%TOTROT(ii(2)),gbuf%TOTROT(ii(3)),tf ,npf ,
225 6 gbuf%POSX,gbuf%POSY,gbuf%POSZ,gbuf%POSXX,gbuf%POSYY,
226 7 gbuf%POSZZ,igeo )
227C-------------------------
228C MATRICE de RIGIDITE
229C-------------------------
230 CALL r13sumg3(jft ,jlt ,al ,kx ,ky ,
231 2 kz ,mx ,my ,mz ,r11 ,
232 3 r12 ,r13 ,r21 ,r22 ,r23 ,
233 4 r31 ,r32 ,r33 ,ke11 ,ke12 ,
234 5 ke22 )
235 IF (ikgeo==1) CALL pkeg3(jft ,jlt ,
236 1 gbuf%FOR(ii(1)) ,al ,ke11 ,ke12 ,ke22 )
237
238
239 IF (neig>0) CALL peoff(
240 1 nixr, jft, jlt, ixr, etag,
241 2 off )
242 CALL assem_p(nixr,6 ,
243 1 ixr ,nel ,iddl ,ndof ,k_diag ,
244 2 k_lt ,iadk ,jdik ,ke11 ,ke12 ,
245 3 ke22 ,off )
246
247C-----------------------------------------------
248 RETURN
249 END
assem_p
subroutine assem_p(nixpl, nd, ixp, nel, iddl, ndof, k_diag, k_lt, iadk, jdik, kc11, kc12, kc22, off)
Definition assem_p.F:40
my_real
#define my_real
Definition cppsort.cpp:32
peoff
subroutine peoff(nixpl, jft, jlt, ixp, etag, off)
Definition peoff.F:35
pkeg3
subroutine pkeg3(jft, jlt, fx, al, ke11, ke12, ke22)
Definition pke3.F:129
r13ke3
subroutine r13ke3(jft, jlt, nel, mtn, pm, geo, ixr, x, elbuf_tab, npf, tf, skew, offset, fr_wave, ikgeo, igeo, etag, iddl, ndof, k_diag, k_lt, iadk, jdik)
Definition r13ke3.F:43
r13mat3
subroutine r13mat3(jft, jlt, geo, kx, ky, kz, mx, my, mz, mgn, al0, fx, fy, fz, dx, dy, dz, xmom, ymom, zmom, rx, ry, rz, tf, npf, posx, posy, posz, posxx, posyy, poszz, igeo)
Definition r13mat3.F:38
r13sumg3
subroutine r13sumg3(jft, jlt, al, kx, ky, kz, mx, my, mz, r11, r12, r13, r21, r22, r23, r31, r32, r33, ke11, ke12, ke22)
Definition r13sumg3.F:33
r4coork3
subroutine r4coork3(jft, jlt, x, ncc, pid, ngl, al, r11, r12, r13, r21, r22, r23, r31, r32, r33, rloc, offg, off)
Definition r4coork3.F:32