a = 4;(* nodes Х (-a,a)*)
b = 4;(* nodes У (-b,b)*)
n = (2 a + 1) ( 2 b + 1);(*all nodes *)
x0 = 2;
y0 = 2;
k = 1/2 // N;
Ax = k/(x0)^2;
Ay = k/(y0)^2;
EE = 2 10^5;
μ = 0.3;
h = 0.8;
Dc = EE*h^3/12/(1 - μ^2) ;
ϕ = {Flatten[Table[Exp[-Ax (x - xi)^2] Exp[-Ay (y - yi)^2], {xi, -a,
a}, {yi, -b, b} ]]};
ϕT = Transpose[ϕ];
solK = Dc/2 Laplacian[(ϕT . ϕ),{x,y}];
K1 = NIntegrate[solK, {x, -a, a}, {y, -b, b},
Method -> {"MultidimensionalRule", "Generators" -> 5,
"SymbolicProcessing" -> 0}, PrecisionGoal -> 3, AccuracyGoal -> 3];
Hi! I want find deformation (strain) energy for plane by pressure .In one of the steps I need to use NIntegrate
function so problem with speed of my code when I take (a=6, b=6) it works too slowly! What can I do to inprove code?
Compile
can i do smth without it function (Because I'm not good at Parallel Computing) to make code faster or this is better way to useCompile
? $\endgroup$ – Arthuros Apr 23 '15 at 18:37solK
is symmetric, so you shouldn't integrate redundant terms. $\endgroup$ – george2079 Apr 23 '15 at 20:38