Question: Help ...Maple will not plot

Hi

II ve managed to build an expression with one variable 
I me trying to plot this expression on a defined range but maple doesn t not let me.

How do I manage to plot this expression ?

If you have any advise on how to improve my code I m open to comments 

 

Thanks a lot in advance
 

restart


#data

`ϵc1` := 2.1*10^(-3);

0.2100000000e-2

(1)

`ϵcu1` := 3.5*10^(-3):

Ecm := 31000:

fcm := 25:

Fy := 500:

Es := 200000:

#geometry

As := 1885:

b := 250:

d := 450:


#coefficients

eta := epsilon/`ϵc1`:

k := 1.05*Ecm*`ϵc1`/fcm:

NULL


#formula

 

NULL

`σc` := proc (epsilon) options operator, arrow; fcm*(k*epsilon/`ϵc1`-epsilon^2/`ϵc1`^2)/(1+(k-2)*epsilon/`ϵc1`) end proc;

proc (epsilon) options operator, arrow; fcm*(k*epsilon/`ϵc1`-epsilon^2/`ϵc1`^2)/(1+(k-2)*epsilon/`ϵc1`) end proc

(2)

plot(`σc`, 0 .. `ϵcu1`);

 

NULL

#Pressure as a function of y (y=`ϵc`*y/x) :

`σc` := proc (y) options operator, arrow; fcm*(k*`ϵc`*y/(x*`ϵc1`)-`ϵc`^2*y^2/(x^2*`ϵc1`^2))/(1+(k-2)*`ϵc`*y/(x*`ϵc1`)) end proc;

proc (y) options operator, arrow; fcm*(k*`ϵc`*y/(x*`ϵc1`)-`ϵc`^2*y^2/(x^2*`ϵc1`^2))/(1+(k-2)*`ϵc`*y/(x*`ϵc1`)) end proc

(3)

`assuming`([int(`&sigma;c`, 0 .. x)], [0 <= x and x <= (1/2)*d, 0 <= `&varepsilon;c` and `&varepsilon;c` <= 0.35e-2]);

int(25*(1302.000000*`&varepsilon;c`*x/x-226757.3696*`&varepsilon;c`^2*x^2/x^2)/(1+349.6190476*`&varepsilon;c`*x/x), x = 0 .. x)

(4)

C := `assuming`([int(fcm*(k*`&varepsilon;c`*s/(x*`&varepsilon;c1`)-(`&varepsilon;c`*s/(x*`&varepsilon;c1`))^2)*b/(1+(k-2)*`&varepsilon;c`*s/(x*`&varepsilon;c1`)), s = 0 .. x)], [0 < x and x <= (1/2)*d, 0 <= `&varepsilon;c` and `&varepsilon;c` <= 0.35e-2]);

-0.7487980799e-19*x*(0.2706771684e26*`&varepsilon;c`^2+0.1331955800e22*ln(874047619.*`&varepsilon;c`+2500000.)-0.1962210817e23-0.4656771182e24*`&varepsilon;c`)/`&varepsilon;c`

(5)

T := `assuming`([(d-x)*`&varepsilon;c`*Es*As/x], [0 < x and x <= (1/2)*d, 0 <= `&varepsilon;c` and `&varepsilon;c` <= 0.35e-2]);

377000000*(450-x)*`&varepsilon;c`/x

(6)

``

x := `assuming`([solve(C = T, x)], [0 < x and x <= (1/2)*d, 0 <= `&varepsilon;c` and `&varepsilon;c` <= 0.35e-2]);

868.0577815*(0.2900000000e12*`&varepsilon;c`+(-0.7297483689e24*`&varepsilon;c`^2+0.1400157114e23*`&varepsilon;c`-0.4004807874e20*ln(874047619.*`&varepsilon;c`+2500000.)+0.5899803379e21)^(1/2))*`&varepsilon;c`/(0.2706771684e13*`&varepsilon;c`^2+133195580.*ln(874047619.*`&varepsilon;c`+2500000.)-0.4656771182e11*`&varepsilon;c`-1962210817.), -868.0577815*(-0.2900000000e12*`&varepsilon;c`+(-0.7297483689e24*`&varepsilon;c`^2+0.1400157114e23*`&varepsilon;c`-0.4004807874e20*ln(874047619.*`&varepsilon;c`+2500000.)+0.5899803379e21)^(1/2))*`&varepsilon;c`/(0.2706771684e13*`&varepsilon;c`^2+133195580.*ln(874047619.*`&varepsilon;c`+2500000.)-0.4656771182e11*`&varepsilon;c`-1962210817.)

(7)

``

``

``

l

(8)

plot(l, `&varepsilon;c` = 0 .. 0.34e-2)

Error, (in plot) expected a range but received `&varepsilon;c` = 0 .. 0.34e-2

 

``

NULL


 

Download HW1_-_EC2_strain-pressure_graph.mw
 

restart


#data

`&varepsilon;c1` := 2.1*10^(-3);

0.2100000000e-2

(1)

`&varepsilon;cu1` := 3.5*10^(-3):

Ecm := 31000:

fcm := 25:

Fy := 500:

Es := 200000:

#geometry

As := 1885:

b := 250:

d := 450:


#coefficients

eta := epsilon/`&varepsilon;c1`:

k := 1.05*Ecm*`&varepsilon;c1`/fcm:

NULL


#formula

 

NULL

`&sigma;c` := proc (epsilon) options operator, arrow; fcm*(k*epsilon/`&varepsilon;c1`-epsilon^2/`&varepsilon;c1`^2)/(1+(k-2)*epsilon/`&varepsilon;c1`) end proc;

proc (epsilon) options operator, arrow; fcm*(k*epsilon/`&varepsilon;c1`-epsilon^2/`&varepsilon;c1`^2)/(1+(k-2)*epsilon/`&varepsilon;c1`) end proc

(2)

plot(`&sigma;c`, 0 .. `&varepsilon;cu1`);

 

NULL

#Pressure as a function of y (y=`&varepsilon;c`*y/x) :

`&sigma;c` := proc (y) options operator, arrow; fcm*(k*`&varepsilon;c`*y/(x*`&varepsilon;c1`)-`&varepsilon;c`^2*y^2/(x^2*`&varepsilon;c1`^2))/(1+(k-2)*`&varepsilon;c`*y/(x*`&varepsilon;c1`)) end proc;

proc (y) options operator, arrow; fcm*(k*`&varepsilon;c`*y/(x*`&varepsilon;c1`)-`&varepsilon;c`^2*y^2/(x^2*`&varepsilon;c1`^2))/(1+(k-2)*`&varepsilon;c`*y/(x*`&varepsilon;c1`)) end proc

(3)

`assuming`([int(`&sigma;c`, 0 .. x)], [0 <= x and x <= (1/2)*d, 0 <= `&varepsilon;c` and `&varepsilon;c` <= 0.35e-2]);

int(25*(1302.000000*`&varepsilon;c`*x/x-226757.3696*`&varepsilon;c`^2*x^2/x^2)/(1+349.6190476*`&varepsilon;c`*x/x), x = 0 .. x)

(4)

C := `assuming`([int(fcm*(k*`&varepsilon;c`*s/(x*`&varepsilon;c1`)-(`&varepsilon;c`*s/(x*`&varepsilon;c1`))^2)*b/(1+(k-2)*`&varepsilon;c`*s/(x*`&varepsilon;c1`)), s = 0 .. x)], [0 < x and x <= (1/2)*d, 0 <= `&varepsilon;c` and `&varepsilon;c` <= 0.35e-2]);

-0.7487980799e-19*x*(0.2706771684e26*`&varepsilon;c`^2+0.1331955800e22*ln(874047619.*`&varepsilon;c`+2500000.)-0.1962210817e23-0.4656771182e24*`&varepsilon;c`)/`&varepsilon;c`

(5)

T := `assuming`([(d-x)*`&varepsilon;c`*Es*As/x], [0 < x and x <= (1/2)*d, 0 <= `&varepsilon;c` and `&varepsilon;c` <= 0.35e-2]);

377000000*(450-x)*`&varepsilon;c`/x

(6)

``

x := `assuming`([solve(C = T, x)], [0 < x and x <= (1/2)*d, 0 <= `&varepsilon;c` and `&varepsilon;c` <= 0.35e-2]);

868.0577815*(0.2900000000e12*`&varepsilon;c`+(-0.7297483689e24*`&varepsilon;c`^2+0.1400157114e23*`&varepsilon;c`-0.4004807874e20*ln(874047619.*`&varepsilon;c`+2500000.)+0.5899803379e21)^(1/2))*`&varepsilon;c`/(0.2706771684e13*`&varepsilon;c`^2+133195580.*ln(874047619.*`&varepsilon;c`+2500000.)-0.4656771182e11*`&varepsilon;c`-1962210817.), -868.0577815*(-0.2900000000e12*`&varepsilon;c`+(-0.7297483689e24*`&varepsilon;c`^2+0.1400157114e23*`&varepsilon;c`-0.4004807874e20*ln(874047619.*`&varepsilon;c`+2500000.)+0.5899803379e21)^(1/2))*`&varepsilon;c`/(0.2706771684e13*`&varepsilon;c`^2+133195580.*ln(874047619.*`&varepsilon;c`+2500000.)-0.4656771182e11*`&varepsilon;c`-1962210817.)

(7)

``

``

``

l

(8)

plot(l, `&varepsilon;c` = 0 .. 0.34e-2)

Error, (in plot) expected a range but received `&varepsilon;c` = 0 .. 0.34e-2

 

``

NULL


 

Download HW1_-_EC2_strain-pressure_graph.mw

 

 

Please Wait...