Code: TMA4500 Industrial Mathematics, Specialization Project

My code for building a reduced-order model for the linear elasticity equation on a square in 2D, solving for plane stress, under the conditons that the body force f, the prescribed traction force(s) h on the Neumann bondary and the prescribed displacment u on the Dirichlet boundary do not depend on my choice of parameters, the Young's module E and the Poisson ration ν.

Please see the Code used in report for code for the examples and the plotting results from the report, as the Patch Test, Example 1 traction forces and Example 2 ``Gravity in 2D´´.

Code not use in report has some unused code (example code) and old plotting/example code.

The solver Class

Please see LinearElasticity2DProblem for more documentation.

Default constants

Please see Default constants for documentation.

Useful helper functions

Please see Helpers for documentation.

Exceptions

Please see Exceptions for documentation.

Triangulation

Please see Get plate for documentation.

Other documentation

High-fidelity assembly

Please see Default constants for documentation.

Gauss quadrature

Please see Gauss quadrature for documentation.

Plotting

Please see Plotting for documentation.

Proper orthogonal decomposition with respect to the energy norm

Please see POD for documentation.

Reduced-order data class

Please see RB data for documentation.

Saving and loading matrices and vectors form files

Please see Save and load for documentation.

Solution Function class

Please see Solution function class for documentation.

Stress recovery

Please see Stress recovery for documentation.