1 : * @FIS_3DS PROCEDUR CHAT 97/09/25 21:15:11 2851 2 : * 3 : ************************************************* 4 : * VERSION 3.0 fournie par EDF/MMM * 5 : * Version du 15/06/96 * 6 : ************************************************* 7 : * 8 : * parametres de la procedure principale : 9 : * --------------------------------------- 10 : * 11 : * c demi grand axe de l'ellipse 12 : * a demi petit axe de l'ellipse 13 : * rc0 rayon du tore 14 : * rc1 coefficient multiplicateur du parametre rc0 15 : * definissant l'epaisseur de la premiere couronne 16 : * de deraffinement 17 : * rc2 coefficient multiplicateur du parametre rc0 18 : * definissant l'epaisseur de la deuxieme couronne 19 : * de deraffinement 20 : * rc3 coefficient multiplicateur du parametre rc0 21 : * definissant l'epaisseur de la troisieme couronne 22 : * de deraffinement (si ndt=2) 23 : * nc nombre de couronnes 24 : * ns nombre de secteurs sur 90 degres 25 : * nt nombre de divisions sur un quart d'ellipse 26 : * eps demi-angle d'ouverture de la fissure (degres) 27 : * lo longueur du bloc 28 : * to largeur du bloc 29 : * ho hauteur du bloc 30 : * ndt nombre de couronnes de deraffinement (1 ou 2) 31 : * nsdt Nombre de secteurs sur 90 degres au niveau des 32 : * couronnes de deraffinement des tranches (2 ou 4) 33 : * beta impose le decoupage le long de la generatrice 34 : * alpha impose l'angle des differentes tranches 35 : * xl impose le nombre d'elements pour la prolongation 36 : * du bloc initial suivant l'axe x (longueur) 37 : * xt impose le nombre d'elements pour la prolongation 38 : * du bloc initial suivant l'axe y (largeur) 39 : * xh impose le nombre d'elements pour la prolongation 40 : * du bloc initial suivant l'axe -z (hauteur) 41 : * d1 'oui' si deplacement des noeuds intermediares du 42 : * boudin 43 : * d2 'oui' si maillage 'rayonnant' sur la face 44 : * exterieure au boudin torique 45 : * On doit avoir d2 = 'non' dans le cas de l'utilisation 46 : * de la procedure plaque epaisse fissuree 47 : debproc @FIS_3DS c*flottant a*flottant lo*flottant 48 : to*flottant ho*flottant 49 : nt*entier nc*entier ns*entier 50 : rc0*flottant rc1*flottant rc2*flottant rc3*flottant 51 : alpha*flottant 52 : ndt*entier nsdt*entier 53 : xl*entier xt*entier xh*entier ; 54 : * 55 : * DEBUT PARAMETRES UTILISATEUR FIXES 56 : * 57 : eps = 0.1 ; 58 : beta = 1. ; 59 : d1 = 'oui' ; 60 : d2 = 'non' ; 61 : * 62 : * FIN PARAMETRES UTILISATEUR FIXES 63 : * 64 : si (c < a) ; 65 : saut ligne ; 66 : mess '================================================== ' ; 67 : mess ' ATTENTION !!! ' ; 68 : saut ligne ; 69 : mess ' Le premier parametre (demi grand axe) doit etre ' ; 70 : mess ' superieur aux deuxieme parametre (demi petit axe)' ; 71 : mess ' Demi grand axe : c =' c ; 72 : mess ' Demi petit axe : a =' a ; 73 : saut ligne ; 74 : mess ' ------ ON SORT DE LA PROCEDURE ----- ' ; 75 : mess '================================================== ' ; 76 : quit @FIS_3DS ; 77 : finsi ; 78 : * 79 : asurc = a / c ; 80 : si ((asurc < .09) ou (asurc > 1.01)) ; 81 : saut ligne ; 82 : mess '================================================== ' ; 83 : mess ' ATTENTION !!! ' ; 84 : saut ligne ; 85 : mess ' L"allongement de l"ellipse : ' ; 86 : mess ' (demi petit axe)/(demi grand axe) doit etre ' ; 87 : mess ' compris entre 0,1 et 1 ' ; 88 : mess ' (a / c) =' asurc ; 89 : saut ligne ; 90 : mess ' ------ ON SORT DE LA PROCEDURE ----- ' ; 91 : mess '================================================== ' ; 92 : quit @FIS_3DS ; 93 : finsi ; 94 : * 95 : asurto = a / to ; 96 : si ((asurto < .0999) ou (asurto > .8001)) ; 97 : saut ligne ; 98 : mess '================================================== ' ; 99 : mess ' ATTENTION !!! ' ; 100 : saut ligne ; 101 : mess ' Le rapport a/to (demi-petit axe / largeur du bloc)' ; 102 : mess ' doit etre compris entre 0,1 et 0,8 ' ; 103 : mess ' (a/to) =' a/to ; 104 : saut ligne ; 105 : mess ' ------ ON SORT DE LA PROCEDURE ----- ' ; 106 : mess '================================================== ' ; 107 : quit @FIS_3DS ; 108 : finsi ; 109 : * 110 : si ((nc < 1) ou (nc > 6)) ; 111 : saut ligne ; 112 : mess '================================================== ' ; 113 : mess ' ATTENTION !!! ' ; 114 : saut ligne ; 115 : mess ' Le nombre de couronnes du "boudin torique" ' ; 116 : mess ' doit etre compris entre 1 et 6 ' ; 117 : mess ' nc =' nc ; 118 : saut ligne ; 119 : mess ' ------ ON SORT DE LA PROCEDURE ----- ' ; 120 : mess '================================================== ' ; 121 : quit @FIS_3DS ; 122 : finsi ; 123 : * 124 : si ((ns < 2) ou (ns > 9)) ; 125 : saut ligne ; 126 : mess '================================================== ' ; 127 : mess ' ATTENTION !!! ' ; 128 : saut ligne ; 129 : mess ' Le nombre de secteurs du "boudin torique" ' ; 130 : mess ' doit etre compris entre 2 et 9 (sur 90 degres) ' ; 131 : mess ' ns =' ns ; 132 : saut ligne ; 133 : mess ' ------ ON SORT DE LA PROCEDURE ----- ' ; 134 : mess '================================================== ' ; 135 : quit @FIS_3DS ; 136 : finsi ; 137 : * 138 : si ((rc1 < .5) ou (rc1 > 3.)) ; 139 : saut ligne ; 140 : mess '================================================== ' ; 141 : mess ' ATTENTION !!! ' ; 142 : saut ligne ; 143 : mess ' Le parametre rc1 doit etre compris entre 0,5 et 3 ' ; 144 : mess ' rc1 =' rc1 ; 145 : saut ligne ; 146 : mess ' ------ ON SORT DE LA PROCEDURE ----- ' ; 147 : mess '================================================== ' ; 148 : quit @FIS_3DS ; 149 : finsi ; 150 : * 151 : si ((rc2 < .5) ou (rc2 > 5.)) ; 152 : saut ligne ; 153 : mess '================================================== ' ; 154 : mess ' ATTENTION !!! ' ; 155 : saut ligne ; 156 : mess ' Le parametre rc2 doit etre compris entre 0,5 et 5 ' ; 157 : mess ' rc2 =' rc2 ; 158 : saut ligne ; 159 : mess ' ------ ON SORT DE LA PROCEDURE ----- ' ; 160 : mess '================================================== ' ; 161 : quit @FIS_3DS ; 162 : finsi ; 163 : * 164 : si ((rc3 < .5) ou (rc3 > 5.)) ; 165 : saut ligne ; 166 : mess '================================================== ' ; 167 : mess ' ATTENTION !!! ' ; 168 : saut ligne ; 169 : mess ' Le parametre rc3 doit etre compris entre 0,5 et 5 ' ; 170 : mess ' rc3 =' rc3 ; 171 : saut ligne ; 172 : mess ' ------ ON SORT DE LA PROCEDURE ----- ' ; 173 : mess '================================================== ' ; 174 : quit @FIS_3DS ; 175 : finsi ; 176 : * 177 : vtot lff levreinf sar slat sinf sav_s ssup_s slat_s boudin ep3 178 : = @FIS_3 c a lo to ho eps nt nc ns rc0 rc1 rc2 rc3 179 : beta alpha ndt nsdt xl xt xh d1 d2 ; 180 : * 181 : finproc vtot lff levreinf sar slat sinf sav_s ssup_s 182 : slat_s boudin ep3 ; 183 :
© Cast3M 2003 - All rights reserved.
Disclaimer