OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
pke3.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!|| pke3 ../engine/source/elements/beam/pke3.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!|| pcoork3 ../engine/source/elements/beam/pcoork3.F
31!|| peoff ../engine/source/elements/beam/peoff.F
32!|| pkeg3 ../engine/source/elements/beam/pke3.F
33!|| pmat3 ../engine/source/elements/beam/pmat3.F
34!|| psumg3 ../engine/source/elements/beam/psumg3.F
35!||--- uses -----------------------------------------------------
36!|| elbufdef_mod ../common_source/modules/mat_elem/elbufdef_mod.F90
37!|| element_mod ../common_source/modules/elements/element_mod.F90
38!||====================================================================
39 SUBROUTINE pke3(JFT ,JLT ,NEL ,MTN ,ISMSTR,
40 1 PM ,NCC ,X ,ELBUF_TAB,GEO ,
41 2 OFFSET, IKGEO,
42 3 ETAG , IDDL ,NDOF ,K_DIAG ,K_LT ,
43 4 IADK ,JDIK )
44C-----------------------------------------------
45C M o d u l e s
46C-----------------------------------------------
47 USE elbufdef_mod
48 use element_mod , only : nixp
49C-----------------------------------------------
50C I m p l i c i t T y p e s
51C-----------------------------------------------
52#include "implicit_f.inc"
53C-----------------------------------------------
54C G l o b a l P a r a m e t e r s
55C-----------------------------------------------
56#include "mvsiz_p.inc"
57#include "com04_c.inc"
58C-----------------------------------------------
59C D u m m y A r g u m e n t s
60C-----------------------------------------------
61 INTEGER NCC(NIXP,*),JFT,JLT,OFFSET,NEL ,MTN,ISMSTR,IKGEO,
62 . ETAG(*),IDDL(*) ,NDOF(*) ,IADK(*) ,JDIK(*)
63 my_real
64 . PM(*), X(*),GEO(*),
65 . ke11(36,mvsiz),ke12(36,mvsiz),ke22(36,mvsiz),off(mvsiz),
66 . k_diag(*) ,k_lt(*)
67C
68 TYPE (ELBUF_STRUCT_), TARGET:: ELBUF_TAB
69C-----------------------------------------------
70C L o c a l V a r i a b l e s
71C-----------------------------------------------
72 INTEGER
73 . I
74 INTEGER MAT(MVSIZ),PID(MVSIZ),NGL(MVSIZ)
75 my_real
76 . ym(mvsiz),g(mvsiz),al(mvsiz),ali(mvsiz),
77 . sh1(mvsiz),sh2(mvsiz),
78 . area(mvsiz),b1(mvsiz),b2(mvsiz),b3(mvsiz),
79 . r11(mvsiz),r12(mvsiz),r13(mvsiz),
80 . r21(mvsiz),r22(mvsiz),r23(mvsiz),
81 . r31(mvsiz),r32(mvsiz),r33(mvsiz),fx(mvsiz)
82 TYPE(g_bufel_) ,POINTER :: GBUF
83C-----------------------------------------------
84C S o u r c e L i n e s
85C-----------------------------------------------
86 gbuf => elbuf_tab%GBUF
87
88 CALL pcoork3(jft ,jlt ,x ,ncc ,gbuf%OFF,
89 2 pid ,mat ,ngl ,al ,ali ,
90 3 gbuf%LENGTH,ismstr ,r11 ,r12 ,r13 ,
91 4 r21 ,r22 ,r23 ,r31 ,r32 ,
92 5 r33 ,gbuf%SKEW ,off )
93 CALL pmat3(jft ,jlt ,pm ,geo ,ym ,
94 2 g ,area ,b1 ,b2 ,b3 ,
95 3 al ,ali ,sh1 ,sh2 ,mtn ,
96 4 mat ,pid )
97C-------------------------
98C STIFFNESS MATRIX
99C-------------------------
100 CALL psumg3(jft ,jlt ,ym ,g ,area ,
101 2 b1 ,b2 ,b3 ,al ,ali ,
102 3 sh1 ,sh2 ,r11 ,r12 ,r13 ,
103 4 r21 ,r22 ,r23 ,r31 ,r32 ,
104 5 r33 ,ke11 ,ke12 ,ke22 )
105 IF (ikgeo>0) THEN
106 DO i =jft ,jlt
107 fx(i)=gbuf%FOR(2*nel+i)
108 ENDDO
109 CALL pkeg3(jft ,jlt ,
110 1 fx ,al ,ke11 ,ke12 ,ke22 )
111 ENDIF
112
113 IF (neig>0) CALL peoff(
114 1 nixp, jft, jlt, ncc, etag,
115 2 off )
116 CALL assem_p(nixp,6 ,
117 1 ncc ,nel ,iddl ,ndof ,k_diag ,
118 2 k_lt ,iadk ,jdik ,ke11 ,ke12 ,
119 3 ke22 ,off )
120C
121 RETURN
122 END
123!||====================================================================
124!|| pkeg3 ../engine/source/elements/beam/pke3.F
125!||--- called by ------------------------------------------------------
126!|| pke3 ../engine/source/elements/beam/pke3.F
127!|| r13ke3 ../engine/source/elements/spring/r13ke3.F
128!||====================================================================
129 SUBROUTINE pkeg3(JFT ,JLT ,
130 1 FX ,AL ,KE11 ,KE12 ,KE22 )
131C-----------------------------------------------
132C I M P L I C I T T Y P E S
133C-----------------------------------------------
134#include "implicit_f.inc"
135#include "mvsiz_p.inc"
136C-----------------------------------------------
137C D U M M Y A R G U M E N T S
138C-----------------------------------------------
139 INTEGER JFT,JLT
140 MY_REAL
141 . FX(*),AL(*),KE11(6,6,*),KE12(6,6,*),KE22(6,6,*)
142C-----------------------------------------------
143C L O C A L V A R I A B L E S
144C-----------------------------------------------
145 INTEGER EP,I
146 MY_REAL
147 . H(MVSIZ)
148C-----------
149 DO EP=jft,jlt
150 h(ep)=fx(ep)/al(ep)
151 ENDDO
152 DO i=1,3
153 DO ep=jft,jlt
154 ke11(i,i,ep) = ke11(i,i,ep)+h(ep)
155 ke12(i,i,ep) = ke12(i,i,ep)-h(ep)
156 ke22(i,i,ep) = ke22(i,i,ep)+h(ep)
157 ENDDO
158 ENDDO
159C
160 RETURN
161 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
area
subroutine area(d1, x, x2, y, y2, eint, stif0)
Definition hm_read_prop32.F:567
pcoork3
subroutine pcoork3(jft, jlt, x, ncc, offg, pid, mat, ngl, al, ali, als, ismstr, r11, r12, r13, r21, r22, r23, r31, r32, r33, rloc, off)
Definition pcoork3.F:33
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:131
pke3
subroutine pke3(jft, jlt, nel, mtn, ismstr, pm, ncc, x, elbuf_tab, geo, offset, ikgeo, etag, iddl, ndof, k_diag, k_lt, iadk, jdik)
Definition pke3.F:44
pmat3
subroutine pmat3(jft, jlt, pm, geo, ym, g, a1, b1, b2, b3, al, ali, sh1, sh2, mtn, mat, pid)
Definition pmat3.F:32
psumg3
subroutine psumg3(jft, jlt, ym, g, area, b1, b2, b3, al, ali, sh1, sh2, r11, r12, r13, r21, r22, r23, r31, r32, r33, ke11, ke12, ke22)
Definition psumg3.F:33