1 : C BSIGMA SOURCE PV 09/03/12 21:15:56 6325
2 : SUBROUTINE BSIGMA
3 : C_______________________________________________________________________
4 : C
5 : C OPERATEUR FORCES INTERNES
6 : C
7 : C FOR1 = BSIGMA MODL1 SIG1 ( CAR1 ) (HOO1) ;
8 : C
9 : C MODL1 objet de type MMODEL
10 : C SIG1 MCHAML de contraintes
11 : C CAR1 MCHAML de caract{ristiques (facultatif)
12 : C HOO1 MCHAML DE MATERIAU OU DE HOOKE (FACULTATIF)
13 : C FOR1 CHPOINT donnant les foces nodales
14 : C___________________________________________________________________
15 : C
16 : IMPLICIT INTEGER(I-N)
17 : IMPLICIT REAL*8(A-H,O-Z)
18 : -INC CCOPTIO
19 : -INC SMCHAML
20 : -INC SMCHPOI
21 : C
22 : ICONV=0
23 : IPCHE1=0
24 : IPCHE2=0
25 : IPCHE3=0
26 : IPCHA2=0
27 : IPCHA3=0
27 : CALL LIROBJ('MMODEL',IPMODL,1,IRT1)
28 : IF (IERR.NE.0) RETURN
29 : C
30 : C 1 ER CHAMP/ELEMENT
31 : C
31 : CALL LIROBJ('MCHAML',IPCHE1,1,IRT1)
32 : IF (IERR.NE.0) RETURN
32 : call reduaf(ipche1,ipmodl,ipch,0,ir,ker)
32 : if (ir.ne.1) call erreur(ker)
33 : IF (IERR.NE.0) RETURN
34 : ipche1=ipch
35 : C
36 : C 2 EME CHAMP/ELEMENT (FACULTATIF)
37 : C
37 : CALL LIROBJ('MCHAML',IPCHA2,0,IRT1)
38 : IF (IERR.NE.0) RETURN
39 : if (irt1.ne.0) then
39 : call reduaf(ipcha2,ipmodl,ipch,0,ir,ker)
39 : if (ir.ne.1) call erreur(ker)
40 : IF (IERR.NE.0) RETURN
41 : ipcha2=ipch
42 : endif
43 : C
44 : C 3 EME CHAMP/ELEMENT (FACULTATIF)
45 : C
45 : CALL LIROBJ('MCHAML',IPCHA3,0,IRT1)
46 : IF (IERR.NE.0) RETURN
47 : if (irt1.ne.0) then
47 : call reduaf(ipcha3,ipmodl,ipch,0,ir,ker)
47 : if (ir.ne.1) call erreur(ker)
48 : IF (IERR.NE.0) RETURN
49 : ipcha3=ipch
50 : endif
51 : C
52 : IF(IPCHA2.NE.0.OR.IPCHA3.NE.0)THEN
52 : CALL RNGCHA(IPCHA2,IPCHA3,'CARACTERISTIQUES',
53 : 1 'MATRICE DE HOOKE',IPCHE2,IPCHE3)
54 : IF(IERR.NE.0) GOTO 666
55 : IF (IPCHE3.EQ.0)THEN
56 : IMAT=1
57 : ELSE
58 : IMAT=2
59 : ENDIF
60 : ELSE
61 : IMAT=0
62 : IPCHE2=0
63 : IPCHE3=0
64 : ENDIF
65 : C
65 : CALL BSIGMP(IPMODL,IPCHE1,IPCHE2,IPCHE3,IMAT,
66 : 1 ICONV,IPCHE4,IRET)
67 : if (ierr.ne.0) return
68 : C
69 : IF (IRET.EQ.1) THEN
70 : C
71 : C ATTRIBUTION D'UNE NATURE DISCRETE AU CHPO QUI SORT
72 : C
73 : MCHPOI = IPCHE4
74 : SEGACT MCHPOI*MOD
75 : JATTRI(1) = 2
76 : SEGDES MCHPOI
77 : C
77 : CALL ECROBJ('CHPOINT',IPCHE4)
78 : ENDIF
79 : C
80 : 666 CONTINUE
81 : RETURN
82 : END
83 :
84 :
85 :
86 :
87 :
© Cast3M 2003 - All rights reserved.
Disclaimer