forked from sandialabs/Albany
-
Notifications
You must be signed in to change notification settings - Fork 0
Commit
This commit does not belong to any branch on this repository, and may belong to a fork outside of the repository.
LCM: rework kinematics and constittive model parameters
Essentially create Evaluator interfaces to reduce the number of evaluators that we are required to register. Also change NeoHookean to Neohookean everywhere, because of reasons.
- Loading branch information
Showing
33 changed files
with
687 additions
and
64 deletions.
There are no files selected for viewing
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| Original file line number | Diff line number | Diff line change |
|---|---|---|
|
|
@@ -19,3 +19,4 @@ phalanx_graph | |
| tags | ||
| *12 | ||
| *13 | ||
| *.save | ||
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| Original file line number | Diff line number | Diff line change |
|---|---|---|
| @@ -0,0 +1,13 @@ | ||
| //*****************************************************************// | ||
| // Albany 2.0: Copyright 2012 Sandia Corporation // | ||
| // This Software is released under the BSD license detailed // | ||
| // in the file "license.txt" in the top-level Albany directory // | ||
| //*****************************************************************// | ||
|
|
||
| #include "PHAL_AlbanyTraits.hpp" | ||
|
|
||
| #include "ConstitutiveModelParameters.hpp" | ||
| #include "ConstitutiveModelParameters_Def.hpp" | ||
|
|
||
| PHAL_INSTANTIATE_TEMPLATE_CLASS(LCM::ConstitutiveModelParameters) | ||
|
|
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| Original file line number | Diff line number | Diff line change |
|---|---|---|
| @@ -0,0 +1,129 @@ | ||
| //*****************************************************************// | ||
| // Albany 2.0: Copyright 2012 Sandia Corporation // | ||
| // This Software is released under the BSD license detailed // | ||
| // in the file "license.txt" in the top-level Albany directory // | ||
| //*****************************************************************// | ||
|
|
||
| #if !defined(LCM_ConstitutiveModelParameters_hpp) | ||
| #define LCM_ConstitutiveModelParameters_hpp | ||
|
|
||
| #include "Phalanx_ConfigDefs.hpp" | ||
| #include "Phalanx_Evaluator_WithBaseImpl.hpp" | ||
| #include "Phalanx_Evaluator_Derived.hpp" | ||
| #include "Phalanx_MDField.hpp" | ||
|
|
||
| #include "Teuchos_ParameterList.hpp" | ||
| #include "Epetra_Vector.h" | ||
| #include "Sacado_ParameterAccessor.hpp" | ||
| #include "Stokhos_KL_ExponentialRandomField.hpp" | ||
| #include "Teuchos_Array.hpp" | ||
| #include "Albany_Layouts.hpp" | ||
|
|
||
| namespace LCM { | ||
| /// | ||
| /// \brief Evaluates a selecltion of Constitutive Model Parameters | ||
| /// Either as a constant or a truncated KL expansion. | ||
| /// | ||
| template<typename EvalT, typename Traits> | ||
| class ConstitutiveModelParameters : | ||
| public PHX::EvaluatorWithBaseImpl<Traits>, | ||
| public PHX::EvaluatorDerived<EvalT, Traits>, | ||
| public Sacado::ParameterAccessor<EvalT, SPL_Traits> { | ||
|
|
||
| public: | ||
| typedef typename EvalT::ScalarT ScalarT; | ||
| typedef typename EvalT::MeshScalarT MeshScalarT; | ||
|
|
||
| /// | ||
| /// Constructor | ||
| /// | ||
| ConstitutiveModelParameters(Teuchos::ParameterList& p, | ||
| const Teuchos::RCP<Albany::Layouts>& dl); | ||
|
|
||
| /// | ||
| /// Phalanx method to allocate space | ||
| /// | ||
| void postRegistrationSetup(typename Traits::SetupData d, | ||
| PHX::FieldManager<Traits>& vm); | ||
|
|
||
| /// | ||
| /// Implementation of physics | ||
| /// | ||
| void evaluateFields(typename Traits::EvalData d); | ||
|
|
||
| /// | ||
| /// Sacado method to access parameter values | ||
| /// | ||
| ScalarT& getValue(const std::string &n); | ||
|
|
||
| /// | ||
| /// Helper method to parse a parameter | ||
| /// | ||
| void parseParameters(const std::string &n, | ||
| Teuchos::ParameterList &pl, | ||
| Teuchos::RCP<ParamLib> paramLib); | ||
|
|
||
| private: | ||
|
|
||
| /// | ||
| /// Number of integration points | ||
| /// | ||
| std::size_t num_pts_; | ||
|
|
||
| /// | ||
| /// Number of spatial dimensions | ||
| /// | ||
| std::size_t num_dims_; | ||
|
|
||
| /// | ||
| /// spatial locations of integration points | ||
| /// | ||
| PHX::MDField<MeshScalarT,Cell,QuadPoint,Dim> coord_vec_; | ||
|
|
||
| /// | ||
| /// Constitutive Model Parameters | ||
| /// | ||
| /// Elastic Moduli | ||
| PHX::MDField<ScalarT,Cell,QuadPoint> elastic_mod_; | ||
| PHX::MDField<ScalarT,Cell,QuadPoint> poissons_ratio_; | ||
| PHX::MDField<ScalarT,Cell,QuadPoint> bulk_mod_; | ||
| PHX::MDField<ScalarT,Cell,QuadPoint> shear_mod_; | ||
| /// Plasticity Parameters | ||
| PHX::MDField<ScalarT,Cell,QuadPoint> yield_strength_; | ||
| PHX::MDField<ScalarT,Cell,QuadPoint> hardening_mod_; | ||
|
|
||
| /// | ||
| /// map of strings to specify parameter names to MDFields | ||
| /// | ||
| std::map<std::string, PHX::MDField<ScalarT,Cell,QuadPoint> > field_map_; | ||
|
|
||
| /// | ||
| /// map of flags to specify if a parameter is constant | ||
| /// | ||
| std::map<std::string, bool> is_constant_map_; | ||
|
|
||
| /// | ||
| /// map of strings to ScalarTs to specify constant values | ||
| /// | ||
| std::map<std::string, ScalarT> constant_value_map_; | ||
|
|
||
| /// | ||
| /// Optional dependence on Temperature | ||
| /// | ||
| bool have_temperature_; | ||
| PHX::MDField<ScalarT,Cell,QuadPoint> temperature_; | ||
| std::map<std::string, RealType> dparam_dtemp_map_; | ||
| std::map<std::string, RealType> ref_temp_map_; | ||
|
|
||
| //! map of strings to exponential random fields | ||
| std::map<std::string, Teuchos::RCP< Stokhos::KL::ExponentialRandomField<MeshScalarT> > > exp_rf_kl_map_; | ||
|
|
||
| //! map of strings to Arrays of values of the random variables | ||
| std::map<std::string, Teuchos::Array<ScalarT> > rv_map_; | ||
|
|
||
| //! storing the DataLayouts | ||
| const Teuchos::RCP<Albany::Layouts>& dl_; | ||
| }; | ||
| } | ||
|
|
||
| #endif |
Oops, something went wrong.