1 : C TUFH1 SOURCE CHAT 05/01/13 03:55:08 5004 2 : SUBROUTINE TUFH1(RAP,XN,AN,XLSLP1,H1,KERR1) 3 : C ----------------------------------------------------------------- 4 : C VALEURS DE H1 (ZAHOOR) POUR THETA/XPI=0.0625D0 R/t=10 POUR 5 : C N=2,5,7,10 EN FONCTION DU RAPPORT LAMDA/(1+LAMDA) 6 : C****************************************************************** 7 : IMPLICIT INTEGER(I-N) 8 : IMPLICIT REAL*8(A-H,O-Z) 9 : DIMENSION XL(21),XN21(21),XN51(21),XN71(21),XN101(21) 10 : DIMENSION XN22(21),XN52(21),XN72(21),XN102(21) 11 : DIMENSION XN23(21),XN53(21),XN73(21),XN103(21) 12 : DIMENSION XN24(21),XN54(21),XN74(21),XN104(21) 13 : DIMENSION XN25(21),XN55(21),XN75(21),XN105(21) 14 : DIMENSION XXN(5,4,21) 15 : DIMENSION TSP(13),TSPI(19) 16 : DIMENSION YN11(13),YN21(13),YN31(13),YN51(13),YN71(13) 17 : DIMENSION YN101(13),YN201(13) 18 : DIMENSION YN12(13),YN22(13),YN32(13),YN52(13),YN72(13) 19 : DIMENSION YN102(13),YN202(13) 20 : DIMENSION YN13(13),YN23(13),YN33(13),YN53(13),YN73(13) 21 : DIMENSION YN103(13),YN203(13) 22 : DIMENSION YN(13,7,3) 23 : DIMENSION ZN11(19),ZN21(19),ZN31(19),ZN51(19),ZN71(19) 24 : DIMENSION ZN12(19),ZN22(19),ZN32(19),ZN52(19),ZN72(19) 25 : DIMENSION ZN13(19),ZN23(19),ZN33(19),ZN53(19),ZN73(19) 26 : DIMENSION ZN(19,5,3) 27 : UN = 1.D0 28 : DEUX=2.D0 29 : TROIS = 3.D0 30 : CINQ=5.D0 31 : SEPT = 7.D0 32 : DIX = 10.D0 33 : VINGT=20.D0 34 : KERR1=0 35 : C 36 : C COMBINED BENDING AND TENSION 37 : C 38 : C LISTE DE LAMDA/(1+LAMDA) 39 : C************************* 40 : DATA XL /0.00D0,0.05D0,0.10D0,0.15D0,0.20D0,0.25D0,0.30D0 41 : # ,0.35D0,0.40D0,0.45D0,0.50D0,0.55D0,0.60D0,0.65D0,0.70D0 42 : # ,0.75D0,0.80D0,0.85D0,0.90D0,0.95D0,1.00D0/ 43 : C VALEURS DE H1 POUR LE CAS 0.0625 44 : C********************************** 45 : 46 : DATA XN21 /3.967D0,4.313D0,4.736D0,5.125D0,5.614D0,6.000D0 47 : # ,6.438D0,6.789D0,7.140D0,7.500D0,7.901D0,8.094D0,8.287D0 48 : # ,8.344D0,8.257D0,8.125D0,7.811D0,7.500D0,7.063D0,6.563D0 49 : # ,6.018D0/ 50 : DATA XN51 /5.567D0,6.500D0,7.375D0,8.250D0,9.080D0,9.750D0 51 : # ,10.501D0,11.000D0,11.457D0,11.875D0,12.150D0,12.313D0 52 : # ,12.236D0,11.938D0,11.642D0,11.125D0,10.617D0,9.875D0 53 : # ,9.190D0,8.500D0,7.620D0/ 54 : DATA XN71 /6.104D0,6.500D0,7.080D0,7.875D0,8.787D0,9.875D0 55 : # ,11.078D0,12.125D0,13.188D0,14.000D0,14.610D0,15.000D0 56 : # ,15.130D0,14.875D0,14.408D0,13.625D0,12.729D0,11.688D0 57 : # ,10.447D0,9.250D0,8.160D0/ 58 : DATA XN101 /6.510D0,7.969D0,9.721D0,11.250D0,12.937D0 59 : # ,14.250D0,15.463D0,16.375D0,17.063D0,17.500D0,17.839D0 60 : # ,17.550D0,17.241D0,16.375D0,15.500D0,14.500D0,13.366D0 61 : # ,12.125D0,10.738D0,9.313D0,7928D0/ 62 : C VALEUR DE H1 POUR LE CAS 0.125 63 : C******************************* 64 : DATA XN22 /4.157D0,4.625D0,5.085D0,5.531D0,5.956D0 65 : # ,6.375D0,6.750D0,7.125D0,7.500D0,7.719D0,7.930D0 66 : # ,8.000D0,8.064D0,8.063D0,7.926D0,7.750D0,7.435D0 67 : # ,7.125D0,6.672D0,6.375D0,5.987D0/ 68 : DATA XN52 /5.163D0,5.750D0,6.331D0,6.875D0,7.323D0 69 : # ,7.875D0,8.309D0,8.688D0,9.008D0,9.313D0,9.500D0 70 : # ,9.625D0,9.566D0,9.406D0,9.165D0,8.813D0,8.390D0 71 : # ,7.938D0,7.391D0,6.875D0,6.311D0/ 72 : DATA XN72 /5.102D0,5.719D0,6.331D0,6.906D0,7.518D0 73 : # ,8.000D0,8.500D0,9.000D0,9.381D0,9.630D0,9.875D0 74 : # ,10.000D0,10.011D0,9.813D0,9.495D0,9.000D0,8.535D0 75 : # ,7.875D0,7.313D0,6.625D0,5.996D0/ 76 : DATA XN102 /4.750D0,5.250D0,5.982D0,6.813D0,7.688D0 77 : # ,8.563D0,9.416D0,10.000D0,10.543D0,10.813D0,10.935D0 78 : # ,10.844D0,10.567D0,10.219D0,9.702D0,9.125D0,8.578D0 79 : # ,7.844D0,7.158D0,6.375D0,5.688D0/ 80 : C VALEURS DE H1 POUR LE CAS 0.25 81 : C******************************* 82 : DATA XN23 /4.159D0,4.522D0,4.885D0,5.250D0,5.614D0 83 : # ,6.000D0,6.347D0,6.650D0,6.900D0,7.150D0,7.363D0 84 : # ,7.500D0,7.508D0,7.450D0,7.307D0,7.100D0,6.885D0 85 : # ,6.519D0,6.152D0,5.700D0,5.312D0/ 86 : DATA XN53 /3.238D0,3.550D0,3.988D0,4.400D0,4.882D0 87 : # ,5.400D0,5.908D0,6.350D0,6.800D0,7.100D0,7.333D0 88 : # ,7.450D0,7.397D0,7.250D0,7.018D0,6.619D0,6.220D0 89 : # ,5.700D0,5.108D0,4.550D0,3.969D0/ 90 : DATA XN73 /2.605D0,3.000D0,3.490D0,4.050D0,4.491D0 91 : # ,5.200D0,5.816D0,6.350D0,6.800D0,7.150D0,7.375D0 92 : # ,7.500D0,7.397D0,7.075D0,6.688D0,6.150D0,5.641D0 93 : # ,5.000D0,4.372D0,3.800D0,3.240D0/ 94 : DATA XN103 /3.000D0,3.225D0,3.589D0,4.000D0,4.540D0 95 : # ,5.050D0,5.700D0,6.400D0,7.264D0,7.707D0,8.150D0 96 : # ,8.150D0,7.900D0,7.500D0,7.018D0,6.300D0,5.670D0 97 : # ,4.950D0,4.353D0,3.700D0,3.125D0/ 98 : C VALEURS DE H1 POUR LE CAS 0.37 99 : C*********** ******************** 100 : DATA XN24 /2.892D0,3.315D0,3.739D0,4.164D0,4.589D0 101 : # ,4.957D0,5.325D0,5.651D0,5.977D0,6.237D0,6.497D0 102 : # ,6.586D0,6.674D0,6.536D0,6.399D0,6.093D0,5.786D0 103 : # ,5.360D0,4.934D0,4.459D0,3.984D0/ 104 : DATA XN54 /1.992D0,2.000D0,2.044D0,2.213D0,2.441D0 105 : # ,2.651D0,2.862D0,3.091D0,3.233D0,3.285D0,3.336D0 106 : # ,3.300D0,3.200D0,3.045D0,2.890D0,2.646D0,2.401D0 107 : # ,2.120D0,1.838D0,1.603D0,1.368D0/ 108 : DATA XN74 /1.496D0,1.700D0,1.894D0,2.094D0,2.294D0 109 : # ,2.484D0,2.675D0,2.791D0,2.906D0,2.945D0,2.985D0 110 : # ,2.883D0,2.781D0,2.629D0,2.477D0,2.251D0,2.025D0 111 : # ,1.767D0,1.509D0,1.301D0,1.092D0/ 112 : DATA XN104 /2.000D0,2.100D0,2.243D0,2.342D0,2.441D0 113 : # ,2.558D0,2.675D0,2.708D0,2.740D0,2.750D0,2.739D0 114 : # ,2.650D0,2.447D0,2.260D0,2.064D0,1.856D0,1.649D0 115 : # ,1.444D0,1.238D0,1.049D0,0.860D0/ 116 : C VALEURS DE H1 POUR LE CAS 0.5 117 : C******************************* 118 : DATA XN25 /2.220D0,2.300D0,2.525D0,2.850D0,3.320D0 119 : # ,3.784D0,4.247D0,4.800D0,5.230D0,5.600D0,5.830D0 120 : # ,5.900D0,5.850D0,5.700D0,5.367D0,5.000D0,4.542D0 121 : # ,4.100D0,3.600D0,3.100D0,2.682D0/ 122 : DATA XN55 /1.137D0,1.300D0,1.496D0,1.700D0,1.953D0 123 : # ,2.200D0,2.475D0,2.700D0,2.900D0,3.100D0,3.160D0 124 : # ,3.200D0,3.170D0,3.050D0,2.890D0,2.700D0,2.430D0 125 : # ,2.183D0,1.935D0,1.700D0,1.424D0/ 126 : DATA XN75 /0.816D0,0.950D0,1.122D0,1.350D0,1.611D0 127 : # ,1.800D0,2.050D0,2.300D0,2.438D0,2.575D0,2.550D0 128 : # ,2.525D0,2.447D0,2.297D0,2.147D0,2.000D0,1.823D0 129 : # ,1.589D0,1.354D0,1.200D0,1.008D0/ 130 : DATA XN105 /1.500D0,1.525D0,1.570D0,1.650D0,1.757D0 131 : # ,1.800D0,1.900D0,2.000D0,2.050D0,2.100D0,2.107D0 132 : # ,2.100D0,2.000D0,1.900D0,1.734D0,1.600D0,1.446D0 133 : # ,1.200D0,1.083D0,0.900D0,0.760D0/ 134 : 135 : 136 : C 137 : C ONLY BENDING MOMENT 138 : C 139 : C************* 140 : C VALEURS DE H1 EN FONCTION DE THETA/XPI POUR R/t=5 ET POUR 141 : C N = 1,2,3,5,7,10,20 142 : C******* 143 : DATA TSP /0.000D0,0.063D0,0.100D0,0.125D0,0.150D0 144 : # ,0.175D0,0.200D0,0.250D0,0.300D0,0.350D0,0.400D0 145 : # ,0.450D0,0.500D0/ 146 : DATA YN11 /0.000D0,0.313D0,0.520D0,0.664D0,0.823D0 147 : # ,0.983D0,1.142D0,1.461D0,1.760D0,2.048D0,2.320D0 148 : # ,2.576D0,2.795D0/ 149 : DATA YN21 /0.000D0,0.369D0,0.560D0,0.708D0,0.836D0 150 : # ,0.960D0,1.080D0,1.333D0,1.526D0,1.720D0,1.872D0 151 : # ,1.990D0,2.059D0/ 152 : DATA YN31 /0.000D0,0.403D0,0.559D0,0.713D0,0.848D0 153 : # ,0.960D0,1.056D0,1.198D0,1.328D0,1.420D0,1.488D0 154 : # ,1.556D0,1.598D0/ 155 : DATA YN51 /0.000D0,0.433D0,0.558D0,0.675D0,0.760D0 156 : # ,0.802D0,0.880D0,0.941D0,0.990D0,1.016D0,1.036D0 157 : # ,1.040D0,1.048D0/ 158 : DATA YN71 /0.000D0,0.432D0,0.557D0,0.611D0,0.638D0 159 : # ,0.664D0,0.691D0,0.744D0,0.754D0,0.760D0,0.774D0 160 : # ,0.780D0,0.794D0/ 161 : DATA YN101 /0.000D0,0.431D0,0.536D0,0.577D0,0.582D0 162 : # ,0.586D0,0.591D0,0.600D0,0.604D0,0.608D0,0.612D0 163 : # ,0.616D0,0.620D0/ 164 : DATA YN201 /0.000D0,0.431D0,0.496D0,0.529D0,0.510D0 165 : # ,0.490D0,0.471D0,0.432D0,0.428D0,0.425D0,0.421D0 166 : # ,0.418D0,0.414D0/ 167 : C************* 168 : C VALEURS DE H1 EN FONCTION DE THETA/XPI POUR R/t=10 ET POUR 169 : C N = 1,2,3,5,7,10,20 170 : C******* 171 : DATA YN12 /0.000D0,0.333D0,0.600D0,0.766D0,0.987D0 172 : # ,1.209D0,1.403D0,1.873D0,2.280D0,2.680D0,3.040D0 173 : # ,3.360D0,3.646D0/ 174 : DATA YN22 /0.000D0,0.401D0,0.760D0,0.855D0,1.040D0 175 : # ,1.202D0,1.400D0,1.771D0,2.003D0,2.260D0,2.440D0 176 : # ,2.572D0,2.682D0/ 177 : DATA YN32 /0.000D0,0.450D0,0.770D0,0.897D0,1.080D0 178 : # ,1.240D0,1.400D0,1.629D0,1.760D0,1.970D0,1.880D0 179 : # ,1.970D0,2.105D0/ 180 : DATA YN52 /0.000D0,0.508D0,0.788D0,0.902D0,1.072D0 181 : # ,1.160D0,1.244D0,1.323D0,1.343D0,1.364D0,1.383D0 182 : # ,1.400D0,1.424D0/ 183 : DATA YN72 /0.000D0,0.531D0,0.750D0,0.857D0,0.930D0 184 : # ,0.988D0,1.020D0,1.080D0,1.100D0,1.100D0,1.080D0 185 : # ,1.052D0,1.035D0/ 186 : DATA YN102 /0.000D0,0.549D0,0.710D0,0.816D0,0.834D0 187 : # ,0.852D0,0.876D0,0.907D0,0.889D0,0.880D0,0.840D0 188 : # ,0.800D0,0.760D0/ 189 : DATA YN202 /0.000D0,0.570D0,0.715D0,0.777D0,0.759D0 190 : # ,0.742D0,0.724D0,0.688D0,0.640D0,0.592D0,0.524D0 191 : # ,0.472D0,0.392D0/ 192 : C************* 193 : C VALEURS DE H1 EN FONCTION DE THETA/XPI POUR R/t=20 ET POUR 194 : C N = 1,2,3,5,7,10,20 195 : C******* 196 : DATA YN13 /0.000D0,0.365D0,0.650D0,0.933D0,1.242D0 197 : # ,1.552D0,1.861D0,2.480D0,2.960D0,3.450D0,3.950D0 198 : # ,4.400D0,4.859D0/ 199 : DATA YN23 /0.000D0,0.456D0,0.850D0,1.091D0,1.350D0 200 : # ,1.600D0,1.863D0,2.391D0,2.700D0,3.020D0,3.260D0 201 : # ,3.445D0,3.571D0/ 202 : DATA YN33 /0.000D0,0.527D0,0.938D0,1.187D0,1.415D0 203 : # ,1.650D0,1.865D0,2.224D0,2.435D0,2.575D0,2.670D0 204 : # ,2.750D0,2.821D0/ 205 : DATA YN53 /0.000D0,0.631D0,0.985D0,1.241D0,1.463D0 206 : # ,1.610D0,1.708D0,1.834D0,1.857D0,1.886D0,1.904D0 207 : # ,1.925D0,1.950D0/ 208 : DATA YN73 /0.000D0,0.696D0,1.025D0,1.208D0,1.271D0 209 : # ,1.334D0,1.397D0,1.524D0,1.539D0,1.560D0,1.570D0 210 : # ,1.595D0,1.600D0/ 211 : DATA YN103 /0.000D0,0.757D0,1.010D0,1.177D0,1.195D0 212 : # ,1.213D0,1.231D0,1.267D0,1.277D0,1.285D0,1.299D0 213 : # ,1.303D0,1.320D0/ 214 : DATA YN203 /0.000D0,0.864D0,1.010D0,1.143D0,1.114D0 215 : # ,1.085D0,1.056D0,1.010D0,1.005D0,1.003D0,1.002D0 216 : # ,1.001D0,1.000D0/ 217 : C 218 : C EFFORT NORMAL 219 : C************* 220 : C VALEURS DE H1 EN FONCTION DE THETA/XPI POUR R/t=5 ET POUR 221 : C N = 1,2,3,5,7 222 : C******* 223 : DATA TSPI /0.000D0,0.063D0,0.100D0,0.125D0,0.150D0,0.175D0 224 : # ,0.200D0,0.225D0,0.250D0,0.275D0,0.300D0,0.325D0,0.350D0 225 : # ,0.375D0,0.400D0,0.425D0,0.450D0,0.475D0,0.500D0/ 226 : DATA ZN11 /0.000D0,0.177D0,0.280D0,0.352D0,0.427D0,0.504D0 227 : # ,0.578D0,0.652D0,0.731D0,0.793D0,0.858D0,0.918D0,0.976D0 228 : # ,1.025D0,1.076D0,1.117D0,1.154D0,1.184D0,1.224D0/ 229 : DATA ZN21 /0.000D0,0.230D0,0.355D0,0.437D0,0.520D0,0.600D0 230 : # ,0.674D0,0.740D0,0.804D0,0.840D0,0.850D0,0.853D0,0.856D0 231 : # ,0.859D0,0.862D0,0.865D0,0.868D0,0.873D0,0.879D0/ 232 : DATA ZN31 /0.000D0,0.265D0,0.400D0,0.479D0,0.560D0,0.630D0 233 : # ,0.695D0,0.742D0,0.765D0,0.786D0,0.795D0,0.793D0,0.786D0 234 : # ,0.778D0,0.760D0,0.743D0,0.720D0,0.696D0,0.669D0/ 235 : DATA ZN51 /0.000D0,0.307D0,0.433D0,0.500D0,0.554D0,0.580D0 236 : # ,0.595D0,0.596D0,0.596D0,0.595D0,0.590D0,0.580D0,0.562D0 237 : # ,0.545D0,0.525D0,0.506D0,0.480D0,0.454D0,0.426D0/ 238 : DATA ZN71 /0.000D0,0.326D0,0.425D0,0.478D0,0.494D0,0.492D0 239 : # ,0.485D0,0.477D0,0.464D0,0.455D0,0.440D0,0.421D0,0.404D0 240 : # ,0.387D0,0.370D0,0.348D0,0.330D0,0.308D0,0.288D0/ 241 : C************* 242 : C VALEURS DE H1 EN FONCTION DE THETA/XPI POUR R/t=10 ET POUR 243 : C N = 1,2,3,5,7 244 : C******* 245 : DATA ZN12 /0.000D0,0.186D0,0.330D0,0.403D0,0.530D0,0.626D0 246 : # ,0.723D0,0.824D0,0.919D0,1.005D0,1.084D0,1.163D0,1.235D0 247 : # ,1.300D0,1.360D0,1.420D0,1.460D0,1.505D0,1.546D0/ 248 : DATA ZN22 /0.000D0,0.248D0,0.415D0,0.520D0,0.620D0,0.720D0 249 : # ,0.822D0,0.930D0,1.040D0,1.073D0,1.094D0,1.100D0,1.105D0 250 : # ,1.107D0,1.110D0,1.110D0,1.110D0,1.110D0,1.110D0/ 251 : DATA ZN32 /0.000D0,0.291D0,0.472D0,0.589D0,0.683D0,0.790D0 252 : # ,0.884D0,0.970D0,1.008D0,1.030D0,1.037D0,1.035D0,1.021D0 253 : # ,1.007D0,0.982D0,0.959D0,0.925D0,0.890D0,0.857D0/ 254 : DATA ZN52 /0.000D0,0.323D0,0.518D0,0.645D0,0.730D0,0.777D0 255 : # ,0.800D0,0.808D0,0.810D0,0.804D0,0.794D0,0.778D0,0.759D0 256 : # ,0.734D0,0.704D0,0.674D0,0.640D0,0.603D0,0.568D0/ 257 : DATA ZN72 /0.000D0,0.382D0,0.595D0,0.638D0,0.660D0,0.666D0 258 : # ,0.670D0,0.661D0,0.652D0,0.635D0,0.615D0,0.597D0,0.573D0 259 : # ,0.550D0,0.524D0,0.496D0,0.465D0,0.439D0,0.408D0/ 260 : C************* 261 : C VALEURS DE H1 EN FONCTION DE THETA/XPI POUR R/t=20 ET POUR 262 : C N = 1,2,3,5,7 263 : C******* 264 : DATA ZN13 /0.000D0,0.204D0,0.368D0,0.489D0,0.630D0,0.770D0 265 : # ,0.920D0,1.065D0,1.202D0,1.340D0,1.448D0,1.548D0,1.640D0 266 : # ,1.728D0,1.800D0,1.868D0,1.930D0,1.970D0,2.020D0/ 267 : DATA ZN23 /0.000D0,0.281D0,0.504D0,0.657D0,0.800D0,0.945D0 268 : # ,1.085D0,1.240D0,1.350D0,1.460D0,1.498D0,1.496D0,1.488D0 269 : # ,1.479D0,1.470D0,1.464D0,1.458D0,1.452D0,1.445D0/ 270 : DATA ZN33 /0.000D0,0.339D0,0.584D0,0.767D0,0.928D0,1.100D0 271 : # ,1.241D0,1.310D0,1.340D0,1.370D0,1.375D0,1.364D0,1.350D0 272 : # ,1.328D0,1.294D0,1.260D0,1.220D0,1.176D0,1.121D0/ 273 : DATA ZN53 /0.000D0,0.426D0,0.684D0,0.873D0,1.000D0,1.050D0 274 : # ,1.080D0,1.099D0,1.094D0,1.088D0,1.075D0,1.049D0,1.019D0 275 : # ,0.984D0,0.945D0,0.904D0,0.855D0,0.804D0,0.758D0/ 276 : DATA ZN73 /0.000D0,0.480D0,0.710D0,0.801D0,0.880D0,0.900D0 277 : # ,0.907D0,0.902D0,0.889D0,0.875D0,0.855D0,0.824D0,0.800D0 278 : # ,0.760D0,0.720D0,0.684D0,0.650D0,0.606D0,0.564D0/ 279 : DO 6 I=1,21 280 : XXN(1,1,I)=XN21(I) 281 : XXN(1,2,I)=XN51(I) 282 : XXN(1,3,I)=XN71(I) 283 : XXN(1,4,I)=XN101(I) 284 : XXN(2,1,I)=XN22(I) 285 : XXN(2,2,I)=XN52(I) 286 : XXN(2,3,I)=XN72(I) 287 : XXN(2,4,I)=XN102(I) 288 : XXN(3,1,I)=XN23(I) 289 : XXN(3,2,I)=XN53(I) 290 : XXN(3,3,I)=XN73(I) 291 : XXN(3,4,I)=XN103(I) 292 : XXN(4,1,I)=XN24(I) 293 : XXN(4,2,I)=XN54(I) 294 : XXN(4,3,I)=XN74(I) 295 : XXN(4,4,I)=XN104(I) 296 : XXN(5,1,I)=XN25(I) 297 : XXN(5,2,I)=XN55(I) 298 : XXN(5,3,I)=XN75(I) 299 : XXN(5,4,I)=XN105(I) 300 : 6 CONTINUE 301 : DO 10 I=1,13 302 : YN(I,1,1) =YN11(I) 303 : YN(I,2,1) =YN21(I) 304 : YN(I,3,1) =YN31(I) 305 : YN(I,4,1) =YN51(I) 306 : YN(I,5,1) =YN71(I) 307 : YN(I,6,1) =YN101(I) 308 : YN(I,7,1) =YN201(I) 309 : 310 : YN(I,1,2) =YN12(I) 311 : YN(I,2,2) =YN22(I) 312 : YN(I,3,2) =YN32(I) 313 : YN(I,4,2) =YN52(I) 314 : YN(I,5,2) =YN72(I) 315 : YN(I,6,2) =YN102(I) 316 : YN(I,7,2) =YN202(I) 317 : 318 : YN(I,1,3) =YN13(I) 319 : YN(I,2,3) =YN23(I) 320 : YN(I,3,3) =YN33(I) 321 : YN(I,4,3) =YN53(I) 322 : YN(I,5,3) =YN73(I) 323 : YN(I,6,3) =YN103(I) 324 : YN(I,7,3) =YN203(I) 325 : 10 CONTINUE 326 : DO 20 I=1,19 327 : ZN(I,1,1) =ZN11(I) 328 : ZN(I,2,1) =ZN21(I) 329 : ZN(I,3,1) =ZN31(I) 330 : ZN(I,4,1) =ZN51(I) 331 : ZN(I,5,1) =ZN71(I) 332 : 333 : ZN(I,1,2) =ZN12(I) 334 : ZN(I,2,2) =ZN22(I) 335 : ZN(I,3,2) =ZN32(I) 336 : ZN(I,4,2) =ZN52(I) 337 : ZN(I,5,2) =ZN72(I) 338 : 339 : ZN(I,1,3) =ZN13(I) 340 : ZN(I,2,3) =ZN23(I) 341 : ZN(I,3,3) =ZN33(I) 342 : ZN(I,4,3) =ZN53(I) 343 : ZN(I,5,3) =ZN73(I) 344 : 20 CONTINUE 345 : 346 : C INTERPOLATION 347 : 348 : IF (XLSLP1.GE.(0.99D0)) THEN 349 : C FLEXION PURE 350 : 351 : IF (XN.GE.UN.AND.XN.LT.DEUX) THEN 352 : XN1=UN 353 : XN2=DEUX 354 : J1=1 355 : J2=2 356 : ENDIF 357 : IF (XN.GE.DEUX.AND.XN.LT.TROIS) THEN 358 : XN1=DEUX 359 : XN2=TROIS 360 : J1=2 361 : J2=3 362 : ENDIF 363 : IF (XN.GE.TROIS.AND.XN.LT.CINQ) THEN 364 : XN1=TROIS 365 : XN2=CINQ 366 : J1=3 367 : J2=4 368 : ENDIF 369 : IF (XN.GE.CINQ.AND.XN.LT.SEPT) THEN 370 : XN1=CINQ 371 : XN2=SEPT 372 : J1=4 373 : J2=5 374 : ENDIF 375 : IF (XN.GE.SEPT.AND.XN.LT.DIX) THEN 376 : XN1=SEPT 377 : XN2=DIX 378 : J1=5 379 : J2=6 380 : ENDIF 381 : IF (XN.GE.DIX.AND.XN.LE.VINGT) THEN 382 : XN1=DIX 383 : XN2=VINGT 384 : J1=6 385 : J2=7 386 : ENDIF 387 : IF (XN.LT.UN.OR.XN.GT.VINGT) THEN 388 : KERR1=683 389 : RETURN 390 : ENDIF 391 : DO 30 I=1,11 392 : IF (AN.GE.TSP(I).AND.AN.LT.TSP(I+1)) THEN 393 : AN1=TSP(I) 394 : AN2=TSP(I+1) 395 : I1=I 396 : I2=I+1 397 : ENDIF 398 : 30 CONTINUE 399 : IF (AN.GE.TSP(11).AND.AN.LE.TSP(12)) THEN 400 : AN1=TSP(11) 401 : AN2=TSP(12) 402 : I1=11 403 : I2=12 404 : ENDIF 405 : IF (AN.GT.TSP(12)) THEN 406 : KERR1=683 407 : RETURN 408 : ENDIF 409 : 410 : IF (RAP.GE.CINQ.AND.RAP.LT.DIX) THEN 411 : RAP1=CINQ 412 : RAP2 =DIX 413 : K1=1 414 : K2=2 415 : ENDIF 416 : IF (RAP.GE.DIX.AND.RAP.LE.VINGT) THEN 417 : RAP1=DIX 418 : RAP2 =VINGT 419 : K1=2 420 : K2=3 421 : ENDIF 422 : IF (RAP.LT.CINQ.OR.RAP.GT.VINGT) THEN 423 : KERR1=683 424 : RETURN 425 : ENDIF 426 : 427 : X11=YN(I1,J1,K1) 428 : X21=YN(I1,J2,K1) 429 : X12=YN(I2,J1,K1) 430 : X22=YN(I2,J2,K1) 431 : CO1=(XN2-XN)/(XN2-XN1) 432 : X1=CO1*X11+(1-CO1)*X21 433 : X2=CO1*X12+(1-CO1)*X22 434 : CO2=(AN2-AN)/(AN2-AN1) 435 : X=CO2*X1+(1-CO2)*X2 436 : Y11=YN(I1,J1,K2) 437 : Y21=YN(I1,J2,K2) 438 : Y12=YN(I2,J1,K2) 439 : Y22=YN(I2,J2,K2) 440 : CO1=(XN2-XN)/(XN2-XN1) 441 : Y1=CO1*Y11+(1-CO1)*Y21 442 : Y2=CO1*Y12+(1-CO1)*Y22 443 : CO2=(AN2-AN)/(AN2-AN1) 444 : Y=CO2*Y1+(1-CO2)*Y2 445 : CO3=(RAP2-RAP)/(RAP2-RAP1) 446 : H1=CO3*X+(1-CO3)*Y 447 : ENDIF 448 : 449 : IF (XLSLP1.LE.(0.01D0)) THEN 450 : C MEMBRANE PURE 451 : 452 : IF (XN.GE.UN.AND.XN.LT.DEUX) THEN 453 : XN1=UN 454 : XN2=DEUX 455 : J1=1 456 : J2=2 457 : ENDIF 458 : IF (XN.GE.DEUX.AND.XN.LT.TROIS) THEN 459 : XN1=DEUX 460 : XN2=TROIS 461 : J1=2 462 : J2=3 463 : ENDIF 464 : IF (XN.GE.TROIS.AND.XN.LT.CINQ) THEN 465 : XN1=TROIS 466 : XN2=CINQ 467 : J1=3 468 : J2=4 469 : ENDIF 470 : IF (XN.GE.CINQ.AND.XN.LE.SEPT) THEN 471 : XN1=CINQ 472 : XN2=SEPT 473 : J1=4 474 : J2=5 475 : ENDIF 476 : IF (XN.LT.UN.OR.XN.GT.SEPT) THEN 477 : KERR1=683 478 : RETURN 479 : ENDIF 480 : 481 : DO 40 I=1,17 482 : IF (AN.GE.TSPI(I).AND.AN.LT.TSPI(I+1)) THEN 483 : AN1=TSPI(I) 484 : AN2=TSPI(I+1) 485 : I1=I 486 : I2=I+1 487 : ENDIF 488 : 40 CONTINUE 489 : IF (AN.GE.TSPI(17).AND.AN.LE.TSPI(18)) THEN 490 : AN1=TSP(11) 491 : AN2=TSP(12) 492 : I1=11 493 : I2=12 494 : ENDIF 495 : IF (AN.GT.TSPI(18)) THEN 496 : KERR1=683 497 : RETURN 498 : ENDIF 499 : 500 : IF (RAP.GE.CINQ.AND.RAP.LT.DIX) THEN 501 : RAP1=CINQ 502 : RAP2 =DIX 503 : K1=1 504 : K2=2 505 : ENDIF 506 : IF (RAP.GE.DIX.AND.RAP.LE.VINGT) THEN 507 : RAP1=DIX 508 : RAP2 =VINGT 509 : K1=2 510 : K2=3 511 : ENDIF 512 : IF (RAP.LT.CINQ.OR.RAP.GT.VINGT) THEN 513 : KERR1=683 514 : RETURN 515 : ENDIF 516 : 517 : X11=ZN(I1,J1,K1) 518 : X21=ZN(I1,J2,K1) 519 : X12=ZN(I2,J1,K1) 520 : X22=ZN(I2,J2,K1) 521 : CO1=(XN2-XN)/(XN2-XN1) 522 : X1=CO1*X11+(1-CO1)*X21 523 : X2=CO1*X12+(1-CO1)*X22 524 : CO2=(AN2-AN)/(AN2-AN1) 525 : X=CO2*X1+(1-CO2)*X2 526 : Y11=ZN(I1,J1,K2) 527 : Y21=ZN(I1,J2,K2) 528 : Y12=ZN(I2,J1,K2) 529 : Y22=ZN(I2,J2,K2) 530 : CO1=(XN2-XN)/(XN2-XN1) 531 : Y1=CO1*Y11+(1-CO1)*Y21 532 : Y2=CO1*Y12+(1-CO1)*Y22 533 : CO2=(AN2-AN)/(AN2-AN1) 534 : Y=CO2*Y1+(1-CO2)*Y2 535 : CO3=(RAP2-RAP)/(RAP2-RAP1) 536 : H1=CO3*X+(1-CO3)*Y 537 : ENDIF 538 : 539 : IF (XLSLP1.GT.(0.01D0).AND.XLSLP1.LT.(0.99D0)) THEN 540 : 541 : IF (XN.GE.DEUX.AND.XN.LT.CINQ) THEN 542 : XN1 =DEUX 543 : XN2 =CINQ 544 : J1=1 545 : J2=2 546 : ENDIF 547 : IF (XN.GE.CINQ.AND.XN.LT.SEPT) THEN 548 : XN1 =CINQ 549 : XN2 =SEPT 550 : J1=2 551 : J2=3 552 : ENDIF 553 : IF (XN.GE.SEPT.AND.XN.LE.DIX) THEN 554 : XN1 =SEPT 555 : XN2 =DIX 556 : J1=3 557 : J2=4 558 : ENDIF 559 : IF (XN.LT.DEUX.OR.XN.GT.DIX) THEN 560 : KERR1=683 561 : RETURN 562 : ENDIF 563 : 564 : 565 : DO 50 I=1,20 566 : IF (XLSLP1.GE.XL(I).AND.XLSLP1.LT.XL(I+1)) THEN 567 : XL1=XL(I) 568 : XL2=XL(I+1) 569 : K1=I 570 : K2=I+1 571 : ENDIF 572 : 50 CONTINUE 573 : 574 : IF (AN.LT.(0.125D0)) THEN 575 : AN1=0.0625D0 576 : AN2 =0.125D0 577 : I1=1 578 : I2=2 579 : ENDIF 580 : IF (AN.GE.(0.125D0).AND.AN.LT.(0.25D0)) THEN 581 : AN1=0.125D0 582 : AN2 =0.25D0 583 : I1=2 584 : I2=3 585 : ENDIF 586 : IF (AN.GE.(0.25D0).AND.AN.LT.(0.37D0)) THEN 587 : AN1=0.25D0 588 : AN2 =0.37D0 589 : I1=3 590 : I2=4 591 : ENDIF 592 : IF (AN.GE.(0.37D0).AND.AN.LE.(0.5D0)) THEN 593 : AN1=0.37D0 594 : AN2 =0.50D0 595 : I1=4 596 : I2=5 597 : ENDIF 598 : IF (AN.GT.(0.5D0)) THEN 599 : KERR1=683 600 : RETURN 601 : ENDIF 602 : 603 : IF (RAP.LT.CINQ.OR.RAP.GT.VINGT) THEN 604 : KERR1=683 605 : RETURN 606 : ENDIF 607 : IF (RAP.LT.(9.D0).OR.RAP.GT.(11.D0)) THEN 608 : KERR1=684 609 : RETURN 610 : ENDIF 611 : 612 : X11=XXN(I1,J1,K1) 613 : X21=XXN(I1,J2,K1) 614 : X12=XXN(I2,J1,K1) 615 : X22=XXN(I2,J2,K1) 616 : CO1=(XN2-XN)/(XN2-XN1) 617 : X1=CO1*X11+(1-CO1)*X21 618 : X2=CO1*X12+(1-CO1)*X22 619 : CO2=(AN2-AN)/(AN2-AN1) 620 : X=CO2*X1+(1-CO2)*X2 621 : Y11=XXN(I1,J1,K2) 622 : Y21=XXN(I1,J2,K2) 623 : Y12=XXN(I2,J1,K2) 624 : Y22=XXN(I2,J2,K2) 625 : CO1=(XN2-XN)/(XN2-XN1) 626 : Y1=CO1*Y11+(1-CO1)*Y21 627 : Y2=CO1*Y12+(1-CO1)*Y22 628 : CO2=(AN2-AN)/(AN2-AN1) 629 : Y=CO2*Y1+(1-CO2)*Y2 630 : CO3=(XL2-XLSLP1)/(XL2-XL1) 631 : H1=CO3*X+(1-CO3)*Y 632 : ENDIF 633 : RETURN 634 : END 635 : 636 : 637 : 638 : 639 :
© Cast3M 2003 - All rights reserved.
Disclaimer