Commit 028bf2b2 authored by Turnhout, M.C. van's avatar Turnhout, M.C. van
Browse files

fix axisymmetric correction

parent 5820eb4c
......@@ -108,6 +108,7 @@ for int = 1:nint
B(3, ii+1) = dndxy(:, 2)';
B(3, ii+2) = dndxy(:, 1)';
if axi == 1 % axi-symmetric correction
x = n(int, :)*nodcoord(:, 1); % ipts x-coordinate (r)
B(4, ii+1) = n(int, :)/x;
end
......
......@@ -108,6 +108,7 @@ for int = 1:nint
B(3, ii+1) = dndxy(:, 2)';
B(3, ii+2) = dndxy(:, 1)';
if axi == 1 % axi-symmetric correction
x = n(int, :)*nodcoord(:, 1); % ipts x-coordinate (r)
B(4, ii+1) = n(int, :)/x;
end
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment