common_assemble_functions_elasticity_3D.cpp File Reference

#include "common_assemble_functions_elasticity_3D.h"

Go to the source code of this file.

Functions
void	set_displaced_border_translation (libMesh::ImplicitSystem &elasticity_system, border_displacement_values &displ, int boundary_id)
	Set a displacement border. More...
void	set_clamped_border (libMesh::ImplicitSystem &elasticity_system, int boundary_id)
	Set a clamped border. More...
libMesh::ExplicitSystem &	add_stress (libMesh::EquationSystems &input_systems)
	Add a stress libMesh::ExplicitSystem to the input libMesh::EquationSystems. More...
libMesh::LinearImplicitSystem &	add_elasticity (libMesh::EquationSystems &input_systems, libMesh::Order order, libMesh::FEFamily family)
	Add a linear elasticity libMesh::LinearImplicitSystem to the input libMesh::EquationSystems& input_systems. More...

Function Documentation

libMesh::LinearImplicitSystem& add_elasticity	(	libMesh::EquationSystems &	input_systems,
		libMesh::Order	order,
		libMesh::FEFamily	family
	)

Add a linear elasticity libMesh::LinearImplicitSystem to the input libMesh::EquationSystems& input_systems.

Definition at line 68 of file common_assemble_functions_elasticity_3D.cpp.

{
    libMesh::ExplicitSystem& physical_variables =
            input_systems.add_system<libMesh::ExplicitSystem> ("PhysicalConstants");

    // Physical constants are set as constant, monomial
    physical_variables.add_variable("E", libMesh::CONSTANT, libMesh::MONOMIAL);
    physical_variables.add_variable("mu", libMesh::CONSTANT, libMesh::MONOMIAL);
    physical_variables.add_variable("Index", libMesh::CONSTANT, libMesh::MONOMIAL);

    libMesh::LinearImplicitSystem& elasticity_system =
            input_systems.add_system<libMesh::LinearImplicitSystem> ("Elasticity");

    elasticity_system.add_variable("u", order, family);
    elasticity_system.add_variable("v", order, family);
    elasticity_system.add_variable("w", order, family);

    return elasticity_system;
}

libMesh::ExplicitSystem& add_stress

(

libMesh::EquationSystems &

input_systems

)

Add a stress libMesh::ExplicitSystem to the input libMesh::EquationSystems.

Definition at line 49 of file common_assemble_functions_elasticity_3D.cpp.

{
    libMesh::ExplicitSystem& stress_system =
            input_systems.add_system<libMesh::ExplicitSystem> ("StressSystem");

    stress_system.add_variable("sigma_00", libMesh::CONSTANT, libMesh::MONOMIAL);
    stress_system.add_variable("sigma_01", libMesh::CONSTANT, libMesh::MONOMIAL);
    stress_system.add_variable("sigma_02", libMesh::CONSTANT, libMesh::MONOMIAL);
    stress_system.add_variable("sigma_10", libMesh::CONSTANT, libMesh::MONOMIAL);
    stress_system.add_variable("sigma_11", libMesh::CONSTANT, libMesh::MONOMIAL);
    stress_system.add_variable("sigma_12", libMesh::CONSTANT, libMesh::MONOMIAL);
    stress_system.add_variable("sigma_20", libMesh::CONSTANT, libMesh::MONOMIAL);
    stress_system.add_variable("sigma_21", libMesh::CONSTANT, libMesh::MONOMIAL);
    stress_system.add_variable("sigma_22", libMesh::CONSTANT, libMesh::MONOMIAL);
    stress_system.add_variable("vonMises", libMesh::CONSTANT, libMesh::MONOMIAL);

    return stress_system;
}

void set_clamped_border	(	libMesh::ImplicitSystem &	elasticity_system,
		int	boundary_id
	)

Set a clamped border.

Definition at line 27 of file common_assemble_functions_elasticity_3D.cpp.

{
    // Defining the boundaries with Dirichlet conditions ...
    std::set<libMesh::boundary_id_type> boundary_id_displacement;

    boundary_id_displacement.insert(boundary_id);

    std::vector<unsigned int> variables(3);
    variables[0] = elasticity_system.variable_number("u");
    variables[1] = elasticity_system.variable_number("v");
    variables[2] = elasticity_system.variable_number("w");

    libMesh::ZeroFunction<> zero_function;

    // ... and set them
    libMesh::DirichletBoundary dirichlet_bc(    boundary_id_displacement,
                                                variables,
                                                &zero_function);

    elasticity_system.get_dof_map().add_dirichlet_boundary(dirichlet_bc);
}

void set_displaced_border_translation	(	libMesh::ImplicitSystem &	elasticity_system,
		border_displacement_values &	displ,
		int	boundary_id
	)

Set a displacement border.

Definition at line 4 of file common_assemble_functions_elasticity_3D.cpp.

{
    // Defining the boundaries with Dirichlet conditions ...
    std::set<libMesh::boundary_id_type> boundary_id_displacement;

    boundary_id_displacement.insert(boundary_id);

    std::vector<unsigned int> variables(3);
    variables[0] = elasticity_system.variable_number("u");
    variables[1] = elasticity_system.variable_number("v");
    variables[2] = elasticity_system.variable_number("w");

    border_displacement_function move_border(   variables[0],variables[1],variables[2],
                                        displ.x_displ,displ.y_displ,displ.z_displ);

    // ... and set them
    libMesh::DirichletBoundary dirichlet_bc(    boundary_id_displacement,
                                                variables,
                                                &move_border);

    elasticity_system.get_dof_map().add_dirichlet_boundary(dirichlet_bc);
}