Télécharger clmbda.eso

Retour à la liste

Numérotation des lignes :

  1. C CLMBDA SOURCE AM 12/12/20 21:15:01 7631
  2. SUBROUTINE CLMBDA(D1,D2,D3,LAMBDA,YOUNG,NU,ALPH,COMP)
  3. IMPLICIT INTEGER(I-N)
  4. implicit real*8 (a-h,o-z)
  5. REAL*8 D1,D2,D3,ALPH
  6. REAL*8 YOUNG,NU
  7. REAL*8 S0302,S0320,S2100
  8. REAL*8 LAMBDA(6)
  9. REAL*8 PI
  10. REAL*8 X1,sinpixs2,sinpiys2,sinpizs2
  11. LOGICAL COMP
  12. ******************************************************************
  13. ******************************************************************
  14. * CALCUL DE LA MATRICE LAMDA DANS LE CAS OU IL N'Y A PAS DE DELTA*
  15. * (PAS DE DECOMPOSITION DE SIGMA-N EN SIGMA-D + SIGMA-V) *
  16. ******************************************************************
  17. ******************************************************************
  18. C
  19. C INITIALISATION DES VARIABLES
  20. C
  21. PI=4.D0*ATAN(1.D0)
  22. IF (COMP) THEN
  23. AF=1.D0-(0.5*(D2+D3))**ALPH
  24. BF=1.D0-(0.5*(D3+D1))**ALPH
  25. CF=1.D0-(0.5*(D1+D2))**ALPH
  26. ELSE
  27. AF=1.D0-D1
  28. BF=1.D0-D2
  29. CF=1.D0-D3
  30. END IF
  31. C print*,'dans callambda: alpha=',alph
  32. C print*,'af,bf,cf',af,bf,cf
  33. C print*,'D3',D3,'d1',d1,'D2',D2
  34. C
  35. C ON A DIVISE EN HAUT ET EN BAS DE X1/DET PAR AF*BF*CF
  36. C -> X1 VAUT 1 (CF FIJ)
  37. C
  38. X1=AF*BF*CF
  39. C
  40. C CALCUL DES INTEGRALES SUR 1/4 DE SPHERE
  41. C
  42. sinpixs2=sin(pi*AF/2.D0)
  43. sinpiys2=sin(pi*BF/2.D0)
  44. sinpizs2=sin(pi*CF/2.D0)
  45. S0302=
  46. & + .49587724947700D+02*sinpixs2 *sinpiys2 *sinpizs2
  47. & -.11826019614403D+03 *sinpixs2 *sinpiys2 *sinpizs2**2
  48. & + .15888740080185D+03*sinpixs2 *sinpiys2 *sinpizs2**3
  49. & -.12090181599425D+03 *sinpixs2 *sinpiys2 *sinpizs2**4
  50. & + .48134780378626D+02*sinpixs2 *sinpiys2 *sinpizs2**5
  51. & -.76097170495877D+01 *sinpixs2 *sinpiys2 *sinpizs2**6
  52. & -.42982186153999D+02 *sinpixs2 *sinpiys2**2 *sinpizs2
  53. & + .92334205518749D+02*sinpixs2 *sinpiys2**2 *sinpizs2**2
  54. & -.94595653059631D+02 *sinpixs2 *sinpiys2**2 *sinpizs2**3
  55. & + .46508202980390D+02*sinpixs2 *sinpiys2**2 *sinpizs2**4
  56. & -.85557270110063D+01 *sinpixs2 *sinpiys2**2 *sinpizs2**5
  57. & -.59795082278790D+01 *sinpixs2 *sinpiys2**3 *sinpizs2
  58. & -.26820923709862D+02 *sinpixs2 *sinpiys2**3 *sinpizs2**2
  59. & + .23967186388666D+02*sinpixs2 *sinpiys2**3 *sinpizs2**3
  60. & -.68917781168937D+01 *sinpixs2 *sinpiys2**3 *sinpizs2**4
  61. & + .49597006251687D+02*sinpixs2 *sinpiys2**4 *sinpizs2
  62. & -.38087148378371D+01 *sinpixs2 *sinpiys2**4 *sinpizs2**2
  63. S0302=S0302
  64. & -.62685158801033D+00 *sinpixs2 *sinpiys2**4 *sinpizs2**3
  65. & -.45758829115125D+02 *sinpixs2 *sinpiys2**5 *sinpizs2
  66. & + .21517303738659D+01*sinpixs2 *sinpiys2**5 *sinpizs2**2
  67. & + .15441307857503D+02*sinpixs2 *sinpiys2**6 *sinpizs2
  68. & -.11825524477172D+03 *sinpixs2**2 *sinpiys2 *sinpizs2
  69. & + .16747889813112D+03*sinpixs2**2 *sinpiys2 *sinpizs2**2
  70. & -.13676189436854D+03 *sinpixs2**2 *sinpiys2 *sinpizs2**3
  71. & + .59638830553943D+02*sinpixs2**2 *sinpiys2 *sinpizs2**4
  72. & -.10868512220306D+02 *sinpixs2**2 *sinpiys2 *sinpizs2**5
  73. & + .92341016512039D+02*sinpixs2**2 *sinpiys2**2 *sinpizs2
  74. & -.82914555064696D+02 *sinpixs2**2 *sinpiys2**2 *sinpizs2**2
  75. & + .41590292781164D+02*sinpixs2**2 *sinpiys2**2 *sinpizs2**3
  76. & -.83827514727064D+01 *sinpixs2**2 *sinpiys2**2 *sinpizs2**4
  77. & -.26813764702232D+02 *sinpixs2**2 *sinpiys2**3 *sinpizs2
  78. & + .14967505013368D+02*sinpixs2**2 *sinpiys2**3 *sinpizs2**2
  79. & -.43826806253910D+01 *sinpixs2**2 *sinpiys2**3 *sinpizs2**3
  80. & -.38031226905467D+01 *sinpixs2**2 *sinpiys2**4 *sinpizs2
  81. & + .10826472143751D+01*sinpixs2**2 *sinpiys2**4 *sinpizs2**2
  82. S0302=S0302
  83. & + .21539434629142D+01*sinpixs2**2 *sinpiys2**5 *sinpizs2
  84. & + .15885885518183D+03*sinpixs2**3 *sinpiys2 *sinpizs2
  85. & -.13675927885166D+03 *sinpixs2**3 *sinpiys2 *sinpizs2**2
  86. & + .62883150455112D+02*sinpixs2**3 *sinpiys2 *sinpizs2**3
  87. & -.11888439526808D+02 *sinpixs2**3 *sinpiys2 *sinpizs2**4
  88. & -.94627040869973D+02 *sinpixs2**3 *sinpiys2**2 *sinpizs2
  89. & + .41591649499595D+02*sinpixs2**3 *sinpiys2**2 *sinpizs2**2
  90. & -.81213735252988D+01 *sinpixs2**3 *sinpiys2**2 *sinpizs2**3
  91. & + .23945386833297D+02*sinpixs2**3 *sinpiys2**3 *sinpizs2
  92. & -.43826106229617D+01 *sinpixs2**3 *sinpiys2**3 *sinpizs2**2
  93. & -.63409891017392D+00 *sinpixs2**3 *sinpiys2**4 *sinpizs2
  94. & -.12083164855111D+03 *sinpixs2**4 *sinpiys2 *sinpizs2
  95. & + .59635568525663D+02*sinpixs2**4 *sinpiys2 *sinpizs2**2
  96. & -.11888504501967D+02 *sinpixs2**4 *sinpiys2 *sinpizs2**3
  97. & + .46560572564192D+02*sinpixs2**4 *sinpiys2**2 *sinpizs2
  98. & -.83833965641336D+01 *sinpixs2**4 *sinpiys2**2 *sinpizs2**2
  99. & -.68748680983310D+01 *sinpixs2**4 *sinpiys2**3 *sinpizs2
  100. & + .48058063063930D+02*sinpixs2**5 *sinpiys2 *sinpizs2
  101. S0302=S0302
  102. & -.10867120963554D+02 *sinpixs2**5 *sinpiys2 *sinpizs2**2
  103. & -.85838947053063D+01 *sinpixs2**5 *sinpiys2**2 *sinpizs2
  104. & -.75792226760522D+01 *sinpixs2**6 *sinpiys2 *sinpizs2
  105. C
  106. C
  107. C
  108. S0320=
  109. & + .49587727033586D+02*sinpixs2 *sinpiys2 *sinpizs2
  110. & -.11826018579995D+03 *sinpixs2 *sinpiys2 *sinpizs2**2
  111. & + .15888737682738D+03*sinpixs2 *sinpiys2 *sinpizs2**3
  112. & -.12090180560884D+03 *sinpixs2 *sinpiys2 *sinpizs2**4
  113. & + .48134782941498D+02*sinpixs2 *sinpiys2 *sinpizs2**5
  114. & -.76097182547511D+01 *sinpixs2 *sinpiys2 *sinpizs2**6
  115. & -.11825524141371D+03 *sinpixs2 *sinpiys2**2 *sinpizs2
  116. & + .16747884821826D+03*sinpixs2 *sinpiys2**2 *sinpizs2**2
  117. & -.13676182945011D+03 *sinpixs2 *sinpiys2**2 *sinpizs2**3
  118. & + .59638796128496D+02*sinpixs2 *sinpiys2**2 *sinpizs2**4
  119. & -.10868505584411D+02 *sinpixs2 *sinpiys2**2 *sinpizs2**5
  120. & + .15885885515815D+03*sinpixs2 *sinpiys2**3 *sinpizs2
  121. & -.13675922155796D+03 *sinpixs2 *sinpiys2**3 *sinpizs2**2
  122. & + .62883118221758D+02*sinpixs2 *sinpiys2**3 *sinpizs2**3
  123. & -.11888433966211D+02 *sinpixs2 *sinpiys2**3 *sinpizs2**4
  124. & -.12083167330411D+03 *sinpixs2 *sinpiys2**4 *sinpizs2
  125. & + .59635549771322D+02*sinpixs2 *sinpiys2**4 *sinpizs2**2
  126. & -.11888500424114D+02 *sinpixs2 *sinpiys2**4 *sinpizs2**3
  127. S0320=S0320
  128. & + .48058086050282D+02*sinpixs2 *sinpiys2**5 *sinpizs2
  129. & -.10867121049526D+02 *sinpixs2 *sinpiys2**5 *sinpizs2**2
  130. & -.75792275647182D+01 *sinpixs2 *sinpiys2**6 *sinpizs2
  131. & -.42982213136941D+02 *sinpixs2**2 *sinpiys2 *sinpizs2
  132. & + .92334241631458D+02*sinpixs2**2 *sinpiys2 *sinpizs2**2
  133. & -.94595651357197D+02 *sinpixs2**2 *sinpiys2 *sinpizs2**3
  134. & + .46508203600717D+02*sinpixs2**2 *sinpiys2 *sinpizs2**4
  135. & -.85557322153729D+01 *sinpixs2**2 *sinpiys2 *sinpizs2**5
  136. & + .92341059860355D+02*sinpixs2**2 *sinpiys2**2 *sinpizs2
  137. & -.82914586469419D+02 *sinpixs2**2 *sinpiys2**2 *sinpizs2**2
  138. & + .41590274527138D+02*sinpixs2**2 *sinpiys2**2 *sinpizs2**3
  139. & -.83827399930577D+01 *sinpixs2**2 *sinpiys2**2 *sinpizs2**4
  140. & -.94627049453367D+02 *sinpixs2**2 *sinpiys2**3 *sinpizs2
  141. & + .41591614676660D+02*sinpixs2**2 *sinpiys2**3 *sinpizs2**2
  142. & -.81213580886453D+01 *sinpixs2**2 *sinpiys2**3 *sinpizs2**3
  143. & + .46560584004622D+02*sinpixs2**2 *sinpiys2**4 *sinpizs2
  144. & -.83833812204707D+01 *sinpixs2**2 *sinpiys2**4 *sinpizs2**2
  145. & -.85839039617849D+01 *sinpixs2**2 *sinpiys2**5 *sinpizs2
  146. S0320=S0320
  147. & -.59794809896161D+01 *sinpixs2**3 *sinpiys2 *sinpizs2
  148. & -.26820979413895D+02 *sinpixs2**3 *sinpiys2 *sinpizs2**2
  149. & + .23967192097949D+02*sinpixs2**3 *sinpiys2 *sinpizs2**3
  150. & -.68917762822240D+01 *sinpixs2**3 *sinpiys2 *sinpizs2**4
  151. & -.26813831452726D+02 *sinpixs2**3 *sinpiys2**2 *sinpizs2
  152. & + .14967600520974D+02*sinpixs2**3 *sinpiys2**2 *sinpizs2**2
  153. & -.43826992963629D+01 *sinpixs2**3 *sinpiys2**2 *sinpizs2**3
  154. & + .23945404982568D+02*sinpixs2**3 *sinpiys2**3 *sinpizs2
  155. & -.43826231635356D+01 *sinpixs2**3 *sinpiys2**3 *sinpizs2**2
  156. & -.68748689102157D+01 *sinpixs2**3 *sinpiys2**4 *sinpizs2
  157. & + .49597019950936D+02*sinpixs2**4 *sinpiys2 *sinpizs2
  158. & -.38087036456114D+01 *sinpixs2**4 *sinpiys2 *sinpizs2**2
  159. & -.62684827232429D+00 *sinpixs2**4 *sinpiys2 *sinpizs2**3
  160. & -.38031088909489D+01 *sinpixs2**4 *sinpiys2**2 *sinpizs2
  161. & + .10826222088093D+01*sinpixs2**4 *sinpiys2**2 *sinpizs2**2
  162. & -.63410176286672D+00 *sinpixs2**4 *sinpiys2**3 *sinpizs2
  163. & -.45758849889839D+02 *sinpixs2**5 *sinpiys2 *sinpizs2
  164. & + .21517325427401D+01*sinpixs2**5 *sinpiys2 *sinpizs2**2
  165. S0320=S0320
  166. & + .21539481683039D+01*sinpixs2**5 *sinpiys2**2 *sinpizs2
  167. & + .15441312720879D+02*sinpixs2**6 *sinpiys2 *sinpizs2
  168. C
  169. C
  170. C
  171. S2100=
  172. & + .49587467973979D+02*sinpixs2 *sinpiys2 *sinpizs2
  173. & -.42980769227454D+02 *sinpixs2 *sinpiys2 *sinpizs2**2
  174. & -.59805404737685D+01 *sinpixs2 *sinpiys2 *sinpizs2**3
  175. & + .49599317855294D+02*sinpixs2 *sinpiys2 *sinpizs2**4
  176. & -.45759239038724D+02 *sinpixs2 *sinpiys2 *sinpizs2**5
  177. & + .15441895510618D+02*sinpixs2 *sinpiys2 *sinpizs2**6
  178. & -.11825825078992D+03 *sinpixs2 *sinpiys2**2 *sinpizs2
  179. & + .92332492764995D+02*sinpixs2 *sinpiys2**2 *sinpizs2**2
  180. & -.26819401505354D+02 *sinpixs2 *sinpiys2**2 *sinpizs2**3
  181. & -.38091597509388D+01 *sinpixs2 *sinpiys2**2 *sinpizs2**4
  182. & + .21520684276939D+01*sinpixs2 *sinpiys2**2 *sinpizs2**5
  183. & + .15888385621676D+03*sinpixs2 *sinpiys2**3 *sinpizs2
  184. & -.94594703123626D+02 *sinpixs2 *sinpiys2**3 *sinpizs2**2
  185. & + .23966025562355D+02*sinpixs2 *sinpiys2**3 *sinpizs2**3
  186. & -.62703651861082D+00 *sinpixs2 *sinpiys2**3 *sinpizs2**4
  187. & -.12089720898763D+03 *sinpixs2 *sinpiys2**4 *sinpizs2
  188. & + .46508836417526D+02*sinpixs2 *sinpiys2**4 *sinpizs2**2
  189. & -.68912765229417D+01 *sinpixs2 *sinpiys2**4 *sinpizs2**3
  190. S2100=S2100
  191. & + .48129967460623D+02*sinpixs2 *sinpiys2**5 *sinpizs2
  192. & -.85566605377921D+01 *sinpixs2 *sinpiys2**5 *sinpizs2**2
  193. & -.76074803611026D+01 *sinpixs2 *sinpiys2**6 *sinpizs2
  194. & -.11825824332720D+03 *sinpixs2**2 *sinpiys2 *sinpizs2
  195. & + .92332489273672D+02*sinpixs2**2 *sinpiys2 *sinpizs2**2
  196. & -.26819394292098D+02 *sinpixs2**2 *sinpiys2 *sinpizs2**3
  197. & -.38091608415685D+01 *sinpixs2**2 *sinpiys2 *sinpizs2**4
  198. & + .21520663179297D+01*sinpixs2**2 *sinpiys2 *sinpizs2**5
  199. & + .16747419796953D+03*sinpixs2**2 *sinpiys2**2 *sinpizs2
  200. & -.82911722127338D+02 *sinpixs2**2 *sinpiys2**2 *sinpizs2**2
  201. & + .14966176247075D+02*sinpixs2**2 *sinpiys2**2 *sinpizs2**3
  202. & + .10825591719184D+01*sinpixs2**2 *sinpiys2**2 *sinpizs2**4
  203. & -.13675642532919D+03 *sinpixs2**2 *sinpiys2**3 *sinpizs2
  204. & + .41589140819666D+02*sinpixs2**2 *sinpiys2**3 *sinpizs2**2
  205. & -.43821050567095D+01 *sinpixs2**2 *sinpiys2**3 *sinpizs2**3
  206. & + .59639580135436D+02*sinpixs2**2 *sinpiys2**4 *sinpizs2
  207. & -.83812695652628D+01 *sinpixs2**2 *sinpiys2**4 *sinpizs2**2
  208. & -.10870415781963D+02 *sinpixs2**2 *sinpiys2**5 *sinpizs2
  209. S2100=S2100
  210. & + .15888383027537D+03*sinpixs2**3 *sinpiys2 *sinpizs2
  211. & -.94594699500019D+02 *sinpixs2**3 *sinpiys2 *sinpizs2**2
  212. & + .23966016553017D+02*sinpixs2**3 *sinpiys2 *sinpizs2**3
  213. & -.62703218218355D+00 *sinpixs2**3 *sinpiys2 *sinpizs2**4
  214. & -.13675640659286D+03 *sinpixs2**3 *sinpiys2**2 *sinpizs2
  215. & + .41589154124408D+02*sinpixs2**3 *sinpiys2**2 *sinpizs2**2
  216. & -.43821101691482D+01 *sinpixs2**3 *sinpiys2**2 *sinpizs2**3
  217. & + .62875012679184D+02*sinpixs2**3 *sinpiys2**3 *sinpizs2
  218. & -.81231387874529D+01 *sinpixs2**3 *sinpiys2**3 *sinpizs2**2
  219. & -.11886187812260D+02 *sinpixs2**3 *sinpiys2**4 *sinpizs2
  220. & -.12089717335754D+03 *sinpixs2**4 *sinpiys2 *sinpizs2
  221. & + .46508832417159D+02*sinpixs2**4 *sinpiys2 *sinpizs2**2
  222. & -.68912742079210D+01 *sinpixs2**4 *sinpiys2 *sinpizs2**3
  223. & + .59639553850800D+02*sinpixs2**4 *sinpiys2**2 *sinpizs2
  224. & -.83812725667904D+01 *sinpixs2**4 *sinpiys2**2 *sinpizs2**2
  225. & -.11886185734126D+02 *sinpixs2**4 *sinpiys2**3 *sinpizs2
  226. & + .48129948086467D+02*sinpixs2**5 *sinpiys2 *sinpizs2
  227. & -.85566592717808D+01 *sinpixs2**5 *sinpiys2 *sinpizs2**2
  228. S2100=S2100
  229. & -.10870405546200D+02 *sinpixs2**5 *sinpiys2**2 *sinpizs2
  230. & -.76074772239863D+01 *sinpixs2**6 *sinpiys2 *sinpizs2
  231. C
  232. C
  233. C
  234. C
  235. C DEFINITION DE LA MATRICE DE RIGIDITE LAMBDA
  236. C ON CALCULE LES INTEGRALES ENTRE 0 ET PI/2 POUR TETA ET 0 ET PI POUR PHI
  237. C -> IL FAUT MULTIPLIER PAR 4 LES INTEGRALES
  238. C
  239. LAMBDA(1)=S0320
  240. LAMBDA(2)=S0302
  241. LAMBDA(3)=S2100
  242. LAMBDA(4)=(S0320+S0302)/2.D0
  243. LAMBDA(5)=(S0320+S2100)/2.D0
  244. LAMBDA(6)=(S0302+S2100)/2.D0
  245. C WRITE (6,200) LAMBDA
  246. C WRITE (6,100) S0320,S0302,S2100,S0300
  247. C WRITE (6,100) S2340,S0322,S2300,S0302,S2100,S0300
  248. C WRITE (6,300) D1,D2,D3
  249. C200 FORMAT ('LAMBDA = ',/,6 ( ' |',6 (E13.6,' , '),' |',/),//)
  250. C100 FORMAT ('SOM = ',/,6 ( ' |',6 (E11.4,' , '),' |',/),//)
  251. C300 FORMAT ('D1=',E11.4,'D2=',E11.4,'D3=',E11.4)
  252. RETURN
  253. END
  254.  
  255.  
  256.  
  257.  
  258.  

© Cast3M 2003 - Tous droits réservés.
Mentions légales