strain_computation.h

Go to the documentation of this file.

 
#ifndef STRAIN_COMPUTATION_H
#define STRAIN_COMPUTATION_H
 
#include <Eigen/Dense>
#include <vector>
 
namespace cppxdic {
 
struct DeformationResult {
    // Director vectors
    std::vector<Eigen::Vector3d> D1;  // Reference director 1
    std::vector<Eigen::Vector3d> D2;  // Reference director 2
    std::vector<Eigen::Vector3d> D3;  // Reference normal
    std::vector<Eigen::Vector3d> d1;  // Current director 1
    std::vector<Eigen::Vector3d> d2;  // Current director 2
    std::vector<Eigen::Vector3d> d3;  // Current normal
    
    // Reciprocal vectors
    std::vector<Eigen::Vector3d> Drec1;
    std::vector<Eigen::Vector3d> Drec2;
    std::vector<double> Dnorm;
    
    // Deformation gradient tensor (3x3 per face)
    std::vector<Eigen::Matrix3d> Fmat;
    
    // Cauchy-Green tensor
    std::vector<Eigen::Matrix3d> Cmat;
    
    // Principal stretches
    std::vector<double> Lamda1;
    std::vector<double> Lamda2;
    
    // Dilatation (area change)
    std::vector<double> J;
    
    // Lagrangian strain tensor
    std::vector<Eigen::Matrix3d> Emat;
    
    // Eulerian-Almansi strain tensor
    std::vector<Eigen::Matrix3d> emat;
    
    // Strain magnitudes
    std::vector<double> Emgn;  // Lagrangian
    std::vector<double> emgn;  // Eulerian
    
    // Planar principal strains (Lagrangian)
    std::vector<double> Epc1;  // Smallest
    std::vector<double> Epc2;  // Largest
    std::vector<Eigen::Vector3d> Epc1vec;     // Direction (reference)
    std::vector<Eigen::Vector3d> Epc2vec;     // Direction (reference)
    std::vector<Eigen::Vector3d> Epc1vecCur;  // Direction (current)
    std::vector<Eigen::Vector3d> Epc2vecCur;  // Direction (current)
    
    // Planar principal strains (Eulerian)
    std::vector<double> epc1;
    std::vector<double> epc2;
    std::vector<Eigen::Vector3d> epc1vec;
    std::vector<Eigen::Vector3d> epc2vec;
    
    // Maximum shear strain
    std::vector<double> EShearMax;  // Lagrangian
    std::vector<double> eShearMax;  // Eulerian
    std::vector<Eigen::Vector3d> EShearMaxVec1;
    std::vector<Eigen::Vector3d> EShearMaxVec2;
    std::vector<Eigen::Vector3d> EShearMaxVecCur1;
    std::vector<Eigen::Vector3d> EShearMaxVecCur2;
    std::vector<Eigen::Vector3d> eShearMaxVec1;
    std::vector<Eigen::Vector3d> eShearMaxVec2;
    
    // Equivalent strain (von Mises)
    std::vector<double> Eeq;  // Lagrangian
    std::vector<double> eeq;  // Eulerian
    
    // Area
    std::vector<double> Area;
};
 
struct FrameDeformationResult {
    size_t n_frames;
    size_t n_faces;
    
    // Each cell corresponds to one time frame
    std::vector<DeformationResult> frames;
};
 
FrameDeformationResult computeTriSurfaceDeformation(
    const std::vector<int>& faces,
    const std::vector<Eigen::Vector3d>& vertices_ref,
    const std::vector<std::vector<Eigen::Vector3d>>& vertices_def,
    bool cumulative = true
);
 
DeformationResult computeSingleFrameDeformation(
    const std::vector<int>& faces,
    const std::vector<Eigen::Vector3d>& vertices_ref,
    const std::vector<Eigen::Vector3d>& vertices_cur
);
 
} // namespace cppxdic
 
#endif // STRAIN_COMPUTATION_H