OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
cepsini.F File Reference
#include "implicit_f.inc"
#include "mvsiz_p.inc"
#include "param_c.inc"
#include "scr03_c.inc"

Go to the source code of this file.

Functions/Subroutines

subroutine cepsini (elbuf_str, mat_param, jft, jlt, ismstr, ilaw, ithk, pm, geo, ixc, x, xrefc, for, thk, eint, gstr, nlay, px1g, px2g, py1g, py2g, x2s, y2s, x3s, y3s, x4s, y4s, uvar, ipm, igeo, imat, skew, nel, dir_a, dir_b, sigi, npf, tf, irep)
subroutine cmlawi (mat_param, jft, jlt, ilaw, ithk, pm, for, thk, eint, gstr, dir1, dir2, uvar, ipm, nel, sigi, npf, tf, area, exx, eyy, exy, imat)
subroutine cpxpyi (jft, jlt, ismstr, area, px1g, px2g, py1g, py2g, px1, px2, py1, py2, x1g, x2g, x3g, x4g, y1g, y2g, y3g, y4g, z1g, z2g, z3g, z4g, e1x, e2x, e3x, e1y, e2y, e3y, e1z, e2z, e3z, x2l, y2l, x3l, y3l, x4l, y4l)
subroutine cdefoi (jft, jlt, vl, gstr, px1, px2, py1, py2, nel, area, exx, eyy, exy, eyz, ezx, e1x, e2x, e3x, e1y, e2y, e3y, e1z, e2z, e3z)
subroutine cm1ini (jft, jlt, ithk, pm, for, thk, eint, nel, imat, area, exx, eyy, exy)
subroutine cm19ini (jft, jlt, ithk, pm, for, thk, eint, gstr, dir, sigi, nel, imat, area, exx, eyy, exy)
subroutine roto (jft, jlt, tab, ltab, dir, nel)
subroutine uroto (jft, jlt, tab, ltab, dir, nel)
subroutine cepschk (jft, jlt, nft, pm, geo, ixc, gstr, thk, nel, cpt_eltens)

Function/Subroutine Documentation

◆ cdefoi()

subroutine cdefoi ( integer jft,
integer jlt,
vl,
gstr,
px1,
px2,
py1,
py2,
integer nel,
area,
exx,
eyy,
exy,
eyz,
ezx,
e1x,
e2x,
e3x,
e1y,
e2y,
e3y,
e1z,
e2z,
e3z )

Definition at line 429 of file cepsini.F.

433C-----------------------------------------------
434C I m p l i c i t T y p e s
435C-----------------------------------------------
436#include "implicit_f.inc"
437C-----------------------------------------------
438C G l o b a l P a r a m e t e r s
439C-----------------------------------------------
440#include "mvsiz_p.inc"
441C-----------------------------------------------
442C D u m m y A r g u m e n t s
443C-----------------------------------------------
444 INTEGER JFT, JLT,NEL
445 my_real
446 . vl(3,4,*) ,gstr(nel,8),
447 . px1(nel) ,px2(nel) ,py1(nel) ,py2(nel),
448 . e1x(mvsiz) , e1y(mvsiz) , e1z(mvsiz) ,
449 . e2x(mvsiz) , e2y(mvsiz) , e2z(mvsiz) ,
450 . e3x(mvsiz) , e3y(mvsiz) , e3z(mvsiz) ,
451 . area(mvsiz),exx(mvsiz),eyy(mvsiz),exy(mvsiz),eyz(mvsiz),ezx(mvsiz)
452C-----------------------------------------------
453C L o c a l V a r i a b l e s
454C-----------------------------------------------
455 INTEGER I
456C REAL
457 my_real
458 . vx1(mvsiz) , vx2(mvsiz) , vx3(mvsiz) , vx4(mvsiz) ,
459 . vy1(mvsiz) , vy2(mvsiz) , vy3(mvsiz) , vy4(mvsiz) ,
460 . vz1(mvsiz) , vz2(mvsiz) , vz3(mvsiz) , vz4(mvsiz) ,
461 . vx13(mvsiz), vx24(mvsiz), vy13(mvsiz), vy24(mvsiz),
462 . vz13(mvsiz), vz24(mvsiz)
463 my_real tmp1, tmp2, tmp3, fac1
464C=======================================================================
465 DO i=jft,jlt
466 vx1(i) = e1x(i)*vl(1,1,i) + e1y(i)*vl(2,1,i) + e1z(i)*vl(3,1,i)
467 vx2(i) = e1x(i)*vl(1,2,i) + e1y(i)*vl(2,2,i) + e1z(i)*vl(3,2,i)
468 vx3(i) = e1x(i)*vl(1,3,i) + e1y(i)*vl(2,3,i) + e1z(i)*vl(3,3,i)
469 vx4(i) = e1x(i)*vl(1,4,i) + e1y(i)*vl(2,4,i) + e1z(i)*vl(3,4,i)
470C
471 vy4(i) = e2x(i)*vl(1,4,i) + e2y(i)*vl(2,4,i) + e2z(i)*vl(3,4,i)
472 vy3(i) = e2x(i)*vl(1,3,i) + e2y(i)*vl(2,3,i) + e2z(i)*vl(3,3,i)
473 vy2(i) = e2x(i)*vl(1,2,i) + e2y(i)*vl(2,2,i) + e2z(i)*vl(3,2,i)
474 vy1(i) = e2x(i)*vl(1,1,i) + e2y(i)*vl(2,1,i) + e2z(i)*vl(3,1,i)
475C
476 vz1(i) = e3x(i)*vl(1,1,i) + e3y(i)*vl(2,1,i) + e3z(i)*vl(3,1,i)
477 vz2(i) = e3x(i)*vl(1,2,i) + e3y(i)*vl(2,2,i) + e3z(i)*vl(3,2,i)
478 vz3(i) = e3x(i)*vl(1,3,i) + e3y(i)*vl(2,3,i) + e3z(i)*vl(3,3,i)
479 vz4(i) = e3x(i)*vl(1,4,i) + e3y(i)*vl(2,4,i) + e3z(i)*vl(3,4,i)
480 ENDDO
481C
482 DO i=jft,jlt
483 vz13(i) = vz1(i)-vz3(i)
484 vz24(i) = vz2(i)-vz4(i)
485 tmp2 = py2(i)+py1(i)
486 tmp3 = sign(max(abs(tmp2),em20),tmp2)
487 tmp1 = fourth*(vz13(i)-vz24(i))**2/tmp3
488 vx13(i) = vx1(i)-vx3(i)
489 vx24(i) = vx2(i)-vx4(i)
490 vx13(i) = vx13(i)-tmp1
491 vx24(i) = vx24(i)+tmp1
492C
493 exx(i) = px1(i)*vx13(i)+px2(i)*vx24(i)
494 exy(i) = py1(i)*vx13(i)+py2(i)*vx24(i)
495C
496 tmp1=px2(i)-px1(i)
497 tmp3=sign(max(abs(tmp1),em20),tmp1)
498 tmp2=fourth*(vz13(i)+vz24(i))**2/tmp3
499 vy13(i)=vy1(i)-vy3(i)
500 vy24(i)=vy2(i)-vy4(i)
501 vy13(i)=vy13(i)+tmp2
502 vy24(i)=vy24(i)+tmp2
503C
504 exy(i)=exy(i)+px1(i)*vy13(i)+px2(i)*vy24(i)
505 eyy(i)=py1(i)*vy13(i)+py2(i)*vy24(i)
506 eyz(i)=py1(i)*vz13(i)+py2(i)*vz24(i)
507 ezx(i)=px1(i)*vz13(i)+px2(i)*vz24(i)
508 ENDDO
509C
510 DO i=jft,jlt
511 fac1 = one/area(i)
512 exx(i)= exx(i)*fac1
513 eyy(i)= eyy(i)*fac1
514 exy(i)= exy(i)*fac1
515 ENDDO
516C
517 DO i=jft,jlt
518 gstr(i,1)=gstr(i,1)+exx(i)
519 gstr(i,2)=gstr(i,2)+eyy(i)
520 gstr(i,3)=gstr(i,3)+exy(i)
521 ENDDO
522 RETURN
my_real
#define my_real
Definition cppsort.cpp:32
area
subroutine area(d1, x, x2, y, y2, eint, stif0)
Definition hm_read_prop32.F:567
max
#define max(a, b)
Definition macros.h:21

◆ cepschk()

subroutine cepschk ( integer jft,
integer jlt,
integer nft,
pm,
geo,
integer, dimension(nixc,*) ixc,
gstr,
thk,
integer nel,
integer cpt_eltens )

Definition at line 827 of file cepsini.F.

828 USE message_mod
829C-----------------------------------------------
830C I m p l i c i t T y p e s
831C-----------------------------------------------
832#include "implicit_f.inc"
833C-----------------------------------------------
834C G l o b a l P a r a m e t e r s
835C-----------------------------------------------
836#include "mvsiz_p.inc"
837C-----------------------------------------------
838C C o m m o n B l o c k s
839C-----------------------------------------------
840#include "param_c.inc"
841#include "scr03_c.inc"
842C-----------------------------------------------
843C D u m m y A r g u m e n t s
844C-----------------------------------------------
845 INTEGER JFT, JLT,NFT ,NEL ,IXC(NIXC,*),CPT_ELTENS
846 my_real
847 . pm(npropm,*), gstr(nel,8),thk(*), geo(npropg,*)
848C------------------------------------------------------
849C L o c a l V a r i a b l e s
850C------------------------------------------------------
851 INTEGER I
852 my_real
853 . delt(mvsiz),x1(mvsiz),x2(mvsiz)
854C-------------------------------------------------------------
855 DO i=jft,jlt
856 delt(i)=((gstr(i,1)+gstr(i,2))*(gstr(i,1)+gstr(i,2)))-four*
857 . (gstr(i,1)*gstr(i,2)-fourth*gstr(i,3)*gstr(i,3))
858 IF(delt(i)>=zero)THEN
859 x1(i) = (gstr(i,1)+gstr(i,2)-sqrt(delt(i)))/two
860 x2(i) = (gstr(i,1)+gstr(i,2)+sqrt(delt(i)))/two
861 IF((x1(i)>zero.AND.(x1(i)>em10)).OR.(x2(i)>zero
862 . .AND.(x2(i)>em10)))THEN
863 IF(ipri == 5) THEN
864 CALL ancmsg(msgid=607,
865 . msgtype=msgwarning,
866 . anmode=aninfo_blind_1,
867 . r2=min(x1(i),x2(i)),
868 . r1=max(x1(i),x2(i)),
869 . i1=ixc(7,i))
870 ELSE
871 cpt_eltens = cpt_eltens + 1
872 ENDIF
873 ENDIF
874 ENDIF
875 ENDDO
min
#define min(a, b)
Definition macros.h:20
message_mod
Definition message_mod.F:1249
ancmsg
subroutine ancmsg(msgid, msgtype, anmode, i1, i2, i3, i4, i5, i6, i7, i8, i9, i10, i11, i12, i13, i14, i15, i16, i17, i18, i19, i20, r1, r2, r3, r4, r5, r6, r7, r8, r9, c1, c2, c3, c4, c5, c6, c7, c8, c9, prmode)
Definition message.F:889

◆ cepsini()

subroutine cepsini ( type (elbuf_struct_), target elbuf_str,
type (matparam_struct_), intent(in) mat_param,
integer jft,
integer jlt,
integer ismstr,
integer ilaw,
integer ithk,
pm,
geo,
integer, dimension(nixc,*) ixc,
x,
xrefc,
for,
thk,
eint,
gstr,
integer nlay,
px1g,
px2g,
py1g,
py2g,
x2s,
y2s,
x3s,
y3s,
x4s,
y4s,
uvar,
integer, dimension(npropmi,*) ipm,
integer, dimension(npropgi,*) igeo,
integer imat,
skew,
integer nel,
dir_a,
dir_b,
sigi,
integer, dimension(*) npf,
tf,
integer irep )

Definition at line 35 of file cepsini.F.

44C-----------------------------------------------
45C M o d u l e s
46C-----------------------------------------------
47 USE mat_elem_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 "param_c.inc"
60#include "scr03_c.inc"
61C-----------------------------------------------
62C D u m m y A r g u m e n t s
63C-----------------------------------------------
64 INTEGER JFT, JLT, ISMSTR, ILAW, ITHK, NEL, NLAY,IREP,IMAT
65 INTEGER IXC(NIXC,*),IGEO(NPROPGI,*),IPM(NPROPMI,*),NPF(*)
66 my_real
67 . pm(npropm,*),geo(npropg,*),x(3,*),xrefc(4,3,*),skew(lskew,*),
68 . for(nel,5) ,thk(*) ,eint(jlt,2),gstr(nel,8),
69 . px1g(*) ,px2g(*) ,py1g(*) ,py2g(*), uvar(nel,*),
70 . x2s(*), y2s(*), x3s(*), y3s(*), x4s(*), y4s(*),
71 . dir_a(nel,*),dir_b(nel,*),sigi(nel,3),tf(*)
72 TYPE (ELBUF_STRUCT_), TARGET :: ELBUF_STR
73 TYPE (MATPARAM_STRUCT_) ,INTENT(IN) :: MAT_PARAM
74C------------------------------------------------------
75C L o c a l V a r i a b l e s
76C-----------------------------------------------
77 INTEGER I, NT, I1, I2, I3, I4,IDRAPE,IGTYP
78 my_real
79 . fac
80 my_real
81 . x1(mvsiz) , x2(mvsiz) , x3(mvsiz) , x4(mvsiz) ,
82 . y1(mvsiz) , y2(mvsiz) , y3(mvsiz) , y4(mvsiz) ,
83 . z1(mvsiz) , z2(mvsiz) , z3(mvsiz) , z4(mvsiz),
84 . px1(mvsiz), px2(mvsiz), py1(mvsiz), py2(mvsiz),
85 . vl(3,4,mvsiz), xl(3,4,mvsiz), area(mvsiz),
86 . e1x(mvsiz), e1y(mvsiz), e1z(mvsiz),
87 . e2x(mvsiz), e2y(mvsiz), e2z(mvsiz),
88 . e3x(mvsiz), e3y(mvsiz), e3z(mvsiz),
89 . x4i(mvsiz), y4i(mvsiz), z4i(mvsiz),
90 . exx(mvsiz),eyy(mvsiz),exy(mvsiz),eyz(mvsiz),ezx(mvsiz)
91 my_real,
92 . DIMENSION(:) , POINTER :: dir1, dir2
93C=======================================================================
94c initial state - rep global
95 idrape = elbuf_str%IDRAPE
96 igtyp = elbuf_str%IGTYP
97 DO i=jft,jlt
98 i1 = ixc(2,i)
99 i2 = ixc(3,i)
100 i3 = ixc(4,i)
101 i4 = ixc(5,i)
102 x1(i) = zero
103 y1(i) = zero
104 z1(i) = zero
105 x2(i) = x(1,i2) - x(1,i1)
106 y2(i) = x(2,i2) - x(2,i1)
107 z2(i) = x(3,i2) - x(3,i1)
108 x3(i) = x(1,i3) - x(1,i1)
109 y3(i) = x(2,i3) - x(2,i1)
110 z3(i) = x(3,i3) - x(3,i1)
111 x4(i) = x(1,i4) - x(1,i1)
112 y4(i) = x(2,i4) - x(2,i1)
113 z4(i) = x(3,i4) - x(3,i1)
114 ENDDO
115c
116 CALL ceveci(jft ,jlt ,area,
117 . x1 ,x2 ,x3 ,x4 ,y1 ,y2 ,
118 . y3 ,y4 ,z1 ,z2 ,z3 ,z4 ,
119 . e1x, e2x, e3x, e1y, e2y, e3y ,e1z, e2z, e3z )
120c
121 DO i=jft,jlt
122 xl(1,2,i) = e1x(i)*x2(i) + e1y(i)*y2(i) + e1z(i)*z2(i)
123 xl(2,2,i) = e2x(i)*x2(i) + e2y(i)*y2(i) + e2z(i)*z2(i)
124 xl(1,3,i) = e1x(i)*x3(i) + e1y(i)*y3(i) + e1z(i)*z3(i)
125 xl(2,3,i) = e2x(i)*x3(i) + e2y(i)*y3(i) + e2z(i)*z3(i)
126 xl(1,4,i) = e1x(i)*x4(i) + e1y(i)*y4(i) + e1z(i)*z4(i)
127 xl(2,4,i) = e2x(i)*x4(i) + e2y(i)*y4(i) + e2z(i)*z4(i)
128 ENDDO
129c
130c ref state - rep global
131 DO i=jft,jlt
132 x2(i) = xrefc(2,1,i) - xrefc(1,1,i)
133 y2(i) = xrefc(2,2,i) - xrefc(1,2,i)
134 z2(i) = xrefc(2,3,i) - xrefc(1,3,i)
135 x3(i) = xrefc(3,1,i) - xrefc(1,1,i)
136 y3(i) = xrefc(3,2,i) - xrefc(1,2,i)
137 z3(i) = xrefc(3,3,i) - xrefc(1,3,i)
138 x4(i) = xrefc(4,1,i) - xrefc(1,1,i)
139 y4(i) = xrefc(4,2,i) - xrefc(1,2,i)
140 z4(i) = xrefc(4,3,i) - xrefc(1,3,i)
141 ENDDO
142c
143 CALL ceveci(jft ,jlt ,area,
144 . x1 ,x2 ,x3 ,x4 ,y1 ,y2 ,
145 . y3 ,y4 ,z1 ,z2 ,z3 ,z4 ,
146 . e1x, e2x, e3x, e1y, e2y, e3y ,e1z, e2z, e3z )
147c
148 fac = one / float(nitrs)
149 DO i= jft,jlt
150 vl(1,1,i) = zero
151 vl(2,1,i) = zero
152 vl(3,1,i) = zero
153 vl(1,2,i) = e1x(i)*xl(1,2,i) + e2x(i)*xl(2,2,i)
154 vl(2,2,i) = e1y(i)*xl(1,2,i) + e2y(i)*xl(2,2,i)
155 vl(3,2,i) = e1z(i)*xl(1,2,i) + e2z(i)*xl(2,2,i)
156 vl(1,3,i) = e1x(i)*xl(1,3,i) + e2x(i)*xl(2,3,i)
157 vl(2,3,i) = e1y(i)*xl(1,3,i) + e2y(i)*xl(2,3,i)
158 vl(3,3,i) = e1z(i)*xl(1,3,i) + e2z(i)*xl(2,3,i)
159 vl(1,4,i) = e1x(i)*xl(1,4,i) + e2x(i)*xl(2,4,i)
160 vl(2,4,i) = e1y(i)*xl(1,4,i) + e2y(i)*xl(2,4,i)
161 vl(3,4,i) = e1z(i)*xl(1,4,i) + e2z(i)*xl(2,4,i)
162C
163 vl(1,2,i) = (vl(1,2,i) - x2(i)) * fac
164 vl(2,2,i) = (vl(2,2,i) - y2(i)) * fac
165 vl(3,2,i) = (vl(3,2,i) - z2(i)) * fac
166 vl(1,3,i) = (vl(1,3,i) - x3(i)) * fac
167 vl(2,3,i) = (vl(2,3,i) - y3(i)) * fac
168 vl(3,3,i) = (vl(3,3,i) - z3(i)) * fac
169 vl(1,4,i) = (vl(1,4,i) - x4(i)) * fac
170 vl(2,4,i) = (vl(2,4,i) - y4(i)) * fac
171 vl(3,4,i) = (vl(3,4,i) - z4(i)) * fac
172 ENDDO
173c--------------------------------------------------------
174 IF (ismstr /= 1 .AND. ismstr /= 11)THEN
175c-----------------
176 DO nt = 1,nitrs
177 fac = float(nt) - one
178 DO i=jft,jlt
179 x1(i) = xrefc(1,1,i) + fac*vl(1,1,i)
180 y1(i) = xrefc(1,2,i) + fac*vl(2,1,i)
181 z1(i) = xrefc(1,3,i) + fac*vl(3,1,i)
182 x2(i) = xrefc(2,1,i) + fac*vl(1,2,i)
183 y2(i) = xrefc(2,2,i) + fac*vl(2,2,i)
184 z2(i) = xrefc(2,3,i) + fac*vl(3,2,i)
185 x3(i) = xrefc(3,1,i) + fac*vl(1,3,i)
186 y3(i) = xrefc(3,2,i) + fac*vl(2,3,i)
187 z3(i) = xrefc(3,3,i) + fac*vl(3,3,i)
188 x4(i) = xrefc(4,1,i) + fac*vl(1,4,i)
189 y4(i) = xrefc(4,2,i) + fac*vl(2,4,i)
190 z4(i) = xrefc(4,3,i) + fac*vl(3,4,i)
191 ENDDO
192c rep intermediaire
193 CALL ceveci(jft ,jlt ,area,
194 . x1 ,x2 ,x3 ,x4 ,y1 ,y2 ,
195 . y3 ,y4 ,z1 ,z2 ,z3 ,z4 ,
196 . e1x, e2x, e3x, e1y, e2y, e3y ,e1z, e2z, e3z )
197c DIRA, DIRB : rep intermediaire
198 CALL corth3(elbuf_str,dir_a,dir_b,jft,jlt,
199 . nlay ,irep ,nel ,
200 . x1 ,x2 ,x3 ,x4 ,y1 ,y2 ,
201 . y3 ,y4 ,z1 ,z2 ,z3 ,z4 ,
202 . e1x, e2x, e3x, e1y, e2y, e3y ,e1z, e2z, e3z ,
203 . idrape, igtyp)
204c
205c PX1..PX3, X2S..X3S : rep intermediaire
206 CALL cpxpyi(jft ,jlt ,ismstr ,area ,
207 . px1g ,px2g ,py1g ,py2g ,
208 . px1 ,px2 ,py1 ,py2 ,
209 . x1 ,x2 ,x3 ,x4 ,y1 ,y2 ,
210 . y3 ,y4 ,z1 ,z2 ,z3 ,z4 ,
211 . e1x ,e2x ,e3x ,e1y ,e2y ,e3y ,e1z ,e2z ,e3z ,
212 . x2s ,y2s ,x3s ,y3s ,x4s ,y4s )
213c
214 CALL cdefoi(jft ,jlt ,vl ,gstr ,
215 . px1 ,px2 ,py1 ,py2 ,nel,
216 . area ,exx ,eyy ,exy ,eyz ,ezx ,
217 . e1x ,e2x ,e3x ,e1y ,e2y ,e3y ,e1z ,e2z ,e3z )
218c
219 CALL cmlawi(mat_param,
220 . jft ,jlt ,ilaw ,ithk ,pm ,
221 . for ,thk ,eint ,gstr ,dir_a ,
222 . dir_b ,uvar ,ipm ,
223 . nel ,sigi ,npf ,tf ,
224 . area ,exx ,eyy ,exy ,imat )
225C
226 ENDDO
227C------
228 ELSE ! SMALL STRAIN FORMULATION (ISMSTR = 1,11)
229C------
230c
231 CALL cpxpyi(jft ,jlt ,ismstr ,area ,
232 . px1g ,px2g ,py1g ,py2g ,
233 . px1 ,px2 ,py1 ,py2 ,
234 . x1 ,x2 ,x3 ,x4 ,y1 ,y2 ,
235 . y3 ,y4 ,z1 ,z2 ,z3 ,z4 ,
236 . e1x, e2x, e3x, e1y, e2y, e3y ,e1z, e2z, e3z ,
237 . x2s ,y2s ,x3s ,y3s ,x4s ,y4s )
238C
239 dir1 => elbuf_str%BUFLY(1)%DIRA
240 dir2 => elbuf_str%BUFLY(1)%DIRB
241c
242 DO nt=1,nitrs
243 CALL cdefoi(jft ,jlt ,vl ,gstr ,
244 . px1 ,px2 ,py1 ,py2 ,nel,
245 . area ,exx ,eyy ,exy ,eyz ,ezx ,
246 . e1x ,e2x ,e3x ,e1y ,e2y ,e3y ,e1z ,e2z ,e3z )
247c
248 CALL cmlawi(mat_param,
249 . jft ,jlt ,ilaw ,ithk ,pm ,
250 . for ,thk ,eint ,gstr ,dir1 ,
251 . dir2 ,uvar ,ipm ,
252 . nel ,sigi ,npf ,tf ,
253 . area ,exx ,eyy ,exy ,imat )
254 ENDDO
255 ENDIF
256c-----------
257 RETURN
cmlawi
subroutine cmlawi(mat_param, jft, jlt, ilaw, ithk, pm, for, thk, eint, gstr, dir1, dir2, uvar, ipm, nel, sigi, npf, tf, area, exx, eyy, exy, imat)
Definition cepsini.F:279
cdefoi
subroutine cdefoi(jft, jlt, vl, gstr, px1, px2, py1, py2, nel, area, exx, eyy, exy, eyz, ezx, e1x, e2x, e3x, e1y, e2y, e3y, e1z, e2z, e3z)
Definition cepsini.F:433
cpxpyi
subroutine cpxpyi(jft, jlt, ismstr, area, px1g, px2g, py1g, py2g, px1, px2, py1, py2, x1g, x2g, x3g, x4g, y1g, y2g, y3g, y4g, z1g, z2g, z3g, z4g, e1x, e2x, e3x, e1y, e2y, e3y, e1z, e2z, e3z, x2l, y2l, x3l, y3l, x4l, y4l)
Definition cepsini.F:350
ceveci
subroutine ceveci(jft, jlt, area, x1, x2, x3, x4, y1, y2, y3, y4, z1, z2, z3, z4, e1x, e2x, e3x, e1y, e2y, e3y, e1z, e2z, e3z)
Definition ceveci.F:37
corth3
subroutine corth3(elbuf_str, dir_a, dir_b, jft, jlt, nlay, irep, nel, x1, x2, x3, x4, y1, y2, y3, y4, z1, z2, z3, z4, e1x, e2x, e3x, e1y, e2y, e3y, e1z, e2z, e3z, idrape, igtyp)
Definition corth3.F:42
for
for(i8=*sizetab-1;i8 >=0;i8--)
Definition mumps_common.c:91

◆ cm19ini()

subroutine cm19ini ( integer jft,
integer jlt,
integer ithk,
pm,
for,
thk,
eint,
gstr,
dir,
sigi,
integer nel,
integer imat,
area,
exx,
eyy,
exy )

Definition at line 596 of file cepsini.F.

599C-----------------------------------------------
600C I m p l i c i t T y p e s
601C-----------------------------------------------
602#include "implicit_f.inc"
603C-----------------------------------------------
604C G l o b a l P a r a m e t e r s
605C-----------------------------------------------
606#include "mvsiz_p.inc"
607C-----------------------------------------------
608C C o m m o n B l o c k s
609C-----------------------------------------------
610#include "param_c.inc"
611C-----------------------------------------------
612C D u m m y A r g u m e n t s
613C-----------------------------------------------
614 INTEGER JFT, JLT, ITHK,NEL,IMAT
615 my_real
616 . pm(npropm,*) ,for(nel,5) ,thk(*),
617 . area(mvsiz),exx(mvsiz) , eyy(mvsiz) , exy(mvsiz) ,
618 . eint(nel,2) ,gstr(nel,8) ,dir(nel,*), sigi(nel,3)
619C-----------------------------------------------
620C L o c a l V a r i a b l e s
621C-----------------------------------------------
622 INTEGER I
623 my_real
624 . degmb(mvsiz), epst(5,mvsiz),
625 . e11(mvsiz), e22(mvsiz), n12(mvsiz), n21(mvsiz),
626 . g12(mvsiz), de(mvsiz) ,
627 . a11(mvsiz), a12(mvsiz), a22(mvsiz)
628 my_real
629 . ezz, scale, zerostress, s, d, r, p1, p2, beta
630C-------------------------------------------------------------
631 DO i=jft,jlt
632 e11(i) =pm(33,imat)
633 e22(i) =pm(34,imat)
634 n12(i) =pm(35,imat)
635 n21(i) =pm(36,imat)
636 g12(i) =pm(37,imat)
637 de(i) =pm(44,imat)
638 ENDDO
639C
640 DO i=jft,jlt
641 degmb(i) = for(i,1)*exx(i)+for(i,2)*eyy(i)+for(i,3)*exy(i)
642 ENDDO
643C
644 DO i=jft,jlt
645 epst(1,i)= gstr(i,1)
646 epst(2,i)= gstr(i,2)
647 epst(3,i)= gstr(i,3)*half
648 epst(4,i)= gstr(i,4)*half
649 epst(5,i)= gstr(i,5)*half
650 ENDDO
651C-------------------------------------------------
652C ROTATIONS DANS LE REPERE D ORTHOTROPIE EPST
653C-------------------------------------------------
654 CALL roto(jft,jlt,epst,5,dir,nel)
655C-------------------------------------------------
656 DO i=jft,jlt
657 a12(i)=one-n12(i)*n21(i)
658 a11(i)=e11(i)/a12(i)
659 a22(i)=e22(i)/a12(i)
660 a12(i)=n21(i)*a11(i)
661C-----------------------------
662C FORCES ELASTIQUES
663C-----------------------------
664 for(i,1)=a11(i)*epst(1,i)+a12(i)*epst(2,i)
665 for(i,2)=a12(i)*epst(1,i)+a22(i)*epst(2,i)
666 for(i,3)=g12(i)*epst(3,i)*two
667C-----------------------------
668C CONTRAINTES reduites en compression
669C-----------------------------
670 s = half*(for(i,1)+for(i,2))
671 d = half*(for(i,1)-for(i,2))
672 r = sqrt(for(i,3)*for(i,3) + d*d)
673 p1 = s - r
674 IF(p1<zero)THEN
675 p2 = s + r
676 IF(p2>zero)THEN
677 beta = half*((one-de(i))*s/r + one+de(i))
678 for(i,1)=beta*(for(i,1)-p2) + p2
679 for(i,2)=beta*(for(i,2)-p2) + p2
680 for(i,3)=beta*for(i,3)
681 ELSE
682 beta = de(i)
683 for(i,1)=beta*for(i,1)
684 for(i,2)=beta*for(i,2)
685 for(i,3)=beta*for(i,3)
686 ENDIF
687 ENDIF
688 ENDDO
689C-----------------------------------------------------------
690C REF-STATE ZEROSTRESS OPTION
691C-----------------------------------------------------------
692 DO i=jft,jlt
693 zerostress =pm(55,imat)
694 IF(zerostress /= zero)THEN
695 sigi(i,1) = for(i,1)
696 sigi(i,2) = for(i,2)
697 sigi(i,3) = for(i,3)
698 ENDIF
699 ENDDO
700C--------------------------------
701C RETOUR DANS LE REPERE COQUE
702C--------------------------------
703 CALL uroto(jft,jlt,for,5,dir,nel)
704C
705 DO i=jft,jlt
706 degmb(i) = degmb(i)+
707 + for(i,1)*exx(i)+for(i,2)*eyy(i)+for(i,3)*exy(i)
708 eint(i,1) = eint(i,1) + degmb(i)*half*thk(i)*area(i)
709 ENDDO
710C IF(ITHK == 1) THEN
711C DO I=JFT,JLT
712C EZZ = -NU(I) * (EXX(I) + EYY(I)) / (ONE-NU(I))
713C THK(I) = THK(I) * (ONE + EZZ)
714C ENDDO
715C ENDIF
716 RETURN
roto
subroutine roto(jft, jlt, tab, ltab, dir, nel)
Definition cepsini.F:724
uroto
subroutine uroto(jft, jlt, tab, ltab, dir, nel)
Definition cepsini.F:774

◆ cm1ini()

subroutine cm1ini ( integer jft,
integer jlt,
integer ithk,
pm,
for,
thk,
eint,
integer nel,
integer imat,
area,
exx,
eyy,
exy )

Definition at line 530 of file cepsini.F.

533C-----------------------------------------------
534C I m p l i c i t T y p e s
535C-----------------------------------------------
536#include "implicit_f.inc"
537C-----------------------------------------------
538C G l o b a l P a r a m e t e r s
539C-----------------------------------------------
540#include "mvsiz_p.inc"
541C-----------------------------------------------
542C C o m m o n B l o c k s
543C-----------------------------------------------
544#include "param_c.inc"
545C-----------------------------------------------
546C D u m m y A r g u m e n t s
547C-----------------------------------------------
548 INTEGER JFT, JLT, ITHK,NEL,IMAT
549 my_real
550 . exx(mvsiz) , eyy(mvsiz) , exy(mvsiz) ,area(mvsiz),
551 . pm(npropm,*) ,for(nel,5) ,thk(*),eint(nel,2)
552C-----------------------------------------------
553C L o c a l V a r i a b l e s
554C-----------------------------------------------
555 INTEGER I
556 my_real ezz
557 my_real nu(mvsiz),g(mvsiz),a1(mvsiz),a2(mvsiz),degmb(mvsiz)
558C-------------------------------------------------------------
559 DO i=jft,jlt
560 nu(i) = pm(21,imat)
561 g(i) = pm(22,imat)
562 a1(i) = pm(24,imat)
563 a2(i) = pm(25,imat)
564 ENDDO
565 DO i=jft,jlt
566 degmb(i) = for(i,1)*exx(i)+for(i,2)*eyy(i)+for(i,3)*exy(i)
567 ENDDO
568 DO i=jft,jlt
569 for(i,1)=for(i,1)+ a1(i)*exx(i)+a2(i)*eyy(i)
570 for(i,2)=for(i,2)+ a1(i)*eyy(i)+a2(i)*exx(i)
571 for(i,3)=for(i,3)+ g(i)*exy(i)
572 ENDDO
573 DO i=jft,jlt
574 degmb(i) = degmb(i)+
575 + for(i,1)*exx(i)+for(i,2)*eyy(i)+for(i,3)*exy(i)
576 eint(i,1) = eint(i,1) + degmb(i)*half*thk(i)*area(i)
577 ENDDO
578 IF(ithk == 1) THEN
579 DO i=jft,jlt
580 ezz = -nu(i) * (exx(i) + eyy(i)) / (one-nu(i))
581 thk(i) = thk(i) * (one + ezz)
582 ENDDO
583 ENDIF
584c-----------
585 RETURN

◆ cmlawi()

subroutine cmlawi ( type (matparam_struct_), intent(in) mat_param,
integer jft,
integer jlt,
integer ilaw,
integer ithk,
pm,
for,
thk,
eint,
gstr,
dir1,
dir2,
uvar,
integer, dimension(npropmi,*) ipm,
integer nel,
sigi,
integer, dimension(*) npf,
tf,
area,
exx,
eyy,
exy,
integer imat )

Definition at line 273 of file cepsini.F.

279C-----------------------------------------------
280C M o d u l e s
281C-----------------------------------------------
282 USE mat_elem_mod
283C-----------------------------------------------
284C I m p l i c i t T y p e s
285C-----------------------------------------------
286#include "implicit_f.inc"
287C-----------------------------------------------
288C C o m m o n B l o c k s
289C-----------------------------------------------
290#include "param_c.inc"
291C-----------------------------------------------
292C D u m m y A r g u m e n t s
293C-----------------------------------------------
294 INTEGER IPM(NPROPMI,*),NPF(*)
295 INTEGER JFT, JLT, ITHK, ILAW,NEL,IMAT
296 my_real
297 . pm(*),for(nel,5),thk(*),eint(nel,2),dir1(nel,*),dir2(nel,*),
298 . uvar(nel,*),sigi(nel,3), tf(*),gstr(*),
299 . area(nel),exx(nel),eyy(nel),exy(nel)
300 TYPE (ELBUF_STRUCT_), TARGET :: ELBUF_STR
301 TYPE (MATPARAM_STRUCT_) ,INTENT(IN) :: MAT_PARAM
302C-----------------------------------------------
303C L o c a l V a r i a b l e s
304C-----------------------------------------------
305 INTEGER I,J,NUVAR,NUPARAM,NIPARAM,NBFUNC
306 INTEGER IFUNC(100)
307C=======================================================================
308C INITIALIZE STRESS,THICKNESS & ENERGY
309C-------------------------------------
310 IF (ilaw == 1) THEN
311 CALL cm1ini(jft,jlt,ithk,pm,for,
312 . thk,eint,nel,imat,area,
313 . exx,eyy ,exy )
314 ELSEIF (ilaw == 19) THEN
315 CALL cm19ini(jft,jlt,ithk,pm,for,
316 . thk,eint,gstr,dir1,sigi,
317 . nel,imat,area,exx,eyy ,exy )
318 ELSEIF (ilaw == 58) THEN
319 nuvar = ipm(8,imat)
320 nbfunc= ipm(10,imat)
321 DO j= 1,nbfunc
322 ifunc(j) = ipm(10+j,imat)
323 ENDDO
324 nuparam = mat_param%NUPARAM
325 niparam = mat_param%NIPARAM
326 CALL cm58_refsta(nel ,nuparam,niparam,mat_param%UPARAM,mat_param%IPARAM,
327 1 for ,eint ,gstr ,
328 2 thk ,dir1 ,dir2 ,nuvar ,uvar ,
329 3 nbfunc ,ifunc ,npf ,tf ,area ,
330 4 exx ,eyy ,exy )
331!
332 ENDIF
333c-----------
334 RETURN
cm1ini
subroutine cm1ini(jft, jlt, ithk, pm, for, thk, eint, nel, imat, area, exx, eyy, exy)
Definition cepsini.F:533
cm19ini
subroutine cm19ini(jft, jlt, ithk, pm, for, thk, eint, gstr, dir, sigi, nel, imat, area, exx, eyy, exy)
Definition cepsini.F:599
cm58_refsta
subroutine cm58_refsta(nel, nuparam, niparam, uparam, iparam, for, eint, gstr, thk, dir1, dir2, nuvar, uvar, nbfunc, ifunc, npf, tf, area, exx, eyy, exy)
Definition cm58_refsta.F:35

◆ cpxpyi()

subroutine cpxpyi ( integer jft,
integer jlt,
integer ismstr,
intent(out) area,
px1g,
px2g,
py1g,
py2g,
px1,
px2,
py1,
py2,
intent(in) x1g,
intent(in) x2g,
intent(in) x3g,
intent(in) x4g,
intent(in) y1g,
intent(in) y2g,
intent(in) y3g,
intent(in) y4g,
intent(in) z1g,
intent(in) z2g,
intent(in) z3g,
intent(in) z4g,
intent(in) e1x,
intent(in) e2x,
intent(in) e3x,
intent(in) e1y,
intent(in) e2y,
intent(in) e3y,
intent(in) e1z,
intent(in) e2z,
intent(in) e3z,
x2l,
y2l,
x3l,
y3l,
x4l,
y4l )

Definition at line 343 of file cepsini.F.

350C-----------------------------------------------
351C I m p l i c i t T y p e s
352C-----------------------------------------------
353#include "implicit_f.inc"
354C-----------------------------------------------
355C G l o b a l P a r a m e t e r s
356C-----------------------------------------------
357#include "mvsiz_p.inc"
358C-----------------------------------------------
359C D u m m y A r g u m e n t s
360C-----------------------------------------------
361 INTEGER JFT, JLT, ISMSTR
362 my_real px1(*), px2(*), py1(*), py2(*),
363 . px1g(*), px2g(*), py1g(*), py2g(*),
364 . x2l(*) ,y2l(*),x3l(*),y3l(*),x4l(*),y4l(*)
365 my_real, DIMENSION(MVSIZ), INTENT(IN) ::
366 . x1g,x2g,x3g,x4g,y1g,y2g,y3g,y4g,z1g,z2g,z3g,z4g,
367 . e1x, e2x, e3x, e1y, e2y, e3y ,e1z, e2z, e3z
368 my_real, DIMENSION(MVSIZ) , INTENT(OUT) :: area
369C-----------------------------------------------
370C L o c a l V a r i a b l e s
371C-----------------------------------------------
372 INTEGER I
373 my_real, DIMENSION(MVSIZ) ::
374 . x21g,y21g,z21g,x31g,y31g,z31g,x41g, y41g, z41g
375C=======================================================================
376 IF (ismstr /= 3)THEN
377 DO i=jft,jlt
378 x21g(i)=x2g(i)-x1g(i)
379 y21g(i)=y2g(i)-y1g(i)
380 z21g(i)=z2g(i)-z1g(i)
381 x31g(i)=x3g(i)-x1g(i)
382 y31g(i)=y3g(i)-y1g(i)
383 z31g(i)=z3g(i)-z1g(i)
384 x41g(i)=x4g(i)-x1g(i)
385 y41g(i)=y4g(i)-y1g(i)
386 z41g(i)=z4g(i)-z1g(i)
387 ENDDO
388C
389 DO i=jft,jlt
390 x2l(i)=e1x(i)*x21g(i)+e1y(i)*y21g(i)+e1z(i)*z21g(i)
391 y2l(i)=e2x(i)*x21g(i)+e2y(i)*y21g(i)+e2z(i)*z21g(i)
392 y3l(i)=e2x(i)*x31g(i)+e2y(i)*y31g(i)+e2z(i)*z31g(i)
393 x3l(i)=e1x(i)*x31g(i)+e1y(i)*y31g(i)+e1z(i)*z31g(i)
394 x4l(i)=e1x(i)*x41g(i)+e1y(i)*y41g(i)+e1z(i)*z41g(i)
395 y4l(i)=e2x(i)*x41g(i)+e2y(i)*y41g(i)+e2z(i)*z41g(i)
396 ENDDO
397C
398 DO i=jft,jlt
399 px1(i)= half*(y2l(i)-y4l(i))
400 py1(i)= half*(x4l(i)-x2l(i))
401 px2(i)= half*y3l(i)
402 py2(i)=-half*x3l(i)
403 ENDDO
404C
405 ELSE
406C
407 DO i=jft,jlt
408 px1(i) = px1g(i)
409 px2(i) = px2g(i)
410 py1(i) = py1g(i)
411 py2(i) = py2g(i)
412 ENDDO
413C
414 ENDIF
415C
416 DO i=jft,jlt
417 area(i)=two*(py2(i)*px1(i)-py1(i)*px2(i))
418 ENDDO
419C
420 RETURN

◆ roto()

subroutine roto ( integer jft,
integer jlt,
tab,
integer ltab,
dir,
integer nel )

Definition at line 723 of file cepsini.F.

724C-----------------------------------------------
725C I m p l i c i t T y p e s
726C-----------------------------------------------
727#include "implicit_f.inc"
728C-----------------------------------------------
729C G l o b a l P a r a m e t e r s
730C-----------------------------------------------
731#include "mvsiz_p.inc"
732C-----------------------------------------------
733C D u m m y A r g u m e n t s
734C-----------------------------------------------
735 INTEGER JFT, JLT,NEL,LTAB
736 my_real tab(ltab,*), dir(nel,*)
737C-----------------------------------------------
738C L o c a l V a r i a b l e s
739C-----------------------------------------------
740 INTEGER I
741 my_real
742 . s1(mvsiz), s2(mvsiz), s3(mvsiz), s4(mvsiz), s5(mvsiz)
743C---+----1----+----2----+----3----+----4----+----5----+----6----+----7--
744 DO i=jft,jlt
745 s1(i) = dir(i,1)*dir(i,1)*tab(1,i)
746 . +dir(i,2)*dir(i,2)*tab(2,i)
747 . +two*dir(i,1)*dir(i,2)*tab(3,i)
748 s2(i) = dir(i,2)*dir(i,2)*tab(1,i)
749 . +dir(i,1)*dir(i,1)*tab(2,i)
750 . -two*dir(i,2)*dir(i,1)*tab(3,i)
751 s3(i) =-dir(i,1)*dir(i,2)*tab(1,i)
752 . +dir(i,2)*dir(i,1)*tab(2,i)
753 . +(dir(i,1)*dir(i,1)-dir(i,2)*dir(i,2))*tab(3,i)
754 s4(i) =-dir(i,2)*tab(5,i)+dir(i,1)*tab(4,i)
755 s5(i) = dir(i,1)*tab(5,i)+dir(i,2)*tab(4,i)
756 ENDDO
757C
758 DO i=jft,jlt
759 tab(1,i)=s1(i)
760 tab(2,i)=s2(i)
761 tab(3,i)=s3(i)
762 tab(4,i)=s4(i)
763 tab(5,i)=s5(i)
764 ENDDO
765C
766 RETURN

◆ uroto()

subroutine uroto ( integer jft,
integer jlt,
tab,
integer ltab,
dir,
integer nel )

Definition at line 773 of file cepsini.F.

774C-----------------------------------------------
775C I m p l i c i t T y p e s
776C-----------------------------------------------
777#include "implicit_f.inc"
778C-----------------------------------------------
779C G l o b a l P a r a m e t e r s
780C-----------------------------------------------
781#include "mvsiz_p.inc"
782C-----------------------------------------------
783C D u m m y A r g u m e n t s
784C-----------------------------------------------
785 INTEGER JFT, JLT,NEL,LTAB
786 my_real tab(nel,ltab), dir(nel,*)
787C-----------------------------------------------
788C L o c a l V a r i a b l e s
789C-----------------------------------------------
790 INTEGER I
791 my_real
792 . s1(mvsiz), s2(mvsiz), s3(mvsiz), s4(mvsiz), s5(mvsiz)
793C-------------------------------------------------------------
794 DO i=jft,jlt
795 s1(i) = dir(i,1)*dir(i,1)*tab(i,1)
796 . +dir(i,2)*dir(i,2)*tab(i,2)
797 . -two*dir(i,1)*dir(i,2)*tab(i,3)
798 s2(i) = dir(i,2)*dir(i,2)*tab(i,1)
799 . +dir(i,1)*dir(i,1)*tab(i,2)
800 . +two*dir(i,2)*dir(i,1)*tab(i,3)
801 s3(i) =+dir(i,1)*dir(i,2)*tab(i,1)
802 . -dir(i,2)*dir(i,1)*tab(i,2)
803 . +(dir(i,1)*dir(i,1)-dir(i,2)*dir(i,2))*tab(i,3)
804 s4(i) = dir(i,2)*tab(i,5)+dir(i,1)*tab(i,4)
805 s5(i) = dir(i,1)*tab(i,5)-dir(i,2)*tab(i,4)
806 ENDDO
807 DO i=jft,jlt
808 tab(i,1)=s1(i)
809 tab(i,2)=s2(i)
810 tab(i,3)=s3(i)
811 tab(i,4)=s4(i)
812 tab(i,5)=s5(i)
813 ENDDO
814C
815 RETURN