Diffusion akantu (features/doc_material_parameters) Browse

doc_material_p…

material_damage_iterative.hh
View Options

			/**
			* @file material_damage_iterative.hh
			*
			* @author Aurelia Isabel Cuba Ramos <aurelia.cubaramos@epfl.ch>
			*
			*
			* @brief Specialization of the class material damage to damage only one gauss
			* point at a time and propagate damage in a linear way. Max principal stress
			* criterion is used as a failure criterion.
			*
			* @section LICENSE
			*
			* Copyright (©) 2010-2012, 2014 EPFL (Ecole Polytechnique Fédérale de Lausanne)
			* Laboratory (LSMS - Laboratoire de Simulation en Mécanique des Solides)
			*
			*/

			/* -------------------------------------------------------------------------- */
			#include "aka_common.hh"
			#include "material_damage.hh"
			#include "material.hh"
			/* -------------------------------------------------------------------------- */

			#ifndef __AKANTU_MATERIAL_DAMAGE_ITERATIVE_HH__
			#define __AKANTU_MATERIAL_DAMAGE_ITERATIVE_HH__

			__BEGIN_AKANTU__

			/**
			* Material damage iterative
			*
			* parameters in the material files :
			* - Sc
			*/
			template<UInt spatial_dimension>
			class MaterialDamageIterative : public MaterialDamage<spatial_dimension> {
			/* ------------------------------------------------------------------------ */
			/* Constructors/Destructors */
			/* ------------------------------------------------------------------------ */
			public:

			MaterialDamageIterative(SolidMechanicsModel & model, const ID & id = "");

			virtual ~MaterialDamageIterative() {};

			/* ------------------------------------------------------------------------ */
			/* Methods */
			/* ------------------------------------------------------------------------ */
			public:

			/// virtual void updateInternalParameters();

			virtual void computeAllStresses(GhostType ghost_type = _not_ghost);

			/// update internal field damage
			virtual UInt updateDamage();

			UInt updateDamage(UInt quad_index, const Real eq_stress, const ElementType & el_type, const GhostType & ghost_type);

			/// update energies after damage has been updated
			virtual void updateEnergiesAfterDamage(ElementType el_type, GhostType ghost_typ);

			virtual void onBeginningSolveStep(const AnalysisMethod & method) { };

			virtual void onEndSolveStep(const AnalysisMethod & method) { };

			///compute the equivalent stress on each Gauss point (i.e. the max prinicpal stress) and normalize it by the tensile strength
			virtual void computeNormalizedEquivalentStress(const Array<Real> & grad_u,
			ElementType el_type, GhostType ghost_type = _not_ghost);

			/// find max normalized equivalent stress
			void findMaxNormalizedEquivalentStress(ElementType el_type, GhostType ghost_type = _not_ghost);

			protected:
			/// constitutive law for all element of a type
			virtual void computeStress(ElementType el_type, GhostType ghost_type = _not_ghost);

			inline void computeDamageAndStressOnQuad(Matrix<Real> & sigma,
			Real & dam);

			/* ------------------------------------------------------------------------ */
			/* DataAccessor inherited members */
			/* ------------------------------------------------------------------------ */

			inline UInt getNbDataForElements(const Array<Element> & elements,
			SynchronizationTag tag) const;

			inline void packElementData(CommunicationBuffer & buffer,
			const Array<Element> & elements,
			SynchronizationTag tag) const;

			inline void unpackElementData(CommunicationBuffer & buffer,
			const Array<Element> & elements,
			SynchronizationTag tag);

			/* ------------------------------------------------------------------------ */
			/* Accessors */
			/* ------------------------------------------------------------------------ */
			public:
			/// get max normalized equivalent stress
			AKANTU_GET_MACRO(NormMaxEquivalentStress, norm_max_equivalent_stress, Real);

			/// get a non-const reference to the max normalized equivalent stress
			AKANTU_GET_MACRO_NOT_CONST(NormMaxEquivalentStress, norm_max_equivalent_stress, Real &);

			/* ------------------------------------------------------------------------ */
			/* Class Members */
			/* ------------------------------------------------------------------------ */
			protected:

			/// resistance to damage
			RandomInternalField<Real> Sc;

			/// the reduction
			InternalField<UInt> reduction_step;

			/// internal field to store equivalent stress on each Gauss point
			InternalField<Real> equivalent_stress;

			/// the number of total reductions steps until complete failure
			UInt max_reductions;

			/// damage increment
			Real prescribed_dam;

			/// maximum equivalent stress
			Real norm_max_equivalent_stress;

			/// deviation from max stress at which Gauss point will still get damaged
			Real dam_tolerance;

			/// define damage threshold at which damage will be set to 1
			Real dam_threshold;

			/// maximum damage value
			Real max_damage;
			};

			/* -------------------------------------------------------------------------- */
			/* inline functions */
			/* -------------------------------------------------------------------------- */

			#include "material_damage_iterative_inline_impl.cc"

			__END_AKANTU__

			#endif /* __AKANTU_MATERIAL_DAMAGE_ITERATIVE_HH__ */

rAKA/extra_packages/extra-materials/src/material_damage/material_damage_iterative.hh13604ff6935ffeatures/doc_material_p…

material_damage_iterative.hhView Options

rAKA/extra_packages/extra-materials/src/material_damage/material_damage_iterative.hh
13604ff6935ffeatures/doc_material_p…

material_damage_iterative.hh
View Options