OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
tforc3.F File Reference
#include "implicit_f.inc"
#include "mvsiz_p.inc"
#include "param_c.inc"
#include "parit_c.inc"
#include "com08_c.inc"

Go to the source code of this file.

Functions/Subroutines

subroutine tforc3 (elbuf_str, jft, jlt, pm, geo, nct, x, f, v, partsav, bufmat, dt2t, neltst, ityptst, stifn, fsky, iadt, offset, ipartt, tani, fx1, fx2, fy1, fy2, fz1, fz2, nel, gresav, grth, igrth, mstr, dmeltr, ioutprt, ipm, npf, tf, itask, h3d_data, nft, mtn, jsms, igre, preld1, stf_f, dt, sensors)

Function/Subroutine Documentation

◆ tforc3()

subroutine tforc3 ( type (elbuf_struct_), target elbuf_str,
integer jft,
integer jlt,
pm,
geo,
integer, dimension(*) nct,
x,
f,
v,
partsav,
bufmat,
dt2t,
integer neltst,
integer ityptst,
stifn,
fsky,
integer, dimension(2,*) iadt,
integer offset,
integer, dimension(*) ipartt,
tani,
fx1,
fx2,
fy1,
fy2,
fz1,
fz2,
integer nel,
gresav,
integer, dimension(*) grth,
integer, dimension(*) igrth,
mstr,
dmeltr,
integer ioutprt,
integer, dimension(npropmi,*) ipm,
integer, dimension(*) npf,
tf,
integer itask,
type(h3d_database) h3d_data,
integer, intent(in) nft,
integer, intent(in) mtn,
integer, intent(in) jsms,
integer, intent(in) igre,
intent(in) preld1,
intent(in) stf_f,
type (dt_), intent(in) dt,
type (sensors_), intent(inout) sensors )

Definition at line 46 of file tforc3.F.

59C-----------------------------------------------
60C M o d u l e s
61C-----------------------------------------------
62 USE elbufdef_mod
63 USE h3d_mod
64 USE preload_axial_mod
65 USE dt_mod
66 USE sensor_mod
67C-----------------------------------------------
68C I m p l i c i t T y p e s
69C-----------------------------------------------
70#include "implicit_f.inc"
71C-----------------------------------------------
72C G l o b a l P a r a m e t e r s
73C-----------------------------------------------
74#include "mvsiz_p.inc"
75C-----------------------------------------------
76C C o m m o n B l o c k s
77C-----------------------------------------------
78#include "param_c.inc"
79#include "parit_c.inc"
80#include "com08_c.inc"
81C-----------------------------------------------
82C D u m m y A r g u m e n t s
83C-----------------------------------------------
84 INTEGER, INTENT(IN) :: IGRE
85 INTEGER, INTENT(IN) :: NFT
86 INTEGER, INTENT(IN) :: MTN
87 INTEGER, INTENT(IN) :: JSMS
88 INTEGER NCT(*),IADT(2,*),IPARTT(*),OFFSET,
89 . JFT,JLT,NELTST,ITYPTST,NEL,GRTH(*),IPM(NPROPMI,*),
90 . IGRTH(*),IOUTPRT,NPF(*),ITASK
91 my_real dt2t,
92 . bufmat(*),pm(*),geo(npropg,*),x(*),f(*),v(*),partsav(*),tf(*),
93 . stifn(*),fsky(*),tani(15,*),fx1(mvsiz),fy1(mvsiz),fz1(mvsiz),
94 . fx2(mvsiz),fy2(mvsiz),fz2(mvsiz),gresav(*),mstr(*),dmeltr(*)
95 my_real, INTENT(IN) :: preld1,stf_f
96C
97 TYPE (ELBUF_STRUCT_), TARGET :: ELBUF_STR
98 TYPE(H3D_DATABASE) :: H3D_DATA
99 TYPE (DT_), INTENT(IN) :: DT
100 type (sensors_),INTENT(INOUT) :: SENSORS
101
102C-----------------------------------------------
103C L o c a l V a r i a b l e s
104C-----------------------------------------------
105 INTEGER LCO,NUVAR,NFUNC,IFUNC(100),I
106 INTEGER MAT(MVSIZ),PID(MVSIZ),NGL(MVSIZ),NC1(MVSIZ),NC2(MVSIZ)
107 my_real
108 . sti(mvsiz),off(mvsiz),
109 . eps(mvsiz),al(mvsiz),vx1(mvsiz),vx2(mvsiz),vy1(mvsiz),
110 . vy2(mvsiz),vz1(mvsiz),vz2(mvsiz),ex(mvsiz),ey(mvsiz),ez(mvsiz),
111 . x1(mvsiz),x2(mvsiz),y1(mvsiz),y2(mvsiz),z1(mvsiz),z2(mvsiz)
112C
113 my_real ,DIMENSION(:) ,POINTER :: uvar
114 TYPE(G_BUFEL_),POINTER :: GBUF
115C-----------------------------------------------
116C S o u r c e L i n e s
117C-----------------------------------------------
118 gbuf => elbuf_str%GBUF
119C
120 lco=1+5*nft
121 CALL tcoor3(
122 1 x, nct(lco),mat, pid,
123 2 ngl, nc1, nc2, x1,
124 3 x2, y1, y2, z1,
125 4 z2, nel)
126 CALL tdefo3(
127 1 v, gbuf%OFF,off, nc1,
128 2 nc2, eps, al, vx1,
129 3 vx2, vy1, vy2, vz1,
130 4 vz2, ex, ey, ez,
131 5 x1, x2, y1, y2,
132 6 z1, z2, nel)
133C------- activation only in compression in laws
134 DO i=1,nel
135 IF (geo(2,pid(i))>zero .AND. gbuf%OFF(i)>zero) off(i)=zero
136 ENDDO
137 CALL tdlen3(
138 1 jft, jlt, pm, gbuf%OFF,
139 2 off, dt2t, neltst, ityptst,
140 3 mstr, dmeltr, gbuf%DT, nel,
141 4 gbuf%G_DT,mat, ngl, al,
142 5 jsms)
143!
144 IF (mtn == 1) THEN
145 CALL m1lawt(
146 1 pm, geo, off, gbuf%FOR,
147 2 gbuf%EINT, gbuf%AREA, gbuf%LENGTH,sti,
148 3 mat, pid, eps, al,
149 4 nel)
150c
151 ELSEIF (mtn == 2) THEN
152 CALL m2lawt(
153 1 pm, geo, off, gbuf%FOR,
154 2 gbuf%EINT, gbuf%AREA, gbuf%LENGTH,gbuf%PLA,
155 3 sti, mat, pid, ngl,
156 4 eps, al, nel)
157c
158 ELSEIF (mtn == 34) THEN
159 nuvar = gbuf%G_NUVAR
160 uvar => gbuf%VAR
161 CALL sigeps34t(nel ,ngl ,mat ,pid ,bufmat ,
162 . ipm ,geo ,off ,gbuf%FOR ,sti ,
163 . gbuf%EINT,gbuf%AREA,gbuf%LENGTH,al ,eps ,
164 . nuvar ,uvar )
165c
166c
167 ELSEIF (mtn == 44) THEN
168 nuvar = gbuf%G_NUVAR
169 uvar => gbuf%VAR
170 nfunc = ipm(10,mat(1))
171 DO i=1,nfunc
172 ifunc(i) = ipm(10+i,mat(1))
173 ENDDO
174 CALL sigeps44t(nel ,ngl ,mat ,pid ,bufmat ,
175 . ipm ,geo ,off ,gbuf%FOR ,sti ,
176 . gbuf%PLA ,gbuf%EINT,gbuf%AREA,gbuf%LENGTH,al ,
177 . eps ,nuvar ,uvar ,npf ,tf ,
178 . nfunc ,ifunc )
179c
180 ENDIF ! IF (MTN)
181C------- activation for GBUF%OFF + save the strain in buffer
182 DO i=1,nel
183 IF (geo(2,pid(i))>zero) THEN
184 IF (gbuf%OFF(i)==zero.AND.off(i)==one) gbuf%OFF(i)=one
185 END IF
186 gbuf%STRA(i) = gbuf%STRA(i) + eps(i)*dt1
187 ENDDO
188C--------------------------
189C BILANS PAR MATERIAU
190C--------------------------
191 IF (ioutprt>0) THEN
192 CALL tbilan(
193 1 pm, v, gbuf%EINT,gbuf%AREA,
194 2 partsav, ipartt, tani, gbuf%FOR,
195 3 gresav, grth, igrth, mat,
196 4 al, vx1, vx2, vy1,
197 5 vy2, vz1, vz2, x1,
198 6 x2, y1, y2, z1,
199 7 z2, itask, h3d_data, nel,
200 8 igre, sensors)
201 END IF !(IOUTPRT>0) THEN
202C-------------------------
203C /PRELOAD/AXIAL
204C-------------------------
205 IF (preld1>zero) THEN
206 CALL preload_axial(nel,preld1,gbuf%BPRELD,eps,stf_f,gbuf%FOR)
207 END IF
208C-------------------------
209C ASSEMBLE
210C-------------------------
211 IF (iparit == 0) THEN
212 CALL tfcum3(
213 1 f, gbuf%FOR,stifn, sti,
214 2 fx1, fx2, fy1, fy2,
215 3 fz1, fz2, gbuf%OFF,off,
216 4 nc1, nc2, ex, ey,
217 5 ez, nel)
218 ELSE
219 CALL tfcum3p(
220 1 gbuf%FOR,sti, fsky, fsky,
221 2 iadt, fx1, fx2, fy1,
222 3 fy2, fz1, fz2, gbuf%OFF,
223 4 off, nc1, nc2, ex,
224 5 ey, ez, nel, nft)
225 ENDIF
226C-----------------------------------------------
227 RETURN
my_real
#define my_real
Definition cppsort.cpp:32
m1lawt
subroutine m1lawt(pm, geo, off, for, eint, area, al0, sti, mat, mgm, eps, al, nel)
Definition m1lawt.F:33
m2lawt
subroutine m2lawt(pm, geo, off, for, eint, area, al0, pla, sti, mat, mgm, ngl, eps, al, nel)
Definition m2lawt.F:33
dt_mod
Definition dt_mod.F:79
h3d_mod
Definition h3d_mod.F:308
sigeps34t
subroutine sigeps34t(nel, ngl, mat, pid, uparam, ipm, geo, off, for, sti, eint, area, al0, al, epsp, nuvar, uvar)
Definition sigeps34t.F:32
sigeps44t
subroutine sigeps44t(nel, ngl, mat, pid, uparam, ipm, geo, off, for, sti, pla, eint, area, al0, al, epsp, nuvar, uvar, npf, tf, nfunc, ifunc)
Definition sigeps44t.F:35
tbilan
subroutine tbilan(pm, v, eint, area, partsav, ipartt, tani, for, gresav, grth, igrth, mat, al, vx1, vx2, vy1, vy2, vz1, vz2, x1, x2, y1, y2, z1, z2, itask, h3d_data, nel, igre, sensors)
Definition tbilan.F:43
tcoor3
subroutine tcoor3(x, ncc, mat, mgm, ngl, nc1, nc2, x1, x2, y1, y2, z1, z2, nel)
Definition tcoor3.F:33
tdefo3
subroutine tdefo3(v, offg, off, nc1, nc2, eps, al, vx1, vx2, vy1, vy2, vz1, vz2, ex, ey, ez, x1, x2, y1, y2, z1, z2, nel)
Definition tdefo3.F:35
tdlen3
subroutine tdlen3(jft, jlt, pm, offg, off, dt2t, neltst, ityptst, mstr, dmeltr, dtel, nel, g_dt, mat, ngl, al, jsms)
Definition tdlen3.F:34
tfcum3
subroutine tfcum3(f, for, stifn, sti, fx1, fx2, fy1, fy2, fz1, fz2, offg, off, nc1, nc2, ex, ey, ez, nel)
Definition tfcum3.F:34
tfcum3p
subroutine tfcum3p(for, sti, fsky, fskyv, iadt, fx1, fx2, fy1, fy2, fz1, fz2, offg, off, nc1, nc2, ex, ey, ez, nel, nft)
Definition tfcum3p.F:34