1 : C KP SOURCE CHAT 06/08/24 21:48:35 5529
2 : SUBROUTINE KP
3 : C_______________________________________________________________________
4 : C
5 : C OPERATEUR KP
6 : C
7 : C Construction de la matrice des pressions associees à la linearisation
8 : C des actions de pression sur la géometrie deformée.
9 : C
10 : C
11 : C Nouvelle syntaxe:
12 : C _________________
13 : *
14 : * 1ere option
15 : C
16 : C RIG1 =KP MODL1 PRES1 ('FLAM') ('ASYM');
17 : C
18 : C MODL1 = objet modele de type MMODEL
19 : C
20 : C PRES1=MCHAML ou CHPOINT de PRESSION
21 : C
22 : C 'FLAM' = MOT-CLE N{CESSAIRE SI L'ON VEUT UTILISER LA MATRICE
23 : C POUR FAIRE UNE ANALYSE DE FLAMBAGE.
24 : C
25 : C RIG1=Matrice KP de sous type RIGIDITE
26 : *
27 : * 2ème option
28 : *
29 : * RIG1=KP MODL1 RG (VEC1) ('ASYM')
30 : *
31 : * MOD1 : MODELE DE CALCUL, TYPE MMODEL
32 : * RG : COEFFICIENT MULTIPLICATIF DE RIGIDITE
33 : * VEC1 : vecteur (type point) qui definit le sens de
34 : * la "pesanTeur"(facultatif). S'il n'est pas
35 : * fourni vec1 est consideré suivant la normale
36 : * à chaque élément.
37 : * ASYM : mot clé pour avoir la matrice dissymetrique
38 : *
39 : * la deuxième version CALCULE LA MATRICE DE RAIDEUR LIEE A LA
40 : * VARIATION DE PRESSION DUE AU MOUVEMENT VERTICAL DANS UN CHAMP DE
41 : * PESANTEUR (on peut l'utiliser dans d'autres cas aussi, p.ex. plaque
42 : * sur sol elastique)
43 : C
44 : C ces matrices sont en general non symetriques, mais leur combinaison peut
45 : c etre symetrique (p.ex. en cas de pression hydrostatique). C'est pourquoi
46 : c nous voulons avoir la possibilité de construire les matrices symetriques
47 : c ou disymetriques.
48 : C
49 : C
50 : C CODE COMBESCURE JANV 87
51 : C
52 : C Passage aux nouveaux CHAMELEMs par P.DOWLATYARI le 5/4/91
53 : C
54 : C Ajout de la 2ème option I. Politopoulos Aout 1995
55 : C_______________________________________________________________________
56 : C
57 : IMPLICIT INTEGER(I-N)
58 : IMPLICIT REAL*8(A-H,O-Z)
59 : C
60 : -INC CCOPTIO
61 : C
62 : C
63 : C
64 : IPMODL=0
65 : IPCHP=0
66 : IPRIG=0
67 : IFLAM=0
68 : ICHA=0
69 : IDISS=0
70 :
71 : C
71 : CALL LIRMOT('FLAM',1,IFLAM,0)
71 : IF (IERR.NE.0) CALL ERREUR(7)
72 : C
72 : CALL LIRMOT('ASYM',1,IDISS,0)
72 : IF (IERR.NE.0) CALL ERREUR(7)
73 :
74 : C
74 : CALL LIROBJ('MMODEL',IPMODL,1,IRT1)
75 : IF (IERR.NE.0) RETURN
76 :
76 : CALL LIROBJ('MCHAML',IPCHP,0,IRT2)
77 : IF (IERR.NE.0) RETURN
78 : IF(IRT2.EQ.1) THEN
79 : ICHA=1
80 : ELSE
80 : CALL LIROBJ('CHPOINT ',IPCHP,0,IRT3)
81 : IF (IERR.NE.0) RETURN
82 : IF(IRT3.EQ.0) GO TO 10
83 : ENDIF
84 : C
84 : CALL KPRES(IPMODL,IPCHP,ICHA,ICONV,IFLAM,IDISS,IPRIG,IRET)
85 : GOTO 20
86 : C
87 : *** matrice associee à un gradient de pression imposée
88 : 10 CONTINUE
88 : CALL LIRREE(XRG,1,IRET)
89 : IF (IRET.EQ.0) GOTO 666
89 : CALL LIROBJ('POINT',IPOIN1,0,IRET2)
90 :
90 : CALL KGRAV1 (IPMODL,XRG,IPOIN1,IPRIG,IRET,IDISS,IFLAM)
91 : *
92 : * ecriture de la rigidite
93 :
94 : 20 CONTINUE
95 : IF (IRET.EQ.1) THEN
95 : CALL ECROBJ('RIGIDITE',IPRIG)
96 : ENDIF
97 :
98 : 666 CONTINUE
99 : RETURN
100 : END
101 :
102 :
103 :
104 :
105 :
106 :
© Cast3M 2003 - All rights reserved.
Disclaimer