Created
December 15, 2016 13:31
-
-
Save anonymous/3d1b979692d460e73acce7de6405e740 to your computer and use it in GitHub Desktop.
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
static void diag_MK_matrix(FLOAT_XX EM, FLOAT_XX mu, FLOAT_XX ro, | |
FLOAT_XX a, FLOAT_XX b, FLOAT_XX xx, | |
FLOAT_XX* poutM, FLOAT_XX* poutK) { | |
*poutM = 2 * a*b*Pi*ro*xx; | |
FLOAT_XX matrixK[8][8] = { List(8 * a*b*Pi*(EM / (32.*Power(b, 2)*(1 + mu)) - (-(EM*(1 - mu)) / (4.*a*(1 - 2 * mu)*(1 + mu)) + (EM*mu) / (4.*(1 - 2 * mu)*(1 + mu)*xx)) / (4.*a) + (-(EM*mu) / (4.*a*(1 - 2 * mu)*(1 + mu)) + (EM*(1 - mu)) / (4.*(1 - 2 * mu)*(1 + mu)*xx)) / (4.*xx))*xx, 8 * a*b*Pi*(EM / (32.*a*b*(1 + mu)) - (-(EM*mu) / (4.*a*(1 - 2 * mu)*(1 + mu)) + (EM*mu) / (4.*(1 - 2 * mu)*(1 + mu)*xx)) / (4.*b))*xx, | |
8 * a*b*Pi*(-EM / (32.*Power(b, 2)*(1 + mu)) - (-(EM*(1 - mu)) / (4.*a*(1 - 2 * mu)*(1 + mu)) + (EM*mu) / (4.*(1 - 2 * mu)*(1 + mu)*xx)) / (4.*a) + (-(EM*mu) / (4.*a*(1 - 2 * mu)*(1 + mu)) + (EM*(1 - mu)) / (4.*(1 - 2 * mu)*(1 + mu)*xx)) / (4.*xx))*xx, 8 * a*b*Pi*(EM / (32.*a*b*(1 + mu)) + (-(EM*mu) / (4.*a*(1 - 2 * mu)*(1 + mu)) + (EM*mu) / (4.*(1 - 2 * mu)*(1 + mu)*xx)) / (4.*b))*xx, | |
8 * a*b*Pi*(EM / (32.*Power(b, 2)*(1 + mu)) + (-(EM*(1 - mu)) / (4.*a*(1 - 2 * mu)*(1 + mu)) + (EM*mu) / (4.*(1 - 2 * mu)*(1 + mu)*xx)) / (4.*a) + (-(EM*mu) / (4.*a*(1 - 2 * mu)*(1 + mu)) + (EM*(1 - mu)) / (4.*(1 - 2 * mu)*(1 + mu)*xx)) / (4.*xx))*xx, 8 * a*b*Pi*(-EM / (32.*a*b*(1 + mu)) - (-(EM*mu) / (4.*a*(1 - 2 * mu)*(1 + mu)) + (EM*mu) / (4.*(1 - 2 * mu)*(1 + mu)*xx)) / (4.*b))*xx, | |
8 * a*b*Pi*(-EM / (32.*Power(b, 2)*(1 + mu)) + (-(EM*(1 - mu)) / (4.*a*(1 - 2 * mu)*(1 + mu)) + (EM*mu) / (4.*(1 - 2 * mu)*(1 + mu)*xx)) / (4.*a) + (-(EM*mu) / (4.*a*(1 - 2 * mu)*(1 + mu)) + (EM*(1 - mu)) / (4.*(1 - 2 * mu)*(1 + mu)*xx)) / (4.*xx))*xx, 8 * a*b*Pi*(-EM / (32.*a*b*(1 + mu)) + (-(EM*mu) / (4.*a*(1 - 2 * mu)*(1 + mu)) + (EM*mu) / (4.*(1 - 2 * mu)*(1 + mu)*xx)) / (4.*b))*xx), | |
List(8 * a*b*Pi*(EM / (32.*a*b*(1 + mu)) + (EM*mu) / (16.*a*b*(1 - 2 * mu)*(1 + mu)) - (EM*mu) / (16.*b*(1 - 2 * mu)*(1 + mu)*xx))*xx, 8 * a*b*(EM / (32.*Power(a, 2)*(1 + mu)) + (EM*(1 - mu)) / (16.*Power(b, 2)*(1 - 2 * mu)*(1 + mu)))*Pi*xx, | |
8 * a*b*Pi*(-EM / (32.*a*b*(1 + mu)) + (EM*mu) / (16.*a*b*(1 - 2 * mu)*(1 + mu)) - (EM*mu) / (16.*b*(1 - 2 * mu)*(1 + mu)*xx))*xx, 8 * a*b*(EM / (32.*Power(a, 2)*(1 + mu)) - (EM*(1 - mu)) / (16.*Power(b, 2)*(1 - 2 * mu)*(1 + mu)))*Pi*xx, | |
8 * a*b*Pi*(EM / (32.*a*b*(1 + mu)) - (EM*mu) / (16.*a*b*(1 - 2 * mu)*(1 + mu)) - (EM*mu) / (16.*b*(1 - 2 * mu)*(1 + mu)*xx))*xx, 8 * a*b*(-EM / (32.*Power(a, 2)*(1 + mu)) + (EM*(1 - mu)) / (16.*Power(b, 2)*(1 - 2 * mu)*(1 + mu)))*Pi*xx, | |
8 * a*b*Pi*(-EM / (32.*a*b*(1 + mu)) - (EM*mu) / (16.*a*b*(1 - 2 * mu)*(1 + mu)) - (EM*mu) / (16.*b*(1 - 2 * mu)*(1 + mu)*xx))*xx, 8 * a*b*(-EM / (32.*Power(a, 2)*(1 + mu)) - (EM*(1 - mu)) / (16.*Power(b, 2)*(1 - 2 * mu)*(1 + mu)))*Pi*xx), | |
List(8 * a*b*Pi*(-EM / (32.*Power(b, 2)*(1 + mu)) - (-(EM*(1 - mu)) / (4.*a*(1 - 2 * mu)*(1 + mu)) + (EM*mu) / (4.*(1 - 2 * mu)*(1 + mu)*xx)) / (4.*a) + (-(EM*mu) / (4.*a*(1 - 2 * mu)*(1 + mu)) + (EM*(1 - mu)) / (4.*(1 - 2 * mu)*(1 + mu)*xx)) / (4.*xx))*xx, 8 * a*b*Pi*(-EM / (32.*a*b*(1 + mu)) - (-(EM*mu) / (4.*a*(1 - 2 * mu)*(1 + mu)) + (EM*mu) / (4.*(1 - 2 * mu)*(1 + mu)*xx)) / (4.*b))*xx, | |
8 * a*b*Pi*(EM / (32.*Power(b, 2)*(1 + mu)) - (-(EM*(1 - mu)) / (4.*a*(1 - 2 * mu)*(1 + mu)) + (EM*mu) / (4.*(1 - 2 * mu)*(1 + mu)*xx)) / (4.*a) + (-(EM*mu) / (4.*a*(1 - 2 * mu)*(1 + mu)) + (EM*(1 - mu)) / (4.*(1 - 2 * mu)*(1 + mu)*xx)) / (4.*xx))*xx, 8 * a*b*Pi*(-EM / (32.*a*b*(1 + mu)) + (-(EM*mu) / (4.*a*(1 - 2 * mu)*(1 + mu)) + (EM*mu) / (4.*(1 - 2 * mu)*(1 + mu)*xx)) / (4.*b))*xx, | |
8 * a*b*Pi*(-EM / (32.*Power(b, 2)*(1 + mu)) + (-(EM*(1 - mu)) / (4.*a*(1 - 2 * mu)*(1 + mu)) + (EM*mu) / (4.*(1 - 2 * mu)*(1 + mu)*xx)) / (4.*a) + (-(EM*mu) / (4.*a*(1 - 2 * mu)*(1 + mu)) + (EM*(1 - mu)) / (4.*(1 - 2 * mu)*(1 + mu)*xx)) / (4.*xx))*xx, 8 * a*b*Pi*(EM / (32.*a*b*(1 + mu)) - (-(EM*mu) / (4.*a*(1 - 2 * mu)*(1 + mu)) + (EM*mu) / (4.*(1 - 2 * mu)*(1 + mu)*xx)) / (4.*b))*xx, | |
8 * a*b*Pi*(EM / (32.*Power(b, 2)*(1 + mu)) + (-(EM*(1 - mu)) / (4.*a*(1 - 2 * mu)*(1 + mu)) + (EM*mu) / (4.*(1 - 2 * mu)*(1 + mu)*xx)) / (4.*a) + (-(EM*mu) / (4.*a*(1 - 2 * mu)*(1 + mu)) + (EM*(1 - mu)) / (4.*(1 - 2 * mu)*(1 + mu)*xx)) / (4.*xx))*xx, 8 * a*b*Pi*(EM / (32.*a*b*(1 + mu)) + (-(EM*mu) / (4.*a*(1 - 2 * mu)*(1 + mu)) + (EM*mu) / (4.*(1 - 2 * mu)*(1 + mu)*xx)) / (4.*b))*xx), | |
List(8 * a*b*Pi*(EM / (32.*a*b*(1 + mu)) - (EM*mu) / (16.*a*b*(1 - 2 * mu)*(1 + mu)) + (EM*mu) / (16.*b*(1 - 2 * mu)*(1 + mu)*xx))*xx, 8 * a*b*(EM / (32.*Power(a, 2)*(1 + mu)) - (EM*(1 - mu)) / (16.*Power(b, 2)*(1 - 2 * mu)*(1 + mu)))*Pi*xx, | |
8 * a*b*Pi*(-EM / (32.*a*b*(1 + mu)) - (EM*mu) / (16.*a*b*(1 - 2 * mu)*(1 + mu)) + (EM*mu) / (16.*b*(1 - 2 * mu)*(1 + mu)*xx))*xx, 8 * a*b*(EM / (32.*Power(a, 2)*(1 + mu)) + (EM*(1 - mu)) / (16.*Power(b, 2)*(1 - 2 * mu)*(1 + mu)))*Pi*xx, | |
8 * a*b*Pi*(EM / (32.*a*b*(1 + mu)) + (EM*mu) / (16.*a*b*(1 - 2 * mu)*(1 + mu)) + (EM*mu) / (16.*b*(1 - 2 * mu)*(1 + mu)*xx))*xx, 8 * a*b*(-EM / (32.*Power(a, 2)*(1 + mu)) - (EM*(1 - mu)) / (16.*Power(b, 2)*(1 - 2 * mu)*(1 + mu)))*Pi*xx, | |
8 * a*b*Pi*(-EM / (32.*a*b*(1 + mu)) + (EM*mu) / (16.*a*b*(1 - 2 * mu)*(1 + mu)) + (EM*mu) / (16.*b*(1 - 2 * mu)*(1 + mu)*xx))*xx, 8 * a*b*(-EM / (32.*Power(a, 2)*(1 + mu)) + (EM*(1 - mu)) / (16.*Power(b, 2)*(1 - 2 * mu)*(1 + mu)))*Pi*xx), | |
List(8 * a*b*Pi*(EM / (32.*Power(b, 2)*(1 + mu)) - ((EM*(1 - mu)) / (4.*a*(1 - 2 * mu)*(1 + mu)) + (EM*mu) / (4.*(1 - 2 * mu)*(1 + mu)*xx)) / (4.*a) + ((EM*mu) / (4.*a*(1 - 2 * mu)*(1 + mu)) + (EM*(1 - mu)) / (4.*(1 - 2 * mu)*(1 + mu)*xx)) / (4.*xx))*xx, 8 * a*b*Pi*(EM / (32.*a*b*(1 + mu)) - ((EM*mu) / (4.*a*(1 - 2 * mu)*(1 + mu)) + (EM*mu) / (4.*(1 - 2 * mu)*(1 + mu)*xx)) / (4.*b))*xx, | |
8 * a*b*Pi*(-EM / (32.*Power(b, 2)*(1 + mu)) - ((EM*(1 - mu)) / (4.*a*(1 - 2 * mu)*(1 + mu)) + (EM*mu) / (4.*(1 - 2 * mu)*(1 + mu)*xx)) / (4.*a) + ((EM*mu) / (4.*a*(1 - 2 * mu)*(1 + mu)) + (EM*(1 - mu)) / (4.*(1 - 2 * mu)*(1 + mu)*xx)) / (4.*xx))*xx, 8 * a*b*Pi*(EM / (32.*a*b*(1 + mu)) + ((EM*mu) / (4.*a*(1 - 2 * mu)*(1 + mu)) + (EM*mu) / (4.*(1 - 2 * mu)*(1 + mu)*xx)) / (4.*b))*xx, | |
8 * a*b*Pi*(EM / (32.*Power(b, 2)*(1 + mu)) + ((EM*(1 - mu)) / (4.*a*(1 - 2 * mu)*(1 + mu)) + (EM*mu) / (4.*(1 - 2 * mu)*(1 + mu)*xx)) / (4.*a) + ((EM*mu) / (4.*a*(1 - 2 * mu)*(1 + mu)) + (EM*(1 - mu)) / (4.*(1 - 2 * mu)*(1 + mu)*xx)) / (4.*xx))*xx, 8 * a*b*Pi*(-EM / (32.*a*b*(1 + mu)) - ((EM*mu) / (4.*a*(1 - 2 * mu)*(1 + mu)) + (EM*mu) / (4.*(1 - 2 * mu)*(1 + mu)*xx)) / (4.*b))*xx, | |
8 * a*b*Pi*(-EM / (32.*Power(b, 2)*(1 + mu)) + ((EM*(1 - mu)) / (4.*a*(1 - 2 * mu)*(1 + mu)) + (EM*mu) / (4.*(1 - 2 * mu)*(1 + mu)*xx)) / (4.*a) + ((EM*mu) / (4.*a*(1 - 2 * mu)*(1 + mu)) + (EM*(1 - mu)) / (4.*(1 - 2 * mu)*(1 + mu)*xx)) / (4.*xx))*xx, 8 * a*b*Pi*(-EM / (32.*a*b*(1 + mu)) + ((EM*mu) / (4.*a*(1 - 2 * mu)*(1 + mu)) + (EM*mu) / (4.*(1 - 2 * mu)*(1 + mu)*xx)) / (4.*b))*xx), | |
List(8 * a*b*Pi*(-EM / (32.*a*b*(1 + mu)) + (EM*mu) / (16.*a*b*(1 - 2 * mu)*(1 + mu)) - (EM*mu) / (16.*b*(1 - 2 * mu)*(1 + mu)*xx))*xx, 8 * a*b*(-EM / (32.*Power(a, 2)*(1 + mu)) + (EM*(1 - mu)) / (16.*Power(b, 2)*(1 - 2 * mu)*(1 + mu)))*Pi*xx, | |
8 * a*b*Pi*(EM / (32.*a*b*(1 + mu)) + (EM*mu) / (16.*a*b*(1 - 2 * mu)*(1 + mu)) - (EM*mu) / (16.*b*(1 - 2 * mu)*(1 + mu)*xx))*xx, 8 * a*b*(-EM / (32.*Power(a, 2)*(1 + mu)) - (EM*(1 - mu)) / (16.*Power(b, 2)*(1 - 2 * mu)*(1 + mu)))*Pi*xx, | |
8 * a*b*Pi*(-EM / (32.*a*b*(1 + mu)) - (EM*mu) / (16.*a*b*(1 - 2 * mu)*(1 + mu)) - (EM*mu) / (16.*b*(1 - 2 * mu)*(1 + mu)*xx))*xx, 8 * a*b*(EM / (32.*Power(a, 2)*(1 + mu)) + (EM*(1 - mu)) / (16.*Power(b, 2)*(1 - 2 * mu)*(1 + mu)))*Pi*xx, | |
8 * a*b*Pi*(EM / (32.*a*b*(1 + mu)) - (EM*mu) / (16.*a*b*(1 - 2 * mu)*(1 + mu)) - (EM*mu) / (16.*b*(1 - 2 * mu)*(1 + mu)*xx))*xx, 8 * a*b*(EM / (32.*Power(a, 2)*(1 + mu)) - (EM*(1 - mu)) / (16.*Power(b, 2)*(1 - 2 * mu)*(1 + mu)))*Pi*xx), | |
List(8 * a*b*Pi*(-EM / (32.*Power(b, 2)*(1 + mu)) - ((EM*(1 - mu)) / (4.*a*(1 - 2 * mu)*(1 + mu)) + (EM*mu) / (4.*(1 - 2 * mu)*(1 + mu)*xx)) / (4.*a) + ((EM*mu) / (4.*a*(1 - 2 * mu)*(1 + mu)) + (EM*(1 - mu)) / (4.*(1 - 2 * mu)*(1 + mu)*xx)) / (4.*xx))*xx, 8 * a*b*Pi*(-EM / (32.*a*b*(1 + mu)) - ((EM*mu) / (4.*a*(1 - 2 * mu)*(1 + mu)) + (EM*mu) / (4.*(1 - 2 * mu)*(1 + mu)*xx)) / (4.*b))*xx, | |
8 * a*b*Pi*(EM / (32.*Power(b, 2)*(1 + mu)) - ((EM*(1 - mu)) / (4.*a*(1 - 2 * mu)*(1 + mu)) + (EM*mu) / (4.*(1 - 2 * mu)*(1 + mu)*xx)) / (4.*a) + ((EM*mu) / (4.*a*(1 - 2 * mu)*(1 + mu)) + (EM*(1 - mu)) / (4.*(1 - 2 * mu)*(1 + mu)*xx)) / (4.*xx))*xx, 8 * a*b*Pi*(-EM / (32.*a*b*(1 + mu)) + ((EM*mu) / (4.*a*(1 - 2 * mu)*(1 + mu)) + (EM*mu) / (4.*(1 - 2 * mu)*(1 + mu)*xx)) / (4.*b))*xx, | |
8 * a*b*Pi*(-EM / (32.*Power(b, 2)*(1 + mu)) + ((EM*(1 - mu)) / (4.*a*(1 - 2 * mu)*(1 + mu)) + (EM*mu) / (4.*(1 - 2 * mu)*(1 + mu)*xx)) / (4.*a) + ((EM*mu) / (4.*a*(1 - 2 * mu)*(1 + mu)) + (EM*(1 - mu)) / (4.*(1 - 2 * mu)*(1 + mu)*xx)) / (4.*xx))*xx, 8 * a*b*Pi*(EM / (32.*a*b*(1 + mu)) - ((EM*mu) / (4.*a*(1 - 2 * mu)*(1 + mu)) + (EM*mu) / (4.*(1 - 2 * mu)*(1 + mu)*xx)) / (4.*b))*xx, | |
8 * a*b*Pi*(EM / (32.*Power(b, 2)*(1 + mu)) + ((EM*(1 - mu)) / (4.*a*(1 - 2 * mu)*(1 + mu)) + (EM*mu) / (4.*(1 - 2 * mu)*(1 + mu)*xx)) / (4.*a) + ((EM*mu) / (4.*a*(1 - 2 * mu)*(1 + mu)) + (EM*(1 - mu)) / (4.*(1 - 2 * mu)*(1 + mu)*xx)) / (4.*xx))*xx, 8 * a*b*Pi*(EM / (32.*a*b*(1 + mu)) + ((EM*mu) / (4.*a*(1 - 2 * mu)*(1 + mu)) + (EM*mu) / (4.*(1 - 2 * mu)*(1 + mu)*xx)) / (4.*b))*xx), | |
List(8 * a*b*Pi*(-EM / (32.*a*b*(1 + mu)) - (EM*mu) / (16.*a*b*(1 - 2 * mu)*(1 + mu)) + (EM*mu) / (16.*b*(1 - 2 * mu)*(1 + mu)*xx))*xx, 8 * a*b*(-EM / (32.*Power(a, 2)*(1 + mu)) - (EM*(1 - mu)) / (16.*Power(b, 2)*(1 - 2 * mu)*(1 + mu)))*Pi*xx, | |
8 * a*b*Pi*(EM / (32.*a*b*(1 + mu)) - (EM*mu) / (16.*a*b*(1 - 2 * mu)*(1 + mu)) + (EM*mu) / (16.*b*(1 - 2 * mu)*(1 + mu)*xx))*xx, 8 * a*b*(-EM / (32.*Power(a, 2)*(1 + mu)) + (EM*(1 - mu)) / (16.*Power(b, 2)*(1 - 2 * mu)*(1 + mu)))*Pi*xx, | |
8 * a*b*Pi*(-EM / (32.*a*b*(1 + mu)) + (EM*mu) / (16.*a*b*(1 - 2 * mu)*(1 + mu)) + (EM*mu) / (16.*b*(1 - 2 * mu)*(1 + mu)*xx))*xx, 8 * a*b*(EM / (32.*Power(a, 2)*(1 + mu)) - (EM*(1 - mu)) / (16.*Power(b, 2)*(1 - 2 * mu)*(1 + mu)))*Pi*xx, | |
8 * a*b*Pi*(EM / (32.*a*b*(1 + mu)) + (EM*mu) / (16.*a*b*(1 - 2 * mu)*(1 + mu)) + (EM*mu) / (16.*b*(1 - 2 * mu)*(1 + mu)*xx))*xx, 8 * a*b*(EM / (32.*Power(a, 2)*(1 + mu)) + (EM*(1 - mu)) / (16.*Power(b, 2)*(1 - 2 * mu)*(1 + mu)))*Pi*xx) | |
}; | |
//memcpy(poutM, &(matrixM[0][0]), sizeof(FLOAT_XX) * 8 * 8); | |
memcpy(poutK, &(matrixK[0][0]), sizeof(FLOAT_XX) * 8 * 8); | |
} | |
int main() | |
{ | |
FLOAT_XX M[8][8],K[8][8],sM; | |
FLOAT_XX t = Power(1.3333, 0.2); | |
//MK_matrix(1.1, 2.2, 3.3, 4.4, 5.5, 6.6, &M[0][0], &K[0][0]); | |
diag_MK_matrix(1.1, 2.2, 3.3, 4.4, 5.5, 6.6, &sM, &K[0][0]); | |
return 0; | |
} |
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment