SPACE ARC-SYSTEM ANALYSIS AS LINEAR STATIC 05-08-2007 (c)
function [Pl,Pg,K_l,K_g,T,Hu,D,Re,R] = arc_system(arc_no)
SPACE ARC-SYSTEM ANALYSIS AS LINEAR STATIC 05-08-2007 (c)
H.C.E. Ali ÖZGÜL
Description:
In this matlab application descriptioned space-arc-system
analysis as linear static form. This application had
selected minimal enegry theorem as produced the flexibility
matrix. This sub-function be possible plane-arc-system
analysis.
This function need more input veriables and not
satify for all input veriables in function procedure.
Also, other input veriables in --arc_system-- function. [1]
Syntax:
Input-variables:
arc_no = Graphical interface for selected arc-element no.
[1] Function input variables:
m_p = Material properties matrix
p_p = Per arc element angles theta and beta matrix
Cor = Global system all nodes cartesian coordinates
Pos = Position matrix
Re = Reology matrix
Output-variables:
Pl = Local axis system node reactions
Pg = Globa axis system node reactions
K_l = Local axis system stiffness matrix
K_g = Global axis system stiffness matrix
T = Transformation matrix
S = Static connection matrix for (j) to (i) node
F = Flexibility matrix
Hu = Modal displacement matrix for per element
D = System displacement
Re = Reology matrix
R = Boundary conditions system
Example:
[Pl,Pg,K_l,K_g,T,Hu,D,Re,R] = arc_system([17:19 25:28]) ;
Pl = 12x32 double
Pg = 12x32 double
K_l = 12x12x32 double
K_g = 12x12x32 double
T = 12x12x32 double
Hu = 32x12 double
D = 58x1 double
Re = 17x6 double
References:
{0001 function [Pl,Pg,K_l,K_g,T,Hu,D,Re,R] = arc_system(arc_no) 0002 % SPACE ARC-SYSTEM ANALYSIS AS LINEAR STATIC 05-08-2007 (c) 0003 % H.C.E. Ali ÖZGÜL 0004 % 0005 % Description: 0006 % 0007 % In this matlab application descriptioned space-arc-system 0008 % analysis as linear static form. This application had 0009 % selected minimal enegry theorem as produced the flexibility 0010 % matrix. This sub-function be possible plane-arc-system 0011 % analysis. 0012 % 0013 % This function need more input veriables and not 0014 % satify for all input veriables in function procedure. 0015 % Also, other input veriables in --arc_system-- function. [1] 0016 % Syntax: 0017 % Input-variables: 0018 % 0019 % arc_no = Graphical interface for selected arc-element no. 0020 % 0021 % [1] Function input variables: 0022 % m_p = Material properties matrix 0023 % p_p = Per arc element angles theta and beta matrix 0024 % Cor = Global system all nodes cartesian coordinates 0025 % Pos = Position matrix 0026 % Re = Reology matrix 0027 % 0028 % Output-variables: 0029 % Pl = Local axis system node reactions 0030 % Pg = Globa axis system node reactions 0031 % K_l = Local axis system stiffness matrix 0032 % K_g = Global axis system stiffness matrix 0033 % T = Transformation matrix 0034 % S = Static connection matrix for (j) to (i) node 0035 % F = Flexibility matrix 0036 % Hu = Modal displacement matrix for per element 0037 % D = System displacement 0038 % Re = Reology matrix 0039 % R = Boundary conditions system 0040 % 0041 % Example: 0042 % 0043 % [Pl,Pg,K_l,K_g,T,Hu,D,Re,R] = arc_system([17:19 25:28]) ; 0044 % 0045 % Pl = 12x32 double 0046 % Pg = 12x32 double 0047 % K_l = 12x12x32 double 0048 % K_g = 12x12x32 double 0049 % T = 12x12x32 double 0050 % Hu = 32x12 double 0051 % D = 58x1 double 0052 % Re = 17x6 double 0053 % 0054 % References: 0055 %{ 0056 0057 [1] Robert D. Cook, David S. Malkus, Michael E. Plesha, Concepts and Applications 0058 of Finite Element Analysis, John-Wiley Sons,1989, no:88-27929. 0059 0060 [2] T.Y. Yang, Finite Element Structural Analysis, Prentice-Hall International 0061 Series in Civil Engineering and Engineering Mechanics,1986, ISBN:0-13-317116-7. 0062 0063 [3] Kasimzade A.A., Finite Element Method : Foundation and Application to Earthquake 0064 Engineering (is included education and finite element analysis programs CD) , 0065 Istanbul, Beta Publication , (First edition 1997) Second edition, 2004, p.827. 0066 (ISBN 975-511-379-7). 0067 0068 [4] Kasimzade A.A., Structural Dynamics : Theory and Application to Earthquake Engineering 0069 (is included education and dynamic analysis programs CD) , Istanbul, Beta Publication , 0070 (First edition 1998) Second edition, 2004, p.527. (ISBN 975-511-381-9). 0071 0072 [5] J.N.Reddy, D.K. Gartling, The Finite Element Method in Heat Transfer and 0073 Fluid Dynamics, Second edition,CRC Press,no:DE-AC04-76DP00789. 0074 0075 [6] Chuen-Yuan Chia, Nonlinear Analysis of Plates, McGraww-Hill Press,1980, 0076 ISBN: 0-07-010746-7. 0077 0078 [7] Harry Kraus, Thin Elastic Shells, John Wiley &Sons inc.,1967,Library 0079 of Congress Catalog Card Number: 67-23328. 0080 0081 [8] William Weaver, Paul R. Johnston, Finite Elementsd for Structural Analysis 0082 Prentice-Hall Inc.,1984, ISBN: 0-13-317099-3. 0083 0084 [9] Mei Duan,Yutaka Miyatoma, Shoji Iwasaki,Hideaki Deto, 5-Node hybrid/mixed 0085 finite element for Reissner-Mindlin Plate, Elsevier Science, 1999, 0086 PII:S0168-874X(99)-00020-7. 0087 0088 [10] A. Cazzani, E. Garusi, A. Tralli, S.N. Alturi, A four node hybrid assumed-strain 0089 finite element for laminated composite plates, CMC, vol-2, no:1,pp:23-28,2005 0090 0091 [11] A. Laulusa,O.A. Bauchau ,J-Y. Choi,V.B.C. Tan c, L. Li, Evaluation of some shear 0092 deformable shell elements, Elsevier Science, 18 October 2005,pp:43-(2006)-50335054 0093 0094 [12] emsettin Temiz, Adnan Özel, Murat Demir Aydn, FE Stress Analysis of Thick 0095 Composite Laminates with a Hole in Bending,Composite Materials,Kluwer Academic 0096 Publishers,28 November 2002,pp:10:103117 0097 0098 [13] Kasimzade A.A., Theory of Elasticity and Structural Analysis of space frame , 0099 membrane, plate, shell, asolid, solid systems (is included education and finite 0100 element analysis programs-2 diskettes ) , Istanbul ,"Beta" Publication, 2000, 0101 p.401 . (ISBN 975-486-866-2). 0102 0103 %} 0104 0105 0106 clc 0107 %Error control-1 0108 if nargin==false || (nargin==true && (isempty(arc_no)==true || isa(arc_no,'double')==false)) 0109 disp('warning: This function analysis need arc_node vector') 0110 disp('Runing automatic example:') 0111 arc_no = [1 2]; 0112 end 0113 0114 0115 0116 0117 0118 % ########################### 0119 % INPUT-VERIABLES [1,2,3,4,5] 0120 % ########################### 0121 0122 0123 0124 %if area-section rectangular than ka=3/2 0125 %if area-section circular than ka=4/3 0126 0127 % MATERIAL PROPERTIES (SI) 0128 % ####################### ##### [1] ##### 0129 m_p =[ 2E+6 % E_x (KN/m**2) 0130 2E+6 % E_z (KN/m**2) 0131 2E+6 % G (KN/m**2) 0132 3.2E-3 % I_x (m**4) 0133 7.2E-3 % I_z (m**4) 0134 10.40E-3 % J_p (m**4) 0135 5.00 % R = radius (m) 0136 0.24 % A = cross sectional area (m**2) 0137 3/2 ]; % ka = shear twisting factor 0138 0139 0140 p_p = [ 90.00 % theta angles 0141 0.00 ]; % beta angles 0142 % ####################### 0143 0144 %if all arc-element material properties is equal than. 0145 for i=1:32 0146 p_p(:,i)=p_p(:,1); 0147 m_p(:,i)=m_p(:,1); 0148 end 0149 0150 0151 0152 0153 %GLOBAL COORDINATES 0154 %###################### #### [2] #### 0155 % {x_g} {y_g} {z_g} 0156 Cor = [ 0.00 -sqrt(41) 0.00 %1 0157 sqrt(41) 0.00 0.00 0158 0.00 sqrt(41) 0.00 0159 -sqrt(41) 0.00 0.00 %4 0160 0.00 -5.00 0.00 0161 5.00 0.00 0.00 0162 0.00 5.00 0.00 0163 -5.00 0.00 0.00 %8 0164 0.00 -5.00 4.00 0165 5.00 0.00 0.00 0166 0.00 5.00 0.00 0167 -5.00 0.00 0.00 %12 0168 0.00 -3.00 4.00 0169 3.00 0.00 4.00 0170 0.00 3.00 4.00 0171 -3.00 0.00 4.00 0172 0.00 0.00 sqrt(41)]; %17 0173 %###################### 0174 0175 Node=size(Cor,1); 0176 0177 0178 %POSITION MATRIX 0179 %###################### ##### [3] ##### 0180 Pos = [1 2 %1 0181 2 3 0182 1 4 0183 4 3 0184 5 6 %5 0185 6 7 0186 5 8 0187 8 7 0188 9 10 0189 10 11 %10 0190 9 12 0191 12 11 0192 13 14 0193 14 15 0194 13 16 %15 0195 16 15 0196 1 9 0197 2 10 0198 3 11 0199 4 12 %20 0200 9 17 0201 10 17 0202 11 17 0203 12 17 0204 5 13 %25 0205 6 14 0206 7 15 0207 8 16 0208 13 17 0209 14 17 %30 0210 15 17 0211 16 17]; 0212 %###################### 0213 %Total arc-element value 0214 No=size(Pos,1); 0215 0216 0217 0218 for i=1:Node; 0219 Re(i,:)=[1 1 1 1 1 1]; 0220 end 0221 %One node's total degree of freedom. 0222 Nom=size(Re,2); 0223 0224 0225 0226 0227 %SYSTEM SUPPORT 0228 %###################### ##### [4] ##### 0229 %Re(Node number,:)=[ u(x) v(y) w(z) Qxx Qyy Qzz] 0230 Re(1,:)=[0 0 0 0 0 0]; 0231 Re(2,:)=[1 0 0 0 0 0]; 0232 Re(3,:)=[0 0 0 0 0 0]; 0233 Re(4,:)=[1 0 0 0 0 0]; 0234 Re(5,:)=[0 0 0 0 0 0]; 0235 Re(6,:)=[1 0 0 0 0 0]; 0236 Re(7,:)=[0 0 0 0 0 0]; 0237 Re(8,:)=[1 0 0 0 0 0]; 0238 %###################### 0239 0240 0241 0242 %ACCUMULATE METHOD 0243 freedom=0; 0244 for i=1:Node; 0245 for j=1:Nom; 0246 if Re(i,j)==1 ; 0247 freedom = freedom +1; 0248 Re(i,j) = freedom; 0249 end 0250 end 0251 end 0252 0253 0254 0255 0256 %TOPOLOGY METHOD 0257 for i=1:No 0258 R(i,:)=[Re(Pos(i,1),:) Re(Pos(i,2),:) ]; 0259 end 0260 0261 0262 0263 % SYSTEM LOAD 0264 %#################### ##### [5] ##### 0265 %Re(Node number, Freedoom number)[fx fy fz Mxx Myy Mzz] 0266 P(freedom)=0; 0267 P(Re(9:12,3)) = -5.00; 0268 P(Re(17,3)) = 5.00; 0269 P(Re(2,1)) = -5.00; 0270 P(Re(4,1)) = 5.00; 0271 P(Re(6,1)) = -5.00; 0272 P(Re(8,1)) = 5.00; 0273 %#################### 0274 0275 0276 0277 % SYSTEM STIFFNESS MATRIX COMPONENTS 0278 % Calculating; 0279 % F = Flexibility matrix 0280 % S = Static equilibrium matrix 0281 % K_l = Local axis system stiffness matrix 0282 % K_g = Global axis system stiffness matrix 0283 % T = Transformation matrix for global to local axis system 0284 0285 for s=1:No; 0286 %Flexibility matrix 0287 F = flexibility( m_p(1,s) , m_p(2,s) , m_p(3,s) , ... 0288 m_p(4,s) , m_p(5,s) , m_p(6,s) , ... 0289 m_p(7,s) , m_p(8,s) , m_p(9,s) , p_p(1,s) ); %>>> 5 <<< 0290 0291 0292 %Static connection matrix 0293 S = equilibrium(p_p(1,s) , m_p(7,s) ); %>>> 6 <<< 0294 0295 0296 %Local axis system stiffness matrix 0297 K_l(:,:,s) = stiffness_local(F,S); %>>> 7 <<< 0298 0299 0300 %Transformation matrix for local to global axis system 0301 T(:,:,s) = transformation(s,Pos,Cor,p_p(1,s),p_p(2,s)); %>>> 1 <<< 0302 0303 0304 %Global axis system stiffness matrix 0305 K_g(:,:,s)=T(:,:,s)'*K_l(:,:,s)*T(:,:,s); 0306 0307 end %|---s=1:No-->| 0308 0309 0310 0311 0312 0313 %System stiffness matrix with topology superposition method. 0314 Ksis(freedom,freedom)=0; 0315 for n=1:No; 0316 for sat=1:12; 0317 for sut=1:12; 0318 if (R(n,sat)~=0) 0319 if (R(n,sut)~=0); 0320 Ksis(R(n,sut),R(n,sat))=Ksis(R(n,sut),R(n,sat)) + K_g(sat,sut,n); 0321 end 0322 end 0323 end 0324 end 0325 end 0326 0327 0328 %Error control-2 0329 %System stiffness matrix singularity control 0330 if size(Ksis,2) == rank(Ksis) 0331 D = Ksis^-1*P'; 0332 else 0333 disp('warning: System stiffness matrix solution error') 0334 disp('control the system boundary conditions') 0335 disp(R); 0336 end 0337 0338 0339 0340 %Moving per arc-element node freedom global system 0341 for v = 1 : No; 0342 for m = 1 :12; 0343 u = R(v, m); 0344 if u ~=0 0345 Hu(v, m,:) = D(u) ; 0346 else 0347 Hu(v,m)=0; 0348 end 0349 end 0350 end 0351 0352 0353 0354 for v=1 :No; 0355 %Global axis system arc-element node reactions 0356 Pg(:,v) = K_g(:,:,v)*Hu(v,:)';%+Q(:,v); 0357 %Local axis system arc-element node reactions. 0358 Pl(:,v) = T(:,:,v)*Pg(:,v); 0359 end 0360 0361 0362 0363 %Plotting arc-element node reactions. 0364 graphical_interface(Pl,m_p,p_p,arc_no); % >>> 8 <<< 0365 0366 0367 0368 0369 0370 0371 0372 0373 0374 %% SUB-FUNCTIONS 0375 0376 function T = transformation(s, Pos ,Cor,teta,beta )%______________||| 1 ||| 0377 % Description: 0378 % 0379 % This function produced base-transformation matrix for 0380 % arc-element direction cosinesses ny==1 or ny<>1 stattions. 0381 % 0382 % Syntax: 0383 % s = Total arc-element value. 0384 % Pos = Position matrix 0385 % Cor = System all node cartesian coordinates matrix. 0386 % teta = Arc-angle 0387 % beta = Arc-surface position angel (right-hand-rule) 0388 % 0389 % function calls: 0390 % | 0391 % |------ny==1|-->arc_transfrom_special(teta,beta,ly,my,ny,Q) 0392 % | | 0393 % | |-->arc_t_matrix(Ti,Tj) 0394 % | 0395 % |-----ny<>0-|-->arc_transfrom (+teta,beta,ly,my,ny,Q); 0396 % | 0397 % |-->arc_t_matrix(Ti,Tj) 0398 % 0399 0400 0401 0402 %[K_global]XYZ= [T]'.[K_local].[T] 0403 0404 %Moving global arc-element (i) and (j) node position global system 0405 px = [Cor(Pos(s,1),1) Cor(Pos(s,2),1)]; 0406 py = [Cor(Pos(s,1),2) Cor(Pos(s,2),2)]; 0407 pz = [Cor(Pos(s,1),3) Cor(Pos(s,2),3)]; 0408 0409 0410 %Cartesian coordinates 0411 x1 = px(1); x2 = px(2); 0412 y1 = py(1); y2 = py(2); 0413 z1 = pz(1); z2 = pz(2); 0414 0415 0416 %Direction cosinesses 0417 L = sqrt( (x2-x1)^2 + (y2-y1)^2 + (z2-z1)^2 ); 0418 0419 ly = (x2-x1)/L ; 0420 my = (y2-y1)/L ; 0421 ny = (z2-z1)/L ; 0422 Q = sqrt(1-ny^2); 0423 0424 0425 %Special status for ny==1 or ny<>1 than 0426 if ny==1 0427 %Special transform 0428 Ti = arc_transfrom_special(+teta,beta,ly,my,ny); %>>> 2 <<< 0429 Tj = arc_transfrom_special(-teta,beta,ly,my,ny); %>>> 2 <<< 0430 T = arc_t_matrix(Ti,Tj); 0431 else 0432 %Normal transform 0433 Ti = arc_transfrom (+teta,beta,ly,my,ny,Q); %>>> 2 <<< 0434 Tj = arc_transfrom (-teta,beta,ly,my,ny,Q); %>>> 2 <<< 0435 T = arc_t_matrix(Ti,Tj) ; 0436 end 0437 %___T = transformation(s, Pos ,Cor,teta,beta )%_____| 0438 0439 0440 0441 0442 function [T] = arc_transfrom (teta,beta,ly,my,ny,Q) %_____________||| 2 ||| 0443 % Description: 0444 % 0445 % In this sub-function produced arc-element transformation 0446 % matrix for local to global axis system. 0447 % 0448 % Syntax: 0449 % teta: Arc-angle 0450 % beta: Arc-surface rotation angel (right-hand-rule) 0451 % ly,my,ny,Q: Drection cosinesses 0452 % 0453 0454 %radii<----degree 0455 teta=teta*pi/180; 0456 beta=beta*pi/180; 0457 0458 0459 %1. tranformation matrix (Multi-axis rotation matrix with theta,beta) 0460 t1=[cos(beta) -sin(beta)*sin(teta/2) -sin(beta)*cos(teta/2) 0461 sin(beta)*sin(teta/2) (1-2*sin(teta/2)^2*sin(beta/2)^2) -1/2*(1-cos(beta))*sin(teta) 0462 sin(beta)*cos(teta/2) -1/2*(1-cos(beta))*sin(teta) (1-2*cos(teta/2)^2*sin(beta/2)^2) ]; 0463 0464 0465 %2. tranformation matrix (Rotation matrix with theta) 0466 t2 = [ 1 0 0 0467 0 cos(teta/2) sin(teta/2) 0468 0 -sin(teta/2) cos(teta/2) ]; 0469 0470 0471 %3. tranformation matrix (General space vector transform form) 0472 t3 = [ my/Q -ly/Q 0 0473 ly my ny 0474 -ly*ny/Q -my*ny/Q Q ]; 0475 0476 0477 T = t1*t2*t3; 0478 0479 0480 0481 % Transformation matrix components; [Maxima-Output] 0482 %{ 0483 %for arc-element's (i) node 0484 T_11 = my / Q * cos(beta) + ly * ny / Q * sin(beta); 0485 T_12 = my / Q * sin(teta / 2) * sin(beta) + ly * (cos(teta / 2) - ny / Q * sin(teta / 2) * cos(beta)) ; 0486 T_13 = my / Q * cos(teta / 2) * sin(beta) - ly * (sin(teta / 2) + ny / Q * cos(teta / 2) * cos(beta)) ; 0487 T_21 =-ly / Q * cos(beta) + my * ny / Q * sin(beta) ; 0488 T_22 =-ly / Q * sin(teta / 2) * sin(beta) + my * (cos(teta / 2) - ny / Q * sin(teta / 2) * cos(beta)) ; 0489 T_23 =-ly / Q * cos(teta / 2) * sin(beta) - my * (sin(teta / 2) + ny / Q * cos(teta / 2) * cos(beta)) ; 0490 T_31 =-Q * sin(beta) ; 0491 T_32 = ny * cos(teta / 2) + Q * sin(teta / 2) * cos(beta) ; 0492 T_33 = Q * cos(teta / 2) * cos(beta) - ny * sin(teta / 2) ; 0493 0494 %sub transformation matrix for [i] node 0495 Ti = [ T_11 T_12 T_13 0496 T_21 T_22 T_23 0497 T_31 T_32 T_33]' 0498 0499 0500 %for arc-element's (j) node 0501 T_77 = my / Q * cos(beta) + ly * ny / Q * sin(beta); 0502 T_78 = my / Q * sin(-teta / 2) * sin(beta) + ly * (cos(-teta / 2) - ny / Q * sin(-teta / 2) * cos(beta)); 0503 T_79 = my / Q * cos(-teta / 2) * sin(beta) - ly * (sin(-teta / 2) + ny / Q * cos(-teta / 2) * cos(beta)); 0504 T_87 =-ly / Q * cos(beta) + my * ny / Q * sin(beta); 0505 T_88 =-ly / Q * sin(-teta / 2) * sin(beta) + my * (cos(-teta / 2) - ny / Q * sin(-teta / 2) * cos(beta)); 0506 T_89 =-ly / Q * cos(-teta / 2) * sin(beta) - my * (sin(-teta / 2) + ny / Q * cos(-teta / 2) * cos(beta)); 0507 T_97 =-Q * sin(beta); 0508 T_98 = ny * cos(-teta / 2) + Q * sin(-teta / 2) * cos(beta); 0509 T_99 = Q * cos(-teta / 2) * cos(beta) - ny * sin(-teta / 2); 0510 0511 %sub transformation matrix for [j] node 0512 Tj = [ T_77 T_78 T_79 0513 T_87 T_88 T_89 0514 T_97 T_98 T_99]' 0515 %} 0516 %___[T] = arc_transfrom (teta,beta,ly,my,ny,Q)____| 0517 0518 0519 0520 0521 function T_matrix = arc_t_matrix(Ti,Tj) %_________________________||| 3 ||| 0522 % Description: 0523 % 0524 % In sub-function produced transformation matrix general 0525 % form with (i) and (j) node transformation matrix. 0526 % Sytnax: 0527 % Ti = Arc-element (i) node transformation matrix 0528 % Tj = Arc-element (j) node transformation matrix 0529 % 0530 % Example: 0531 % {P}global = [C].{P}lokal 0532 % {P}lokal = [T].{P}global 0533 0534 0535 0536 zn = zeros(3,3); 0537 %Tranformation matrix 0538 T_matrix = [Ti zn zn zn 0539 zn Ti zn zn 0540 zn zn Tj zn 0541 zn zn zn Tj]; 0542 %____T_matrix = arc_t_matrix(Ti,Tj)____| 0543 0544 0545 0546 0547 function T = arc_transfrom_special(teta,beta,ly,my,ny) %________||| 4 ||| 0548 % Description: 0549 % 0550 % Where is calculate special tranformation matrix for 0551 % "z" axis direction cosines of ny = 1 and Q=sqrt(1-ny^2) = 0. 0552 % 0553 % Syntax: 0554 % teta = Arc-angle 0555 % beta = Arc-surface rotation angel (right-hand-rule) 0556 % ly,my,ny,Q = Drection cosinesses 0557 % 0558 0559 %radii<-----degree 0560 teta=teta*pi/180; 0561 beta=beta*pi/180; 0562 0563 0564 t_11 = abs(sin(beta)) + ny*cos(beta); 0565 t_21 = abs(sin(beta)) * sin(teta/2) + ly * cos(teta/2) - ny * cos(beta) * sin(teta/2); 0566 t_31 = abs(sin(beta)) * cos(teta/2) - ly * sin(teta/2) - ny * cos(beta) * cos(teta/2); 0567 t_12 = -cos(beta) + ny * abs(sin(beta)); 0568 t_22 = -cos(beta) * sin(teta/2) + my * cos(teta/2) - ny * abs(sin(beta)) * sin(teta/2); 0569 t_32 = -cos(beta) * cos(teta/2) - my * sin(teta/2) - ny * abs(sin(beta)) * cos(teta/2); 0570 t_13 = 0; 0571 t_23 = ny * cos(teta/2); 0572 t_33 =-ny * sin(teta/2); 0573 0574 0575 T =[t_11 t_12 t_13 0576 t_21 t_22 t_23 0577 t_31 t_32 t_33]; 0578 %____T = arc_transfrom_special(teta,beta,ly,my,ny,Q) __| 0579 0580 0581 0582 0583 0584 0585 0586 function F = flexibility(E_x,E_z,G,I_x,I_z,J_p,R,A,ka,teta)%______||| 5 ||| 0587 % Description: 0588 % 0589 % In this sub-function calculated element flexibility matrix. 0590 % 0591 % 0592 % Syntax: 0593 % E_x = x axis elasticity module (KN/m**2) 0594 % E_z = z axis elasticity module (KN/m**2) 0595 % G = y axis shear elasticity modules (KN/m**2) 0596 % I_x = x axis inertia module 0597 % I_z = z axis inertia module 0598 % J_p = Cross sectional area polar inertia module. 0599 % R = Arc-radius 0600 % A = Arc cross sectional area 0601 % ka = corss sectional area-shear correction factor for 0602 % twisting motion 0603 % teta = arc angle 0604 % 0605 0606 0607 %radian <----degree 0608 Q = teta*pi/180; 0609 0610 0611 %Flexibility matrix constants. (Chain-rule) 0612 a = 0.5*(Q - 0.5*sin(2*Q)) ; 0613 b = Q - sin(Q) ; 0614 c = 0.5*(Q + 0.5*sin(2*Q)) ; 0615 d = 1 - cos(Q) ; 0616 e = 0.5*((sin(Q))^2) ; 0617 0618 0619 %flexibilite matrix components 0620 %f_12 = 0;%f_13 = 0;%f_14 = 0; 0621 f_11 = (R^3)*a / (E_z * I_z) + (R^3) * (2*b-a) / (G*J_p) + ka*R*Q/(G*A); 0622 f_15 = (R^2)*a / (E_z * I_z) - (R^2) * (sin(Q) - c) / (G*J_p); 0623 f_16 = (R^2)*e / (E_z * I_z) + (R^2) * (d - e) / (G*J_p) ; 0624 0625 0626 %f_21 = 0;%f_25 = 0; %f_26 = 0 0627 f_22 = (R^3) * (2*b - a) / (E_x * I_x) + R*c / (E_x * A) + ka*R*a / (G*A); 0628 f_23 = (R^3) * (d - e) / (E_x * I_x) - R*e / (E_x * A) + ka*e / (G*A) ; 0629 f_24 =-(R^2) * b / (E_x * I_x) ; 0630 0631 0632 %f_31 = 0;%f_35 = 0;%f_36 = 0 0633 f_32 = f_23 ; 0634 f_33 = (R^3) * a / (E_x * I_x) + R*a / (E_x * A) + ka*R*c / (G*A); 0635 f_34 =-(R^2) * d / (E_x * I_x) ; 0636 0637 0638 %f_41 = 0; %f_45 = 0;%f_46 = 0 0639 f_42 = f_24; 0640 f_43 = f_34; 0641 f_44 = R*Q / (E_x*I_x); 0642 0643 0644 %f_52 = 0;%f_53 = 0;%f_54 = 0 0645 f_51 = f_15; 0646 f_55 = R*a / (E_z * I_z) + R*c / (G*J_p); 0647 f_56 = R*e / (E_z * I_z) - R*e / (G*J_p); 0648 0649 0650 %f_62 = 0;%f_63 = 0;%f_64 = 0 0651 f_61 = f_16; 0652 f_65 = f_56; 0653 f_66 = R*c / (E_z * I_z) + R*a / (G*J_p); 0654 0655 0656 %flexibilite matrix 0657 F= [ f_11 0 0 0 f_15 f_16 0658 0 f_22 f_23 f_24 0 0 0659 0 f_32 f_33 f_34 0 0 0660 0 f_42 f_43 f_44 0 0 0661 f_51 0 0 0 f_55 f_56 0662 f_61 0 0 0 f_65 f_66]; 0663 0664 0665 %______F = flexibility(E_x,E_z,G,I_x,I_z,J_p,R,A,ka,teta)___| 0666 0667 0668 0669 0670 0671 function S = equilibrium(teta,R) %________________________________||| 6 ||| 0672 % Description: 0673 % 0674 % In this sub-function produced static-connection matrix 0675 % for arc element (j) node. This matrix used element 0676 % local axis system stiffness matrix, graphical plotting 0677 % in arc-element (i) and (j) node reactions. 0678 % 0679 % Syntax: 0680 % teta: Arc-element angle. 0681 % R : Arc-element radius. 0682 % 0683 0684 0685 0686 % radian <----degree 0687 Q = teta*pi/180; 0688 0689 %Static connection matrix for arc element's (i) to (j) node 0690 %{x_j } [-1.00 0 0 0 0 0 ] {x_i } 0691 %{y_j } [ 0 -cos(Q) sin(Q) 0 0 0 ] {y_i } 0692 %{z_j }=[ 0 -sin(Q) -cos(Q) 0 0 0 ]*{z_i } 0693 %{Mx_j} [ 0 R*(1-cos(Q)) R*sin(Q) -1 0 0 ] {Mx_i} 0694 %{My_j} [ R*(1-cos(Q)) 0 0 0 -cos(Q) sin(Q) ] {My_i} 0695 %{Mz_j} [ -R*sin(Q) 0 0 0 -sin(Q) -cos(Q) ] {Mz_i} 0696 0697 %Statical connection matrix for arc element's (i) to (j) node 0698 S = [-1.00 0 0 0 0 0 0699 0 -cos(Q) sin(Q) 0 0 0 0700 0 -sin(Q) -cos(Q) 0 0 0 0701 0 R*(1-cos(Q)) R*sin(Q) -1 0 0 0702 R*(1-cos(Q)) 0 0 0 -cos(Q) sin(Q) 0703 -R*sin(Q) 0 0 0 -sin(Q) -cos(Q) ]; 0704 0705 %_____________________________S = equilibrium(teta,R)__| 0706 0707 0708 0709 0710 0711 function K_local= stiffness_local (F,S) %_________________________||| 7 ||| 0712 % Description: 0713 % 0714 % In this sub-function produced arc-element local axis system 0715 % stiffness matrix with arc-element flexibility matrix. 0716 % 0717 % Syntax: 0718 % F = Flexibility matrix 0719 % S = Static connection matrix for arc element's (i) to (j) node 0720 0721 0722 0723 %Local axis system element stiffness matrix with Flexibility matrix 0724 K_local = [ F^-1 (F^-1)*S' ; 0725 S*(F^-1) S*(F^-1)*S']; 0726 0727 % Local axis system stiffness matrix components [Maxima-output] 0728 %{ 0729 a1 = (f(1,1) * f(5,5) - f(1,5)^2) 0730 a2 = (f(1,1) * f(5,6) - f(1,6)^2) 0731 a3 = (f(1,1) * f(5,6) - f(1,5)*f(1,6))^2 0732 a4 = (f(2,2) * f(3,3) - f(2,3)^2) 0733 a5 = (f(2,2) * f(4,4) - f(2,4)^2) 0734 a6 = (f(2,2) * f(3,4) - f(2,3)*f(2,4))^2 0735 0736 W = ((a1 * a2) - a3) / f(1, 1) 0737 U = ((a4 * a5) - a6) / f(2, 2) 0738 0739 K(1,1) = (f(5, 5) * f(6, 6) - f(5, 6) ^ 2) / W 0740 K(1,2) = 0 0741 K(1,3) = 0 0742 K(1,4) = 0 0743 K(1,5) = (f(6, 1) * f(5, 6) - f(5, 1) * f(6, 6)) / W 0744 K(1,6) = (f(5, 1) * f(5, 6) - f(5, 5) * f(1, 6)) / W 0745 0746 K(1,7) = -K(1, 1) 0747 K(1,8) = 0 0748 K(1,9) = 0 0749 K(1,10) = 0 0750 K(1,11) = -K(1, 5) 0751 K(1,12) = K(1, 6) 0752 0753 K(2,1) = 0 0754 K(2,2) = (f(3, 3) * f(4, 4) - f(3, 4) ^ 2) / U 0755 K(2,3) = (f(4, 2) * f(4, 3) - f(3, 2) * f(4, 4)) / U 0756 K(2,4) = (f(3, 2) * f(4, 3) - f(3, 3) * f(4, 2)) / U 0757 K(2,5) = 0 0758 K(2,6) = 0 0759 K(2,7) = 0 0760 K(2,8) = -K(2, 2) * cos(teta) + K(2, 3) * sin(teta) 0761 0762 K(3,3) = (f(2, 2) * f(4, 4) - f(4, 2) ^ 2) / U 0763 K(3,4) = (f(3, 2) * f(4, 2) - f(2, 2) * f(4, 3)) / U 0764 K(2,9) =-(-K(2, 3) * cos(teta) + K(3, 3) * sin(teta)) 0765 K(2,10) = -K(2, 4) * cos(teta) + K(3, 4) * sin(teta) 0766 K(2,11) = 0 0767 K(2,12) = 0 0768 0769 K(3,1) = 0 0770 K(3,2) = K(2, 3) 0771 K(3,3) = (f(2, 2) * f(4, 4) - f(4, 2) ^ 2) / U 0772 K(3,4) = (f(3, 2) * f(4, 2) - f(2, 2) * f(4, 3)) / U 0773 K(3,5) = 0 0774 K(3,6) = 0 0775 K(3,7) = 0 0776 K(3,8) =-K(2, 9) 0777 K(3,9) = -K(2, 3) * sin(teta) - K(3, 3) * cos(teta) 0778 K(3,10) = K(2, 4) * sin(teta) + K(3, 4) * cos(teta) 0779 K(3,11) = 0 0780 K(3,12) = 0 0781 0782 K(4,1) = 0 0783 K(4,2) = K(2, 4) 0784 K(4,3) = K(3, 4) 0785 K(4,4) = (f(2, 2) * f(3, 3) - f(2, 3) ^ 2) / U 0786 K(4,5) = 0 0787 K(4,6) = 0 0788 K(4,7) = 0 0789 K(4,8) = K(2, 10) 0790 K(4,9) = -K(3, 10) 0791 K(4,10) = K(2, 4) * R * (1 - cos(teta)) + K(3, 4) * R * sin(teta) - K(4, 4) 0792 K(4,11) = 0 0793 K(4,12) = 0 0794 0795 K(5,1) = K(1, 5) 0796 K(5,2) = 0 0797 K(5,3) = 0 0798 K(5,4) = 0 0799 K(5,5) = (f(1, 1) * f(6, 6) - f(6, 1) ^ 2) / W 0800 K(5,6) = (f(5, 1) * f(6, 1) - f(1, 1) * f(6, 5)) / W 0801 K(5,7) = -K(1, 5) 0802 K(5,8) = 0 0803 K(5,9) = 0 0804 K(5,10) = 0 0805 K(5,11) = K(1, 5) * R * (1 - cos(teta)) - K(5, 5) * cos(teta) + K(5, 6) * sin(teta) 0806 0807 K(6,6) = (f(1,1) * f(5, 5) - f(5, 1) ^ 2) / W 0808 K(5,12) = -((K(1, 6) * R * (1 - cos(teta)) - K(5, 6) * cos(teta) + K(6, 6) * sin(teta))) 0809 0810 K(6,1) = K(1, 6) 0811 K(6,2) = 0 0812 K(6,3) = 0 0813 K(6,4) = 0 0814 K(6,5) = K(5, 6) 0815 K(6,6) = (f(1, 1) * f(5, 5) - f(5, 1) ^ 2) / W 0816 K(6,7) = -K(1, 6) 0817 K(6,8) = 0 0818 K(6,9) = 0 0819 K(6,10) = 0 0820 K(6,11) = -K(5, 12) 0821 K(6,12) = -K(1, 6) * R * sin(teta) - K(5, 6) * sin(teta) - K(6, 6) * cos(teta) 0822 0823 K(7,1) = -K(1, 1) 0824 K(7,2) = 0 0825 K(7,3) = 0 0826 K(7,4) = 0 0827 K(7,5) = -K(1, 5) 0828 K(7,6) = -K(1, 6) 0829 K(7,7) = K(1, 1) 0830 K(7,8) = 0 0831 K(7,9) = 0 0832 K(7,10) = 0 0833 K(7,11) = K(1, 5) 0834 K(7,12) = -K(1, 6) 0835 0836 0837 K(8,1) = 0 0838 K(8,2) = -K(2, 8) 0839 K(8,3) = K(3, 8) 0840 K(8,4) = K(4, 8) 0841 K(8,5) = 0 0842 K(8,6) = 0 0843 K(8,7) = 0 0844 K(8,8) = K(2, 2) 0845 K(8,9) = -K(2, 3) 0846 K(8,10) = K(2, 4) 0847 K(8,11) = 0 0848 K(8,12) = 0 0849 0850 K(9,1) = 0 0851 K(9,2) = -K(3, 8) 0852 K(9,3) = K(3, 9) 0853 K(9,4) = K(4, 9) 0854 K(9,5) = 0 0855 K(9,6) = 0 0856 K(9,7) = 0 0857 K(9,8) = -K(2, 3) 0858 K(9,9) = K(3, 3) 0859 K(9,10) = -K(3, 4) 0860 K(9,11) = 0 0861 K(9,12) = 0 0862 0863 K(10,1) = 0 0864 K(10,2) = K(4, 8) 0865 K(10,3) = -K(4, 9) 0866 K(10,4) = K(4, 10) 0867 K(10,5) = 0 0868 K(10,6) = 0 0869 K(10,7) = 0 0870 K(10,8) = K(2, 4) 0871 K(10,9) = -K(3, 4) 0872 K(10,10) = K(4, 4) 0873 K(10,11) = 0 0874 K(10,12) = 0 0875 0876 K(11,1) = -K(1, 5) 0877 K(11,2) = 0 0878 K(11,3) = 0 0879 K(11,4) = 0 0880 K(11,5) = K(5, 11) 0881 K(11,6) = K(6, 11) 0882 K(11,7) = K(1, 5) 0883 K(11,8) = 0 0884 K(11,9) = 0 0885 K(11,10) = 0 0886 K(11,11) = K(5, 5) 0887 K(11,12) = -K(5, 6) 0888 0889 K(12,1) = K(1, 6) 0890 K(12,2) = 0 0891 K(12,3) = 0 0892 K(12,4) = 0 0893 K(12,5) = -K(6, 11) 0894 K(12,6) = K(6, 12) 0895 K(12,7) = -K(1, 6) 0896 K(12,8) = 0 0897 K(12,9) = 0 0898 K(12,10) = 0 0899 K(12,11) = -K(5, 6) 0900 K(12,12) = K(6, 6) 0901 0902 K_local=K; 0903 %} 0904 %_______________________________| 0905 0906 0907 0908 0909 0910 function graphical_interface(Pl,m_p,p_p,arc_no) %_________________||| 8 ||| 0911 % Description: 0912 % 0913 % In this function produced graphical-interface for per unit 0914 % arc-element node reactions as normal,shear force,bending 0915 % and twisting moments components. 0916 % Where per unit arc-element (j) node reactions was find 0917 % (i)-node as known reactions with static connection matrix[S]. 0918 % 0919 % Syntax: 0920 % Pg = Global axis system arc-element node reactions. 0921 % Pl = Local axis systen arc-element node reactions. 0922 % m_p = Material properties matrix. 0923 % p_p = Per arc element angle of theta and beta matrix 0924 % arc_no = Active graphical plotting arc-element reactions. 0925 % 0926 % Example: 0927 % 0928 %Static connection matrix for arc element's (i) to (j) node 0929 %{x_j } [-1.00 0 0 0 0 0 ] {x_i } 0930 %{y_j } [ 0 -cos(Q) sin(Q) 0 0 0 ] {y_i } 0931 %{z_j }=[ 0 -sin(Q) -cos(Q) 0 0 0 ]*{z_i } 0932 %{Mx_j} [ 0 R*(1-cos(Q)) R*sin(Q) -1 0 0 ] {Mx_i} 0933 %{My_j} [ R*(1-cos(Q)) 0 0 0 -cos(Q) sin(Q) ] {My_i} 0934 %{Mz_j} [ -R*sin(Q) 0 0 0 -sin(Q) -cos(Q) ] {Mz_i} 0935 0936 0937 0938 0939 for i=1:length(arc_no) 0940 %Active graphical ploting 0941 s = arc_no(i); 0942 %Local axis system (i) node reactions 0943 x_i = Pl(1,s) ; %fx(i) 0944 y_i = Pl(2,s) ; %fy(i) 0945 z_i = Pl(3,s) ; %fz(i) 0946 Mx_i = Pl(4,s) ; %Mx-x(i) 0947 My_i = Pl(5,s) ; %My-y(i) 0948 Mz_i = Pl(6,s) ; %Mz-z(i) 0949 0950 %Produce theta angles-vector 0951 teta = 0 :p_p(1,s)/50 : p_p(1,s); 0952 %radian<-----degree 0953 teta = teta*pi/180; 0954 %Arc-radius 0955 R = m_p(7,s); 0956 0957 %Calculate the (j) node reaction with matlab-vector applications. 0958 x_j = -x_i.*(teta>=0) ;%fx(j) 0959 y_j = -y_i.*cos(teta) + z_i.*sin(teta) ;%fy(j) 0960 z_j = -(y_i.*sin(teta) + z_i.*cos(teta)) ;%fz(j) 0961 Mx_j = -Mx_i + R*y_i.*(1-cos(teta)) + R*z_i.*sin(teta) ;%Mx-x(j) 0962 My_j = R*x_i.*(1-cos(teta)) - My_i.*cos(teta) + Mz_i.*sin(teta) ;%My-y(j) 0963 Mz_j = -R*x_i.*sin(teta) - My_i.*sin(teta) - Mz_i.*cos(teta) ;%Mz-z(j) 0964 0965 0966 %Open new figure 0967 figure(s); 0968 0969 subplot(6,1,1); 0970 plot(teta,x_j,'LineWidth',8); 0971 graphical_text(teta,x_j) 0972 title ('Nx: x local axis direction normal reactions', ... 0973 'FontSize',8,'horizontalAlignment','right','FontWeight','Bold') 0974 0975 subplot(6,1,2); 0976 plot(teta,y_j,'LineWidth',8); 0977 graphical_text(teta,y_j) 0978 title ('Ny: y local axis direction axial reactions', ... 0979 'FontSize',8,'horizontalAlignment','right','FontWeight','Bold' ) 0980 0981 subplot(6,1,3); 0982 plot(teta,z_j,'LineWidth',8); 0983 graphical_text(teta,z_j) 0984 title ('Nz: z local axis direction shear reactions',... 0985 'FontSize',8,'horizontalAlignment','right','FontWeight','Bold' ) 0986 0987 subplot(6,1,4); 0988 plot(teta,Mx_j,'LineWidth',8); 0989 graphical_text(teta,Mx_j) 0990 title ('Mx-x: x local axis bending monent reactions', ... 0991 'FontSize',8,'horizontalAlignment','right','FontWeight','Bold') 0992 0993 subplot(6,1,5); 0994 plot(teta,My_j,'LineWidth',8); 0995 graphical_text(teta,My_j) 0996 title ('My-y: y local axis bending monent reactions', ... 0997 'FontSize',8,'horizontalAlignment','right','FontWeight','Bold') 0998 0999 subplot(6,1,6); 1000 plot(teta,Mz_j,'LineWidth',8); 1001 graphical_text(teta,Mz_j) 1002 title ('Mz-z: z local axis twisting monent reactions', ... 1003 'FontSize',8,'horizontalAlignment','right','FontWeight','Bold') 1004 1005 end % |--- i=1:length(arc_no)--->>| 1006 %_______________________________graphical_interface(Pg,Pl,m_p,p_p,arc_no)__| 1007 1008 1009 1010 1011 1012 function graphical_text(x,y) %______________________________________||||||| 1013 % Description: 1014 % 1015 % In this sub-function produced graphical plotting for per 1016 % arc-element (i) and (j) node's local reactions. 1017 % 1018 % Sytnax: 1019 % x = plot data as x-vector 1020 % y = plot data as y-vector 1021 % 1022 %Graphical text-plotting. 1023 for i = 1:3:size(x,2) 1024 text(x(1,i),y(1,i),num2str(y(1,i)),... 1025 'FontSize',6,'FontWeight','Bold','Color',[1 1 1]) 1026 1027 end 1028 1029 %Graphical stem structure 1030 line([0 max(x)],[0 0]) ; 1031 1032 for i=1:size(x,2) 1033 line([x(1,i) x(1,i)],[0 y(1,i)]) ; 1034 end 1035 1036 %__________________________graphical_text(x,y)___|
Pl =
Columns 1 through 13
0.3336 -0.3336 -0.3336 0.3336 0.2421 -0.2421 -0.2421 0.2421 0.4254 -0.4270 -0.4254 0.4270 -0.0059
2.1505 2.1505 2.1505 2.1505 1.5603 1.5603 1.5603 1.5603 1.2701 1.4419 1.2701 1.4419 -0.5281
-2.1505 -2.1505 -2.1505 -2.1505 -1.5603 -1.5603 -1.5603 -1.5603 -0.7000 -1.4523 -0.7000 -1.4523 0.3758
-2.9380 -2.9380 -2.9380 -2.9380 -2.1317 -2.1317 -2.1317 -2.1317 -0.3739 -2.6208 -0.3739 -2.6208 0.1283
-0.6855 0.6855 0.6855 -0.6855 -0.4974 0.4974 0.4974 -0.4974 -0.8336 1.1902 0.8336 -1.1902 0.1661
-0.9826 0.9826 0.9826 -0.9826 -0.7129 0.7129 0.7129 -0.7129 -1.2590 0.9329 1.2590 -0.9329 -0.0823
-0.3336 0.3336 0.3336 -0.3336 -0.2421 0.2421 0.2421 -0.2421 -0.4254 0.4270 0.4254 -0.4270 0.0059
-2.1505 -2.1505 -2.1505 -2.1505 -1.5603 -1.5603 -1.5603 -1.5603 -0.7000 -1.4523 -0.7000 -1.4523 0.3758
-2.1505 -2.1505 -2.1505 -2.1505 -1.5603 -1.5603 -1.5603 -1.5603 -1.2701 -1.4419 -1.2701 -1.4419 0.5281
2.9380 2.9380 2.9380 2.9380 2.1317 2.1317 2.1317 2.1317 3.2240 2.5687 3.2240 2.5687 -0.8900
0.6855 -0.6855 -0.6855 0.6855 0.4974 -0.4974 -0.4974 0.4974 0.8682 -1.2022 -0.8682 1.2022 -0.1117
-0.9826 0.9826 0.9826 -0.9826 -0.7129 0.7129 0.7129 -0.7129 -1.2935 0.9448 1.2935 -0.9448 -0.1368
Columns 14 through 26
-0.0126 0.0059 0.0126 0.0000 0.2008 -0.0000 -0.2008 0.0000 -0.2475 0.0000 0.2475 -0.0000 -0.2254
-0.8458 -0.5281 -0.8458 2.8193 0.6126 0.2738 0.6126 -0.4362 -4.3091 -3.9731 -4.3091 0.5876 1.9350
0.7966 0.3758 0.7966 -2.9945 0.2913 0.7052 0.2913 -0.2419 3.3719 3.1039 3.3719 -0.4599 -2.2652
0.7752 0.1283 0.7752 -3.9170 2.3709 3.1029 2.3709 -1.8519 2.5731 2.0977 2.5731 -0.2717 -3.5760
0.1177 -0.1661 -0.1177 -0.0000 -0.7942 0.0000 0.7942 -0.0000 0.4993 -0.0000 -0.4993 0.0000 0.5564
-0.1088 0.0823 0.1088 0.0000 -0.1980 0.0000 0.1980 -0.0000 0.7858 -0.0000 -0.7858 0.0000 0.8678
0.0126 -0.0059 -0.0126 -0.0000 -0.2008 0.0000 0.2008 -0.0000 0.2475 -0.0000 -0.2475 0.0000 0.2254
0.7966 0.3758 0.7966 -2.9945 0.2913 0.7052 0.2913 -0.2419 3.3719 3.1039 3.3719 -0.4599 -2.2652
0.8458 0.5281 0.8458 -2.8193 -0.6126 -0.2738 -0.6126 0.4362 4.3091 3.9731 4.3091 -0.5876 -1.9350
-1.0210 -0.8900 -1.0210 3.0409 2.1483 1.7924 2.1483 -1.5385 -7.2591 -6.4437 -7.2591 0.9099 1.9249
-0.1721 0.1117 0.1721 0.0000 0.8062 -0.0000 -0.8062 0.0000 -0.4515 0.0000 0.4515 0.0000 -0.2595
-0.0544 0.1368 0.0544 -0.0000 -0.2100 0.0000 0.2100 -0.0000 0.7381 -0.0000 -0.7381 -0.0000 0.5708
Columns 27 through 32
0.0000 0.2254 0.0000 0.1569 -0.0000 -0.1569
2.5833 1.9350 1.0386 2.7721 3.6949 2.7721
-2.9524 -2.2652 -1.0383 -2.5162 -3.2583 -2.5162
-4.5970 -3.5760 -1.3398 -3.0914 -3.9686 -3.0914
-0.0000 -0.5564 -0.0000 -0.4334 0.0000 0.4334
-0.0000 -0.8678 0.0000 -0.1565 -0.0000 0.1565
-0.0000 -0.2254 -0.0000 -0.1569 0.0000 0.1569
-2.9524 -2.2652 -1.0383 -2.5162 -3.2583 -2.5162
-2.5833 -1.9350 -1.0386 -2.7721 -3.6949 -2.7721
2.7511 1.9249 1.3412 4.3710 6.1516 4.3710
0.0000 0.2595 0.0000 0.6281 -0.0000 -0.6281
-0.0000 -0.5708 -0.0000 -0.3512 -0.0000 0.3512
Pg =
Columns 1 through 13
2.3864 -2.3864 -2.3864 2.3864 1.7315 -1.7315 -1.7315 1.7315 1.2985 -1.7491 -1.2985 1.7491 -0.4561
1.9146 1.9146 1.9146 1.9146 1.3891 1.3891 1.3891 1.3891 0.6969 1.1452 0.6969 1.1452 -0.4478
-0.0000 -0.0000 -0.0000 -0.0000 -0.0000 -0.0000 -0.0000 -0.0000 -0.3351 -0.0074 -0.3351 -0.0074 -0.1077
-1.9289 -1.9289 -1.9289 -1.9289 -1.3995 -1.3995 -1.3995 -1.3995 -0.5945 -1.9818 -0.5945 -1.9818 0.2150
2.2260 -2.2260 -2.2260 2.2260 1.6151 -1.6151 -1.6151 1.6151 -0.0657 -1.7245 0.0657 1.7245 0.0335
-1.1795 1.1795 1.1795 -1.1795 -0.8558 0.8558 0.8558 -0.8558 -1.4360 1.5013 1.4360 -1.5013 0.0592
-2.3864 2.3864 2.3864 -2.3864 -1.7315 1.7315 1.7315 -1.7315 -1.2985 1.7491 1.2985 -1.7491 0.4561
-1.9146 -1.9146 -1.9146 -1.9146 -1.3891 -1.3891 -1.3891 -1.3891 -0.6969 -1.1452 -0.6969 -1.1452 0.4478
-0.0000 -0.0000 -0.0000 -0.0000 -0.0000 -0.0000 -0.0000 -0.0000 0.3351 0.0074 0.3351 0.0074 0.1077
1.9289 1.9289 1.9289 1.9289 1.3995 1.3995 1.3995 1.3995 1.5625 1.9451 1.5625 1.9451 -0.7536
-2.2260 2.2260 2.2260 -2.2260 -1.6151 1.6151 1.6151 -1.6151 -2.9970 1.6877 2.9970 -1.6877 0.5051
-1.1795 1.1795 1.1795 -1.1795 -0.8558 0.8558 0.8558 -0.8558 -1.1823 1.5181 1.1823 -1.5181 -0.0177
Columns 14 through 26
0.8123 0.4561 -0.8123 0.0000 -0.2272 0.0000 0.2272 0.0000 2.8204 -0.0000 -2.8204 -0.0000 -1.5371
-0.8301 -0.4478 -0.8301 1.4777 0.2008 0.3050 0.2008 0.0839 -0.2475 2.5954 -0.2475 0.2505 -0.2254
-0.0348 -0.1077 -0.0348 3.8382 0.6391 0.6923 0.6391 -0.4917 -4.6886 -4.3224 -4.6886 0.7029 2.5520
0.4349 0.2150 0.4349 -3.9170 0.4216 -3.1029 0.4216 -1.8519 0.8409 -2.0977 0.8409 -0.2717 0.9992
0.6614 -0.0335 -0.6614 -0.0000 2.3709 -0.0000 -2.3709 0.0000 2.5731 -0.0000 -2.5731 -0.0000 -3.5760
0.0062 -0.0592 -0.0062 -0.0000 -0.7016 0.0000 0.7016 -0.0000 0.3996 -0.0000 -0.3996 0.0000 0.2535
-0.8123 -0.4561 0.8123 -0.0000 0.2272 -0.0000 -0.2272 -0.0000 -2.8204 0.0000 2.8204 0.0000 1.5371
0.8301 0.4478 0.8301 -1.4777 -0.2008 -0.3050 -0.2008 -0.0839 0.2475 -2.5954 0.2475 -0.2505 0.2254
0.0348 0.1077 0.0348 -3.8382 -0.6391 -0.6923 -0.6391 0.4917 4.6886 4.3224 4.6886 -0.7029 -2.5520
-0.6087 -0.7536 -0.6087 3.0409 -0.4216 -1.7924 -0.4216 -1.5385 0.5382 6.4437 0.5382 0.9099 0.4267
-0.8352 -0.5051 0.8352 0.0000 2.1483 0.0000 -2.1483 0.0000 -7.2591 -0.0000 7.2591 0.0000 1.9249
0.0832 0.0177 -0.0832 -0.0000 -0.7185 0.0000 0.7185 -0.0000 0.6774 -0.0000 -0.6774 -0.0000 0.4595
Columns 27 through 32
-0.0000 1.5371 0.0000 -2.8054 0.0000 2.8054
-1.9840 -0.2254 1.1461 0.1569 -3.6443 0.1569
3.3843 2.5520 0.9183 2.4791 3.3148 2.4791
4.5970 0.9992 -1.3398 -0.1080 3.9686 -0.1080
-0.0000 3.5760 -0.0000 -3.0914 -0.0000 3.0914
-0.0000 -0.2535 -0.0000 -0.4479 0.0000 0.4479
0.0000 -1.5371 -0.0000 2.8054 -0.0000 -2.8054
1.9840 0.2254 -1.1461 -0.1569 3.6443 -0.1569
-3.3843 -2.5520 -0.9183 -2.4791 -3.3148 -2.4791
-2.7511 0.4267 1.3412 -0.5857 -6.1516 -0.5857
-0.0000 -1.9249 0.0000 4.3710 0.0000 -4.3710
-0.0000 -0.4595 -0.0000 -0.4181 0.0000 0.4181
K_l
K_l(:,:,1) =
1.0e+003 *
0.4194 0 0 0 -0.8617 -1.2352 -0.4194 0 0 0 0.8617 -1.2352
0 2.0009 -1.8221 -2.1644 0 0 0 -1.8221 -2.0009 3.0587 0 0
0 -1.8221 2.0009 3.0587 0 0 0 2.0009 1.8221 -2.1644 0 0
0 -2.1644 3.0587 6.6185 0 0 0 3.0587 2.1644 -2.1470 0 0
-0.8617 0 0 0 3.9669 2.2838 0.8617 0 0 0 -2.0248 0.3417
-1.2352 0 0 0 2.2838 5.8342 1.2352 0 0 0 -0.3417 3.8922
-0.4194 0 0 0 0.8617 1.2352 0.4194 0 0 0 -0.8617 1.2352
0 -1.8221 2.0009 3.0587 0 0 0 2.0009 1.8221 -2.1644 0 0
0 -2.0009 1.8221 2.1644 0 0 0 1.8221 2.0009 -3.0587 0 0
0 3.0587 -2.1644 -2.1470 0 0 0 -2.1644 -3.0587 6.6185 0 0
0.8617 0 0 0 -2.0248 -0.3417 -0.8617 0 0 0 3.9669 -2.2838
-1.2352 0 0 0 0.3417 3.8922 1.2352 0 0 0 -2.2838 5.8342
K_l(:,:,2) =
1.0e+003 *
0.4194 0 0 0 -0.8617 -1.2352 -0.4194 0 0 0 0.8617 -1.2352
0 2.0009 -1.8221 -2.1644 0 0 0 -1.8221 -2.0009 3.0587 0 0
0 -1.8221 2.0009 3.0587 0 0 0 2.0009 1.8221 -2.1644 0 0
0 -2.1644 3.0587 6.6185 0 0 0 3.0587 2.1644 -2.1470 0 0
-0.8617 0 0 0 3.9669 2.2838 0.8617 0 0 0 -2.0248 0.3417
-1.2352 0 0 0 2.2838 5.8342 1.2352 0 0 0 -0.3417 3.8922
-0.4194 0 0 0 0.8617 1.2352 0.4194 0 0 0 -0.8617 1.2352
0 -1.8221 2.0009 3.0587 0 0 0 2.0009 1.8221 -2.1644 0 0
0 -2.0009 1.8221 2.1644 0 0 0 1.8221 2.0009 -3.0587 0 0
0 3.0587 -2.1644 -2.1470 0 0 0 -2.1644 -3.0587 6.6185 0 0
0.8617 0 0 0 -2.0248 -0.3417 -0.8617 0 0 0 3.9669 -2.2838
-1.2352 0 0 0 0.3417 3.8922 1.2352 0 0 0 -2.2838 5.8342
K_l(:,:,3) =
1.0e+003 *
0.4194 0 0 0 -0.8617 -1.2352 -0.4194 0 0 0 0.8617 -1.2352
0 2.0009 -1.8221 -2.1644 0 0 0 -1.8221 -2.0009 3.0587 0 0
0 -1.8221 2.0009 3.0587 0 0 0 2.0009 1.8221 -2.1644 0 0
0 -2.1644 3.0587 6.6185 0 0 0 3.0587 2.1644 -2.1470 0 0
-0.8617 0 0 0 3.9669 2.2838 0.8617 0 0 0 -2.0248 0.3417
-1.2352 0 0 0 2.2838 5.8342 1.2352 0 0 0 -0.3417 3.8922
-0.4194 0 0 0 0.8617 1.2352 0.4194 0 0 0 -0.8617 1.2352
0 -1.8221 2.0009 3.0587 0 0 0 2.0009 1.8221 -2.1644 0 0
0 -2.0009 1.8221 2.1644 0 0 0 1.8221 2.0009 -3.0587 0 0
0 3.0587 -2.1644 -2.1470 0 0 0 -2.1644 -3.0587 6.6185 0 0
0.8617 0 0 0 -2.0248 -0.3417 -0.8617 0 0 0 3.9669 -2.2838
-1.2352 0 0 0 0.3417 3.8922 1.2352 0 0 0 -2.2838 5.8342
K_l(:,:,4) =
1.0e+003 *
0.4194 0 0 0 -0.8617 -1.2352 -0.4194 0 0 0 0.8617 -1.2352
0 2.0009 -1.8221 -2.1644 0 0 0 -1.8221 -2.0009 3.0587 0 0
0 -1.8221 2.0009 3.0587 0 0 0 2.0009 1.8221 -2.1644 0 0
0 -2.1644 3.0587 6.6185 0 0 0 3.0587 2.1644 -2.1470 0 0
-0.8617 0 0 0 3.9669 2.2838 0.8617 0 0 0 -2.0248 0.3417
-1.2352 0 0 0 2.2838 5.8342 1.2352 0 0 0 -0.3417 3.8922
-0.4194 0 0 0 0.8617 1.2352 0.4194 0 0 0 -0.8617 1.2352
0 -1.8221 2.0009 3.0587 0 0 0 2.0009 1.8221 -2.1644 0 0
0 -2.0009 1.8221 2.1644 0 0 0 1.8221 2.0009 -3.0587 0 0
0 3.0587 -2.1644 -2.1470 0 0 0 -2.1644 -3.0587 6.6185 0 0
0.8617 0 0 0 -2.0248 -0.3417 -0.8617 0 0 0 3.9669 -2.2838
-1.2352 0 0 0 0.3417 3.8922 1.2352 0 0 0 -2.2838 5.8342
K_l(:,:,5) =
1.0e+003 *
0.4194 0 0 0 -0.8617 -1.2352 -0.4194 0 0 0 0.8617 -1.2352
0 2.0009 -1.8221 -2.1644 0 0 0 -1.8221 -2.0009 3.0587 0 0
0 -1.8221 2.0009 3.0587 0 0 0 2.0009 1.8221 -2.1644 0 0
0 -2.1644 3.0587 6.6185 0 0 0 3.0587 2.1644 -2.1470 0 0
-0.8617 0 0 0 3.9669 2.2838 0.8617 0 0 0 -2.0248 0.3417
-1.2352 0 0 0 2.2838 5.8342 1.2352 0 0 0 -0.3417 3.8922
-0.4194 0 0 0 0.8617 1.2352 0.4194 0 0 0 -0.8617 1.2352
0 -1.8221 2.0009 3.0587 0 0 0 2.0009 1.8221 -2.1644 0 0
0 -2.0009 1.8221 2.1644 0 0 0 1.8221 2.0009 -3.0587 0 0
0 3.0587 -2.1644 -2.1470 0 0 0 -2.1644 -3.0587 6.6185 0 0
0.8617 0 0 0 -2.0248 -0.3417 -0.8617 0 0 0 3.9669 -2.2838
-1.2352 0 0 0 0.3417 3.8922 1.2352 0 0 0 -2.2838 5.8342
K_l(:,:,6) =
1.0e+003 *
0.4194 0 0 0 -0.8617 -1.2352 -0.4194 0 0 0 0.8617 -1.2352
0 2.0009 -1.8221 -2.1644 0 0 0 -1.8221 -2.0009 3.0587 0 0
0 -1.8221 2.0009 3.0587 0 0 0 2.0009 1.8221 -2.1644 0 0
0 -2.1644 3.0587 6.6185 0 0 0 3.0587 2.1644 -2.1470 0 0
-0.8617 0 0 0 3.9669 2.2838 0.8617 0 0 0 -2.0248 0.3417
-1.2352 0 0 0 2.2838 5.8342 1.2352 0 0 0 -0.3417 3.8922
-0.4194 0 0 0 0.8617 1.2352 0.4194 0 0 0 -0.8617 1.2352
0 -1.8221 2.0009 3.0587 0 0 0 2.0009 1.8221 -2.1644 0 0
0 -2.0009 1.8221 2.1644 0 0 0 1.8221 2.0009 -3.0587 0 0
0 3.0587 -2.1644 -2.1470 0 0 0 -2.1644 -3.0587 6.6185 0 0
0.8617 0 0 0 -2.0248 -0.3417 -0.8617 0 0 0 3.9669 -2.2838
-1.2352 0 0 0 0.3417 3.8922 1.2352 0 0 0 -2.2838 5.8342
K_l(:,:,7) =
1.0e+003 *
0.4194 0 0 0 -0.8617 -1.2352 -0.4194 0 0 0 0.8617 -1.2352
0 2.0009 -1.8221 -2.1644 0 0 0 -1.8221 -2.0009 3.0587 0 0
0 -1.8221 2.0009 3.0587 0 0 0 2.0009 1.8221 -2.1644 0 0
0 -2.1644 3.0587 6.6185 0 0 0 3.0587 2.1644 -2.1470 0 0
-0.8617 0 0 0 3.9669 2.2838 0.8617 0 0 0 -2.0248 0.3417
-1.2352 0 0 0 2.2838 5.8342 1.2352 0 0 0 -0.3417 3.8922
-0.4194 0 0 0 0.8617 1.2352 0.4194 0 0 0 -0.8617 1.2352
0 -1.8221 2.0009 3.0587 0 0 0 2.0009 1.8221 -2.1644 0 0
0 -2.0009 1.8221 2.1644 0 0 0 1.8221 2.0009 -3.0587 0 0
0 3.0587 -2.1644 -2.1470 0 0 0 -2.1644 -3.0587 6.6185 0 0
0.8617 0 0 0 -2.0248 -0.3417 -0.8617 0 0 0 3.9669 -2.2838
-1.2352 0 0 0 0.3417 3.8922 1.2352 0 0 0 -2.2838 5.8342
K_l(:,:,8) =
1.0e+003 *
0.4194 0 0 0 -0.8617 -1.2352 -0.4194 0 0 0 0.8617 -1.2352
0 2.0009 -1.8221 -2.1644 0 0 0 -1.8221 -2.0009 3.0587 0 0
0 -1.8221 2.0009 3.0587 0 0 0 2.0009 1.8221 -2.1644 0 0
0 -2.1644 3.0587 6.6185 0 0 0 3.0587 2.1644 -2.1470 0 0
-0.8617 0 0 0 3.9669 2.2838 0.8617 0 0 0 -2.0248 0.3417
-1.2352 0 0 0 2.2838 5.8342 1.2352 0 0 0 -0.3417 3.8922
-0.4194 0 0 0 0.8617 1.2352 0.4194 0 0 0 -0.8617 1.2352
0 -1.8221 2.0009 3.0587 0 0 0 2.0009 1.8221 -2.1644 0 0
0 -2.0009 1.8221 2.1644 0 0 0 1.8221 2.0009 -3.0587 0 0
0 3.0587 -2.1644 -2.1470 0 0 0 -2.1644 -3.0587 6.6185 0 0
0.8617 0 0 0 -2.0248 -0.3417 -0.8617 0 0 0 3.9669 -2.2838
-1.2352 0 0 0 0.3417 3.8922 1.2352 0 0 0 -2.2838 5.8342
K_l(:,:,9) =
1.0e+003 *
0.4194 0 0 0 -0.8617 -1.2352 -0.4194 0 0 0 0.8617 -1.2352
0 2.0009 -1.8221 -2.1644 0 0 0 -1.8221 -2.0009 3.0587 0 0
0 -1.8221 2.0009 3.0587 0 0 0 2.0009 1.8221 -2.1644 0 0
0 -2.1644 3.0587 6.6185 0 0 0 3.0587 2.1644 -2.1470 0 0
-0.8617 0 0 0 3.9669 2.2838 0.8617 0 0 0 -2.0248 0.3417
-1.2352 0 0 0 2.2838 5.8342 1.2352 0 0 0 -0.3417 3.8922
-0.4194 0 0 0 0.8617 1.2352 0.4194 0 0 0 -0.8617 1.2352
0 -1.8221 2.0009 3.0587 0 0 0 2.0009 1.8221 -2.1644 0 0
0 -2.0009 1.8221 2.1644 0 0 0 1.8221 2.0009 -3.0587 0 0
0 3.0587 -2.1644 -2.1470 0 0 0 -2.1644 -3.0587 6.6185 0 0
0.8617 0 0 0 -2.0248 -0.3417 -0.8617 0 0 0 3.9669 -2.2838
-1.2352 0 0 0 0.3417 3.8922 1.2352 0 0 0 -2.2838 5.8342
K_l(:,:,10) =
1.0e+003 *
0.4194 0 0 0 -0.8617 -1.2352 -0.4194 0 0 0 0.8617 -1.2352
0 2.0009 -1.8221 -2.1644 0 0 0 -1.8221 -2.0009 3.0587 0 0
0 -1.8221 2.0009 3.0587 0 0 0 2.0009 1.8221 -2.1644 0 0
0 -2.1644 3.0587 6.6185 0 0 0 3.0587 2.1644 -2.1470 0 0
-0.8617 0 0 0 3.9669 2.2838 0.8617 0 0 0 -2.0248 0.3417
-1.2352 0 0 0 2.2838 5.8342 1.2352 0 0 0 -0.3417 3.8922
-0.4194 0 0 0 0.8617 1.2352 0.4194 0 0 0 -0.8617 1.2352
0 -1.8221 2.0009 3.0587 0 0 0 2.0009 1.8221 -2.1644 0 0
0 -2.0009 1.8221 2.1644 0 0 0 1.8221 2.0009 -3.0587 0 0
0 3.0587 -2.1644 -2.1470 0 0 0 -2.1644 -3.0587 6.6185 0 0
0.8617 0 0 0 -2.0248 -0.3417 -0.8617 0 0 0 3.9669 -2.2838
-1.2352 0 0 0 0.3417 3.8922 1.2352 0 0 0 -2.2838 5.8342
K_l(:,:,11) =
1.0e+003 *
0.4194 0 0 0 -0.8617 -1.2352 -0.4194 0 0 0 0.8617 -1.2352
0 2.0009 -1.8221 -2.1644 0 0 0 -1.8221 -2.0009 3.0587 0 0
0 -1.8221 2.0009 3.0587 0 0 0 2.0009 1.8221 -2.1644 0 0
0 -2.1644 3.0587 6.6185 0 0 0 3.0587 2.1644 -2.1470 0 0
-0.8617 0 0 0 3.9669 2.2838 0.8617 0 0 0 -2.0248 0.3417
-1.2352 0 0 0 2.2838 5.8342 1.2352 0 0 0 -0.3417 3.8922
-0.4194 0 0 0 0.8617 1.2352 0.4194 0 0 0 -0.8617 1.2352
0 -1.8221 2.0009 3.0587 0 0 0 2.0009 1.8221 -2.1644 0 0
0 -2.0009 1.8221 2.1644 0 0 0 1.8221 2.0009 -3.0587 0 0
0 3.0587 -2.1644 -2.1470 0 0 0 -2.1644 -3.0587 6.6185 0 0
0.8617 0 0 0 -2.0248 -0.3417 -0.8617 0 0 0 3.9669 -2.2838
-1.2352 0 0 0 0.3417 3.8922 1.2352 0 0 0 -2.2838 5.8342
K_l(:,:,12) =
1.0e+003 *
0.4194 0 0 0 -0.8617 -1.2352 -0.4194 0 0 0 0.8617 -1.2352
0 2.0009 -1.8221 -2.1644 0 0 0 -1.8221 -2.0009 3.0587 0 0
0 -1.8221 2.0009 3.0587 0 0 0 2.0009 1.8221 -2.1644 0 0
0 -2.1644 3.0587 6.6185 0 0 0 3.0587 2.1644 -2.1470 0 0
-0.8617 0 0 0 3.9669 2.2838 0.8617 0 0 0 -2.0248 0.3417
-1.2352 0 0 0 2.2838 5.8342 1.2352 0 0 0 -0.3417 3.8922
-0.4194 0 0 0 0.8617 1.2352 0.4194 0 0 0 -0.8617 1.2352
0 -1.8221 2.0009 3.0587 0 0 0 2.0009 1.8221 -2.1644 0 0
0 -2.0009 1.8221 2.1644 0 0 0 1.8221 2.0009 -3.0587 0 0
0 3.0587 -2.1644 -2.1470 0 0 0 -2.1644 -3.0587 6.6185 0 0
0.8617 0 0 0 -2.0248 -0.3417 -0.8617 0 0 0 3.9669 -2.2838
-1.2352 0 0 0 0.3417 3.8922 1.2352 0 0 0 -2.2838 5.8342
K_l(:,:,13) =
1.0e+003 *
0.4194 0 0 0 -0.8617 -1.2352 -0.4194 0 0 0 0.8617 -1.2352
0 2.0009 -1.8221 -2.1644 0 0 0 -1.8221 -2.0009 3.0587 0 0
0 -1.8221 2.0009 3.0587 0 0 0 2.0009 1.8221 -2.1644 0 0
0 -2.1644 3.0587 6.6185 0 0 0 3.0587 2.1644 -2.1470 0 0
-0.8617 0 0 0 3.9669 2.2838 0.8617 0 0 0 -2.0248 0.3417
-1.2352 0 0 0 2.2838 5.8342 1.2352 0 0 0 -0.3417 3.8922
-0.4194 0 0 0 0.8617 1.2352 0.4194 0 0 0 -0.8617 1.2352
0 -1.8221 2.0009 3.0587 0 0 0 2.0009 1.8221 -2.1644 0 0
0 -2.0009 1.8221 2.1644 0 0 0 1.8221 2.0009 -3.0587 0 0
0 3.0587 -2.1644 -2.1470 0 0 0 -2.1644 -3.0587 6.6185 0 0
0.8617 0 0 0 -2.0248 -0.3417 -0.8617 0 0 0 3.9669 -2.2838
-1.2352 0 0 0 0.3417 3.8922 1.2352 0 0 0 -2.2838 5.8342
K_l(:,:,14) =
1.0e+003 *
0.4194 0 0 0 -0.8617 -1.2352 -0.4194 0 0 0 0.8617 -1.2352
0 2.0009 -1.8221 -2.1644 0 0 0 -1.8221 -2.0009 3.0587 0 0
0 -1.8221 2.0009 3.0587 0 0 0 2.0009 1.8221 -2.1644 0 0
0 -2.1644 3.0587 6.6185 0 0 0 3.0587 2.1644 -2.1470 0 0
-0.8617 0 0 0 3.9669 2.2838 0.8617 0 0 0 -2.0248 0.3417
-1.2352 0 0 0 2.2838 5.8342 1.2352 0 0 0 -0.3417 3.8922
-0.4194 0 0 0 0.8617 1.2352 0.4194 0 0 0 -0.8617 1.2352
0 -1.8221 2.0009 3.0587 0 0 0 2.0009 1.8221 -2.1644 0 0
0 -2.0009 1.8221 2.1644 0 0 0 1.8221 2.0009 -3.0587 0 0
0 3.0587 -2.1644 -2.1470 0 0 0 -2.1644 -3.0587 6.6185 0 0
0.8617 0 0 0 -2.0248 -0.3417 -0.8617 0 0 0 3.9669 -2.2838
-1.2352 0 0 0 0.3417 3.8922 1.2352 0 0 0 -2.2838 5.8342
K_l(:,:,15) =
1.0e+003 *
0.4194 0 0 0 -0.8617 -1.2352 -0.4194 0 0 0 0.8617 -1.2352
0 2.0009 -1.8221 -2.1644 0 0 0 -1.8221 -2.0009 3.0587 0 0
0 -1.8221 2.0009 3.0587 0 0 0 2.0009 1.8221 -2.1644 0 0
0 -2.1644 3.0587 6.6185 0 0 0 3.0587 2.1644 -2.1470 0 0
-0.8617 0 0 0 3.9669 2.2838 0.8617 0 0 0 -2.0248 0.3417
-1.2352 0 0 0 2.2838 5.8342 1.2352 0 0 0 -0.3417 3.8922
-0.4194 0 0 0 0.8617 1.2352 0.4194 0 0 0 -0.8617 1.2352
0 -1.8221 2.0009 3.0587 0 0 0 2.0009 1.8221 -2.1644 0 0
0 -2.0009 1.8221 2.1644 0 0 0 1.8221 2.0009 -3.0587 0 0
0 3.0587 -2.1644 -2.1470 0 0 0 -2.1644 -3.0587 6.6185 0 0
0.8617 0 0 0 -2.0248 -0.3417 -0.8617 0 0 0 3.9669 -2.2838
-1.2352 0 0 0 0.3417 3.8922 1.2352 0 0 0 -2.2838 5.8342
K_l(:,:,16) =
1.0e+003 *
0.4194 0 0 0 -0.8617 -1.2352 -0.4194 0 0 0 0.8617 -1.2352
0 2.0009 -1.8221 -2.1644 0 0 0 -1.8221 -2.0009 3.0587 0 0
0 -1.8221 2.0009 3.0587 0 0 0 2.0009 1.8221 -2.1644 0 0
0 -2.1644 3.0587 6.6185 0 0 0 3.0587 2.1644 -2.1470 0 0
-0.8617 0 0 0 3.9669 2.2838 0.8617 0 0 0 -2.0248 0.3417
-1.2352 0 0 0 2.2838 5.8342 1.2352 0 0 0 -0.3417 3.8922
-0.4194 0 0 0 0.8617 1.2352 0.4194 0 0 0 -0.8617 1.2352
0 -1.8221 2.0009 3.0587 0 0 0 2.0009 1.8221 -2.1644 0 0
0 -2.0009 1.8221 2.1644 0 0 0 1.8221 2.0009 -3.0587 0 0
0 3.0587 -2.1644 -2.1470 0 0 0 -2.1644 -3.0587 6.6185 0 0
0.8617 0 0 0 -2.0248 -0.3417 -0.8617 0 0 0 3.9669 -2.2838
-1.2352 0 0 0 0.3417 3.8922 1.2352 0 0 0 -2.2838 5.8342
K_l(:,:,17) =
1.0e+003 *
0.4194 0 0 0 -0.8617 -1.2352 -0.4194 0 0 0 0.8617 -1.2352
0 2.0009 -1.8221 -2.1644 0 0 0 -1.8221 -2.0009 3.0587 0 0
0 -1.8221 2.0009 3.0587 0 0 0 2.0009 1.8221 -2.1644 0 0
0 -2.1644 3.0587 6.6185 0 0 0 3.0587 2.1644 -2.1470 0 0
-0.8617 0 0 0 3.9669 2.2838 0.8617 0 0 0 -2.0248 0.3417
-1.2352 0 0 0 2.2838 5.8342 1.2352 0 0 0 -0.3417 3.8922
-0.4194 0 0 0 0.8617 1.2352 0.4194 0 0 0 -0.8617 1.2352
0 -1.8221 2.0009 3.0587 0 0 0 2.0009 1.8221 -2.1644 0 0
0 -2.0009 1.8221 2.1644 0 0 0 1.8221 2.0009 -3.0587 0 0
0 3.0587 -2.1644 -2.1470 0 0 0 -2.1644 -3.0587 6.6185 0 0
0.8617 0 0 0 -2.0248 -0.3417 -0.8617 0 0 0 3.9669 -2.2838
-1.2352 0 0 0 0.3417 3.8922 1.2352 0 0 0 -2.2838 5.8342
K_l(:,:,18) =
1.0e+003 *
0.4194 0 0 0 -0.8617 -1.2352 -0.4194 0 0 0 0.8617 -1.2352
0 2.0009 -1.8221 -2.1644 0 0 0 -1.8221 -2.0009 3.0587 0 0
0 -1.8221 2.0009 3.0587 0 0 0 2.0009 1.8221 -2.1644 0 0
0 -2.1644 3.0587 6.6185 0 0 0 3.0587 2.1644 -2.1470 0 0
-0.8617 0 0 0 3.9669 2.2838 0.8617 0 0 0 -2.0248 0.3417
-1.2352 0 0 0 2.2838 5.8342 1.2352 0 0 0 -0.3417 3.8922
-0.4194 0 0 0 0.8617 1.2352 0.4194 0 0 0 -0.8617 1.2352
0 -1.8221 2.0009 3.0587 0 0 0 2.0009 1.8221 -2.1644 0 0
0 -2.0009 1.8221 2.1644 0 0 0 1.8221 2.0009 -3.0587 0 0
0 3.0587 -2.1644 -2.1470 0 0 0 -2.1644 -3.0587 6.6185 0 0
0.8617 0 0 0 -2.0248 -0.3417 -0.8617 0 0 0 3.9669 -2.2838
-1.2352 0 0 0 0.3417 3.8922 1.2352 0 0 0 -2.2838 5.8342
K_l(:,:,19) =
1.0e+003 *
0.4194 0 0 0 -0.8617 -1.2352 -0.4194 0 0 0 0.8617 -1.2352
0 2.0009 -1.8221 -2.1644 0 0 0 -1.8221 -2.0009 3.0587 0 0
0 -1.8221 2.0009 3.0587 0 0 0 2.0009 1.8221 -2.1644 0 0
0 -2.1644 3.0587 6.6185 0 0 0 3.0587 2.1644 -2.1470 0 0
-0.8617 0 0 0 3.9669 2.2838 0.8617 0 0 0 -2.0248 0.3417
-1.2352 0 0 0 2.2838 5.8342 1.2352 0 0 0 -0.3417 3.8922
-0.4194 0 0 0 0.8617 1.2352 0.4194 0 0 0 -0.8617 1.2352
0 -1.8221 2.0009 3.0587 0 0 0 2.0009 1.8221 -2.1644 0 0
0 -2.0009 1.8221 2.1644 0 0 0 1.8221 2.0009 -3.0587 0 0
0 3.0587 -2.1644 -2.1470 0 0 0 -2.1644 -3.0587 6.6185 0 0
0.8617 0 0 0 -2.0248 -0.3417 -0.8617 0 0 0 3.9669 -2.2838
-1.2352 0 0 0 0.3417 3.8922 1.2352 0 0 0 -2.2838 5.8342
K_l(:,:,20) =
1.0e+003 *
0.4194 0 0 0 -0.8617 -1.2352 -0.4194 0 0 0 0.8617 -1.2352
0 2.0009 -1.8221 -2.1644 0 0 0 -1.8221 -2.0009 3.0587 0 0
0 -1.8221 2.0009 3.0587 0 0 0 2.0009 1.8221 -2.1644 0 0
0 -2.1644 3.0587 6.6185 0 0 0 3.0587 2.1644 -2.1470 0 0
-0.8617 0 0 0 3.9669 2.2838 0.8617 0 0 0 -2.0248 0.3417
-1.2352 0 0 0 2.2838 5.8342 1.2352 0 0 0 -0.3417 3.8922
-0.4194 0 0 0 0.8617 1.2352 0.4194 0 0 0 -0.8617 1.2352
0 -1.8221 2.0009 3.0587 0 0 0 2.0009 1.8221 -2.1644 0 0
0 -2.0009 1.8221 2.1644 0 0 0 1.8221 2.0009 -3.0587 0 0
0 3.0587 -2.1644 -2.1470 0 0 0 -2.1644 -3.0587 6.6185 0 0
0.8617 0 0 0 -2.0248 -0.3417 -0.8617 0 0 0 3.9669 -2.2838
-1.2352 0 0 0 0.3417 3.8922 1.2352 0 0 0 -2.2838 5.8342
K_l(:,:,21) =
1.0e+003 *
0.4194 0 0 0 -0.8617 -1.2352 -0.4194 0 0 0 0.8617 -1.2352
0 2.0009 -1.8221 -2.1644 0 0 0 -1.8221 -2.0009 3.0587 0 0
0 -1.8221 2.0009 3.0587 0 0 0 2.0009 1.8221 -2.1644 0 0
0 -2.1644 3.0587 6.6185 0 0 0 3.0587 2.1644 -2.1470 0 0
-0.8617 0 0 0 3.9669 2.2838 0.8617 0 0 0 -2.0248 0.3417
-1.2352 0 0 0 2.2838 5.8342 1.2352 0 0 0 -0.3417 3.8922
-0.4194 0 0 0 0.8617 1.2352 0.4194 0 0 0 -0.8617 1.2352
0 -1.8221 2.0009 3.0587 0 0 0 2.0009 1.8221 -2.1644 0 0
0 -2.0009 1.8221 2.1644 0 0 0 1.8221 2.0009 -3.0587 0 0
0 3.0587 -2.1644 -2.1470 0 0 0 -2.1644 -3.0587 6.6185 0 0
0.8617 0 0 0 -2.0248 -0.3417 -0.8617 0 0 0 3.9669 -2.2838
-1.2352 0 0 0 0.3417 3.8922 1.2352 0 0 0 -2.2838 5.8342
K_l(:,:,22) =
1.0e+003 *
0.4194 0 0 0 -0.8617 -1.2352 -0.4194 0 0 0 0.8617 -1.2352
0 2.0009 -1.8221 -2.1644 0 0 0 -1.8221 -2.0009 3.0587 0 0
0 -1.8221 2.0009 3.0587 0 0 0 2.0009 1.8221 -2.1644 0 0
0 -2.1644 3.0587 6.6185 0 0 0 3.0587 2.1644 -2.1470 0 0
-0.8617 0 0 0 3.9669 2.2838 0.8617 0 0 0 -2.0248 0.3417
-1.2352 0 0 0 2.2838 5.8342 1.2352 0 0 0 -0.3417 3.8922
-0.4194 0 0 0 0.8617 1.2352 0.4194 0 0 0 -0.8617 1.2352
0 -1.8221 2.0009 3.0587 0 0 0 2.0009 1.8221 -2.1644 0 0
0 -2.0009 1.8221 2.1644 0 0 0 1.8221 2.0009 -3.0587 0 0
0 3.0587 -2.1644 -2.1470 0 0 0 -2.1644 -3.0587 6.6185 0 0
0.8617 0 0 0 -2.0248 -0.3417 -0.8617 0 0 0 3.9669 -2.2838
-1.2352 0 0 0 0.3417 3.8922 1.2352 0 0 0 -2.2838 5.8342
K_l(:,:,23) =
1.0e+003 *
0.4194 0 0 0 -0.8617 -1.2352 -0.4194 0 0 0 0.8617 -1.2352
0 2.0009 -1.8221 -2.1644 0 0 0 -1.8221 -2.0009 3.0587 0 0
0 -1.8221 2.0009 3.0587 0 0 0 2.0009 1.8221 -2.1644 0 0
0 -2.1644 3.0587 6.6185 0 0 0 3.0587 2.1644 -2.1470 0 0
-0.8617 0 0 0 3.9669 2.2838 0.8617 0 0 0 -2.0248 0.3417
-1.2352 0 0 0 2.2838 5.8342 1.2352 0 0 0 -0.3417 3.8922
-0.4194 0 0 0 0.8617 1.2352 0.4194 0 0 0 -0.8617 1.2352
0 -1.8221 2.0009 3.0587 0 0 0 2.0009 1.8221 -2.1644 0 0
0 -2.0009 1.8221 2.1644 0 0 0 1.8221 2.0009 -3.0587 0 0
0 3.0587 -2.1644 -2.1470 0 0 0 -2.1644 -3.0587 6.6185 0 0
0.8617 0 0 0 -2.0248 -0.3417 -0.8617 0 0 0 3.9669 -2.2838
-1.2352 0 0 0 0.3417 3.8922 1.2352 0 0 0 -2.2838 5.8342
K_l(:,:,24) =
1.0e+003 *
0.4194 0 0 0 -0.8617 -1.2352 -0.4194 0 0 0 0.8617 -1.2352
0 2.0009 -1.8221 -2.1644 0 0 0 -1.8221 -2.0009 3.0587 0 0
0 -1.8221 2.0009 3.0587 0 0 0 2.0009 1.8221 -2.1644 0 0
0 -2.1644 3.0587 6.6185 0 0 0 3.0587 2.1644 -2.1470 0 0
-0.8617 0 0 0 3.9669 2.2838 0.8617 0 0 0 -2.0248 0.3417
-1.2352 0 0 0 2.2838 5.8342 1.2352 0 0 0 -0.3417 3.8922
-0.4194 0 0 0 0.8617 1.2352 0.4194 0 0 0 -0.8617 1.2352
0 -1.8221 2.0009 3.0587 0 0 0 2.0009 1.8221 -2.1644 0 0
0 -2.0009 1.8221 2.1644 0 0 0 1.8221 2.0009 -3.0587 0 0
0 3.0587 -2.1644 -2.1470 0 0 0 -2.1644 -3.0587 6.6185 0 0
0.8617 0 0 0 -2.0248 -0.3417 -0.8617 0 0 0 3.9669 -2.2838
-1.2352 0 0 0 0.3417 3.8922 1.2352 0 0 0 -2.2838 5.8342
K_l(:,:,25) =
1.0e+003 *
0.4194 0 0 0 -0.8617 -1.2352 -0.4194 0 0 0 0.8617 -1.2352
0 2.0009 -1.8221 -2.1644 0 0 0 -1.8221 -2.0009 3.0587 0 0
0 -1.8221 2.0009 3.0587 0 0 0 2.0009 1.8221 -2.1644 0 0
0 -2.1644 3.0587 6.6185 0 0 0 3.0587 2.1644 -2.1470 0 0
-0.8617 0 0 0 3.9669 2.2838 0.8617 0 0 0 -2.0248 0.3417
-1.2352 0 0 0 2.2838 5.8342 1.2352 0 0 0 -0.3417 3.8922
-0.4194 0 0 0 0.8617 1.2352 0.4194 0 0 0 -0.8617 1.2352
0 -1.8221 2.0009 3.0587 0 0 0 2.0009 1.8221 -2.1644 0 0
0 -2.0009 1.8221 2.1644 0 0 0 1.8221 2.0009 -3.0587 0 0
0 3.0587 -2.1644 -2.1470 0 0 0 -2.1644 -3.0587 6.6185 0 0
0.8617 0 0 0 -2.0248 -0.3417 -0.8617 0 0 0 3.9669 -2.2838
-1.2352 0 0 0 0.3417 3.8922 1.2352 0 0 0 -2.2838 5.8342
K_l(:,:,26) =
1.0e+003 *
0.4194 0 0 0 -0.8617 -1.2352 -0.4194 0 0 0 0.8617 -1.2352
0 2.0009 -1.8221 -2.1644 0 0 0 -1.8221 -2.0009 3.0587 0 0
0 -1.8221 2.0009 3.0587 0 0 0 2.0009 1.8221 -2.1644 0 0
0 -2.1644 3.0587 6.6185 0 0 0 3.0587 2.1644 -2.1470 0 0
-0.8617 0 0 0 3.9669 2.2838 0.8617 0 0 0 -2.0248 0.3417
-1.2352 0 0 0 2.2838 5.8342 1.2352 0 0 0 -0.3417 3.8922
-0.4194 0 0 0 0.8617 1.2352 0.4194 0 0 0 -0.8617 1.2352
0 -1.8221 2.0009 3.0587 0 0 0 2.0009 1.8221 -2.1644 0 0
0 -2.0009 1.8221 2.1644 0 0 0 1.8221 2.0009 -3.0587 0 0
0 3.0587 -2.1644 -2.1470 0 0 0 -2.1644 -3.0587 6.6185 0 0
0.8617 0 0 0 -2.0248 -0.3417 -0.8617 0 0 0 3.9669 -2.2838
-1.2352 0 0 0 0.3417 3.8922 1.2352 0 0 0 -2.2838 5.8342
K_l(:,:,27) =
1.0e+003 *
0.4194 0 0 0 -0.8617 -1.2352 -0.4194 0 0 0 0.8617 -1.2352
0 2.0009 -1.8221 -2.1644 0 0 0 -1.8221 -2.0009 3.0587 0 0
0 -1.8221 2.0009 3.0587 0 0 0 2.0009 1.8221 -2.1644 0 0
0 -2.1644 3.0587 6.6185 0 0 0 3.0587 2.1644 -2.1470 0 0
-0.8617 0 0 0 3.9669 2.2838 0.8617 0 0 0 -2.0248 0.3417
-1.2352 0 0 0 2.2838 5.8342 1.2352 0 0 0 -0.3417 3.8922
-0.4194 0 0 0 0.8617 1.2352 0.4194 0 0 0 -0.8617 1.2352
0 -1.8221 2.0009 3.0587 0 0 0 2.0009 1.8221 -2.1644 0 0
0 -2.0009 1.8221 2.1644 0 0 0 1.8221 2.0009 -3.0587 0 0
0 3.0587 -2.1644 -2.1470 0 0 0 -2.1644 -3.0587 6.6185 0 0
0.8617 0 0 0 -2.0248 -0.3417 -0.8617 0 0 0 3.9669 -2.2838
-1.2352 0 0 0 0.3417 3.8922 1.2352 0 0 0 -2.2838 5.8342
K_l(:,:,28) =
1.0e+003 *
0.4194 0 0 0 -0.8617 -1.2352 -0.4194 0 0 0 0.8617 -1.2352
0 2.0009 -1.8221 -2.1644 0 0 0 -1.8221 -2.0009 3.0587 0 0
0 -1.8221 2.0009 3.0587 0 0 0 2.0009 1.8221 -2.1644 0 0
0 -2.1644 3.0587 6.6185 0 0 0 3.0587 2.1644 -2.1470 0 0
-0.8617 0 0 0 3.9669 2.2838 0.8617 0 0 0 -2.0248 0.3417
-1.2352 0 0 0 2.2838 5.8342 1.2352 0 0 0 -0.3417 3.8922
-0.4194 0 0 0 0.8617 1.2352 0.4194 0 0 0 -0.8617 1.2352
0 -1.8221 2.0009 3.0587 0 0 0 2.0009 1.8221 -2.1644 0 0
0 -2.0009 1.8221 2.1644 0 0 0 1.8221 2.0009 -3.0587 0 0
0 3.0587 -2.1644 -2.1470 0 0 0 -2.1644 -3.0587 6.6185 0 0
0.8617 0 0 0 -2.0248 -0.3417 -0.8617 0 0 0 3.9669 -2.2838
-1.2352 0 0 0 0.3417 3.8922 1.2352 0 0 0 -2.2838 5.8342
K_l(:,:,29) =
1.0e+003 *
0.4194 0 0 0 -0.8617 -1.2352 -0.4194 0 0 0 0.8617 -1.2352
0 2.0009 -1.8221 -2.1644 0 0 0 -1.8221 -2.0009 3.0587 0 0
0 -1.8221 2.0009 3.0587 0 0 0 2.0009 1.8221 -2.1644 0 0
0 -2.1644 3.0587 6.6185 0 0 0 3.0587 2.1644 -2.1470 0 0
-0.8617 0 0 0 3.9669 2.2838 0.8617 0 0 0 -2.0248 0.3417
-1.2352 0 0 0 2.2838 5.8342 1.2352 0 0 0 -0.3417 3.8922
-0.4194 0 0 0 0.8617 1.2352 0.4194 0 0 0 -0.8617 1.2352
0 -1.8221 2.0009 3.0587 0 0 0 2.0009 1.8221 -2.1644 0 0
0 -2.0009 1.8221 2.1644 0 0 0 1.8221 2.0009 -3.0587 0 0
0 3.0587 -2.1644 -2.1470 0 0 0 -2.1644 -3.0587 6.6185 0 0
0.8617 0 0 0 -2.0248 -0.3417 -0.8617 0 0 0 3.9669 -2.2838
-1.2352 0 0 0 0.3417 3.8922 1.2352 0 0 0 -2.2838 5.8342
K_l(:,:,30) =
1.0e+003 *
0.4194 0 0 0 -0.8617 -1.2352 -0.4194 0 0 0 0.8617 -1.2352
0 2.0009 -1.8221 -2.1644 0 0 0 -1.8221 -2.0009 3.0587 0 0
0 -1.8221 2.0009 3.0587 0 0 0 2.0009 1.8221 -2.1644 0 0
0 -2.1644 3.0587 6.6185 0 0 0 3.0587 2.1644 -2.1470 0 0
-0.8617 0 0 0 3.9669 2.2838 0.8617 0 0 0 -2.0248 0.3417
-1.2352 0 0 0 2.2838 5.8342 1.2352 0 0 0 -0.3417 3.8922
-0.4194 0 0 0 0.8617 1.2352 0.4194 0 0 0 -0.8617 1.2352
0 -1.8221 2.0009 3.0587 0 0 0 2.0009 1.8221 -2.1644 0 0
0 -2.0009 1.8221 2.1644 0 0 0 1.8221 2.0009 -3.0587 0 0
0 3.0587 -2.1644 -2.1470 0 0 0 -2.1644 -3.0587 6.6185 0 0
0.8617 0 0 0 -2.0248 -0.3417 -0.8617 0 0 0 3.9669 -2.2838
-1.2352 0 0 0 0.3417 3.8922 1.2352 0 0 0 -2.2838 5.8342
K_l(:,:,31) =
1.0e+003 *
0.4194 0 0 0 -0.8617 -1.2352 -0.4194 0 0 0 0.8617 -1.2352
0 2.0009 -1.8221 -2.1644 0 0 0 -1.8221 -2.0009 3.0587 0 0
0 -1.8221 2.0009 3.0587 0 0 0 2.0009 1.8221 -2.1644 0 0
0 -2.1644 3.0587 6.6185 0 0 0 3.0587 2.1644 -2.1470 0 0
-0.8617 0 0 0 3.9669 2.2838 0.8617 0 0 0 -2.0248 0.3417
-1.2352 0 0 0 2.2838 5.8342 1.2352 0 0 0 -0.3417 3.8922
-0.4194 0 0 0 0.8617 1.2352 0.4194 0 0 0 -0.8617 1.2352
0 -1.8221 2.0009 3.0587 0 0 0 2.0009 1.8221 -2.1644 0 0
0 -2.0009 1.8221 2.1644 0 0 0 1.8221 2.0009 -3.0587 0 0
0 3.0587 -2.1644 -2.1470 0 0 0 -2.1644 -3.0587 6.6185 0 0
0.8617 0 0 0 -2.0248 -0.3417 -0.8617 0 0 0 3.9669 -2.2838
-1.2352 0 0 0 0.3417 3.8922 1.2352 0 0 0 -2.2838 5.8342
K_l(:,:,32) =
1.0e+003 *
0.4194 0 0 0 -0.8617 -1.2352 -0.4194 0 0 0 0.8617 -1.2352
0 2.0009 -1.8221 -2.1644 0 0 0 -1.8221 -2.0009 3.0587 0 0
0 -1.8221 2.0009 3.0587 0 0 0 2.0009 1.8221 -2.1644 0 0
0 -2.1644 3.0587 6.6185 0 0 0 3.0587 2.1644 -2.1470 0 0
-0.8617 0 0 0 3.9669 2.2838 0.8617 0 0 0 -2.0248 0.3417
-1.2352 0 0 0 2.2838 5.8342 1.2352 0 0 0 -0.3417 3.8922
-0.4194 0 0 0 0.8617 1.2352 0.4194 0 0 0 -0.8617 1.2352
0 -1.8221 2.0009 3.0587 0 0 0 2.0009 1.8221 -2.1644 0 0
0 -2.0009 1.8221 2.1644 0 0 0 1.8221 2.0009 -3.0587 0 0
0 3.0587 -2.1644 -2.1470 0 0 0 -2.1644 -3.0587 6.6185 0 0
0.8617 0 0 0 -2.0248 -0.3417 -0.8617 0 0 0 3.9669 -2.2838
-1.2352 0 0 0 0.3417 3.8922 1.2352 0 0 0 -2.2838 5.8342
K_g
K_g(:,:,1) =
1.0e+003 *
2.1212 1.7018 -0.0000 -1.7146 1.9787 -1.0485 -2.1212 -1.7018 -0.0000 1.7146 -1.9787 -1.0485
1.7018 2.1212 -0.0000 -1.9787 1.7146 1.0485 -1.7018 -2.1212 -0.0000 1.9787 -1.7146 1.0485
-0.0000 -0.0000 0.1789 0.4471 -0.4471 0 0.0000 0.0000 -0.1789 0.4471 -0.4471 0
-1.7146 -1.9787 0.4471 4.6176 -2.0009 -0.6602 1.7146 1.9787 -0.4471 -2.3819 -0.2349 -0.6602
1.9787 1.7146 -0.4471 -2.0009 4.6176 -0.6602 -1.9787 -1.7146 0.4471 -0.2349 -2.3819 -0.6602
-1.0485 1.0485 0 -0.6602 -0.6602 7.1843 1.0485 -1.0485 0 0.6602 0.6602 3.3002
-2.1212 -1.7018 0.0000 1.7146 -1.9787 1.0485 2.1212 1.7018 0.0000 -1.7146 1.9787 1.0485
-1.7018 -2.1212 0.0000 1.9787 -1.7146 -1.0485 1.7018 2.1212 0.0000 -1.9787 1.7146 -1.0485
-0.0000 -0.0000 -0.1789 -0.4471 0.4471 0 0.0000 0.0000 0.1789 -0.4471 0.4471 0
1.7146 1.9787 0.4471 -2.3819 -0.2349 0.6602 -1.7146 -1.9787 -0.4471 4.6176 -2.0009 0.6602
-1.9787 -1.7146 -0.4471 -0.2349 -2.3819 0.6602 1.9787 1.7146 0.4471 -2.0009 4.6176 0.6602
-1.0485 1.0485 0 -0.6602 -0.6602 3.3002 1.0485 -1.0485 0 0.6602 0.6602 7.1843
K_g(:,:,2) =
1.0e+003 *
2.1212 -1.7018 0.0000 1.7146 1.9787 -1.0485 -2.1212 1.7018 0.0000 -1.7146 -1.9787 -1.0485
-1.7018 2.1212 -0.0000 -1.9787 -1.7146 -1.0485 1.7018 -2.1212 -0.0000 1.9787 1.7146 -1.0485
0.0000 -0.0000 0.1789 0.4471 0.4471 0 -0.0000 0.0000 -0.1789 0.4471 0.4471 0
1.7146 -1.9787 0.4471 4.6176 2.0009 0.6602 -1.7146 1.9787 -0.4471 -2.3819 0.2349 0.6602
1.9787 -1.7146 0.4471 2.0009 4.6176 -0.6602 -1.9787 1.7146 -0.4471 0.2349 -2.3819 -0.6602
-1.0485 -1.0485 0 0.6602 -0.6602 7.1843 1.0485 1.0485 0 -0.6602 0.6602 3.3002
-2.1212 1.7018 -0.0000 -1.7146 -1.9787 1.0485 2.1212 -1.7018 -0.0000 1.7146 1.9787 1.0485
1.7018 -2.1212 0.0000 1.9787 1.7146 1.0485 -1.7018 2.1212 0.0000 -1.9787 -1.7146 1.0485
0.0000 -0.0000 -0.1789 -0.4471 -0.4471 0 -0.0000 0.0000 0.1789 -0.4471 -0.4471 0
-1.7146 1.9787 0.4471 -2.3819 0.2349 -0.6602 1.7146 -1.9787 -0.4471 4.6176 2.0009 -0.6602
-1.9787 1.7146 0.4471 0.2349 -2.3819 0.6602 1.9787 -1.7146 -0.4471 2.0009 4.6176 0.6602
-1.0485 -1.0485 0 0.6602 -0.6602 3.3002 1.0485 1.0485 0 -0.6602 0.6602 7.1843
K_g(:,:,3) =
1.0e+003 *
2.1212 -1.7018 0.0000 1.7146 1.9787 -1.0485 -2.1212 1.7018 0.0000 -1.7146 -1.9787 -1.0485
-1.7018 2.1212 -0.0000 -1.9787 -1.7146 -1.0485 1.7018 -2.1212 -0.0000 1.9787 1.7146 -1.0485
0.0000 -0.0000 0.1789 0.4471 0.4471 0 -0.0000 0.0000 -0.1789 0.4471 0.4471 0
1.7146 -1.9787 0.4471 4.6176 2.0009 0.6602 -1.7146 1.9787 -0.4471 -2.3819 0.2349 0.6602
1.9787 -1.7146 0.4471 2.0009 4.6176 -0.6602 -1.9787 1.7146 -0.4471 0.2349 -2.3819 -0.6602
-1.0485 -1.0485 0 0.6602 -0.6602 7.1843 1.0485 1.0485 0 -0.6602 0.6602 3.3002
-2.1212 1.7018 -0.0000 -1.7146 -1.9787 1.0485 2.1212 -1.7018 -0.0000 1.7146 1.9787 1.0485
1.7018 -2.1212 0.0000 1.9787 1.7146 1.0485 -1.7018 2.1212 0.0000 -1.9787 -1.7146 1.0485
0.0000 -0.0000 -0.1789 -0.4471 -0.4471 0 -0.0000 0.0000 0.1789 -0.4471 -0.4471 0
-1.7146 1.9787 0.4471 -2.3819 0.2349 -0.6602 1.7146 -1.9787 -0.4471 4.6176 2.0009 -0.6602
-1.9787 1.7146 0.4471 0.2349 -2.3819 0.6602 1.9787 -1.7146 -0.4471 2.0009 4.6176 0.6602
-1.0485 -1.0485 0 0.6602 -0.6602 3.3002 1.0485 1.0485 0 -0.6602 0.6602 7.1843
K_g(:,:,4) =
1.0e+003 *
2.1212 1.7018 -0.0000 -1.7146 1.9787 -1.0485 -2.1212 -1.7018 -0.0000 1.7146 -1.9787 -1.0485
1.7018 2.1212 -0.0000 -1.9787 1.7146 1.0485 -1.7018 -2.1212 -0.0000 1.9787 -1.7146 1.0485
-0.0000 -0.0000 0.1789 0.4471 -0.4471 0 0.0000 0.0000 -0.1789 0.4471 -0.4471 0
-1.7146 -1.9787 0.4471 4.6176 -2.0009 -0.6602 1.7146 1.9787 -0.4471 -2.3819 -0.2349 -0.6602
1.9787 1.7146 -0.4471 -2.0009 4.6176 -0.6602 -1.9787 -1.7146 0.4471 -0.2349 -2.3819 -0.6602
-1.0485 1.0485 0 -0.6602 -0.6602 7.1843 1.0485 -1.0485 0 0.6602 0.6602 3.3002
-2.1212 -1.7018 0.0000 1.7146 -1.9787 1.0485 2.1212 1.7018 0.0000 -1.7146 1.9787 1.0485
-1.7018 -2.1212 0.0000 1.9787 -1.7146 -1.0485 1.7018 2.1212 0.0000 -1.9787 1.7146 -1.0485
-0.0000 -0.0000 -0.1789 -0.4471 0.4471 0 0.0000 0.0000 0.1789 -0.4471 0.4471 0
1.7146 1.9787 0.4471 -2.3819 -0.2349 0.6602 -1.7146 -1.9787 -0.4471 4.6176 -2.0009 0.6602
-1.9787 -1.7146 -0.4471 -0.2349 -2.3819 0.6602 1.9787 1.7146 0.4471 -2.0009 4.6176 0.6602
-1.0485 1.0485 0 -0.6602 -0.6602 3.3002 1.0485 -1.0485 0 0.6602 0.6602 7.1843
K_g(:,:,5) =
1.0e+003 *
2.1212 1.7018 -0.0000 -1.7146 1.9787 -1.0485 -2.1212 -1.7018 -0.0000 1.7146 -1.9787 -1.0485
1.7018 2.1212 -0.0000 -1.9787 1.7146 1.0485 -1.7018 -2.1212 -0.0000 1.9787 -1.7146 1.0485
-0.0000 -0.0000 0.1789 0.4471 -0.4471 0 0.0000 0.0000 -0.1789 0.4471 -0.4471 0
-1.7146 -1.9787 0.4471 4.6176 -2.0009 -0.6602 1.7146 1.9787 -0.4471 -2.3819 -0.2349 -0.6602
1.9787 1.7146 -0.4471 -2.0009 4.6176 -0.6602 -1.9787 -1.7146 0.4471 -0.2349 -2.3819 -0.6602
-1.0485 1.0485 0 -0.6602 -0.6602 7.1843 1.0485 -1.0485 0 0.6602 0.6602 3.3002
-2.1212 -1.7018 0.0000 1.7146 -1.9787 1.0485 2.1212 1.7018 0.0000 -1.7146 1.9787 1.0485
-1.7018 -2.1212 0.0000 1.9787 -1.7146 -1.0485 1.7018 2.1212 0.0000 -1.9787 1.7146 -1.0485
-0.0000 -0.0000 -0.1789 -0.4471 0.4471 0 0.0000 0.0000 0.1789 -0.4471 0.4471 0
1.7146 1.9787 0.4471 -2.3819 -0.2349 0.6602 -1.7146 -1.9787 -0.4471 4.6176 -2.0009 0.6602
-1.9787 -1.7146 -0.4471 -0.2349 -2.3819 0.6602 1.9787 1.7146 0.4471 -2.0009 4.6176 0.6602
-1.0485 1.0485 0 -0.6602 -0.6602 3.3002 1.0485 -1.0485 0 0.6602 0.6602 7.1843
K_g(:,:,6) =
1.0e+003 *
2.1212 -1.7018 0.0000 1.7146 1.9787 -1.0485 -2.1212 1.7018 0.0000 -1.7146 -1.9787 -1.0485
-1.7018 2.1212 -0.0000 -1.9787 -1.7146 -1.0485 1.7018 -2.1212 -0.0000 1.9787 1.7146 -1.0485
0.0000 -0.0000 0.1789 0.4471 0.4471 0 -0.0000 0.0000 -0.1789 0.4471 0.4471 0
1.7146 -1.9787 0.4471 4.6176 2.0009 0.6602 -1.7146 1.9787 -0.4471 -2.3819 0.2349 0.6602
1.9787 -1.7146 0.4471 2.0009 4.6176 -0.6602 -1.9787 1.7146 -0.4471 0.2349 -2.3819 -0.6602
-1.0485 -1.0485 0 0.6602 -0.6602 7.1843 1.0485 1.0485 0 -0.6602 0.6602 3.3002
-2.1212 1.7018 -0.0000 -1.7146 -1.9787 1.0485 2.1212 -1.7018 -0.0000 1.7146 1.9787 1.0485
1.7018 -2.1212 0.0000 1.9787 1.7146 1.0485 -1.7018 2.1212 0.0000 -1.9787 -1.7146 1.0485
0.0000 -0.0000 -0.1789 -0.4471 -0.4471 0 -0.0000 0.0000 0.1789 -0.4471 -0.4471 0
-1.7146 1.9787 0.4471 -2.3819 0.2349 -0.6602 1.7146 -1.9787 -0.4471 4.6176 2.0009 -0.6602
-1.9787 1.7146 0.4471 0.2349 -2.3819 0.6602 1.9787 -1.7146 -0.4471 2.0009 4.6176 0.6602
-1.0485 -1.0485 0 0.6602 -0.6602 3.3002 1.0485 1.0485 0 -0.6602 0.6602 7.1843
K_g(:,:,7) =
1.0e+003 *
2.1212 -1.7018 0.0000 1.7146 1.9787 -1.0485 -2.1212 1.7018 0.0000 -1.7146 -1.9787 -1.0485
-1.7018 2.1212 -0.0000 -1.9787 -1.7146 -1.0485 1.7018 -2.1212 -0.0000 1.9787 1.7146 -1.0485
0.0000 -0.0000 0.1789 0.4471 0.4471 0 -0.0000 0.0000 -0.1789 0.4471 0.4471 0
1.7146 -1.9787 0.4471 4.6176 2.0009 0.6602 -1.7146 1.9787 -0.4471 -2.3819 0.2349 0.6602
1.9787 -1.7146 0.4471 2.0009 4.6176 -0.6602 -1.9787 1.7146 -0.4471 0.2349 -2.3819 -0.6602
-1.0485 -1.0485 0 0.6602 -0.6602 7.1843 1.0485 1.0485 0 -0.6602 0.6602 3.3002
-2.1212 1.7018 -0.0000 -1.7146 -1.9787 1.0485 2.1212 -1.7018 -0.0000 1.7146 1.9787 1.0485
1.7018 -2.1212 0.0000 1.9787 1.7146 1.0485 -1.7018 2.1212 0.0000 -1.9787 -1.7146 1.0485
0.0000 -0.0000 -0.1789 -0.4471 -0.4471 0 -0.0000 0.0000 0.1789 -0.4471 -0.4471 0
-1.7146 1.9787 0.4471 -2.3819 0.2349 -0.6602 1.7146 -1.9787 -0.4471 4.6176 2.0009 -0.6602
-1.9787 1.7146 0.4471 0.2349 -2.3819 0.6602 1.9787 -1.7146 -0.4471 2.0009 4.6176 0.6602
-1.0485 -1.0485 0 0.6602 -0.6602 3.3002 1.0485 1.0485 0 -0.6602 0.6602 7.1843
K_g(:,:,8) =
1.0e+003 *
2.1212 1.7018 -0.0000 -1.7146 1.9787 -1.0485 -2.1212 -1.7018 -0.0000 1.7146 -1.9787 -1.0485
1.7018 2.1212 -0.0000 -1.9787 1.7146 1.0485 -1.7018 -2.1212 -0.0000 1.9787 -1.7146 1.0485
-0.0000 -0.0000 0.1789 0.4471 -0.4471 0 0.0000 0.0000 -0.1789 0.4471 -0.4471 0
-1.7146 -1.9787 0.4471 4.6176 -2.0009 -0.6602 1.7146 1.9787 -0.4471 -2.3819 -0.2349 -0.6602
1.9787 1.7146 -0.4471 -2.0009 4.6176 -0.6602 -1.9787 -1.7146 0.4471 -0.2349 -2.3819 -0.6602
-1.0485 1.0485 0 -0.6602 -0.6602 7.1843 1.0485 -1.0485 0 0.6602 0.6602 3.3002
-2.1212 -1.7018 0.0000 1.7146 -1.9787 1.0485 2.1212 1.7018 0.0000 -1.7146 1.9787 1.0485
-1.7018 -2.1212 0.0000 1.9787 -1.7146 -1.0485 1.7018 2.1212 0.0000 -1.9787 1.7146 -1.0485
-0.0000 -0.0000 -0.1789 -0.4471 0.4471 0 0.0000 0.0000 0.1789 -0.4471 0.4471 0
1.7146 1.9787 0.4471 -2.3819 -0.2349 0.6602 -1.7146 -1.9787 -0.4471 4.6176 -2.0009 0.6602
-1.9787 -1.7146 -0.4471 -0.2349 -2.3819 0.6602 1.9787 1.7146 0.4471 -2.0009 4.6176 0.6602
-1.0485 1.0485 0 -0.6602 -0.6602 3.3002 1.0485 -1.0485 0 0.6602 0.6602 7.1843
K_g(:,:,9) =
1.0e+003 *
1.6795 1.2601 -1.1043 -1.7017 1.2015 -1.0045 -1.6795 -1.2601 1.1043 1.2830 -2.2429 -0.8206
1.2601 1.6795 -1.1043 -1.2015 1.7017 1.0045 -1.2601 -1.6795 1.1043 2.2429 -1.2830 0.8206
-1.1043 -1.1043 1.0623 1.6750 -1.6750 0 1.1043 1.1043 -1.0623 -0.8966 0.8966 0
-1.7017 -1.2015 1.6750 4.7711 -1.8473 1.0440 1.7017 1.2015 -1.6750 -1.6647 0.4823 1.1328
1.2015 1.7017 -1.6750 -1.8473 4.7711 1.0440 -1.2015 -1.7017 1.6750 0.4823 -1.6647 1.1328
-1.0045 1.0045 0 1.0440 1.0440 6.8773 1.0045 -1.0045 0 2.4532 2.4532 1.8658
-1.6795 -1.2601 1.1043 1.7017 -1.2015 1.0045 1.6795 1.2601 -1.1043 -1.2830 2.2429 0.8206
-1.2601 -1.6795 1.1043 1.2015 -1.7017 -1.0045 1.2601 1.6795 -1.1043 -2.2429 1.2830 -0.8206
1.1043 1.1043 -1.0623 -1.6750 1.6750 0 -1.1043 -1.1043 1.0623 0.8966 -0.8966 0
1.2830 2.2429 -0.8966 -1.6647 0.4823 2.4532 -1.2830 -2.2429 0.8966 5.5714 -1.0471 1.7242
-2.2429 -1.2830 0.8966 0.4823 -1.6647 2.4532 2.2429 1.2830 -0.8966 -1.0471 5.5714 1.7242
-0.8206 0.8206 0 1.1328 1.1328 1.8658 0.8206 -0.8206 0 1.7242 1.7242 5.2768
K_g(:,:,10) =
1.0e+003 *
2.1212 -1.7018 0.0000 1.7146 1.9787 -1.0485 -2.1212 1.7018 0.0000 -1.7146 -1.9787 -1.0485
-1.7018 2.1212 -0.0000 -1.9787 -1.7146 -1.0485 1.7018 -2.1212 -0.0000 1.9787 1.7146 -1.0485
0.0000 -0.0000 0.1789 0.4471 0.4471 0 -0.0000 0.0000 -0.1789 0.4471 0.4471 0
1.7146 -1.9787 0.4471 4.6176 2.0009 0.6602 -1.7146 1.9787 -0.4471 -2.3819 0.2349 0.6602
1.9787 -1.7146 0.4471 2.0009 4.6176 -0.6602 -1.9787 1.7146 -0.4471 0.2349 -2.3819 -0.6602
-1.0485 -1.0485 0 0.6602 -0.6602 7.1843 1.0485 1.0485 0 -0.6602 0.6602 3.3002
-2.1212 1.7018 -0.0000 -1.7146 -1.9787 1.0485 2.1212 -1.7018 -0.0000 1.7146 1.9787 1.0485
1.7018 -2.1212 0.0000 1.9787 1.7146 1.0485 -1.7018 2.1212 0.0000 -1.9787 -1.7146 1.0485
0.0000 -0.0000 -0.1789 -0.4471 -0.4471 0 -0.0000 0.0000 0.1789 -0.4471 -0.4471 0
-1.7146 1.9787 0.4471 -2.3819 0.2349 -0.6602 1.7146 -1.9787 -0.4471 4.6176 2.0009 -0.6602
-1.9787 1.7146 0.4471 0.2349 -2.3819 0.6602 1.9787 -1.7146 -0.4471 2.0009 4.6176 0.6602
-1.0485 -1.0485 0 0.6602 -0.6602 3.3002 1.0485 1.0485 0 -0.6602 0.6602 7.1843
K_g(:,:,11) =
1.0e+003 *
1.6795 -1.2601 1.1043 1.7017 1.2015 -1.0045 -1.6795 1.2601 -1.1043 -1.2830 -2.2429 -0.8206
-1.2601 1.6795 -1.1043 -1.2015 -1.7017 -1.0045 1.2601 -1.6795 1.1043 2.2429 1.2830 -0.8206
1.1043 -1.1043 1.0623 1.6750 1.6750 0 -1.1043 1.1043 -1.0623 -0.8966 -0.8966 0
1.7017 -1.2015 1.6750 4.7711 1.8473 -1.0440 -1.7017 1.2015 -1.6750 -1.6647 -0.4823 -1.1328
1.2015 -1.7017 1.6750 1.8473 4.7711 1.0440 -1.2015 1.7017 -1.6750 -0.4823 -1.6647 1.1328
-1.0045 -1.0045 0 -1.0440 1.0440 6.8773 1.0045 1.0045 0 -2.4532 2.4532 1.8658
-1.6795 1.2601 -1.1043 -1.7017 -1.2015 1.0045 1.6795 -1.2601 1.1043 1.2830 2.2429 0.8206
1.2601 -1.6795 1.1043 1.2015 1.7017 1.0045 -1.2601 1.6795 -1.1043 -2.2429 -1.2830 0.8206
-1.1043 1.1043 -1.0623 -1.6750 -1.6750 0 1.1043 -1.1043 1.0623 0.8966 0.8966 0
-1.2830 2.2429 -0.8966 -1.6647 -0.4823 -2.4532 1.2830 -2.2429 0.8966 5.5714 1.0471 -1.7242
-2.2429 1.2830 -0.8966 -0.4823 -1.6647 2.4532 2.2429 -1.2830 0.8966 1.0471 5.5714 1.7242
-0.8206 -0.8206 0 -1.1328 1.1328 1.8658 0.8206 0.8206 0 -1.7242 1.7242 5.2768
K_g(:,:,12) =
1.0e+003 *
2.1212 1.7018 -0.0000 -1.7146 1.9787 -1.0485 -2.1212 -1.7018 -0.0000 1.7146 -1.9787 -1.0485
1.7018 2.1212 -0.0000 -1.9787 1.7146 1.0485 -1.7018 -2.1212 -0.0000 1.9787 -1.7146 1.0485
-0.0000 -0.0000 0.1789 0.4471 -0.4471 0 0.0000 0.0000 -0.1789 0.4471 -0.4471 0
-1.7146 -1.9787 0.4471 4.6176 -2.0009 -0.6602 1.7146 1.9787 -0.4471 -2.3819 -0.2349 -0.6602
1.9787 1.7146 -0.4471 -2.0009 4.6176 -0.6602 -1.9787 -1.7146 0.4471 -0.2349 -2.3819 -0.6602
-1.0485 1.0485 0 -0.6602 -0.6602 7.1843 1.0485 -1.0485 0 0.6602 0.6602 3.3002
-2.1212 -1.7018 0.0000 1.7146 -1.9787 1.0485 2.1212 1.7018 0.0000 -1.7146 1.9787 1.0485
-1.7018 -2.1212 0.0000 1.9787 -1.7146 -1.0485 1.7018 2.1212 0.0000 -1.9787 1.7146 -1.0485
-0.0000 -0.0000 -0.1789 -0.4471 0.4471 0 0.0000 0.0000 0.1789 -0.4471 0.4471 0
1.7146 1.9787 0.4471 -2.3819 -0.2349 0.6602 -1.7146 -1.9787 -0.4471 4.6176 -2.0009 0.6602
-1.9787 -1.7146 -0.4471 -0.2349 -2.3819 0.6602 1.9787 1.7146 0.4471 -2.0009 4.6176 0.6602
-1.0485 1.0485 0 -0.6602 -0.6602 3.3002 1.0485 -1.0485 0 0.6602 0.6602 7.1843
K_g(:,:,13) =
1.0e+003 *
2.1212 1.7018 -0.0000 -1.7146 1.9787 -1.0485 -2.1212 -1.7018 -0.0000 1.7146 -1.9787 -1.0485
1.7018 2.1212 -0.0000 -1.9787 1.7146 1.0485 -1.7018 -2.1212 -0.0000 1.9787 -1.7146 1.0485
-0.0000 -0.0000 0.1789 0.4471 -0.4471 0 0.0000 0.0000 -0.1789 0.4471 -0.4471 0
-1.7146 -1.9787 0.4471 4.6176 -2.0009 -0.6602 1.7146 1.9787 -0.4471 -2.3819 -0.2349 -0.6602
1.9787 1.7146 -0.4471 -2.0009 4.6176 -0.6602 -1.9787 -1.7146 0.4471 -0.2349 -2.3819 -0.6602
-1.0485 1.0485 0 -0.6602 -0.6602 7.1843 1.0485 -1.0485 0 0.6602 0.6602 3.3002
-2.1212 -1.7018 0.0000 1.7146 -1.9787 1.0485 2.1212 1.7018 0.0000 -1.7146 1.9787 1.0485
-1.7018 -2.1212 0.0000 1.9787 -1.7146 -1.0485 1.7018 2.1212 0.0000 -1.9787 1.7146 -1.0485
-0.0000 -0.0000 -0.1789 -0.4471 0.4471 0 0.0000 0.0000 0.1789 -0.4471 0.4471 0
1.7146 1.9787 0.4471 -2.3819 -0.2349 0.6602 -1.7146 -1.9787 -0.4471 4.6176 -2.0009 0.6602
-1.9787 -1.7146 -0.4471 -0.2349 -2.3819 0.6602 1.9787 1.7146 0.4471 -2.0009 4.6176 0.6602
-1.0485 1.0485 0 -0.6602 -0.6602 3.3002 1.0485 -1.0485 0 0.6602 0.6602 7.1843
K_g(:,:,14) =
1.0e+003 *
2.1212 -1.7018 0.0000 1.7146 1.9787 -1.0485 -2.1212 1.7018 0.0000 -1.7146 -1.9787 -1.0485
-1.7018 2.1212 -0.0000 -1.9787 -1.7146 -1.0485 1.7018 -2.1212 -0.0000 1.9787 1.7146 -1.0485
0.0000 -0.0000 0.1789 0.4471 0.4471 0 -0.0000 0.0000 -0.1789 0.4471 0.4471 0
1.7146 -1.9787 0.4471 4.6176 2.0009 0.6602 -1.7146 1.9787 -0.4471 -2.3819 0.2349 0.6602
1.9787 -1.7146 0.4471 2.0009 4.6176 -0.6602 -1.9787 1.7146 -0.4471 0.2349 -2.3819 -0.6602
-1.0485 -1.0485 0 0.6602 -0.6602 7.1843 1.0485 1.0485 0 -0.6602 0.6602 3.3002
-2.1212 1.7018 -0.0000 -1.7146 -1.9787 1.0485 2.1212 -1.7018 -0.0000 1.7146 1.9787 1.0485
1.7018 -2.1212 0.0000 1.9787 1.7146 1.0485 -1.7018 2.1212 0.0000 -1.9787 -1.7146 1.0485
0.0000 -0.0000 -0.1789 -0.4471 -0.4471 0 -0.0000 0.0000 0.1789 -0.4471 -0.4471 0
-1.7146 1.9787 0.4471 -2.3819 0.2349 -0.6602 1.7146 -1.9787 -0.4471 4.6176 2.0009 -0.6602
-1.9787 1.7146 0.4471 0.2349 -2.3819 0.6602 1.9787 -1.7146 -0.4471 2.0009 4.6176 0.6602
-1.0485 -1.0485 0 0.6602 -0.6602 3.3002 1.0485 1.0485 0 -0.6602 0.6602 7.1843
K_g(:,:,15) =
1.0e+003 *
2.1212 -1.7018 0.0000 1.7146 1.9787 -1.0485 -2.1212 1.7018 0.0000 -1.7146 -1.9787 -1.0485
-1.7018 2.1212 -0.0000 -1.9787 -1.7146 -1.0485 1.7018 -2.1212 -0.0000 1.9787 1.7146 -1.0485
0.0000 -0.0000 0.1789 0.4471 0.4471 0 -0.0000 0.0000 -0.1789 0.4471 0.4471 0
1.7146 -1.9787 0.4471 4.6176 2.0009 0.6602 -1.7146 1.9787 -0.4471 -2.3819 0.2349 0.6602
1.9787 -1.7146 0.4471 2.0009 4.6176 -0.6602 -1.9787 1.7146 -0.4471 0.2349 -2.3819 -0.6602
-1.0485 -1.0485 0 0.6602 -0.6602 7.1843 1.0485 1.0485 0 -0.6602 0.6602 3.3002
-2.1212 1.7018 -0.0000 -1.7146 -1.9787 1.0485 2.1212 -1.7018 -0.0000 1.7146 1.9787 1.0485
1.7018 -2.1212 0.0000 1.9787 1.7146 1.0485 -1.7018 2.1212 0.0000 -1.9787 -1.7146 1.0485
0.0000 -0.0000 -0.1789 -0.4471 -0.4471 0 -0.0000 0.0000 0.1789 -0.4471 -0.4471 0
-1.7146 1.9787 0.4471 -2.3819 0.2349 -0.6602 1.7146 -1.9787 -0.4471 4.6176 2.0009 -0.6602
-1.9787 1.7146 0.4471 0.2349 -2.3819 0.6602 1.9787 -1.7146 -0.4471 2.0009 4.6176 0.6602
-1.0485 -1.0485 0 0.6602 -0.6602 3.3002 1.0485 1.0485 0 -0.6602 0.6602 7.1843
K_g(:,:,16) =
1.0e+003 *
2.1212 1.7018 -0.0000 -1.7146 1.9787 -1.0485 -2.1212 -1.7018 -0.0000 1.7146 -1.9787 -1.0485
1.7018 2.1212 -0.0000 -1.9787 1.7146 1.0485 -1.7018 -2.1212 -0.0000 1.9787 -1.7146 1.0485
-0.0000 -0.0000 0.1789 0.4471 -0.4471 0 0.0000 0.0000 -0.1789 0.4471 -0.4471 0
-1.7146 -1.9787 0.4471 4.6176 -2.0009 -0.6602 1.7146 1.9787 -0.4471 -2.3819 -0.2349 -0.6602
1.9787 1.7146 -0.4471 -2.0009 4.6176 -0.6602 -1.9787 -1.7146 0.4471 -0.2349 -2.3819 -0.6602
-1.0485 1.0485 0 -0.6602 -0.6602 7.1843 1.0485 -1.0485 0 0.6602 0.6602 3.3002
-2.1212 -1.7018 0.0000 1.7146 -1.9787 1.0485 2.1212 1.7018 0.0000 -1.7146 1.9787 1.0485
-1.7018 -2.1212 0.0000 1.9787 -1.7146 -1.0485 1.7018 2.1212 0.0000 -1.9787 1.7146 -1.0485
-0.0000 -0.0000 -0.1789 -0.4471 0.4471 0 0.0000 0.0000 0.1789 -0.4471 0.4471 0
1.7146 1.9787 0.4471 -2.3819 -0.2349 0.6602 -1.7146 -1.9787 -0.4471 4.6176 -2.0009 0.6602
-1.9787 -1.7146 -0.4471 -0.2349 -2.3819 0.6602 1.9787 1.7146 0.4471 -2.0009 4.6176 0.6602
-1.0485 1.0485 0 -0.6602 -0.6602 3.3002 1.0485 -1.0485 0 0.6602 0.6602 7.1843
K_g(:,:,17) =
1.0e+003 *
0.4194 0 0 0 1.4866 -0.2416 -0.4194 0 0 0 1.3117 -0.7400
0 0.5781 1.1383 -1.8192 0 0 0 -0.5781 -1.1383 0.6258 0 0
0 1.1383 3.4238 -3.2757 0 0 0 -1.1383 -3.4238 3.6944 0 0
0 -1.8192 -3.2757 6.6185 0 0 0 1.8192 3.2757 -2.1470 0 0
1.4866 0 0 0 7.2671 -0.6976 -1.4866 0 0 0 2.6519 -2.7818
-0.2416 0 0 0 -0.6976 2.5339 0.2416 0 0 0 -0.9145 -1.9685
-0.4194 0 0 0 -1.4866 0.2416 0.4194 0 0 0 -1.3117 0.7400
0 -0.5781 -1.1383 1.8192 0 0 0 0.5781 1.1383 -0.6258 0 0
0 -1.1383 -3.4238 3.2757 0 0 0 1.1383 3.4238 -3.6944 0 0
0 0.6258 3.6944 -2.1470 0 0 0 -0.6258 -3.6944 6.6185 0 0
1.3117 0 0 0 2.6519 -0.9145 -1.3117 0 0 0 6.1006 -2.1557
-0.7400 0 0 0 -2.7818 -1.9685 0.7400 0 0 0 -2.1557 3.7005
K_g(:,:,18) =
1.0e+003 *
3.8230 0 0.0000 0 3.6932 0 -3.8230 0 0.0000 0 -3.6932 0
0 0.4194 0 -0.2641 0 -1.4827 0 -0.4194 0 0.2641 0 -1.4827
0.0000 0 0.1789 0 0.6324 0 -0.0000 0 -0.1789 0 0.6324 0
0 -0.2641 0 2.6168 0 0.9337 0 0.2641 0 -2.6168 0 0.9337
3.6932 0 0.6324 0 6.6185 0 -3.6932 0 -0.6324 0 -2.1470 0
0 -1.4827 0 0.9337 0 7.1843 0 1.4827 0 -0.9337 0 3.3002
-3.8230 0 -0.0000 0 -3.6932 0 3.8230 0 -0.0000 0 3.6932 0
0 -0.4194 0 0.2641 0 1.4827 0 0.4194 0 -0.2641 0 1.4827
0.0000 0 -0.1789 0 -0.6324 0 -0.0000 0 0.1789 0 -0.6324 0
0 0.2641 0 -2.6168 0 -0.9337 0 -0.2641 0 2.6168 0 -0.9337
-3.6932 0 0.6324 0 -2.1470 0 3.6932 0 -0.6324 0 6.6185 0
0 -1.4827 0 0.9337 0 3.3002 0 1.4827 0 -0.9337 0 7.1843
K_g(:,:,19) =
1.0e+003 *
0.4194 0 0 0 0.2641 1.4827 -0.4194 0 0 0 -0.2641 1.4827
0 3.8230 0.0000 -3.6932 0 0 0 -3.8230 0.0000 3.6932 0 0
0 0.0000 0.1789 -0.6324 0 0 0 -0.0000 -0.1789 -0.6324 0 0
0 -3.6932 -0.6324 6.6185 0 0 0 3.6932 0.6324 -2.1470 0 0
0.2641 0 0 0 2.6168 0.9337 -0.2641 0 0 0 -2.6168 0.9337
1.4827 0 0 0 0.9337 7.1843 -1.4827 0 0 0 -0.9337 3.3002
-0.4194 0 0 0 -0.2641 -1.4827 0.4194 0 0 0 0.2641 -1.4827
0 -3.8230 -0.0000 3.6932 0 0 0 3.8230 -0.0000 -3.6932 0 0
0 0.0000 -0.1789 0.6324 0 0 0 -0.0000 0.1789 0.6324 0 0
0 3.6932 -0.6324 -2.1470 0 0 0 -3.6932 0.6324 6.6185 0 0
-0.2641 0 0 0 -2.6168 -0.9337 0.2641 0 0 0 2.6168 -0.9337
1.4827 0 0 0 0.9337 3.3002 -1.4827 0 0 0 -0.9337 7.1843
K_g(:,:,20) =
1.0e+003 *
3.8230 0 -0.0000 0 3.6932 0 -3.8230 0 -0.0000 0 -3.6932 0
0 0.4194 0 -0.2641 0 1.4827 0 -0.4194 0 0.2641 0 1.4827
-0.0000 0 0.1789 0 -0.6324 0 0.0000 0 -0.1789 0 -0.6324 0
0 -0.2641 0 2.6168 0 -0.9337 0 0.2641 0 -2.6168 0 -0.9337
3.6932 0 -0.6324 0 6.6185 0 -3.6932 0 0.6324 0 -2.1470 0
0 1.4827 0 -0.9337 0 7.1843 0 -1.4827 0 0.9337 0 3.3002
-3.8230 0 0.0000 0 -3.6932 0 3.8230 0 0.0000 0 3.6932 0
0 -0.4194 0 0.2641 0 -1.4827 0 0.4194 0 -0.2641 0 -1.4827
-0.0000 0 -0.1789 0 0.6324 0 0.0000 0 0.1789 0 0.6324 0
0 0.2641 0 -2.6168 0 0.9337 0 -0.2641 0 2.6168 0 0.9337
-3.6932 0 -0.6324 0 -2.1470 0 3.6932 0 0.6324 0 6.6185 0
0 1.4827 0 -0.9337 0 3.3002 0 -1.4827 0 0.9337 0 7.1843
K_g(:,:,21) =
1.0e+003 *
0.4194 0 0 0 0.8803 -1.2220 -0.4194 0 0 0 0.4043 -1.4508
0 3.1392 1.4228 -3.6027 0 0 0 -3.1392 -1.4228 3.0548 0 0
0 1.4228 0.8627 -1.0299 0 0 0 -1.4228 -0.8627 2.1698 0 0
0 -3.6027 -1.0299 6.6185 0 0 0 3.6027 1.0299 -2.1470 0 0
0.8803 0 0 0 4.2030 -2.3666 -0.8803 0 0 0 -1.5064 -3.2439
-1.2220 0 0 0 -2.3666 5.5981 1.2220 0 0 0 -1.3765 2.1899
-0.4194 0 0 0 -0.8803 1.2220 0.4194 0 0 0 -0.4043 1.4508
0 -3.1392 -1.4228 3.6027 0 0 0 3.1392 1.4228 -3.0548 0 0
0 -1.4228 -0.8627 1.0299 0 0 0 1.4228 0.8627 -2.1698 0 0
0 3.0548 2.1698 -2.1470 0 0 0 -3.0548 -2.1698 6.6185 0 0
0.4043 0 0 0 -1.5064 -1.3765 -0.4043 0 0 0 2.7448 -1.2001
-1.4508 0 0 0 -3.2439 2.1899 1.4508 0 0 0 -1.2001 7.0563
K_g(:,:,22) =
1.0e+003 *
1.5592 0 -1.7677 0 2.7714 0 -1.5592 0 1.7677 0 -1.7746 0
0 0.4194 0 -1.3312 0 -0.7044 0 -0.4194 0 -1.0061 0 -1.1207
-1.7677 0 2.4427 0 -2.5217 0 1.7677 0 -2.4427 0 3.3001 0
0 -1.3312 0 6.3600 0 1.9893 0 1.3312 0 1.0590 0 3.8039
2.7714 0 -2.5217 0 6.6185 0 -2.7714 0 2.5217 0 -2.1470 0
0 -0.7044 0 1.9893 0 3.4411 0 0.7044 0 1.9366 0 -0.3755
-1.5592 0 1.7677 0 -2.7714 0 1.5592 0 -1.7677 0 1.7746 0
0 -0.4194 0 1.3312 0 0.7044 0 0.4194 0 1.0061 0 1.1207
1.7677 0 -2.4427 0 2.5217 0 -1.7677 0 2.4427 0 -3.3001 0
0 -1.0061 0 1.0590 0 1.9366 0 1.0061 0 4.5483 0 2.4420
-1.7746 0 3.3001 0 -2.1470 0 1.7746 0 -3.3001 0 6.6185 0
0 -1.1207 0 3.8039 0 -0.3755 0 1.1207 0 2.4420 0 5.2527
K_g(:,:,23) =
1.0e+003 *
0.4194 0 0 0 1.3312 0.7044 -0.4194 0 0 0 1.0061 1.1207
0 1.5592 -1.7677 -2.7714 0 0 0 -1.5592 1.7677 1.7746 0 0
0 -1.7677 2.4427 2.5217 0 0 0 1.7677 -2.4427 -3.3001 0 0
0 -2.7714 2.5217 6.6185 0 0 0 2.7714 -2.5217 -2.1470 0 0
1.3312 0 0 0 6.3600 1.9893 -1.3312 0 0 0 1.0590 3.8039
0.7044 0 0 0 1.9893 3.4411 -0.7044 0 0 0 1.9366 -0.3755
-0.4194 0 0 0 -1.3312 -0.7044 0.4194 0 0 0 -1.0061 -1.1207
0 -1.5592 1.7677 2.7714 0 0 0 1.5592 -1.7677 -1.7746 0 0
0 1.7677 -2.4427 -2.5217 0 0 0 -1.7677 2.4427 3.3001 0 0
0 1.7746 -3.3001 -2.1470 0 0 0 -1.7746 3.3001 6.6185 0 0
1.0061 0 0 0 1.0590 1.9366 -1.0061 0 0 0 4.5483 2.4420
1.1207 0 0 0 3.8039 -0.3755 -1.1207 0 0 0 2.4420 5.2527
K_g(:,:,24) =
1.0e+003 *
1.5592 0 1.7677 0 2.7714 0 -1.5592 0 -1.7677 0 -1.7746 0
0 0.4194 0 -1.3312 0 0.7044 0 -0.4194 0 -1.0061 0 1.1207
1.7677 0 2.4427 0 2.5217 0 -1.7677 0 -2.4427 0 -3.3001 0
0 -1.3312 0 6.3600 0 -1.9893 0 1.3312 0 1.0590 0 -3.8039
2.7714 0 2.5217 0 6.6185 0 -2.7714 0 -2.5217 0 -2.1470 0
0 0.7044 0 -1.9893 0 3.4411 0 -0.7044 0 -1.9366 0 -0.3755
-1.5592 0 -1.7677 0 -2.7714 0 1.5592 0 1.7677 0 1.7746 0
0 -0.4194 0 1.3312 0 -0.7044 0 0.4194 0 1.0061 0 -1.1207
-1.7677 0 -2.4427 0 -2.5217 0 1.7677 0 2.4427 0 3.3001 0
0 -1.0061 0 1.0590 0 -1.9366 0 1.0061 0 4.5483 0 -2.4420
-1.7746 0 -3.3001 0 -2.1470 0 1.7746 0 3.3001 0 6.6185 0
0 1.1207 0 -3.8039 0 -0.3755 0 -1.1207 0 -2.4420 0 5.2527
K_g(:,:,25) =
1.0e+003 *
0.4194 0 0 0 1.4443 -0.4269 -0.4194 0 0 0 1.2081 -0.8993
0 0.9077 1.4577 -2.2173 0 0 0 -0.9077 -1.4577 1.0861 0 0
0 1.4577 3.0942 -3.0205 0 0 0 -1.4577 -3.0942 3.5861 0 0
0 -2.2173 -3.0205 6.6185 0 0 0 2.2173 3.0205 -2.1470 0 0
1.4443 0 0 0 7.0177 -1.2668 -1.4443 0 0 0 2.1168 -3.3005
-0.4269 0 0 0 -1.2668 2.7833 0.4269 0 0 0 -1.4331 -1.4334
-0.4194 0 0 0 -1.4443 0.4269 0.4194 0 0 0 -1.2081 0.8993
0 -0.9077 -1.4577 2.2173 0 0 0 0.9077 1.4577 -1.0861 0 0
0 -1.4577 -3.0942 3.0205 0 0 0 1.4577 3.0942 -3.5861 0 0
0 1.0861 3.5861 -2.1470 0 0 0 -1.0861 -3.5861 6.6185 0 0
1.2081 0 0 0 2.1168 -1.4331 -1.2081 0 0 0 5.5239 -2.3872
-0.8993 0 0 0 -3.3005 -1.4334 0.8993 0 0 0 -2.3872 4.2772
K_g(:,:,26) =
1.0e+003 *
0.9077 0 -1.4577 0 2.2173 0 -0.9077 0 1.4577 0 -1.0861 0
0 0.4194 0 -1.4443 0 -0.4269 0 -0.4194 0 -1.2081 0 -0.8993
-1.4577 0 3.0942 0 -3.0205 0 1.4577 0 -3.0942 0 3.5861 0
0 -1.4443 0 7.0177 0 1.2668 0 1.4443 0 2.1168 0 3.3005
2.2173 0 -3.0205 0 6.6185 0 -2.2173 0 3.0205 0 -2.1470 0
0 -0.4269 0 1.2668 0 2.7833 0 0.4269 0 1.4331 0 -1.4334
-0.9077 0 1.4577 0 -2.2173 0 0.9077 0 -1.4577 0 1.0861 0
0 -0.4194 0 1.4443 0 0.4269 0 0.4194 0 1.2081 0 0.8993
1.4577 0 -3.0942 0 3.0205 0 -1.4577 0 3.0942 0 -3.5861 0
0 -1.2081 0 2.1168 0 1.4331 0 1.2081 0 5.5239 0 2.3872
-1.0861 0 3.5861 0 -2.1470 0 1.0861 0 -3.5861 0 6.6185 0
0 -0.8993 0 3.3005 0 -1.4334 0 0.8993 0 2.3872 0 4.2772
K_g(:,:,27) =
1.0e+003 *
0.4194 0 0 0 1.4443 0.4269 -0.4194 0 0 0 1.2081 0.8993
0 0.9077 -1.4577 -2.2173 0 0 0 -0.9077 1.4577 1.0861 0 0
0 -1.4577 3.0942 3.0205 0 0 0 1.4577 -3.0942 -3.5861 0 0
0 -2.2173 3.0205 6.6185 0 0 0 2.2173 -3.0205 -2.1470 0 0
1.4443 0 0 0 7.0177 1.2668 -1.4443 0 0 0 2.1168 3.3005
0.4269 0 0 0 1.2668 2.7833 -0.4269 0 0 0 1.4331 -1.4334
-0.4194 0 0 0 -1.4443 -0.4269 0.4194 0 0 0 -1.2081 -0.8993
0 -0.9077 1.4577 2.2173 0 0 0 0.9077 -1.4577 -1.0861 0 0
0 1.4577 -3.0942 -3.0205 0 0 0 -1.4577 3.0942 3.5861 0 0
0 1.0861 -3.5861 -2.1470 0 0 0 -1.0861 3.5861 6.6185 0 0
1.2081 0 0 0 2.1168 1.4331 -1.2081 0 0 0 5.5239 2.3872
0.8993 0 0 0 3.3005 -1.4334 -0.8993 0 0 0 2.3872 4.2772
K_g(:,:,28) =
1.0e+003 *
0.9077 0 1.4577 0 2.2173 0 -0.9077 0 -1.4577 0 -1.0861 0
0 0.4194 0 -1.4443 0 0.4269 0 -0.4194 0 -1.2081 0 0.8993
1.4577 0 3.0942 0 3.0205 0 -1.4577 0 -3.0942 0 -3.5861 0
0 -1.4443 0 7.0177 0 -1.2668 0 1.4443 0 2.1168 0 -3.3005
2.2173 0 3.0205 0 6.6185 0 -2.2173 0 -3.0205 0 -2.1470 0
0 0.4269 0 -1.2668 0 2.7833 0 -0.4269 0 -1.4331 0 -1.4334
-0.9077 0 -1.4577 0 -2.2173 0 0.9077 0 1.4577 0 1.0861 0
0 -0.4194 0 1.4443 0 -0.4269 0 0.4194 0 1.2081 0 -0.8993
-1.4577 0 -3.0942 0 -3.0205 0 1.4577 0 3.0942 0 3.5861 0
0 -1.2081 0 2.1168 0 -1.4331 0 1.2081 0 5.5239 0 -2.3872
-1.0861 0 -3.5861 0 -2.1470 0 1.0861 0 3.5861 0 6.6185 0
0 0.8993 0 -3.3005 0 -1.4334 0 -0.8993 0 -2.3872 0 4.2772
K_g(:,:,29) =
1.0e+003 *
0.4194 0 0 0 1.1331 -0.9921 -0.4194 0 0 0 0.7209 -1.3223
0 2.3987 1.7782 -3.2778 0 0 0 -2.3987 -1.7782 2.4871 0 0
0 1.7782 1.6032 -1.8154 0 0 0 -1.7782 -1.6032 2.8025 0 0
0 -3.2778 -1.8154 6.6185 0 0 0 3.2778 1.8154 -2.1470 0 0
1.1331 0 0 0 5.3132 -2.4325 -1.1331 0 0 0 -0.3040 -3.8208
-0.9921 0 0 0 -2.4325 4.4879 0.9921 0 0 0 -1.9535 0.9875
-0.4194 0 0 0 -1.1331 0.9921 0.4194 0 0 0 -0.7209 1.3223
0 -2.3987 -1.7782 3.2778 0 0 0 2.3987 1.7782 -2.4871 0 0
0 -1.7782 -1.6032 1.8154 0 0 0 1.7782 1.6032 -2.8025 0 0
0 2.4871 2.8025 -2.1470 0 0 0 -2.4871 -2.8025 6.6185 0 0
0.7209 0 0 0 -0.3040 -1.9535 -0.7209 0 0 0 3.4909 -2.0249
-1.3223 0 0 0 -3.8208 0.9875 1.3223 0 0 0 -2.0249 6.3102
K_g(:,:,30) =
1.0e+003 *
2.3987 0 -1.7782 0 3.2778 0 -2.3987 0 1.7782 0 -2.4871 0
0 0.4194 0 -1.1331 0 -0.9921 0 -0.4194 0 -0.7209 0 -1.3223
-1.7782 0 1.6032 0 -1.8154 0 1.7782 0 -1.6032 0 2.8025 0
0 -1.1331 0 5.3132 0 2.4325 0 1.1331 0 -0.3040 0 3.8208
3.2778 0 -1.8154 0 6.6185 0 -3.2778 0 1.8154 0 -2.1470 0
0 -0.9921 0 2.4325 0 4.4879 0 0.9921 0 1.9535 0 0.9875
-2.3987 0 1.7782 0 -3.2778 0 2.3987 0 -1.7782 0 2.4871 0
0 -0.4194 0 1.1331 0 0.9921 0 0.4194 0 0.7209 0 1.3223
1.7782 0 -1.6032 0 1.8154 0 -1.7782 0 1.6032 0 -2.8025 0
0 -0.7209 0 -0.3040 0 1.9535 0 0.7209 0 3.4909 0 2.0249
-2.4871 0 2.8025 0 -2.1470 0 2.4871 0 -2.8025 0 6.6185 0
0 -1.3223 0 3.8208 0 0.9875 0 1.3223 0 2.0249 0 6.3102
K_g(:,:,31) =
1.0e+003 *
0.4194 0 0 0 1.1331 0.9921 -0.4194 0 0 0 0.7209 1.3223
0 2.3987 -1.7782 -3.2778 0 0 0 -2.3987 1.7782 2.4871 0 0
0 -1.7782 1.6032 1.8154 0 0 0 1.7782 -1.6032 -2.8025 0 0
0 -3.2778 1.8154 6.6185 0 0 0 3.2778 -1.8154 -2.1470 0 0
1.1331 0 0 0 5.3132 2.4325 -1.1331 0 0 0 -0.3040 3.8208
0.9921 0 0 0 2.4325 4.4879 -0.9921 0 0 0 1.9535 0.9875
-0.4194 0 0 0 -1.1331 -0.9921 0.4194 0 0 0 -0.7209 -1.3223
0 -2.3987 1.7782 3.2778 0 0 0 2.3987 -1.7782 -2.4871 0 0
0 1.7782 -1.6032 -1.8154 0 0 0 -1.7782 1.6032 2.8025 0 0
0 2.4871 -2.8025 -2.1470 0 0 0 -2.4871 2.8025 6.6185 0 0
0.7209 0 0 0 -0.3040 1.9535 -0.7209 0 0 0 3.4909 2.0249
1.3223 0 0 0 3.8208 0.9875 -1.3223 0 0 0 2.0249 6.3102
K_g(:,:,32) =
1.0e+003 *
2.3987 0 1.7782 0 3.2778 0 -2.3987 0 -1.7782 0 -2.4871 0
0 0.4194 0 -1.1331 0 0.9921 0 -0.4194 0 -0.7209 0 1.3223
1.7782 0 1.6032 0 1.8154 0 -1.7782 0 -1.6032 0 -2.8025 0
0 -1.1331 0 5.3132 0 -2.4325 0 1.1331 0 -0.3040 0 -3.8208
3.2778 0 1.8154 0 6.6185 0 -3.2778 0 -1.8154 0 -2.1470 0
0 0.9921 0 -2.4325 0 4.4879 0 -0.9921 0 -1.9535 0 0.9875
-2.3987 0 -1.7782 0 -3.2778 0 2.3987 0 1.7782 0 2.4871 0
0 -0.4194 0 1.1331 0 -0.9921 0 0.4194 0 0.7209 0 -1.3223
-1.7782 0 -1.6032 0 -1.8154 0 1.7782 0 1.6032 0 2.8025 0
0 -0.7209 0 -0.3040 0 -1.9535 0 0.7209 0 3.4909 0 -2.0249
-2.4871 0 -2.8025 0 -2.1470 0 2.4871 0 2.8025 0 6.6185 0
0 1.3223 0 -3.8208 0 0.9875 0 -1.3223 0 -2.0249 0 6.3102
T
T(:,:,1) =
0.7071 -0.7071 0 0 0 0 0 0 0 0 0 0
0.5000 0.5000 0.7071 0 0 0 0 0 0 0 0 0
-0.5000 -0.5000 0.7071 0 0 0 0 0 0 0 0 0
0 0 0 0.7071 -0.7071 0 0 0 0 0 0 0
0 0 0 0.5000 0.5000 0.7071 0 0 0 0 0 0
0 0 0 -0.5000 -0.5000 0.7071 0 0 0 0 0 0
0 0 0 0 0 0 0.7071 -0.7071 0 0 0 0
0 0 0 0 0 0 0.5000 0.5000 -0.7071 0 0 0
0 0 0 0 0 0 0.5000 0.5000 0.7071 0 0 0
0 0 0 0 0 0 0 0 0 0.7071 -0.7071 0
0 0 0 0 0 0 0 0 0 0.5000 0.5000 -0.7071
0 0 0 0 0 0 0 0 0 0.5000 0.5000 0.7071
T(:,:,2) =
0.7071 0.7071 0 0 0 0 0 0 0 0 0 0
-0.5000 0.5000 0.7071 0 0 0 0 0 0 0 0 0
0.5000 -0.5000 0.7071 0 0 0 0 0 0 0 0 0
0 0 0 0.7071 0.7071 0 0 0 0 0 0 0
0 0 0 -0.5000 0.5000 0.7071 0 0 0 0 0 0
0 0 0 0.5000 -0.5000 0.7071 0 0 0 0 0 0
0 0 0 0 0 0 0.7071 0.7071 0 0 0 0
0 0 0 0 0 0 -0.5000 0.5000 -0.7071 0 0 0
0 0 0 0 0 0 -0.5000 0.5000 0.7071 0 0 0
0 0 0 0 0 0 0 0 0 0.7071 0.7071 0
0 0 0 0 0 0 0 0 0 -0.5000 0.5000 -0.7071
0 0 0 0 0 0 0 0 0 -0.5000 0.5000 0.7071
T(:,:,3) =
0.7071 0.7071 0 0 0 0 0 0 0 0 0 0
-0.5000 0.5000 0.7071 0 0 0 0 0 0 0 0 0
0.5000 -0.5000 0.7071 0 0 0 0 0 0 0 0 0
0 0 0 0.7071 0.7071 0 0 0 0 0 0 0
0 0 0 -0.5000 0.5000 0.7071 0 0 0 0 0 0
0 0 0 0.5000 -0.5000 0.7071 0 0 0 0 0 0
0 0 0 0 0 0 0.7071 0.7071 0 0 0 0
0 0 0 0 0 0 -0.5000 0.5000 -0.7071 0 0 0
0 0 0 0 0 0 -0.5000 0.5000 0.7071 0 0 0
0 0 0 0 0 0 0 0 0 0.7071 0.7071 0
0 0 0 0 0 0 0 0 0 -0.5000 0.5000 -0.7071
0 0 0 0 0 0 0 0 0 -0.5000 0.5000 0.7071
T(:,:,4) =
0.7071 -0.7071 0 0 0 0 0 0 0 0 0 0
0.5000 0.5000 0.7071 0 0 0 0 0 0 0 0 0
-0.5000 -0.5000 0.7071 0 0 0 0 0 0 0 0 0
0 0 0 0.7071 -0.7071 0 0 0 0 0 0 0
0 0 0 0.5000 0.5000 0.7071 0 0 0 0 0 0
0 0 0 -0.5000 -0.5000 0.7071 0 0 0 0 0 0
0 0 0 0 0 0 0.7071 -0.7071 0 0 0 0
0 0 0 0 0 0 0.5000 0.5000 -0.7071 0 0 0
0 0 0 0 0 0 0.5000 0.5000 0.7071 0 0 0
0 0 0 0 0 0 0 0 0 0.7071 -0.7071 0
0 0 0 0 0 0 0 0 0 0.5000 0.5000 -0.7071
0 0 0 0 0 0 0 0 0 0.5000 0.5000 0.7071
T(:,:,5) =
0.7071 -0.7071 0 0 0 0 0 0 0 0 0 0
0.5000 0.5000 0.7071 0 0 0 0 0 0 0 0 0
-0.5000 -0.5000 0.7071 0 0 0 0 0 0 0 0 0
0 0 0 0.7071 -0.7071 0 0 0 0 0 0 0
0 0 0 0.5000 0.5000 0.7071 0 0 0 0 0 0
0 0 0 -0.5000 -0.5000 0.7071 0 0 0 0 0 0
0 0 0 0 0 0 0.7071 -0.7071 0 0 0 0
0 0 0 0 0 0 0.5000 0.5000 -0.7071 0 0 0
0 0 0 0 0 0 0.5000 0.5000 0.7071 0 0 0
0 0 0 0 0 0 0 0 0 0.7071 -0.7071 0
0 0 0 0 0 0 0 0 0 0.5000 0.5000 -0.7071
0 0 0 0 0 0 0 0 0 0.5000 0.5000 0.7071
T(:,:,6) =
0.7071 0.7071 0 0 0 0 0 0 0 0 0 0
-0.5000 0.5000 0.7071 0 0 0 0 0 0 0 0 0
0.5000 -0.5000 0.7071 0 0 0 0 0 0 0 0 0
0 0 0 0.7071 0.7071 0 0 0 0 0 0 0
0 0 0 -0.5000 0.5000 0.7071 0 0 0 0 0 0
0 0 0 0.5000 -0.5000 0.7071 0 0 0 0 0 0
0 0 0 0 0 0 0.7071 0.7071 0 0 0 0
0 0 0 0 0 0 -0.5000 0.5000 -0.7071 0 0 0
0 0 0 0 0 0 -0.5000 0.5000 0.7071 0 0 0
0 0 0 0 0 0 0 0 0 0.7071 0.7071 0
0 0 0 0 0 0 0 0 0 -0.5000 0.5000 -0.7071
0 0 0 0 0 0 0 0 0 -0.5000 0.5000 0.7071
T(:,:,7) =
0.7071 0.7071 0 0 0 0 0 0 0 0 0 0
-0.5000 0.5000 0.7071 0 0 0 0 0 0 0 0 0
0.5000 -0.5000 0.7071 0 0 0 0 0 0 0 0 0
0 0 0 0.7071 0.7071 0 0 0 0 0 0 0
0 0 0 -0.5000 0.5000 0.7071 0 0 0 0 0 0
0 0 0 0.5000 -0.5000 0.7071 0 0 0 0 0 0
0 0 0 0 0 0 0.7071 0.7071 0 0 0 0
0 0 0 0 0 0 -0.5000 0.5000 -0.7071 0 0 0
0 0 0 0 0 0 -0.5000 0.5000 0.7071 0 0 0
0 0 0 0 0 0 0 0 0 0.7071 0.7071 0
0 0 0 0 0 0 0 0 0 -0.5000 0.5000 -0.7071
0 0 0 0 0 0 0 0 0 -0.5000 0.5000 0.7071
T(:,:,8) =
0.7071 -0.7071 0 0 0 0 0 0 0 0 0 0
0.5000 0.5000 0.7071 0 0 0 0 0 0 0 0 0
-0.5000 -0.5000 0.7071 0 0 0 0 0 0 0 0 0
0 0 0 0.7071 -0.7071 0 0 0 0 0 0 0
0 0 0 0.5000 0.5000 0.7071 0 0 0 0 0 0
0 0 0 -0.5000 -0.5000 0.7071 0 0 0 0 0 0
0 0 0 0 0 0 0.7071 -0.7071 0 0 0 0
0 0 0 0 0 0 0.5000 0.5000 -0.7071 0 0 0
0 0 0 0 0 0 0.5000 0.5000 0.7071 0 0 0
0 0 0 0 0 0 0 0 0 0.7071 -0.7071 0
0 0 0 0 0 0 0 0 0 0.5000 0.5000 -0.7071
0 0 0 0 0 0 0 0 0 0.5000 0.5000 0.7071
T(:,:,9) =
0.7071 -0.7071 0 0 0 0 0 0 0 0 0 0
0.6814 0.6814 0.2673 0 0 0 0 0 0 0 0 0
-0.1890 -0.1890 0.9636 0 0 0 0 0 0 0 0 0
0 0 0 0.7071 -0.7071 0 0 0 0 0 0 0
0 0 0 0.6814 0.6814 0.2673 0 0 0 0 0 0
0 0 0 -0.1890 -0.1890 0.9636 0 0 0 0 0 0
0 0 0 0 0 0 0.7071 -0.7071 0 0 0 0
0 0 0 0 0 0 0.1890 0.1890 -0.9636 0 0 0
0 0 0 0 0 0 0.6814 0.6814 0.2673 0 0 0
0 0 0 0 0 0 0 0 0 0.7071 -0.7071 0
0 0 0 0 0 0 0 0 0 0.1890 0.1890 -0.9636
0 0 0 0 0 0 0 0 0 0.6814 0.6814 0.2673
T(:,:,10) =
0.7071 0.7071 0 0 0 0 0 0 0 0 0 0
-0.5000 0.5000 0.7071 0 0 0 0 0 0 0 0 0
0.5000 -0.5000 0.7071 0 0 0 0 0 0 0 0 0
0 0 0 0.7071 0.7071 0 0 0 0 0 0 0
0 0 0 -0.5000 0.5000 0.7071 0 0 0 0 0 0
0 0 0 0.5000 -0.5000 0.7071 0 0 0 0 0 0
0 0 0 0 0 0 0.7071 0.7071 0 0 0 0
0 0 0 0 0 0 -0.5000 0.5000 -0.7071 0 0 0
0 0 0 0 0 0 -0.5000 0.5000 0.7071 0 0 0
0 0 0 0 0 0 0 0 0 0.7071 0.7071 0
0 0 0 0 0 0 0 0 0 -0.5000 0.5000 -0.7071
0 0 0 0 0 0 0 0 0 -0.5000 0.5000 0.7071
T(:,:,11) =
0.7071 0.7071 0 0 0 0 0 0 0 0 0 0
-0.6814 0.6814 0.2673 0 0 0 0 0 0 0 0 0
0.1890 -0.1890 0.9636 0 0 0 0 0 0 0 0 0
0 0 0 0.7071 0.7071 0 0 0 0 0 0 0
0 0 0 -0.6814 0.6814 0.2673 0 0 0 0 0 0
0 0 0 0.1890 -0.1890 0.9636 0 0 0 0 0 0
0 0 0 0 0 0 0.7071 0.7071 0 0 0 0
0 0 0 0 0 0 -0.1890 0.1890 -0.9636 0 0 0
0 0 0 0 0 0 -0.6814 0.6814 0.2673 0 0 0
0 0 0 0 0 0 0 0 0 0.7071 0.7071 0
0 0 0 0 0 0 0 0 0 -0.1890 0.1890 -0.9636
0 0 0 0 0 0 0 0 0 -0.6814 0.6814 0.2673
T(:,:,12) =
0.7071 -0.7071 0 0 0 0 0 0 0 0 0 0
0.5000 0.5000 0.7071 0 0 0 0 0 0 0 0 0
-0.5000 -0.5000 0.7071 0 0 0 0 0 0 0 0 0
0 0 0 0.7071 -0.7071 0 0 0 0 0 0 0
0 0 0 0.5000 0.5000 0.7071 0 0 0 0 0 0
0 0 0 -0.5000 -0.5000 0.7071 0 0 0 0 0 0
0 0 0 0 0 0 0.7071 -0.7071 0 0 0 0
0 0 0 0 0 0 0.5000 0.5000 -0.7071 0 0 0
0 0 0 0 0 0 0.5000 0.5000 0.7071 0 0 0
0 0 0 0 0 0 0 0 0 0.7071 -0.7071 0
0 0 0 0 0 0 0 0 0 0.5000 0.5000 -0.7071
0 0 0 0 0 0 0 0 0 0.5000 0.5000 0.7071
T(:,:,13) =
0.7071 -0.7071 0 0 0 0 0 0 0 0 0 0
0.5000 0.5000 0.7071 0 0 0 0 0 0 0 0 0
-0.5000 -0.5000 0.7071 0 0 0 0 0 0 0 0 0
0 0 0 0.7071 -0.7071 0 0 0 0 0 0 0
0 0 0 0.5000 0.5000 0.7071 0 0 0 0 0 0
0 0 0 -0.5000 -0.5000 0.7071 0 0 0 0 0 0
0 0 0 0 0 0 0.7071 -0.7071 0 0 0 0
0 0 0 0 0 0 0.5000 0.5000 -0.7071 0 0 0
0 0 0 0 0 0 0.5000 0.5000 0.7071 0 0 0
0 0 0 0 0 0 0 0 0 0.7071 -0.7071 0
0 0 0 0 0 0 0 0 0 0.5000 0.5000 -0.7071
0 0 0 0 0 0 0 0 0 0.5000 0.5000 0.7071
T(:,:,14) =
0.7071 0.7071 0 0 0 0 0 0 0 0 0 0
-0.5000 0.5000 0.7071 0 0 0 0 0 0 0 0 0
0.5000 -0.5000 0.7071 0 0 0 0 0 0 0 0 0
0 0 0 0.7071 0.7071 0 0 0 0 0 0 0
0 0 0 -0.5000 0.5000 0.7071 0 0 0 0 0 0
0 0 0 0.5000 -0.5000 0.7071 0 0 0 0 0 0
0 0 0 0 0 0 0.7071 0.7071 0 0 0 0
0 0 0 0 0 0 -0.5000 0.5000 -0.7071 0 0 0
0 0 0 0 0 0 -0.5000 0.5000 0.7071 0 0 0
0 0 0 0 0 0 0 0 0 0.7071 0.7071 0
0 0 0 0 0 0 0 0 0 -0.5000 0.5000 -0.7071
0 0 0 0 0 0 0 0 0 -0.5000 0.5000 0.7071
T(:,:,15) =
0.7071 0.7071 0 0 0 0 0 0 0 0 0 0
-0.5000 0.5000 0.7071 0 0 0 0 0 0 0 0 0
0.5000 -0.5000 0.7071 0 0 0 0 0 0 0 0 0
0 0 0 0.7071 0.7071 0 0 0 0 0 0 0
0 0 0 -0.5000 0.5000 0.7071 0 0 0 0 0 0
0 0 0 0.5000 -0.5000 0.7071 0 0 0 0 0 0
0 0 0 0 0 0 0.7071 0.7071 0 0 0 0
0 0 0 0 0 0 -0.5000 0.5000 -0.7071 0 0 0
0 0 0 0 0 0 -0.5000 0.5000 0.7071 0 0 0
0 0 0 0 0 0 0 0 0 0.7071 0.7071 0
0 0 0 0 0 0 0 0 0 -0.5000 0.5000 -0.7071
0 0 0 0 0 0 0 0 0 -0.5000 0.5000 0.7071
T(:,:,16) =
0.7071 -0.7071 0 0 0 0 0 0 0 0 0 0
0.5000 0.5000 0.7071 0 0 0 0 0 0 0 0 0
-0.5000 -0.5000 0.7071 0 0 0 0 0 0 0 0 0
0 0 0 0.7071 -0.7071 0 0 0 0 0 0 0
0 0 0 0.5000 0.5000 0.7071 0 0 0 0 0 0
0 0 0 -0.5000 -0.5000 0.7071 0 0 0 0 0 0
0 0 0 0 0 0 0.7071 -0.7071 0 0 0 0
0 0 0 0 0 0 0.5000 0.5000 -0.7071 0 0 0
0 0 0 0 0 0 0.5000 0.5000 0.7071 0 0 0
0 0 0 0 0 0 0 0 0 0.7071 -0.7071 0
0 0 0 0 0 0 0 0 0 0.5000 0.5000 -0.7071
0 0 0 0 0 0 0 0 0 0.5000 0.5000 0.7071
T(:,:,17) =
1.0000 0 0 0 0 0 0 0 0 0 0 0
0 -0.4332 0.9013 0 0 0 0 0 0 0 0 0
0 -0.9013 -0.4332 0 0 0 0 0 0 0 0 0
0 0 0 1.0000 0 0 0 0 0 0 0 0
0 0 0 0 -0.4332 0.9013 0 0 0 0 0 0
0 0 0 0 -0.9013 -0.4332 0 0 0 0 0 0
0 0 0 0 0 0 1.0000 0 0 0 0 0
0 0 0 0 0 0 0 0.9013 0.4332 0 0 0
0 0 0 0 0 0 0 -0.4332 0.9013 0 0 0
0 0 0 0 0 0 0 0 0 1.0000 0 0
0 0 0 0 0 0 0 0 0 0 0.9013 0.4332
0 0 0 0 0 0 0 0 0 0 -0.4332 0.9013
T(:,:,18) =
0 1.0000 0 0 0 0 0 0 0 0 0 0
-0.7071 0 0.7071 0 0 0 0 0 0 0 0 0
0.7071 0 0.7071 0 0 0 0 0 0 0 0 0
0 0 0 0 1.0000 0 0 0 0 0 0 0
0 0 0 -0.7071 0 0.7071 0 0 0 0 0 0
0 0 0 0.7071 0 0.7071 0 0 0 0 0 0
0 0 0 0 0 0 0 1.0000 0 0 0 0
0 0 0 0 0 0 -0.7071 0 -0.7071 0 0 0
0 0 0 0 0 0 -0.7071 0 0.7071 0 0 0
0 0 0 0 0 0 0 0 0 0 1.0000 0
0 0 0 0 0 0 0 0 0 -0.7071 0 -0.7071
0 0 0 0 0 0 0 0 0 -0.7071 0 0.7071
T(:,:,19) =
-1.0000 0 0 0 0 0 0 0 0 0 0 0
0 -0.7071 0.7071 0 0 0 0 0 0 0 0 0
0 0.7071 0.7071 0 0 0 0 0 0 0 0 0
0 0 0 -1.0000 0 0 0 0 0 0 0 0
0 0 0 0 -0.7071 0.7071 0 0 0 0 0 0
0 0 0 0 0.7071 0.7071 0 0 0 0 0 0
0 0 0 0 0 0 -1.0000 0 0 0 0 0
0 0 0 0 0 0 0 -0.7071 -0.7071 0 0 0
0 0 0 0 0 0 0 -0.7071 0.7071 0 0 0
0 0 0 0 0 0 0 0 0 -1.0000 0 0
0 0 0 0 0 0 0 0 0 0 -0.7071 -0.7071
0 0 0 0 0 0 0 0 0 0 -0.7071 0.7071
T(:,:,20) =
0 -1.0000 0 0 0 0 0 0 0 0 0 0
0.7071 0 0.7071 0 0 0 0 0 0 0 0 0
-0.7071 0 0.7071 0 0 0 0 0 0 0 0 0
0 0 0 0 -1.0000 0 0 0 0 0 0 0
0 0 0 0.7071 0 0.7071 0 0 0 0 0 0
0 0 0 -0.7071 0 0.7071 0 0 0 0 0 0
0 0 0 0 0 0 0 -1.0000 0 0 0 0
0 0 0 0 0 0 0.7071 0 -0.7071 0 0 0
0 0 0 0 0 0 0.7071 0 0.7071 0 0 0
0 0 0 0 0 0 0 0 0 0 -1.0000 0
0 0 0 0 0 0 0 0 0 0.7071 0 -0.7071
0 0 0 0 0 0 0 0 0 0.7071 0 0.7071
T(:,:,21) =
1.0000 0 0 0 0 0 0 0 0 0 0 0
0 0.3310 0.9436 0 0 0 0 0 0 0 0 0
0 -0.9436 0.3310 0 0 0 0 0 0 0 0 0
0 0 0 1.0000 0 0 0 0 0 0 0 0
0 0 0 0 0.3310 0.9436 0 0 0 0 0 0
0 0 0 0 -0.9436 0.3310 0 0 0 0 0 0
0 0 0 0 0 0 1.0000 0 0 0 0 0
0 0 0 0 0 0 0 0.9436 -0.3310 0 0 0
0 0 0 0 0 0 0 0.3310 0.9436 0 0 0
0 0 0 0 0 0 0 0 0 1.0000 0 0
0 0 0 0 0 0 0 0 0 0 0.9436 -0.3310
0 0 0 0 0 0 0 0 0 0 0.3310 0.9436
T(:,:,22) =
0 1.0000 0 0 0 0 0 0 0 0 0 0
0.1221 0 0.9925 0 0 0 0 0 0 0 0 0
0.9925 0 -0.1221 0 0 0 0 0 0 0 0 0
0 0 0 0 1.0000 0 0 0 0 0 0 0
0 0 0 0.1221 0 0.9925 0 0 0 0 0 0
0 0 0 0.9925 0 -0.1221 0 0 0 0 0 0
0 0 0 0 0 0 0 1.0000 0 0 0 0
0 0 0 0 0 0 -0.9925 0 0.1221 0 0 0
0 0 0 0 0 0 0.1221 0 0.9925 0 0 0
0 0 0 0 0 0 0 0 0 0 1.0000 0
0 0 0 0 0 0 0 0 0 -0.9925 0 0.1221
0 0 0 0 0 0 0 0 0 0.1221 0 0.9925
T(:,:,23) =
-1.0000 0 0 0 0 0 0 0 0 0 0 0
0 0.1221 0.9925 0 0 0 0 0 0 0 0 0
0 0.9925 -0.1221 0 0 0 0 0 0 0 0 0
0 0 0 -1.0000 0 0 0 0 0 0 0 0
0 0 0 0 0.1221 0.9925 0 0 0 0 0 0
0 0 0 0 0.9925 -0.1221 0 0 0 0 0 0
0 0 0 0 0 0 -1.0000 0 0 0 0 0
0 0 0 0 0 0 0 -0.9925 0.1221 0 0 0
0 0 0 0 0 0 0 0.1221 0.9925 0 0 0
0 0 0 0 0 0 0 0 0 -1.0000 0 0
0 0 0 0 0 0 0 0 0 0 -0.9925 0.1221
0 0 0 0 0 0 0 0 0 0 0.1221 0.9925
T(:,:,24) =
0 -1.0000 0 0 0 0 0 0 0 0 0 0
-0.1221 0 0.9925 0 0 0 0 0 0 0 0 0
-0.9925 0 -0.1221 0 0 0 0 0 0 0 0 0
0 0 0 0 -1.0000 0 0 0 0 0 0 0
0 0 0 -0.1221 0 0.9925 0 0 0 0 0 0
0 0 0 -0.9925 0 -0.1221 0 0 0 0 0 0
0 0 0 0 0 0 0 -1.0000 0 0 0 0
0 0 0 0 0 0 0.9925 0 0.1221 0 0 0
0 0 0 0 0 0 -0.1221 0 0.9925 0 0 0
0 0 0 0 0 0 0 0 0 0 -1.0000 0
0 0 0 0 0 0 0 0 0 0.9925 0 0.1221
0 0 0 0 0 0 0 0 0 -0.1221 0 0.9925
T(:,:,25) =
1.0000 0 0 0 0 0 0 0 0 0 0 0
0 -0.3162 0.9487 0 0 0 0 0 0 0 0 0
0 -0.9487 -0.3162 0 0 0 0 0 0 0 0 0
0 0 0 1.0000 0 0 0 0 0 0 0 0
0 0 0 0 -0.3162 0.9487 0 0 0 0 0 0
0 0 0 0 -0.9487 -0.3162 0 0 0 0 0 0
0 0 0 0 0 0 1.0000 0 0 0 0 0
0 0 0 0 0 0 0 0.9487 0.3162 0 0 0
0 0 0 0 0 0 0 -0.3162 0.9487 0 0 0
0 0 0 0 0 0 0 0 0 1.0000 0 0
0 0 0 0 0 0 0 0 0 0 0.9487 0.3162
0 0 0 0 0 0 0 0 0 0 -0.3162 0.9487
T(:,:,26) =
0 1.0000 0 0 0 0 0 0 0 0 0 0
0.3162 0 0.9487 0 0 0 0 0 0 0 0 0
0.9487 0 -0.3162 0 0 0 0 0 0 0 0 0
0 0 0 0 1.0000 0 0 0 0 0 0 0
0 0 0 0.3162 0 0.9487 0 0 0 0 0 0
0 0 0 0.9487 0 -0.3162 0 0 0 0 0 0
0 0 0 0 0 0 0 1.0000 0 0 0 0
0 0 0 0 0 0 -0.9487 0 0.3162 0 0 0
0 0 0 0 0 0 0.3162 0 0.9487 0 0 0
0 0 0 0 0 0 0 0 0 0 1.0000 0
0 0 0 0 0 0 0 0 0 -0.9487 0 0.3162
0 0 0 0 0 0 0 0 0 0.3162 0 0.9487
T(:,:,27) =
-1.0000 0 0 0 0 0 0 0 0 0 0 0
0 0.3162 0.9487 0 0 0 0 0 0 0 0 0
0 0.9487 -0.3162 0 0 0 0 0 0 0 0 0
0 0 0 -1.0000 0 0 0 0 0 0 0 0
0 0 0 0 0.3162 0.9487 0 0 0 0 0 0
0 0 0 0 0.9487 -0.3162 0 0 0 0 0 0
0 0 0 0 0 0 -1.0000 0 0 0 0 0
0 0 0 0 0 0 0 -0.9487 0.3162 0 0 0
0 0 0 0 0 0 0 0.3162 0.9487 0 0 0
0 0 0 0 0 0 0 0 0 -1.0000 0 0
0 0 0 0 0 0 0 0 0 0 -0.9487 0.3162
0 0 0 0 0 0 0 0 0 0 0.3162 0.9487
T(:,:,28) =
0 -1.0000 0 0 0 0 0 0 0 0 0 0
-0.3162 0 0.9487 0 0 0 0 0 0 0 0 0
-0.9487 0 -0.3162 0 0 0 0 0 0 0 0 0
0 0 0 0 -1.0000 0 0 0 0 0 0 0
0 0 0 -0.3162 0 0.9487 0 0 0 0 0 0
0 0 0 -0.9487 0 -0.3162 0 0 0 0 0 0
0 0 0 0 0 0 0 -1.0000 0 0 0 0
0 0 0 0 0 0 0.9487 0 0.3162 0 0 0
0 0 0 0 0 0 -0.3162 0 0.9487 0 0 0
0 0 0 0 0 0 0 0 0 0 -1.0000 0
0 0 0 0 0 0 0 0 0 0.9487 0 0.3162
0 0 0 0 0 0 0 0 0 -0.3162 0 0.9487
T(:,:,29) =
1.0000 0 0 0 0 0 0 0 0 0 0 0
0 0.1098 0.9940 0 0 0 0 0 0 0 0 0
0 -0.9940 0.1098 0 0 0 0 0 0 0 0 0
0 0 0 1.0000 0 0 0 0 0 0 0 0
0 0 0 0 0.1098 0.9940 0 0 0 0 0 0
0 0 0 0 -0.9940 0.1098 0 0 0 0 0 0
0 0 0 0 0 0 1.0000 0 0 0 0 0
0 0 0 0 0 0 0 0.9940 -0.1098 0 0 0
0 0 0 0 0 0 0 0.1098 0.9940 0 0 0
0 0 0 0 0 0 0 0 0 1.0000 0 0
0 0 0 0 0 0 0 0 0 0 0.9940 -0.1098
0 0 0 0 0 0 0 0 0 0 0.1098 0.9940
T(:,:,30) =
0 1.0000 0 0 0 0 0 0 0 0 0 0
-0.1098 0 0.9940 0 0 0 0 0 0 0 0 0
0.9940 0 0.1098 0 0 0 0 0 0 0 0 0
0 0 0 0 1.0000 0 0 0 0 0 0 0
0 0 0 -0.1098 0 0.9940 0 0 0 0 0 0
0 0 0 0.9940 0 0.1098 0 0 0 0 0 0
0 0 0 0 0 0 0 1.0000 0 0 0 0
0 0 0 0 0 0 -0.9940 0 -0.1098 0 0 0
0 0 0 0 0 0 -0.1098 0 0.9940 0 0 0
0 0 0 0 0 0 0 0 0 0 1.0000 0
0 0 0 0 0 0 0 0 0 -0.9940 0 -0.1098
0 0 0 0 0 0 0 0 0 -0.1098 0 0.9940
T(:,:,31) =
-1.0000 0 0 0 0 0 0 0 0 0 0 0
0 -0.1098 0.9940 0 0 0 0 0 0 0 0 0
0 0.9940 0.1098 0 0 0 0 0 0 0 0 0
0 0 0 -1.0000 0 0 0 0 0 0 0 0
0 0 0 0 -0.1098 0.9940 0 0 0 0 0 0
0 0 0 0 0.9940 0.1098 0 0 0 0 0 0
0 0 0 0 0 0 -1.0000 0 0 0 0 0
0 0 0 0 0 0 0 -0.9940 -0.1098 0 0 0
0 0 0 0 0 0 0 -0.1098 0.9940 0 0 0
0 0 0 0 0 0 0 0 0 -1.0000 0 0
0 0 0 0 0 0 0 0 0 0 -0.9940 -0.1098
0 0 0 0 0 0 0 0 0 0 -0.1098 0.9940
T(:,:,32) =
0 -1.0000 0 0 0 0 0 0 0 0 0 0
0.1098 0 0.9940 0 0 0 0 0 0 0 0 0
-0.9940 0 0.1098 0 0 0 0 0 0 0 0 0
0 0 0 0 -1.0000 0 0 0 0 0 0 0
0 0 0 0.1098 0 0.9940 0 0 0 0 0 0
0 0 0 -0.9940 0 0.1098 0 0 0 0 0 0
0 0 0 0 0 0 0 -1.0000 0 0 0 0
0 0 0 0 0 0 0.9940 0 -0.1098 0 0 0
0 0 0 0 0 0 0.1098 0 0.9940 0 0 0
0 0 0 0 0 0 0 0 0 0 -1.0000 0
0 0 0 0 0 0 0 0 0 0.9940 0 -0.1098
0 0 0 0 0 0 0 0 0 0.1098 0 0.9940
Hu'
ans =
Columns 1 through 13
0 -0.0011 0 0.0011 0 -0.0008 0 0.0008 0.0000 -0.0007 0.0000 0.0007 0.0000
0 0 0 0 0 0 0 0 0.0008 -0.0006 0.0008 -0.0006 0.0008
0 0 0 0 0 0 0 0 -0.0020 -0.0050 -0.0020 -0.0050 -0.0008
0 0 0 0 0 0 0 0 -0.0006 -0.0002 -0.0006 -0.0002 -0.0002
0 0 0 0 0 0 0 0 0.0000 -0.0004 0.0000 0.0004 0.0000
0 0 0 0 0 0 0 0 0.0000 -0.0000 0.0000 0.0000 -0.0000
-0.0011 0 0.0011 0 -0.0008 0 0.0008 0 -0.0007 -0.0000 0.0007 -0.0000 0.0003
0 0 0 0 0 0 0 0 -0.0006 0.0003 -0.0006 0.0003 0.0009
0 0 0 0 0 0 0 0 -0.0050 -0.0054 -0.0050 -0.0054 -0.0005
0 0 0 0 0 0 0 0 -0.0002 0.0004 -0.0002 0.0004 -0.0001
0 0 0 0 0 0 0 0 -0.0004 0.0000 0.0004 0.0000 -0.0001
0 0 0 0 0 0 0 0 -0.0000 -0.0000 0.0000 -0.0000 -0.0000
Columns 14 through 26
0.0003 0.0000 -0.0003 0 -0.0011 0 0.0011 0.0000 -0.0007 -0.0000 0.0007 0 -0.0008
0.0009 0.0008 0.0009 0 0 0 0 0.0008 -0.0006 0.0003 -0.0006 0 0
-0.0005 -0.0008 -0.0005 0 0 0 0 -0.0020 -0.0050 -0.0054 -0.0050 0 0
-0.0001 -0.0002 -0.0001 0 0 0 0 -0.0006 -0.0002 0.0004 -0.0002 0 0
-0.0001 0.0000 0.0001 0 0 0 0 0.0000 -0.0004 0.0000 0.0004 0 0
-0.0000 -0.0000 0.0000 0 0 0 0 0.0000 -0.0000 -0.0000 0.0000 0 0
0.0000 -0.0003 0.0000 0.0000 -0.0007 -0.0000 0.0007 0.0000 0.0000 0.0000 0.0000 0.0000 0.0003
0.0014 0.0009 0.0014 0.0008 -0.0006 0.0003 -0.0006 0.0013 0.0013 0.0013 0.0013 0.0008 0.0009
-0.0006 -0.0005 -0.0006 -0.0020 -0.0050 -0.0054 -0.0050 -0.0022 -0.0022 -0.0022 -0.0022 -0.0008 -0.0005
0.0001 -0.0001 0.0001 -0.0006 -0.0002 0.0004 -0.0002 -0.0002 -0.0002 -0.0002 -0.0002 -0.0002 -0.0001
0.0000 0.0001 0.0000 0.0000 -0.0004 0.0000 0.0004 0.0000 0.0000 0.0000 0.0000 0.0000 -0.0001
-0.0000 0.0000 -0.0000 0.0000 -0.0000 -0.0000 0.0000 -0.0000 -0.0000 -0.0000 -0.0000 -0.0000 -0.0000
Columns 27 through 32
0 0.0008 0.0000 0.0003 0.0000 -0.0003
0 0 0.0008 0.0009 0.0014 0.0009
0 0 -0.0008 -0.0005 -0.0006 -0.0005
0 0 -0.0002 -0.0001 0.0001 -0.0001
0 0 0.0000 -0.0001 0.0000 0.0001
0 0 -0.0000 -0.0000 -0.0000 0.0000
0.0000 -0.0003 0.0000 0.0000 0.0000 0.0000
0.0014 0.0009 0.0013 0.0013 0.0013 0.0013
-0.0006 -0.0005 -0.0022 -0.0022 -0.0022 -0.0022
0.0001 -0.0001 -0.0002 -0.0002 -0.0002 -0.0002
0.0000 0.0001 0.0000 0.0000 0.0000 0.0000
-0.0000 0.0000 -0.0000 -0.0000 -0.0000 -0.0000
D'
ans =
Columns 1 through 13
-0.0011 0.0011 -0.0008 0.0008 0.0000 0.0008 -0.0020 -0.0006 0.0000 0.0000 -0.0007 -0.0006 -0.0050
Columns 14 through 26
-0.0002 -0.0004 -0.0000 -0.0000 0.0003 -0.0054 0.0004 0.0000 -0.0000 0.0007 -0.0006 -0.0050 -0.0002
Columns 27 through 39
0.0004 0.0000 0.0000 0.0008 -0.0008 -0.0002 0.0000 -0.0000 0.0003 0.0009 -0.0005 -0.0001 -0.0001
Columns 40 through 52
-0.0000 0.0000 0.0014 -0.0006 0.0001 0.0000 -0.0000 -0.0003 0.0009 -0.0005 -0.0001 0.0001 0.0000
Columns 53 through 58
0.0000 0.0013 -0.0022 -0.0002 0.0000 -0.0000
Re
Re =
0 0 0 0 0 0
1 0 0 0 0 0
0 0 0 0 0 0
2 0 0 0 0 0
0 0 0 0 0 0
3 0 0 0 0 0
0 0 0 0 0 0
4 0 0 0 0 0
5 6 7 8 9 10
11 12 13 14 15 16
17 18 19 20 21 22
23 24 25 26 27 28
29 30 31 32 33 34
35 36 37 38 39 40
41 42 43 44 45 46
47 48 49 50 51 52
53 54 55 56 57 58
R =
0 0 0 0 0 0 1 0 0 0 0 0
1 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 2 0 0 0 0 0
2 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 3 0 0 0 0 0
3 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 4 0 0 0 0 0
4 0 0 0 0 0 0 0 0 0 0 0
5 6 7 8 9 10 11 12 13 14 15 16
11 12 13 14 15 16 17 18 19 20 21 22
5 6 7 8 9 10 23 24 25 26 27 28
23 24 25 26 27 28 17 18 19 20 21 22
29 30 31 32 33 34 35 36 37 38 39 40
35 36 37 38 39 40 41 42 43 44 45 46
29 30 31 32 33 34 47 48 49 50 51 52
47 48 49 50 51 52 41 42 43 44 45 46
0 0 0 0 0 0 5 6 7 8 9 10
1 0 0 0 0 0 11 12 13 14 15 16
0 0 0 0 0 0 17 18 19 20 21 22
2 0 0 0 0 0 23 24 25 26 27 28
5 6 7 8 9 10 53 54 55 56 57 58
11 12 13 14 15 16 53 54 55 56 57 58
17 18 19 20 21 22 53 54 55 56 57 58
23 24 25 26 27 28 53 54 55 56 57 58
0 0 0 0 0 0 29 30 31 32 33 34
3 0 0 0 0 0 35 36 37 38 39 40
0 0 0 0 0 0 41 42 43 44 45 46
4 0 0 0 0 0 47 48 49 50 51 52
29 30 31 32 33 34 53 54 55 56 57 58
35 36 37 38 39 40 53 54 55 56 57 58
41 42 43 44 45 46 53 54 55 56 57 58
47 48 49 50 51 52 53 54 55 56 57 58
Published
with MATLAB® 7.1