Télécharger la plaquette Cast3M


fr en

Télécharger clmbda.eso

Retour à la liste

Numérotation des lignes : 

clmbda

  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