SolidBeamDemo.cpp File Reference

#include "Oomph/SolidProblem.h"

Functions
int	main ()
	[SolidBeamDemo] More...

Function Documentation

main()

int main

(

)

[SolidBeamDemo]

Simulates a cubic rod (20mm x 2mm x 2mm) of density 2500 kg/m^3, elastic modulus 100 MPa, fixed on left side, bent by gravity

{
    // Defines a SolidProblem of element type RefineableQDPVDElement<3,2>.
    SolidProblem<RefineableQDPVDElement<3, 2>> problem;
    // Sets up the problem
    problem.setName("SolidBeamUnitTest");
    problem.setElasticModulus(1e8);
    problem.setDensity(2500);
    problem.setSolidCubicMesh(20, 2, 2, 0, 0.2, 0, 0.02, 0, 0.02);
    problem.pinBoundary(problem.Boundary::X_MIN);
    problem.setBodyForceAsGravity();
    problem.setNewtonSolverTolerance(3e-8);
    problem.prepareForSolve();
    // Solve the problem
    problem.solveSteady();
    return 0;
}

References SolidProblem< ELEMENT_TYPE >::pinBoundary(), SolidProblem< ELEMENT_TYPE >::prepareForSolve(), SolidProblem< ELEMENT_TYPE >::setBodyForceAsGravity(), SolidProblem< ELEMENT_TYPE >::setDensity(), SolidProblem< ELEMENT_TYPE >::setElasticModulus(), SolidProblem< ELEMENT_TYPE >::setName(), SolidProblem< ELEMENT_TYPE >::setNewtonSolverTolerance(), SolidProblem< ELEMENT_TYPE >::setSolidCubicMesh(), and SolidProblem< ELEMENT_TYPE >::solveSteady().