b2material_stress_3d_mixed_linear.H

//------------------------------------------------------------------------
// b2material_stress_3d_mixed_linear.H --
//
//     Implements a linear orthotropic material model for 3D elements
//     using displacement/pressure formulation.
//
// written by Thomas Ludwig
//            Neda Ebrahimi Pour <neda.ebrahimipour@dlr.de>
//
// (c) 2009-2012,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_3D_MIXED_LINEAR_H_
#define B2MATERIAL_STRESS_3D_MIXED_LINEAR_H_
 
#include "elements/properties/b2material_stress_3d_mixed.H"
#include "io/b2000_db/b2000_db_v3/b2fortran_element_property.H"
 
namespace b2000 {
 
class MaterialStress3DMixedLinear : virtual public MaterialStress3DMixed {
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,
          b2linalg::Vector<double, b2linalg::Vdense_constref> nodes_interpolation,
          const double bg_coordinates[3], const double covariant_base[3][3], const double volume,
          const double deformation_gradient[3][3], const double velocity[3],
          const double acceleration[3], const double interpolated_pressure, double stress[6],
          double& pressure, double CUU[21], double CUP[6], double& CPP, double inertia_force[3],
          double& density) override;
 
private:
    struct Layer {
        std::string not_supported;
        bool isotropic;
        double density;
        double c[6][6];
        double cuu[21];
        double cup[6];
        double cpp;
        double p[6];
    };
 
    typedef std::vector<Layer> Layers;
 
    Layers layers;
};
 
}  // namespace b2000
 
#endif  // B2_MATERIAL_STRESS_3D_MIXED_LINEAR_H_