%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Phase 2. Formulation 
% (recuperation des resultats du cours)

syms Nsc Lc Ec Sc roSc   Nsp Lp Ep Sp roSp Ip   

% - construction des elements de reference de poutre Eb 2D
KpP=(-Nsp/(30*Lp))*[0 0 0 0 0 0; 0 36 3*Lp 0 -36 3*Lp ; 0 3*Lp 4*Lp^2 0 -3*Lp -Lp^2 ; 0 0 0 0 0 0 ; 0 -36 -3*Lp 0 36 -3*Lp ; 0 3*Lp -Lp^2 0 -3*Lp 4*Lp^2];
KeP=(Ep*Ip/Lp^3)*[Sp*Lp^2/Ip 0 0 -Sp*Lp^2/Ip 0 0 ; 0 12 6*Lp 0 -12 6*Lp ; 0 6*Lp 4*Lp^2 0 -6*Lp 2*Lp^2 ; -Sp*Lp^2/Ip 0 0 Sp*Lp^2/Ip 0  0;  0 -12 -6*Lp 0 12 -6*Lp ; 0 6*Lp 2*Lp^2 0 -6*Lp 4*Lp^2];
MP=(roSp*Lp/420)*[140 0 0 70 0 0 ; 0 156 22*Lp 0 54 -13*Lp ; 0 22*Lp 4*Lp^2 0 13*Lp -3*Lp^2 ; 70 0 0 140 0 0 ; 0 54 13*Lp 0 156 -22*Lp ; 0 -13*Lp -3*Lp^2 0 -22*Lp 4*Lp^2];

% - construction des elements de reference de cable tendu ( = barre 2D)
KpC=(-Nsc/Lc)*[0 0 0 0 ; 0 1 0 -1 ; 0 0 0 0 ; 0 -1 0 1];
KeC=(Ec*Sc/Lc)*[1 0 -1 0 ; 0 0 0 0 ; -1 0 1 0 ; 0 0 0 0];
MC=(roSc*Lc/6)*[2 0 1 0 ; 0 2 0 1 ; 1 0 2 0 ; 0 1 0 2];
