theta = linspace (0,90,18);
Vxy = zeros(size(theta));
[Ex(i)] = E_longitudinal (E1, E2, G12, v12, theta(i));
[Vxy(i)] = Vxy_longitudinal (E1, E2, G12, v12, Ex(i), theta(i));
plot(theta, Ex/E2, '-+r');
xlabel ('\theta (degree)');
ylabel ({'E_x/E_2', '|Vxy|'});
legend ('E_x/E_2', 'Vxy');
pLeft = plot(theta, Ex/E2, 'bs', theta, Gxy/G12, 'ko', theta, -nxy, 'g+');
Unrecognized function or variable 'Gxy'.
pRight = plot(theta, Vxy, 'rs');
xlabel('\theta (degree)');
ylabel('Ex/E_2, Gxy/G12, -\etaxy, x');
legend(h, 'Ex/E2', 'Gxy/G12', '-\etaxy,x', '\nuxy');
function [Ex] = E_longitudinal (E1, E2, G12, v12, theta)
E_Inverse = (1/E1)*(cosd(theta)^4) + (1/G12 - 2*v12/E1)*(sind(theta)^2)*(cosd(theta)^2) + (1/E2)*(sind(theta)^4);
function [Vxy] = Vxy_longitudinal (E1, E2, G12, V12, Ex, theta)
Vxy = Ex*((V12/E1)*((sind(theta))^4+(cosd(theta))^4)-((1/E1)+(1/E2)-(1/G12))*(sind(theta))^2*(cosd(theta))^2);
I cannot fix the Gxy problem as you have no code anywhere that defiens Gxy.
Best Answer