Advanced Search

pooyan1990

35 Reputation

3 Badges

5 years, 129 days
Contact pooyan1990

MaplePrimes Activity


These are questions asked by pooyan1990

Problem in solving a differential equati...

pooyan1990 35 Maple 2019
syntax dsolve differential-equations
April 05 2023
1 1

I receive the following error when I run maple to solve a diff equation.
Thank you in advance for your help.

> 

eq1:=A__n=[X__n+[(C-z)*(tan(theta1)+tan(theta2))]] * [B + (2*(C-z)*tan(theta3))];
eq2:=P__n=[2*(X__n+B)] + 2 * (C-z) * [(tan(theta1)+tan(theta2)) + (2*tan(theta3))];
eq3:=gamma - [cu*(P__n/A__n)] = diff(sigma__v(z),z);
eq4:=subs([eq1,eq2],eq3);
ics:=sigma__v(0)=sigma__s;
sol:=simplify(dsolve(([eq4,ics]), sigma__v(z)));
NULL

A__n = [X__n+[(C-z)*(tan(theta1)+tan(theta2))]]*[B+2*(C-z)*tan(theta3)]

  

P__n = [2*X__n+2*B]+2*(C-z)*[tan(theta1)+tan(theta2)+2*tan(theta3)]

  

gamma+[-cu*P__n/A__n] = diff(sigma__v(z), z)

  

gamma+[-cu*([2*X__n+2*B]+2*(C-z)*[tan(theta1)+tan(theta2)+2*tan(theta3)])/([X__n+[(C-z)*(tan(theta1)+tan(theta2))]]*[B+2*(C-z)*tan(theta3)])] = diff(sigma__v(z), z)

  

sigma__v(0) = sigma__s

  

Error, (in PDEtools/NumerDenom) invalid input: `PDEtools/NumerDenom` expects its 1st argument, ee, to be of type algebraic, but received [tan(theta1)+tan(theta2)+2*tan(theta3)]

  
> 

``

Download diffeq.mw

Can't solve a differential equation...

pooyan1990 35 Maple 2019
syntax dsolve differential-equations
November 03 2019
1 8

I have an ordinary differential equation to be solved. when I give numeric values to B,H and L__1 (for example B=10,H=10,L__1=4), everything works fine and equations eq37a will be solved.

but when B,H and L__1 are symbolic parameters, dsolve give me this error " Error, (in dsolve) give the main variable as a second argument". and I don't know what it means.

is there any way to solve the ODE with symbolic B,H and L__1??

 

> 

restart;

#B:=10;
#H:=10;
#L__1:=4;

eq0:=R__i=(1/2)*B+L__1;
eq01:=R__f=C+(1/2)*B+L__1;
eq02:=theta__1=arctan((H/2)/(B/2 + L__1));
eq03:=theta__2=Pi - arctan((H/2)/(B/2 - L__1));
eq43:= v__beta(r,theta,beta)= v__m(beta)*(1-(r^2/((R__max(theta,beta))^2)));
eq39:=R(beta)=R__i+2*beta*(R__f-R__i)/Pi;
eq39a:=simplify(subs([eq0,eq01], eq39));
eq38:=R__max(theta,beta)=R__max(theta,0)*(R(beta)/R__i);
eq41:=R__max(theta,0)=piecewise(theta<=rhs(eq02),(B/2+L__1)/cos(theta),theta<=rhs(eq03),(H/2)/sin(theta),rhs(eq03)<theta,(B/2-L__1)/cos(Pi-theta));
eq38a:= simplify(subs([eq41,eq39a,eq0], eq38));
eq44 := v__m(beta) = 1/(2*int((1 - r^2/R__max(theta, beta)^2)*r, [r = 0 .. R__max(theta, beta), theta = 0 .. Pi]));
eq44a:=simplify(subs(eq38a,eq44)) assuming B/L__1 >= 2 , H > 0 , B > 0 , L__1 > 0;
eq44b:= v__m(beta)= rhs(eq44a) * B * H;
eq47:=subs([eq44b,eq38a],eq43);

eq0 := R__i = (1/2)*B+L__1

  

eq01 := R__f = C+(1/2)*B+L__1

  

eq02 := `&theta;__1` = arctan(H/(2*((1/2)*B+L__1)))

  

eq03 := `&theta;__2` = Pi-arctan(H/(2*((1/2)*B-L__1)))

  

eq43 := `v__&beta;`(r, theta, beta) = v__m(beta)*(1-r^2/R__max(theta, beta)^2)

  

eq39 := R(beta) = R__i+2*beta*(R__f-R__i)/Pi

  

eq39a := R(beta) = (1/2)*B+L__1+2*beta*C/Pi

  

eq38 := R__max(theta, beta) = R__max(theta, 0)*R(beta)/R__i

  

eq41 := R__max(theta, 0) = piecewise(theta <= arctan(H/(2*((1/2)*B+L__1))), ((1/2)*B+L__1)/cos(theta), theta <= Pi-arctan(H/(2*((1/2)*B-L__1))), H/(2*sin(theta)), Pi-arctan(H/(2*((1/2)*B-L__1))) < theta, -((1/2)*B-L__1)/cos(theta))

  

eq38a := R__max(theta, beta) = (Pi*(B+2*L__1)+4*beta*C)*piecewise(theta <= arctan(H/(B+2*L__1)), (B+2*L__1)/cos(theta), theta <= Pi-arctan(H/(B-2*L__1)), H/sin(theta), Pi-arctan(H/(B-2*L__1)) < theta, (-B+2*L__1)/cos(theta))/(2*Pi*(B+2*L__1))

  

eq44 := v__m(beta) = 2/(int(R__max(theta, beta)^2, theta = 0 .. Pi))

  

eq44a := v__m(beta) = 2*Pi^2*(B+2*L__1)^2/(H*(B*Pi+4*C*beta+2*L__1*Pi)^2*B)

  

eq44b := v__m(beta) = 2*Pi^2*(B+2*L__1)^2/(B*Pi+4*C*beta+2*L__1*Pi)^2

  

v__beta(r, theta, beta) = 2*Pi^2*(B+2*L__1)^2*(1-4*r^2*Pi^2*(B+2*L__1)^2/((Pi*(B+2*L__1)+4*beta*C)^2*piecewise(theta <= arctan(H/(B+2*L__1)), (B+2*L__1)/cos(theta), theta <= Pi-arctan(H/(B-2*L__1)), H/sin(theta), Pi-arctan(H/(B-2*L__1)) < theta, (-B+2*L__1)/cos(theta))^2))/(B*Pi+4*C*beta+2*L__1*Pi)^2

 (1)
> 

eq37:=diff(v__r(r, theta, beta),r)*r*(R__f-r*cos(theta))+v__r(r,theta,beta)*(R__f-2*r*cos(theta))+r*diff(v__beta(r,theta,beta),beta)=0;
eq37a:=simplify(expand(subs(eq47, eq37)));

eq37 := (diff(v__r(r, theta, beta), r))*r*(R__f-r*cos(theta))+v__r(r, theta, beta)*(R__f-2*r*cos(theta))+r*(diff(`v__&beta;`(r, theta, beta), beta)) = 0

  

(-2*(Pi*(B+2*L__1)+4*beta*C)^2*((1/2)*(Pi*(B+2*L__1)+4*beta*C)^3*(r*cos(theta)-R__f)*r*(diff(v__r(r, theta, beta), r))+(Pi*(B+2*L__1)+4*beta*C)^3*(r*cos(theta)-(1/2)*R__f)*v__r(r, theta, beta)+8*C*r*Pi^2*(B+2*L__1)^2)*piecewise(theta <= arctan(H/(B+2*L__1)), (B+2*L__1)/cos(theta), theta <= Pi-arctan(H/(B-2*L__1)), H/sin(theta), Pi-arctan(H/(B-2*L__1)) < theta, (-B+2*L__1)/cos(theta))^2+128*C*r^3*Pi^4*(B+2*L__1)^4)/((Pi*(B+2*L__1)+4*beta*C)^5*piecewise(theta <= arctan(H/(B+2*L__1)), (B+2*L__1)/cos(theta), theta <= Pi-arctan(H/(B-2*L__1)), H/sin(theta), Pi-arctan(H/(B-2*L__1)) < theta, (-B+2*L__1)/cos(theta))^2) = 0

 (2)
> 

bc:=v__r(rhs(eq38a),theta,beta)=0;

bc := v__r((Pi*(B+2*L__1)+4*beta*C)*piecewise(theta <= arctan(H/(B+2*L__1)), (B+2*L__1)/cos(theta), theta <= Pi-arctan(H/(B-2*L__1)), H/sin(theta), Pi-arctan(H/(B-2*L__1)) < theta, (-B+2*L__1)/cos(theta))/(2*Pi*(B+2*L__1)), theta, beta) = 0

 (3)
> 

sol3:=simplify(dsolve([eq37a,bc], v__r(r, theta,beta)));

Error, (in dsolve) give the main variable as a second argument

  
> 

NULL


 

Download Mechrect.mw

Can't solve an integral...

pooyan1990 35 Maple 2019
integration
October 29 2019
1 1

Maple can't solve the following symbolic integration for me and it remains in the state of evaluating when I use simplify and nothing happens after hours.. any solution? 

I was thinking maybe there is something wrong with my maple !!?

 
> 

eq0:=R__i=(1/2)*d+L__1;
eq00:=R__f=c+(1/2)*d+L__1;
eq39:=R(beta)=R__i+2*beta*(R__f-R__i)/Pi;
eq39a:=simplify(subs([eq0,eq00], eq39));
eq38:=R__max(theta,beta)=R__max(theta,0)*(R(beta)/R__i);
eq41:=R__max(theta,0)=((2*L__1*cos(theta)) + sqrt(d^2 - 2*(L__1)^2 + 2*(L__1)^2*cos(2*theta)))/2;
eq38a:= simplify(subs([eq41,eq39a,eq0], eq38));

R__i = L__1+(1/2)*d

  

R__f = c+(1/2)*d+L__1

  

R(beta) = R__i+2*beta*(R__f-R__i)/Pi

  

R(beta) = L__1+(1/2)*d+2*beta*c/Pi

  

R__max(theta, beta) = R__max(theta, 0)*R(beta)/R__i

  

R__max(theta, 0) = L__1*cos(theta)+(1/2)*(d^2-2*L__1^2+2*L__1^2*cos(2*theta))^(1/2)

  

R__max(theta, beta) = (1/2)*(2*L__1*cos(theta)+(d^2-2*L__1^2+2*L__1^2*cos(2*theta))^(1/2))*(2*L__1*Pi+Pi*d+4*beta*c)/(Pi*(2*L__1+d))

 (1)
> 

eq44 := v__m(beta) = 1/int((1 - r^2/R__max(theta, beta)^2)*r, [r = 0 .. R__max(theta, beta), theta = 0 .. 2*Pi]);

v__m(beta) = 4/(int(R__max(theta, beta)^2, theta = 0 .. 2*Pi))

 (2)
> 

eq44a:=simplify(subs(eq38a,eq44));
> 

``

Download int.mw

 

Problem in Solving a first-order linear ...

pooyan1990 35 Maple 2017
October 17 2019
2 8

Hello,

I have a problem in solving a differential equation in maple. I got an error which is "unexpected occurence of the variables..." .  I have uploaded my full problem as a pdf and maple worksheet (mw) files.

Untitled_(4).pdf

Prob.mw

Problem in Solving integral in Maple...

pooyan1990 35 Maple 2017
integration syntax numeric
September 11 2019
1 2

Hello,

I have a problem in solving an integral in maple. I can't solve the below integral in maple and it returns the integral itself to me. I also attach an image from the integral if here is not clearly shown. I want maple to return me just a number. can anyone help me in this?

Thank you

int(sin(beta)*(-0.4447569104e-1*beta(10)^3+1.846983291*beta(10)^2+78.88888890*beta(10)+620.4645491)/(9.+.6366197724*beta(10))^2, beta = 0 .. (1/2)*Pi)

 

 

 

 

Page 1 of 1