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