b2material_stress_2d_mesh_deform.H

//------------------------------------------------------------------------
// b2material_stress_2d_mesh_deform.H --
//
//     Material used for 2D Mesh Deformation (Aero-Elasticity).
//
// 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_2D_MESH_DEFORM_H_
#define B2MATERIAL_STRESS_2D_MESH_DEFORM_H_
 
#include "elements/properties/b2solid_mechanics.H"
#include "utils/b2tensor_calculus.H"
 
namespace b2000 {
 
class MaterialStress2DMeshDeform : public SolidMechanics2D<double> {
public:
    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[2],
          const b2linalg::Vector<double, b2linalg::Vdense_constref> nodes_interpolation,
          const int layer_id, const double bg_coordinates[3], const double covariant_base[2][2],
          const bool use_covariant_base, const double volume,
          const double displacement_gradient[2][2], const double strain[3],
          const double strain_dot[3], const double velocity[2], const double acceleration[2],
          double stress[3], double d_stress_d_strain[6], double d_stress_d_strain_dot[6],
          double d_stress_d_time[3], double inertia_force[2], double& density);
 
    inline void s_eq_c_mul_e(double s[3], const double c[6], const double e[3]) {
        const int ind[3][3] = {{0, 1, 2}, {1, 3, 4}, {2, 4, 5}};
        for (int i = 0; i < 3; ++i) {
            double value = 0.0;
 
            const int* iind = ind[i];
            for (int j = 0; j < 3; ++j) { value += c[iind[j]] * e[j]; }
            s[i] = value;
        }
    }
 
    // Transform the 2nd Piola-Kirchhoff stress to the Cauchy
    // stress.  The stresses are in the Voight notation [sxx, syy,
    // sxy].
    static inline void piola_2_stress_to_cauchy_stress(
          const double F[2][2],  // Deformation gradient, column-major.
          const double S[3],     // 2nd Piola-Kirchhoff stress (Voight not.).
          double T[3]            // Cauchy stress (Voight notation).
    ) {
        const double f00 = F[0][0];
        const double f01 = F[0][1];
        const double f10 = F[1][0];
        const double f11 = F[1][1];
 
        // Get the determinant of the deformation gradient tensor..
        const double det = determinant_2_2(F);
        const double det_inv = 1.0 / det;
 
        T[0] = f00 * f00 * S[0] + f10 * f10 * S[1] + 2 * f00 * f10 * S[2];
 
        T[1] = f01 * f01 * S[0] + f11 * f11 * S[1] + 2 * f01 * f11 * S[2];
 
        T[2] = f00 * f01 * S[0] + f10 * f11 * S[1] + (f00 * f11 + f10 * f01) * S[2];
 
        for (int i = 0; i < 3; ++i) { T[i] *= det_inv; }
    }
 
public:
    double mat_e;
    double mat_p;
    double mat_density;
};
 
}  // namespace b2000
 
#endif  // B2_MATERIAL_STRESS_2D_MESH_DEFORM_H_