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