1 : C SIGMA SOURCE CHAT 11/05/24 21:16:12 6978
2 : SUBROUTINE SIGMA
3 : *_______________________________________________________________________
4 : *
5 : * OPERATEUR CONTRAINTES LINEAIRES
6 : *
7 : * Ancienne syntaxe :
8 : * __________________
9 : *
10 : * SIG1=SIGMA ( II ) MAT1 ( CAR1 ) CHP1 ;
11 : *
12 : * II mot cle indiquant que l'on veut EPS2
13 : * MAT1 champ par {l{ment de mat{riau ou de HOOKE
14 : * CAR1 champ par {l{ment de caract{ristiques (facultatif)
15 : * CHP1 CHPOINT de d{placement
16 : *
17 : * JUILLET 84 EBERSOLT
18 : * FEVRIER 87 BROCHARD :
19 : * - LE CALCUL EFFECTIF EST FAIT DANS SIGMAP
20 : *
21 : *
22 : * Nouvelle syntaxe :
23 : * __________________
24 : *
25 : * SIG1=SIGMA ( II ) MODL1 | CAR1 CHP1 ;
26 : * | HOO1 (CAR1)
27 : *
28 : * II mot cle indiquant que l'on veut EPS2
29 : * MODL1 mod}le de calcul, type MMODEL
30 : * CAR1 champ par {l{ment de caract{ristiques, type MCHAML
31 : * HOO1 champ par {l{ment de matrice de HOOKE, type MCHAML
32 : * CHP1 CHPOINT de d{placement
33 : *
34 : * Passage aux nouveaux CHAMELEMs par I.Monnier le 6 6 1990
35 : *
36 : *_______________________________________________________________________
37 : IMPLICIT INTEGER(I-N)
38 : IMPLICIT REAL*8(A-H,O-Z)
39 : PARAMETER(XZER=0.D0,UNDEMI=.5D0,UN=1.D0)
40 : -INC CCOPTIO
41 : -INC SMCHAML
42 : C
43 : CHARACTER*4 MOEPS2(1)
44 : C
45 : DATA MOEPS2/'LINE'/
46 : C
47 : C IPMODL MODELE MMODEL
48 : C ICHP1 CHAMPOINT DEPLACEMENTS
49 : C IPCHA1 CHAMELEM MATERIAU OU DE HOOKE
50 : C IPCHA2 CHAMELEM CARACTERISTIQUE
51 : C IPSTRS CHAMELEM CONTRAINTES
52 : C 1 SI SIGAM DIRECTIVE 0 SINON
53 : IPCHE1=0
54 : IPCHE2=0
55 : IPCHA1=0
56 : IPCHA2=0
57 : IPCHB1=0
58 : IPCHB2=0
59 : ICONV =0
60 : IRT9=0
61 : C
62 : C VEUT ON CALCULER LES GRANDES DEFORMATIONS ?
63 : C
63 : CALL LIRMOT(MOEPS2,1,IREPS2,0)
64 : C
65 : C ON LIT UN CHAMP/POINT
66 : C
66 : CALL LIROBJ('CHPOINT ',ICHP1,1,IRT1)
67 : IF(IERR.NE.0) RETURN
68 : C
69 : C LECTURE DU MODELE
70 : C
70 : CALL LIROBJ('MMODEL',IPMODL,1,IRT1)
71 : IF(IERR.NE.0) RETURN
72 : C
73 : C LECTURE DU 1ER MCHAML
74 : C
74 : CALL LIROBJ('MCHAML',IPCHA1,1,IRT1)
75 : IF(IERR.NE.0) RETURN
76 : * reduction du mchaml
76 : call reduaf(ipcha1,ipmodl,ipcha10,0,ier,ker)
76 : if (ier.ne.1) call erreur(ker)
77 : ipcha1=ipcha10
78 : C
79 : C LECTURE DU 2EME MCHAML
80 : C
80 : CALL LIROBJ('MCHAML',IPCHA2,0,IRT1)
81 : if (irt1.ne.0) then
82 : * reduction du mchaml
82 : call reduaf(ipcha2,ipmodl,ipcha20,0,ier,ker)
82 : if (ier.ne.1) call erreur(ker)
83 : ipcha2=ipcha20
84 : endif
85 : if (ierr.ne.0) return
86 : C
86 : CALL RNGCHA(IPCHA1,IPCHA2,'CARACTERISTIQUES',
87 : 1 'MATRICE DE HOOKE',IPCHE1,IPCHE2)
88 : IF(IERR.NE.0) RETURN
89 : IF (IPCHE2.EQ.0) THEN
90 : IMAT=1
91 : ELSE
92 : IMAT=2
93 : ENDIF
94 : C
95 : C CALCUL DES CONTRAINTES
96 : C
97 : IRT8=0
97 : CALL SIGMAP(IREPS2,IPMODL,ICHP1,IPCHE1,IPCHE2,ICONV,IMAT,IPSTRS,
98 : 1 IRET)
99 : C
100 : C ECRITURE DU RESULTAT
101 : C
102 : IF(IRT9.EQ.0.AND.IRET.EQ.1) THEN
102 : CALL ECROBJ('MCHAML',IPSTRS)
103 : ENDIF
104 : RETURN
105 : END
106 :
107 :
108 :
109 :
110 :
111 :
© Cast3M 2003 - All rights reserved.
Disclaimer