# Question:Fsolve

## Question:Fsolve

Maple

Good Morning,

Would anybody in the community please help me  to trouble shoot the problem with fsolve not producing result.

Regards,

> restart;
> e[y] := f[y]/E[s];
> f[s] := proc (t) options operator, arrow; piecewise(t <= -e[y], -f[y], t <= e[y], E[s]*t, f[y]) end proc;
> A[tot] := proc (e1, e2, As1, As2) options operator, arrow; As1+As2 end proc;
> M[max] := proc (As1, As2) options operator, arrow; alpha[1]*phi[c]*f[c]*b*(h-c[top])*(h-c[top])+A[s2]*(h-c[bot]-c[top])*f[y] end proc;
> c[neu] := proc (e2, e1) options operator, arrow; h*e2/(e2-e1) end proc;
> a := proc (e2, e1) options operator, arrow; beta[1]*c[neu](e2, e1) end proc;
> F[c2] := proc (e1, e2) options operator, arrow; piecewise(c[top] <= a(e2, e1), alpha[1]*phi[c]*f[c], 0) end proc;
> F[s1] := proc (e1, e2) options operator, arrow; phi[s]*f[s](e1+(e2-e1)*c[bot]/h) end proc;
> F[s2] := proc (e1, e2) options operator, arrow; phi[s]*f[s](e1+(e2-e1)*(h-c[top])/h) end proc;
> Steel[force] := proc (e1, e2, As1, As2) options operator, arrow; As1*F[s1](e1, e2)+As2*(F[s2](e1, e2)-F[c2](e1, e2)) end proc;
> Concrete[force] := proc (e1, e2) options operator, arrow; b*a(e2, e1)*alpha[1]*phi[c]*f[c] end proc;
> Steel[moment] := proc (e1, e2, As1, As2) options operator, arrow; c[bot]*As1*F[s1](e1, e2)+(h-c[top])*As2*(F[s2](e1, e2)-F[c2](e1, e2)) end proc;
> Concrete[moment] := proc (e2, e1) options operator, arrow; (h-(1/2)*a(e2, e1))*b*a(e2, e1)*alpha[1]*phi[c]*f[c] end proc;
> Total[force] := proc (e1, e2, As1, As2) options operator, arrow; Steel[force](e1, e2, As1, As2)+Concrete[force](e1, e2) = 0 end proc;
> Total[moment] := proc (e1, e2, As1, As2) options operator, arrow; Steel[moment](e1, e2, As1, As2)+Concrete[moment](e2, e1) = M[f] end proc;

Input Variables:

> b := 250.;
> d := 500.;
> As1 := 3000; As2 := 3000;
> f[c] := 25.;
> f[y] := 400.;
> epsilon[c] := 0.35e-2;
> phi[c] := .6;
> phi[s] := .85;
> E[s] := 200000.;
> c[bot] := 50;
> c[top] := 50;
> h := 575;
> NULL;
> NULL;
> M[f] := 0.290e9;
> beta[1] := eval(max(.97-0.25e-2*f[c], .67));
> alpha[1] := eval(max(.85-0.15e-2*f[c], .67));

> Sols := fsolve({Total[force], Total[moment]}, {As1, As2, e1 = -0.2e-2, e2 <= 0.35e-2});
> assign(Sols);
>

﻿