Skip to content

Nemus-Project/VKPlate

Repository files navigation

VKPlate

DOI

Implementation of a Föppl–von Kármán (VK) plate.

16-mode modal plate doing with initial excitation
16-mode modal plate doing with initial excitation

About

Simulation of a von Kármán plate written in modal formalism, using scalar auxiliary variable (SAV) method.

Getting started

Cloning this Repository

This project uses the magpie-matlab project as a submodule. To correctly clone this repository yo will need to use the command

git clone --recurse-submodules https://github.com/Nemus-Project/VKPlate

Alternatively, clone via the GitHub Desktop application or CLI Client.

NOTE: Downloading the zip will skip the download of the submodules.

Directory setup

There are some assets that need to be generated or downloaded before the project is useable.

The main.m script needs loads a .mat file containing a specific set of variables.

You can either

  1. generate it with genparams.m
  2. download the param.zip

The contents of param.zip should be put in a directory named param/ at the top level of the repository. The data is quite weighty (~2GB).

The genparams.m will create the param/ directory by default.

At the end the repository directory should look like:

VKPlate
├── README.md
├── genparams.m
├── magpie
├── main.m
├── img
│   └── ...
├── param
│   └── ...
├── private
│   ├── DxBuild.m
│   ├── DxxBuild.m
│   ├── DxyBuild.m
│   ├── DyBuild.m
│   ├── DyyBuild.m
│   ├── eigensign.m
│   ├── eigenMAC.m
│   ├── trapzIntcalc.m
│   └── vkOperator.m
└── test
    ├── Hcalc.m
    ├── ...
    └── SplitDuffing.m

Using main.m

Once the param/ is populates, run the main.m.

Parameter File Format

To avoid needless and time consuming recalculation, the central main.m script loads a .mat with all the require variables.

name description
rho Density
E Young modulus
nu Poisson's ratio
Lz Thickness of the plate
Lx x length of the plate
Ly y length of the plate
Nmodes Number of modes
Phi Eigenvectors of displacement
Om Angular frequency of the Phi modes
Psi Eigenvectors of stress
Om2 Angular frequency of the Psi modes
Nx Number of points in the x direction
Ny Number of points in the y direction
h Space step
X x space
Y y space
zetafourth Eigenvalues of the stress
BCsPhi Displacement boundary conditions
BCsPsi Stress boundary condition
Hv Coupling coefficients

Use the genparams.m script to create the .mat files or download from the latest release (See Directory setup)

Further Goals

Final checks on SAV, different boundary condition, plate under tension

References

see ./references.md

Troubleshooting

TBA

About

Föppl-von Kármán codes

Resources

License

Stars

Watchers

Forks

Contributors