1 : $$$$ SIF NOTICE FD218221 18/03/20 21:15:05 9778 2 : DATE 18/03/20 3 : 4 : Procedure SIF Voir aussi : G_THETA 5 : ------------- 6 : SIF MAT1 DEP1 TAB1 ; 7 : 8 : TAB1.'FRTFISS' 9 : .'LEVRE_1' 10 : .'LEVRE_2' 11 : .'MODMIXTE' 12 : .'MAILLAGE' 13 : .'DEBOUCH' 14 : .'PDEBOUCH' 15 : .'EPAI' 16 : .'FLEXION' 17 : .'MEMBRANE' 18 : 19 : 20 : Objet : 21 : _______ 22 : 23 : Cette procedure calcule le facteur d'intensite de contraintes 24 : en mode I (eventuellement en mode II ), a partir des deplacements 25 : des levres de la fissure. 26 : La valeur calculee est une moyenne sur les trois points les 27 : plus proches de la pointe de fissure. 28 : La procedure est applicable aux cas bidimensionnels, tridimen- 29 : sionnels massifs et coques minces. 30 : - En 2D, le probleme d'un chargement en mode mixte peut etre traite. 31 : - En 3D massif , pour chaque noeud du front de fissure, le facteur 32 : d'intensite de contraintes est calcule a partir des deplacements 33 : des points situes dans le plan normal au front de fissure au noeud 34 : considere (le maillage doit etre elabore de maniere a prevoir 35 : l'existence de ces plans normaux au front de fissure). Si la 36 : fissure est debouchante, la formule de contraintes planes est 37 : appliquee pour les points situes en surface. 38 : - En 3D coques minces, les facteurs d'intensite de contraintes de 39 : membrane et de flexion peuvent etre calculees. Le calcul du facteur 40 : d'intensite de contrainte en flexion est realise grace aux 41 : rotations aux noeuds induisant une deformation de peau. Le calcul 42 : en mode mixte est supporte. 43 : 44 : Commentaire : 45 : _____________ 46 : 47 : En entree : 48 : __________ 49 : 50 : MAT1 : Champ de caracteristiques materielles 51 : 52 : DEP1 : Champ de deplacements 53 : 54 : TAB1 : Objet de type TABLE ,indice par des mots, servant a 55 : definir les options et les parametres du calcul : 56 : 57 : Arguments pour un probleme bidimensionnel 58 : _________________________________________ 59 : 60 : indice type objet commentaires 61 : pointe 62 : 63 : FRTFISS POINT pointe de la fissure 64 : 65 : LEVRE_1 MAILLAGE ligne decrivant les levres de la fissure 66 : si chargement en mode mixte, 1ere levre 67 : 68 : MODMIXTE LOGIQUE VRAI si chargement en mode mixte 69 : 70 : LEVRE_2 MAILLAGE si chargement en mode mixte, ligne 71 : decrivant la 2ieme levre de la fissure 72 : 73 : 74 : Arguments pour un probleme tridimensionnel massif 75 : _________________________________________________ 76 : 77 : indice type objet commentaires 78 : pointe 79 : 80 : FRTFISS MAILLAGE ligne decrivant le front de fissure 81 : 82 : LEVRE_1 MAILLAGE surface decrivant les levres de la fissure 83 : si chargement en mode mixte, 1ere levre 84 : 85 : MODMIXTE LOGIQUE VRAI si chargement en mode mixte 86 : 87 : LEVRE_2 MAILLAGE si chargement en mode mixte, surface 88 : decrivant la 2ieme levre de la fissure 89 : 90 : DEBOUCH LOGIQUE VRAI si la fissure est debouchante 91 : 92 : PDEBOUCH POINT/MAILLAGE points du front situes en surface 93 : 94 : 95 : 96 : Arguments pour un probleme tridimensionnel coques minces 97 : ________________________________________________________ 98 : 99 : indice type objet commentaires 100 : pointe 101 : 102 : FRTFISS POINT pointe de la fissure 103 : 104 : LEVRE_1 MAILLAGE ligne decrivant les levres de la fissure 105 : si chargement en mode mixte, 1ere levre 106 : 107 : MODMIXTE LOGIQUE VRAI si chargement en mode mixte 108 : 109 : LEVRE_2 MAILLAGE si chargement en mode mixte, ligne 110 : decrivant la 2ieme levre de la fissure 111 : 112 : EPAI FLOTTANT epaisseur des coques 113 : 114 : MEMBRANE LOGIQUE VRAI pour le calcul du terme de 115 : membrane 116 : 117 : FLEXION LOGIQUE VRAI pour le calcul du terme de 118 : flexion 119 : 120 : 121 : En sortie : 122 : __________ 123 : 124 : En sortie, TAB1 permet de retrouver les valeurs du facteur 125 : d'intensite de contraintes. 126 : 127 : indice type objet commentaires 128 : pointe 129 : 130 : K1 FLOTTANT/TABLE en 2D, flottant: valeur de K1, 131 : en 3D massif, table contenant 132 : les valeurs de K1 a chaque noeud 133 : du front, 134 : en 3D coques minces, table 135 : contenant 3 flottants : 136 : MEMBRANE : terme de membrane 137 : FLEXION : terme de flexion 138 : TOTAL : somme des deux 139 : 140 : K2 FLOTTANT/TABLE si chargement en mode mixte : 141 : en 2D, flottant : valeur de K2 142 : en 3D coques minces, table 143 : contenant 3 flottants : 144 : MEMBRANE : terme de membrane 145 : FLEXION : terme de flexion 146 : TOTAL : somme des deux 147 : 148 : 149 : Exemple : pour lister le facteur K calcule au noeud P15 du 150 : front de fissure ,il faudra coder : LIST (TAB1.K1.P15 ) 151 : 152 : 153 : Remarque : 154 : _________ 155 : 156 : Si on utilise l'option impi 1 ,le trace de la droite V**2=f(r) 157 : ajustee par la methode des moindres carres est superpose aux 158 : valeurs de V**2 aux 3 premiers noeuds de la fissure ;le coef- 159 : ficient R1 qui qualifie l'ajustement de la droite est egale- 160 : ment indique (R1=1 : ajustement parfait). 161 : D'autre part, l'ouverture de la fissure (courbe V=f(r)) est 162 : tracee . 163 : Le signe de K en mode II est donne de façon conventionnelle : 164 : un deplacement positif dans le repere de la fissure de la levre 165 : superieure (donc negatif de la levre inferieure) implique un 166 : K II positif. Si les levres superieures et inferieures sont 167 : respectivement introduite dans la table sous 'LEVRE_1' et 168 : 'LEVRE_2', la convention de signe est respectee. Sinon les levres 169 : seront echangees et le message suivant sera affiche : 170 : ATTENTION! Echange des lignes de fissure! 171 :
© Cast3M 2003 - Tous droits réservés.
Mentions légales