b2material_stress_2d_mixed_isotropic_linear.H

//------------------------------------------------------------------------
// b2material_stress_2d_mixed_isotropic_linear.H --
//
//     Implements a linear isotropic material model for 2D elements
//     using displacement/pressure formulation.
//
// written by Thomas Ludwig
//            Neda Ebrahimi Pour <neda.ebrahimipour@dlr.de>
//
// (c) 2009-2015,2016 SMR Engineering & Development SA
//               2502 Bienne, Switzerland
//
// (c) 2023 Deutsches Zentrum für Luft- und Raumfahrt (DLR) e.V.
//          Linder Höhe, 51147 Köln
//
// All Rights Reserved.  Proprietary source code.  The contents of
// this file may not be disclosed to third parties, copied or
// duplicated in any form, in whole or in part, without the prior
// written permission of SMR.
//------------------------------------------------------------------------
 
#ifndef B2MATERIAL_STRESS_2D_MIXED_ISOTROPIC_LINEAR_H_
#define B2MATERIAL_STRESS_2D_MIXED_ISOTROPIC_LINEAR_H_
 
#include "elements/properties/b2material_stress_2d_mixed.H"
#include "io/b2000_db/b2000_db_v3/b2fortran_element_property.H"
 
namespace b2000 {
 
class MaterialStress2DMixedIsotropicLinear : virtual public MaterialStress2DMixed {
public:
    void init(b2dbv3::MaterialProperty* material);
 
    void get_stress(
          Model* model, const bool linear, const EquilibriumSolution equilibrium_solution,
          const double time, const double delta_time, GradientContainer* gradient_container,
          SolverHints* solver_hints, const Element* element, const double el_coordinates[3],
          const int layer_id,
          const b2linalg::Vector<double, b2linalg::Vdense_constref> nodes_interpolation,
          const double bg_coordinates[3], const double covariant_base[2][2], const double volume,
          const double deformation_gradient[2][2], const double velocity[2],
          const double acceleration[2], const double interpolated_pressure, double stress[3],
          double& pressure, double CUU[6], double CUP[3], double& CPP, double inertia_force[2],
          double& density) override;
 
private:
    struct Layer {
        std::string not_supported;
        bool planestrain;
        double E;
        double nu;
        double mju;
        double lambda;
        double kappa;
        double cuu[6];
        double cup[3];
        double cpp;
        double density;
    };
 
    Layer d;
};
 
}  // namespace b2000
 
#endif  // B2_MATERIAL_STRESS_2D_MIXED_ISOTROPIC_LINEAR_H_