Implementation of a Föppl–von Kármán (VK) plate.
![]() |
|---|
| 16-mode modal plate doing with initial excitation |
Simulation of a von Kármán plate written in modal formalism, using scalar auxiliary variable (SAV) method.
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/VKPlateAlternatively, clone via the GitHub Desktop application or CLI Client.
NOTE: Downloading the zip will skip the download of the submodules.
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
- generate it with
genparams.m - 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
Once the param/ is populates, run the main.m.
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)
Final checks on SAV, different boundary condition, plate under tension
see ./references.md
TBA
