kill(all)$globalsolve:true$

eq:[ε_x=σ_x/E-ν*σ_y/E-ν*σ_z/E,
ε_y=σ_y/E-ν*σ_z/E-ν*σ_x/E,
ε_z=σ_z/E-ν*σ_x/E-ν*σ_y/E];

sol:linsolve(eq,[ε_x,σ_y,σ_z]);

sol1:ratsimp(ev(sol,ε_y:0,ε_z:0));

ev(sol1,ν:3/10);

σ_x*1/(sol1[1]);