Advanced Search

Question: SOLVING PDE EQUATION

Question:SOLVING PDE EQUATION

9009134 110 Maple 2019
homework pde numeric differential-equations galerkin finite-element
July 21 2019
0

Hello , Anybody can help me to write teh codes for teh PDE to solve by Galerkin finite element method or any other methods can be able to gain results? the parameter omega is unknown and should be determined.

I attached a pdf file for more detais.

Thanks so much

fem2.mw
 

> 

"restart: rho:=7850: E:=0.193e12: n:=1: AD:=10: upsilon:=0.291: mu:=E/(2*(1+upsilon)): l:=0: lambda:=E*upsilon/((1+upsilon)*(1-2*upsilon)): R:=2.5: ii:=2: J:=2: m:=1: `u__theta`(r,theta,phi):= ( V(r,theta))*cos(m*phi): `u__r`(r,theta,phi):= ( U(r,theta))*cos(m*phi): `u__phi`(r,theta,phi):= ( W(r,theta))*sin(m*phi): : eq1:=(r (R+r cos(theta))^2 (mu+lambda) (((∂)^2)/(∂r∂theta) `u__theta`(r,theta,phi))+2 r^2 (mu+lambda/2) (R+r cos(theta))^2 (((∂)^2)/(∂r^2) `u__r`(r,theta,phi))+r^2 (mu+lambda) (R+r cos(theta)) (((∂)^2)/(∂phi∂r) `u__phi`(r,theta,phi))+mu (R+r cos(theta))^2 (((∂)^2)/(∂theta^2) `u__r`(r,theta,phi))+(((∂)^2)/(∂phi^2) `u__r`(r,theta,phi)) mu r^2-3 (R+r cos(theta))^2 (mu+lambda/3) ((∂)/(∂theta) `u__theta`(r,theta,phi))+2 r (mu+lambda/2) (R+2 r cos(theta)) (R+r cos(theta)) ((∂)/(∂r) `u__r`(r,theta,phi))-r^2 sin(theta) (mu+lambda) (R+r cos(theta)) ((∂)/(∂r) `u__theta`(r,theta,phi))-3 r^2 cos(theta) (mu+lambda/3) ((∂)/(∂phi) `u__phi`(r,theta,phi))-r mu sin(theta) (R+r cos(theta)) ((∂)/(∂theta) `u__r`(r,theta,phi))-2 (mu+lambda/2) (2 (cos(theta))^2 r^2+2 cos(theta) R r+R^2) `u__r`(r,theta,phi)+r `u__theta`(r,theta,phi) sin(theta) (3 r (mu+lambda/3) cos(theta)+R mu))/(r^2 (R+r cos(theta))^2): eq2:=(2 (mu+lambda/2) (R+r cos(theta))^2 (((∂)^2)/(∂theta^2) `u__theta`(r,theta,phi))+r (R+r cos(theta))^2 (mu+lambda) (((∂)^2)/(∂r∂theta) `u__r`(r,theta,phi))+r (mu+lambda) (R+r cos(theta)) (((∂)^2)/(∂phi∂theta) `u__phi`(r,theta,phi))+r^2 mu (R+r cos(theta))^2 (((∂)^2)/(∂r^2) `u__theta`(r,theta,phi))+(((∂)^2)/(∂phi^2) `u__theta`(r,theta,phi)) mu r^2+3 (R+r cos(theta)) ((4 r (mu+lambda/2) cos(theta))/3+R (mu+lambda/3)) ((∂)/(∂theta) `u__r`(r,theta,phi))-2 r (mu+lambda/2) sin(theta) (R+r cos(theta)) ((∂)/(∂theta) `u__theta`(r,theta,phi))+r mu (R+2 r cos(theta)) (R+r cos(theta)) ((∂)/(∂r) `u__theta`(r,theta,phi))+3 r^2 sin(theta) (mu+lambda/3) ((∂)/(∂phi) `u__phi`(r,theta,phi))+(-3 r R (mu+lambda/3) cos(theta)+(-lambda-2 mu) r^2-R^2 mu) `u__theta`(r,theta,phi)-2 r (mu+lambda/2) sin(theta) R `u__r`(r,theta,phi))/(r^2 (R+r cos(theta))^2): eq3:=(r (mu+lambda) (R+r cos(theta)) (((∂)^2)/(∂phi∂theta) `u__theta`(r,theta,phi))+r^2 (mu+lambda) (R+r cos(theta)) (((∂)^2)/(∂phi∂r) `u__r`(r,theta,phi))+mu (R+r cos(theta))^2 (((∂)^2)/(∂theta^2) `u__phi`(r,theta,phi))+r (r mu (R+r cos(theta))^2 (((∂)^2)/(∂r^2) `u__phi`(r,theta,phi))+2 r (mu+lambda/2) (((∂)^2)/(∂phi^2) `u__phi`(r,theta,phi))+(4 r (mu+lambda/2) cos(theta)+R (mu+lambda)) ((∂)/(∂phi) `u__r`(r,theta,phi))+mu (R+2 r cos(theta)) (R+r cos(theta)) ((∂)/(∂r) `u__phi`(r,theta,phi))-mu sin(theta) (R+r cos(theta)) ((∂)/(∂theta) `u__phi`(r,theta,phi))-r (3 sin(theta) (mu+lambda/3) ((∂)/(∂phi) `u__theta`(r,theta,phi))+`u__phi`(r,theta,phi) mu)))/(r^2 (R+r cos(theta))^2): "
> 

EQ1 := collect(eq1, cos(m*phi))/cos(m*phi)+rho*omega^2; EQ2 := collect(eq2, cos(m*phi))/cos(m*phi)+rho*omega^2; EQ3 := collect(eq3, sin(m*phi))/sin(m*phi)+rho*omega^2

(0.1788235818e12*r*(2.5+r*cos(theta))^2*(diff(diff(V(r, theta), r), theta))+0.2535718390e12*r^2*(2.5+r*cos(theta))^2*(diff(diff(U(r, theta), r), r))+0.1788235818e12*r^2*(2.5+r*cos(theta))*(diff(W(r, theta), r))+0.7474825716e11*(2.5+r*cos(theta))^2*(diff(diff(U(r, theta), theta), theta))-0.7474825716e11*U(r, theta)*r^2-0.3283200960e12*(2.5+r*cos(theta))^2*(diff(V(r, theta), theta))+0.2535718390e12*r*(2.5+2.*r*cos(theta))*(2.5+r*cos(theta))*(diff(U(r, theta), r))-0.1788235818e12*r^2*sin(theta)*(2.5+r*cos(theta))*(diff(V(r, theta), r))-0.3283200960e12*r^2*cos(theta)*W(r, theta)-0.7474825716e11*r*sin(theta)*(2.5+r*cos(theta))*(diff(U(r, theta), theta))-0.2535718390e12*(2.*cos(theta)^2*r^2+5.0*r*cos(theta)+6.25)*U(r, theta)+r*V(r, theta)*sin(theta)*(0.3283200960e12*r*cos(theta)+0.1868706429e12))/(r^2*(2.5+r*cos(theta))^2)+7850*omega^2

  

(0.2535718390e12*(2.5+r*cos(theta))^2*(diff(diff(V(r, theta), theta), theta))+0.1788235818e12*r*(2.5+r*cos(theta))^2*(diff(diff(U(r, theta), r), theta))+0.1788235818e12*r*(2.5+r*cos(theta))*(diff(W(r, theta), theta))+0.7474825716e11*r^2*(2.5+r*cos(theta))^2*(diff(diff(V(r, theta), r), r))-0.7474825716e11*V(r, theta)*r^2+3.*(2.5+r*cos(theta))*(0.1690478927e12*r*cos(theta)+0.2736000800e12)*(diff(U(r, theta), theta))-0.2535718390e12*r*sin(theta)*(2.5+r*cos(theta))*(diff(V(r, theta), theta))+0.7474825716e11*r*(2.5+2.*r*cos(theta))*(2.5+r*cos(theta))*(diff(V(r, theta), r))+0.3283200960e12*r^2*sin(theta)*W(r, theta)+(-0.8208002400e12*r*cos(theta)-0.2535718389e12*r^2-0.4671766072e12)*V(r, theta)-0.6339295976e12*r*sin(theta)*U(r, theta))/(r^2*(2.5+r*cos(theta))^2)+7850*omega^2

  

(-0.1788235818e12*r*(2.5+r*cos(theta))*(diff(V(r, theta), theta))-0.1788235818e12*r^2*(2.5+r*cos(theta))*(diff(U(r, theta), r))+0.7474825716e11*(2.5+r*cos(theta))^2*(diff(diff(W(r, theta), theta), theta))+r*(0.7474825716e11*r*(2.5+r*cos(theta))^2*(diff(diff(W(r, theta), r), r))-0.2535718390e12*r*W(r, theta)-1.*(0.5071436780e12*r*cos(theta)+0.4470589545e12)*U(r, theta)+0.7474825716e11*(2.5+2.*r*cos(theta))*(2.5+r*cos(theta))*(diff(W(r, theta), r))-0.7474825716e11*sin(theta)*(2.5+r*cos(theta))*(diff(W(r, theta), theta))-1.*r*(-0.3283200960e12*sin(theta)*V(r, theta)+0.7474825716e11*W(r, theta))))/(r^2*(2.5+r*cos(theta))^2)+7850*omega^2

 (1)
> 

#BCs can be selcted from following Bcs
``
U(0, theta) = 0, (D[1](U))(0, theta) = 0, U(1, theta) = 0, (D[1](U))(1, theta) = 0

U(0, theta) = 0, (D[1](U))(0, theta) = 0, U(1, theta) = 0, (D[1](U))(1, theta) = 0

 (2)
> 

NULL
V(0, theta) = 0, (D[1](V))(0, theta) = 0, V(1, theta) = 0, (D[1](V))(1, theta) = 0
NULL
W(0, theta) = 0, (D[1](W))(0, theta) = 0, W(1, theta) = 0, (D[1](W))(1, theta) = 0
``

V(0, theta) = 0, (D[1](V))(0, theta) = 0, V(1, theta) = 0, (D[1](V))(1, theta) = 0

  

W(0, theta) = 0, (D[1](W))(0, theta) = 0, W(1, theta) = 0, (D[1](W))(1, theta) = 0

 (3)
> 

``


 

Download fem2.mw

buchanan2005.pdf

 

 
Please Wait...