w01=42;
w02=49;

R1=1.35;
Kem1=0.1;
Kme1=Kem1;
N1=35;
l1=0.42;
m1=3.312;
J1=122e-6;

R2=8;
Kem2=48/(38*69);
Kme2=Kem2;
N2=69;
l2=0.28;
m2=0.42;
%J2=2.0964e-6;
J2=1.1261e-6;

D11max=1/3*m1*l1^2+1/3*m2*l2^2+m2*l2*l1+m2*l1^2;
D22max=1/3*m2*l2^2;

Kp1=(0.5*w01)^2*R1*(D11max+J1*N1^2)/(N1*Kem1);
Kp2=(0.5*w02)^2*R2*(D22max+J2*N2^2)/(N2*Kem2);

Kd1=2*sqrt(Kp1*(D11max+J1*N1^2)*R1/(Kem1*N1))-Kme1*N1;
Kd2=2*sqrt(Kp2*(D22max+J2*N2^2)*R2/(Kem2*N2))-Kme2*N2;

disp('Kp1, Kd1');
[Kp1, Kd1]

disp('Kp2, Kd2');
[Kp2, Kd2]
