SolidBeamUnitTest.cpp File Reference

#include "Oomph/SolidProblem.h"
#include "Math/Helpers.h"

Classes
class	Beam
	Defines a SolidProblem of element type RefineableQDPVDElement<3,2>. Add functionality to write output. More...

Functions
int	main (int argc, char *argv[])

Function Documentation

main()

int main	(	int argc	,
		char *	argv[]
	)

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

{
    #ifdef OOMPH_HAS_MPI
        MPI_Helpers::init(argc,argv);
    #endif
    // Solve the problem
    Beam 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(Beam::Boundary::X_MIN);
    problem.setBodyForceAsGravity();
    problem.setNewtonSolverTolerance(3e-8);
    problem.prepareForSolve();
    #ifdef OOMPH_HAS_MPI
        // Distribute the problem
        problem.distribute();
        // Check halo schemes (optional)
        //problem.check_halo_schemes();
    #endif
    problem.linear_solver_pt()->disable_doc_time();
    //problem.disable_info_in_newton_solve();
    problem.solveSteady();
    problem.checkBeamDeflection();
    #ifdef OOMPH_HAS_MPI
        MPI_Helpers::finalize();
    #endif
 
    return 0;
}

References Beam::checkBeamDeflection(), 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().