Skip to main content
P
payman1-Visitor
1-Visitor
July 3, 2019
Solved

Rectangular plate vibration in Lame mode

  • July 3, 2019
  • 7 replies
  • 5658 views

Hello everyone,

I'm trying to solve a system of partial differential equations corresponding to in-plane vibration of a thin plate in Lame mode. This is one of a papular resonance mode in microelectromechanical systems (MEMS) devices. An example of the displacement profile of Lame mode is shown in the attached worksheet.

 

For the boundary conditions, the plate is fixed at the four corners, so the displacement U(x,y) and V(x,y) along x and y-axes are zero (in some references, also the center is fixed). On the sides, the normal components of stresses are also zero. 

Gathering all the boundary conditions, I'm still unable to proceed with "Pdesolve"  routine and get the displacement and natural frequency.

 

Could anyone make a suggestion for resolving this problem? I highly appreciate your assistance.

Payman

Best answer by terryhendicott

Hi,

 

This has taken a long time to verify against a commercial software suite Strand 7.

Enclosed is 3 zip files zipped into an over arching archive.

PlateVibration.zip are Mathcad worksheets.

References.zip are PDF references that show how finite elements work and some examples

Strand7_Files.zip is the commercial package files.  A free viewer is available at 

The files with extension lsl and nfl are log files and are text.  nfl standing for natural frequency log shows the results.

 

In summary the Mathcad worksheets perform natural frequency analysis of a rectangular plate with varying support conditions and can plot mode shapes for each natural frequency.

Capture.PNG

It has reasonable accuracy.

Capture-2.PNG

Enjoy

Terry

 

 

 

 

 

 

 

7 replies

T
21-Topaz II
terryhendicott21-Topaz II
21-Topaz II
July 4, 2019

Hi,

" I'm still unable to proceed with "Pdesolve"  routine and get the displacement and natural frequency."

 

What you describe is two problems.  Displacement under load.   Vibration at natural frequency where you solve for the frequency but the displacements values at that frequency are unknown but the displacement pattern is known.  I notice you define no "loads" for displacement under load.

 

I have the solution for you for natural frequency supported at the four corners by finite element analysis but it will take a day or so to transcribe from a worksheet that solves for displacement under load.

 

In terms of your Mathcad solution you need to define "U" and "V".

Capture.PNG

Capture-2.PNG V is defined as default units of Volt so you get this message.

Cheers

Terry

P
payman1-VisitorAuthor
1-Visitor
July 4, 2019

Thank you Terry, I appreciate your help and I'm eager to receive your hints.

 

All the best,

Payman.

PS: I modified the worksheet and posted the new one.

T
21-Topaz II
terryhendicott21-Topaz II
21-Topaz II
July 4, 2019

Hi,

 

Does the reference you are using dimension the thickness of the plate as well as the plate plan dimensions 700 x 10^-6?

 

Cheers

Terry

-
-MFra-21-Topaz II
21-Topaz II
July 5, 2019

 

 

T
21-Topaz II
terryhendicott21-Topaz II
21-Topaz II
July 7, 2019

Hi,

 

I need a little more time with this problem.

Finite element analysis was developed for analysis of airplane structures at an American airforce base in the 50's and 60's.  It may not be applicable to a 0.7 x 0.7 x 0.017 mm chip.

It works to known solutions when solving for static problems of load and displacement even at such relatively small scale..

It does not seem to work when solving for natural frequencies of vibration against known equation solutions.

I think it is ill conditioned matrices when solving the eigen value and eigen vector solutions.

I am having the same problem on a commercial finite element software package as I am having with the Mathcad implementation.  So a little more time investigating is required.

Terry

P
payman1-VisitorAuthor
1-Visitor
July 9, 2019
Hi Terry, Thank you very much for your efforts. I'm a bit confused with lots of examples and articles, yet can't find my way. Most of the people I know in this field simply use COMSOL for such problems, but I wanted to have an insight into the method for further manipulations. All the best, Payman
-
-MFra-21-Topaz II
21-Topaz II
July 12, 2019

Hi,

You could use the finite difference method instead of pdesolve. In this way the problem is clearer and, then,  maybe you could even try using pdesolve. A suggestion how to do it is in the following photo(how to define omega?) :

PDE system FDM.jpg

... the most challenging part comes later.....

P
payman1-VisitorAuthor
1-Visitor
July 12, 2019

Hi Topaz,

I truly appreciate you. This is a great help to start. I'll consider this approach in the following days.

Thanks again,

Payman

T
21-Topaz II
terryhendicott21-Topaz IIAnswer
21-Topaz II
July 17, 2019

Hi,

 

This has taken a long time to verify against a commercial software suite Strand 7.

Enclosed is 3 zip files zipped into an over arching archive.

PlateVibration.zip are Mathcad worksheets.

References.zip are PDF references that show how finite elements work and some examples

Strand7_Files.zip is the commercial package files.  A free viewer is available at 

The files with extension lsl and nfl are log files and are text.  nfl standing for natural frequency log shows the results.

 

In summary the Mathcad worksheets perform natural frequency analysis of a rectangular plate with varying support conditions and can plot mode shapes for each natural frequency.

Capture.PNG

It has reasonable accuracy.

Capture-2.PNG

Enjoy

Terry

 

 

 

 

 

 

 

    P
    payman1-VisitorAuthor
    1-Visitor
    July 18, 2019

    This is a great job, thank you, Terry. 

    It will take some time to digest it.

     

    All the best,

    Payman

    T
    21-Topaz II
    terryhendicott21-Topaz II
    21-Topaz II
    July 22, 2019

    Hi,

    This file shows two new things

    a) How to clamp the top,bottom,left,and right edges of the plate as part of the finite elements.in terms of dx,dy at top of sheet.  By altering these boundary restraint codes all the support cases in RICE0523.pdf tables can be achieved.

    b) Comparison to RICE0523.pdf tables.

     

    The file uses a poissons ratio of 0.3 so comparison to RICE0523.pdf tables is possible.

    Capture.PNG

    b) Comparison to RICE0523.pdf tables is good.

    Capture-2.PNG

    Cheers

    Terry

      T
      21-Topaz II
      terryhendicott21-Topaz II
      21-Topaz II
      July 17, 2019

      Hi

      Next zip file

      Cheers 

      Terry

        T
        21-Topaz II
        terryhendicott21-Topaz II
        21-Topaz II
        July 17, 2019

        Hi,

        Final and biggest zip file

          Terms & ConditionsCookie settingsAccessibility statement