damage equation for uniaxial compression

(4αβ^2 E^(0^2 ) ε_1^2)/〖((β+ γ ) E^0-(1+α)k)〗^2 + ε_1^2-(ε_u^2 k)/(w(E^0-k)) =0

given values for above equation β= 0.2777 γ= 1.111 α=1.13885 E^0=30768.7 MPa START FROM 0.00000 TO 2.00000 WITH INCREAMENT O.00005 ALWAYS ε_u=0.0003 W=0.0049 K=????? AT DIFFERENT VALUES OF ε_1

I need to solve the above equation. Thank you for ur support

Answer 1

[Roger Stafford]

By multiplying by the product of the two denominators, the equation becomes that of a cubic polynomial in k which will in general possess three roots. Which one do you want? You can use 'root's to find all three.

The variable you designate by "E^0" would necessarily have a value of 1 in matlab, since it represents the zero-th power of E. You need to use some other designation. The notation "E^(0^2 )" is even more dubious.