Fourth Order Hyperbolic PDE's

Question

David Koenig 2012-8-17

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/46243-fourth-order-hyperbolic-pde-s

评论： HAMED AKRAMI 2021-9-1

I have been struggling with trying to solve the vibrating plate equation in polar coordinates with angular dependence (non symmetrical initial condition). That is,

c d2y/dt2 = d4y/dr4

except that d4y/dr4 is just short cut symbolism for the double Laplacian in polar coordinates. I have been successful using a relatively simple explicit finite difference approach for symmetrical initial conditions but the asymmetrical IC's cause a problem. I have been holding off trying to reformulate an implicit finite difference method because there is so much algebraic bookkeeping.

I see that the partial differential solver toolbox does not deal with this problem. Does anyone know about Matlab scripts that deal with this problem?

Thanks.

3 个评论
显示 1更早的评论隐藏 1更早的评论

HAMED AKRAMI 2021-9-1

Hi

I have trouble solving 4th degree PDE equations in MATLAB

Please email me if you can help

hamed.akrami.abar@gmail.com

HAMED AKRAMI 2021-9-1

hi alex ; Email me the solution of this equation in MATLAB with full details

hamed.akrami.abar@gmail.com

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Deepak Ramaswamy 2012-10-17

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/46243-fourth-order-hyperbolic-pde-s#answer_62304

Hi

Don't know if you're still looking to solve this problem. The link below shows a static plate example of converting a 4th order PDE to a 2nd order PDE in the PDE Toolbox

http://www.mathworks.com/matlabcentral/fileexchange/34870-deflection-of-a-square-plate-calculated-with-pde-toolbox/content/html/clampedSquarePlateExample.html

The trick shown in the example can be extended to the hyperbolic problem as well. You might also want to consider using the eigensolver in the PDE Toolbox if you're interested in looking at modes.

Deepak

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

David Koenig 2014-10-15

Deepak, I am taking the brute force approach of replacing the partials with finite differences and then developing an explicit numerical method from the result. It seems to work OK but, of course, there is some serious bookkeeping required.

On the other hand, your suggestion of converting to two 2nd order pdes is interesting and I will pursue it. Thanks.

Dave

请先登录，再进行评论。