1 : C HS456 SOURCE CHAT 05/01/13 00:26:29 5004
2 : SUBROUTINE HS456(XE,NSTRS,DDHOMU,HS4,HS5,HS6)
3 : *
4 : * CALCUL DES FONCTIONS HS4,HS5 ET HS6
5 : * POUR L'ELEMENT COQUE DST
6 : *
7 : ************************************************************************
8 : *
9 : * ENTREES :
10 : *
11 : * XE : COORDONNEES DES NOEUDS
12 : * NSTRS :NOMBRE DE COMPOSANTES DE CONTRAINTES
13 : * DDHOMU :MATRICE DE HOOKE
14 : *
15 : * SORTIES :
16 : *
17 : * HS4,HS5,HS6 : LES FONCTIONS QUI RELIE LES TERMES (BETA)SK A {U}
18 : *
19 : * AUTEUR : X.Z. SUO FEVRIER 96
20 : ************************************************************************
21 : IMPLICIT INTEGER(I-N)
22 : IMPLICIT REAL*8(A-H,O-Z)
23 : DIMENSION XE(3,*),DDHOMU(NSTRS,*),HS4(*),HS5(*),HS6(*)
24 : DIMENSION SX(3),SY(3),RL(3),EXX(3),EYY(3),HFX(9),HFY(9)
25 : DIMENSION ASS(6),BSS(6),AXX(3),BXX(3),AYY(3),BYY(3)
26 : *
27 : DO 10 K=4,6
28 : IF(K.EQ.4)THEN
29 : IJ=1
30 : I=2
31 : J=3
32 : ELSEIF(K.EQ.5)THEN
33 : IJ=2
34 : I=3
35 : J=1
36 : ELSE
37 : IJ=3
38 : I=1
39 : J=2
40 : ENDIF
41 : SX(IJ)=XE(1,I)-XE(1,J)
42 : SY(IJ)=XE(2,I)-XE(2,J)
43 : RL(IJ)=SQRT(SX(IJ)*SX(IJ)+SY(IJ)*SY(IJ))
44 : EXX(IJ)=-SX(IJ)/RL(IJ)
45 : EYY(IJ)=-SY(IJ)/RL(IJ)
46 : 10 CONTINUE
47 : AIR=ABS(0.5D0*(SX(1)*SY(2)-SX(2)*SY(1)))
48 : SOMM=SX(IJ)+SY(IJ)+RL(IJ)+EXX(IJ)+EYY(IJ)
49 : *
49 : CALL CISAF(XE,EXX,EYY,DDHOMU,AIR,ASS,AXX,
50 : 1 AYY,BSS,BXX,BYY)
51 : *
51 : CALL CISAR(ASS,BSS,AXX,AYY,BXX,BYY,RL,EXX,EYY,
52 : 1 HS4,HS5,HS6,HFX,HFY)
53 : *
54 : RETURN
55 : END
56 :
57 :
© Cast3M 2003 - All rights reserved.
Disclaimer