Notices | Cast3M

   1 : $$$$ EFFMARTI NOTICE  AF221230  16/11/03    21:15:02     9139           
   2 :                                              DATE     16/11/03
   3 : 
   4 :   procedure EFF_MARTI                       Voir aussi : 
   5 :     ----------------  
   6 :    SIG2 =  EFF_MARTI SIG1 MOD1 MAT1 VECT1 VECT2
   7 :                      H1 ENRE1 ENRI1 COT1
   8 : 
   9 : 
  10 : 
  11 : Objet :
  12 : -------
  13 : 
  14 : Cette procedure transforme les effots generalisés sur un element 
  15 : bidimensionel(coque, voile) sur la base du modele à  trois couche de
  16 : Marti.
  17 : Les efforts (N11, N22 et N12) et les moments appliqués (M11, M22, M12)
  18 : sont projetés sur deux couches externe et interne d epaisseur
  19 : egale à deux fois l enrobage. Les efforts hors-plan de cisaillment 
  20 : (V1 et V2) sont appliqués à la couche intermediaire:
  21 : 
  22 :  /       ---------------------------------------------------- /
  23 :  | 2c_ext     o        o        o        o        o        o  | cou ext
  24 :  /       .................................................... |
  25 :  |                    n^ /t2                                  |
  26 :  |                     |/                                     |
  27 :  | h-2c_int-2c_ext     |---> t1                               |  h
  28 :  |                                                            |
  29 :  |                                                            |
  30 :  /       .................................................... |
  31 :  | 2c_int     o        o        o        o        o        o  | cou int
  32 :  /       ---------------------------------------------------- /
  33 : 
  34 : Les couches externe et interne sont chargées par une triplet des 
  35 : efforts de membrane.
  36 : 
  37 : N11_ext = N11/2 + M11/D + V1^2/(2V0*cotg(theta))
  38 : N22_ext = N22/2 + M22/D + V2^2/(2V0*cotg(theta))
  39 : N12_ext = N12/2 + M12/D + V1*V2/(2V0*cotg(theta))
  40 : 
  41 : N11_int = N11/2 - M11/D + V1^2/(2V0*cotg(theta))
  42 : N22_int = N22/2 - M22/D + V2^2/(2V0*cotg(theta))
  43 : N12_int = N12/2 - M12/D + V1*V2/(2V0*cotg(theta))
  44 : 
  45 : avec 
  46 : 
  47 : h = epaisseur totale
  48 : c_ext = enrobage couche externe
  49 : c_int = enrobage couche interne
  50 : V0 = (V1^2 + V2^2)^(0.5)
  51 : D = h - c_int - c_ext
  52 : 
  53 : Commentaire :
  54 : -------------
  55 : Les couches externe et interne sont definies selon la direction du 
  56 : vecteur normal n=t1xt2
  57 : 
  58 : L angle theta est compris entre 25 et 45 degrés.
  59 : 
  60 : En entree :
  61 : -----------
  62 : 
  63 :     SIG1:          MCHAML des contraintes
  64 :     MOD1:          MMODEL associé au SIG1
  65 :     MAT1:          MCHAML associé au SIG1
  66 :     VEC1:          direction 1 pour le calcul des contraintes
  67 :     VEC2:          direction 2 pour le calcul des contraintes
  68 :     H1:            epaisseur de l element plaque
  69 :     ENRE1:         enrobage couche externe
  70 :     ENRi1:         enrobage couche interne
  71 :     COT1:          terme cotg(theta)
  72 : 
  73 : 
  74 : En sortie :
  75 : -----------
  76 : 
  77 :     SIG2:          MCHAML avec les composents des efforts projeés 
  78 :                    sur les couches externe et interne et les efforts 
  79 :                    globaux qui agissent sur l element coque
  80 :                    Component N11 external layer    N11E 
  81 :                    Component N22 external layer    N22E 
  82 :                    Component N12 external layer    N12E 
  83 :                    Component N11 internal layer    N11I
  84 :                    Component N22 internal layer    N11I
  85 :                    Component N12 internal layer    N11I
  86 :                    Component M11 global            M11T
  87 :                    Component M22 global            M22T
  88 :                    Component M12 global            M12T
  89 :                    Component V1 global             V1T
  90 :                    Component V2 global             V2T
  91 :                    Resultant Vr                    VR
  92 : 
  93 : Remarques :
  94 : -----------
  95 : 
  96 : Seulement en 3D le caLcul est possible
  97 : 
  98 :