Question: Remove zeros after decimal point while retaining some numbers like r1=0.000065


 

restart; eq2 := ef*epsilon+x1*sinh(y1*L/r1)+x2*cosh(y1*L/r1)-.1*(ef*epsilon+x2) = 0; eq3 := ef*epsilon+x1*sinh(-y1*L/r1)+x2*cosh(-y1*L/r1)-.1*(ef*epsilon+x2) = 0; sol1 := solve({eq2, eq3}, {x1, x2})

.9*ef*epsilon+x1*sinh(y1*L/r1)+x2*cosh(y1*L/r1)-.1*x2 = 0

 

.9*ef*epsilon-x1*sinh(y1*L/r1)+x2*cosh(y1*L/r1)-.1*x2 = 0

 

{x1 = 0., x2 = -9.*ef*epsilon/(10.*cosh(y1*L/r1)-1.)}

(1)

assign(sol1); 1; x1, x2

0., -9.*ef*epsilon/(10.*cosh(y1*L/r1)-1.)

(2)

x1;

0.

(3)

x2;

-9.*ef*epsilon/(10.*cosh(y1*L/r1)-1.)

(4)

eq4 := ef*epsilon+x1*sinh(y1*z/r1)+x2*cosh(y1*z/r1);

ef*epsilon-9.*ef*epsilon*cosh(y1*z/r1)/(10.*cosh(y1*L/r1)-1.)

(5)

avg_eq4 := Typesetting:-delayDotProduct(1/L, int(eq4, z));

1/L.(ef*epsilon*z-9.*ef*epsilon*r1*sinh(y1*z/r1)/((10.*cosh(y1*L/r1)-1.)*y1))

(6)

Ecomp := Typesetting:-delayDotProduct(simplify((eval(subs(z = L, avg_eq4)-subs(z = 0, avg_eq4)))/epsilon), Vf)+(1-Vf)*Em;

(1/L.(ef*epsilon*(10.*L*y1*cosh(y1*L/r1)-L*y1-9.*r1*sinh(y1*L/r1))/((10.*cosh(y1*L/r1)-1.)*y1)))/epsilon.Vf+(1-Vf)*Em

(7)

tau := .5*(diff(eq4, z))*r1;

-4.5*ef*epsilon*sinh(y1*z/r1)*y1/(10.*cosh(y1*L/r1)-1.)

(8)

L := lfact*r1;

lfact*r1

 

0.65e-5

(9)

ef := 0.75e11;

0.75e11

 

0.326e10

 

.3

 

.5

(10)

Theo_Ecomp := Vf*ef+(1-Vf)*Em;

0.39130e11

(11)

y1 := 2*Em/((1+nu)*ef*ln(1/Vf))

0.9647560684e-1

(12)

Ecomp

0.5979989095e17*(1/lfact.(epsilon*(0.6270914445e-5*lfact*cosh(0.9647560685e-1*lfact)-0.6270914445e-6*lfact-0.585e-4*sinh(0.9647560685e-1*lfact))/(10.*cosh(0.9647560685e-1*lfact)-1.)))/epsilon+0.1630e10

(13)

epsilon := 0.1e-2;

0.1e-2

(14)

plot(Ecomp, lfact = 1 .. 500)

 

lfact := 90.81506;

90.81506

(15)

eq4

0.75e8-21149.64968*cosh(14842.40105*z)

(16)

plot(eq4, z = 0 .. L)

 

plot(tau, z = 0 .. L)

 

NULL


 

Download Fuly_bonded_updated.mw

Please Wait...