Why does my DSolve command does not solve the Non-linear ODEs?
Clear["Global`*"]
(*Parameters*)
length = 1.5;
width = 0.25;
height = 0.25;
area = width*height;
cforce = -400;
dforce = (500*(X1^2));
stress = 10*displacement'[X1] + 100000*displacement'[X1]^3;
DE1 = D[stress*area, X1] + dforce;
DE2 = 10*displacement'[X1] + 100000*displacement'[X1]^3 + 6400;
(* Solving for the Exact Displacement Function *)
solution =
DSolve[{DE1 == 0, displacement[0] == 0, DE2 == 0 /. X1 -> length},
displacement[X1], X1];
displacement = displacement[X1] /. solution[[1]];
Print["The Exact Displacement Function of the Bar is"]
Print[displacement];
Print[Plot[{displacement}, {X1, 0, 1.5}, AspectRatio -> 0.4,
AxesLabel -> {"Bar Length (m)", "Axial Displacement (m)"},
PlotLegends -> {"Exact Displacement of the Bar"}]];
Print[" "]