Loads.h

Go to the documentation of this file.

// SPDX-License-Identifier: MIT
// Copyright (c) 2021-2025 MPM-Geomechanics Development Team
 
#ifndef LOADS_H_
#define LOADS_H_
 
#include "Body/Body.h"
 
#include <vector>
using std::vector;
 
#include "Eigen/Core"
using Eigen::Vector3d;
 
namespace Loads {
 
    struct LoadDistributedBox
    {
        Vector3d pointP1; 
        Vector3d pointP2; 
        Vector3d load; 
    };
 
    struct NodalPointLoadData
    {
        std::vector<Vector3d> points; 
        std::vector <Vector3d> loads; 
        std::vector<int> nodal_ids; 
    };
 
    struct PressureBox
    {
        Vector3d pointP1; 
        Vector3d pointP2; 
        double pressure; 
    };
 
    struct PressureBoundaryForceBox
    {
        Vector3d pointP1; 
        Vector3d pointP2; 
        Vector3d pressureForce; 
    };
 
    struct PressureMaterial
    {
        int materialId; 
        double pressure; 
    };
 
    struct PrescribedPorePressure{
 
        int bodyIndex; 
        int particleIndex; 
        double pressure; 
 
        PrescribedPorePressure(int body_index, int particle_index, double pressure_value) {
            
            bodyIndex=body_index;
            particleIndex=particle_index;
            pressure=pressure_value;
        }
    };
 
    struct SeismicData
    {
        std::vector<double> time;
        std::vector<Eigen::Vector3d> acceleration;
        std::vector<Eigen::Vector3d> velocity;
    };
 
    void setSeismicData();
 
    SeismicData& getSeismicData();
 
    void setGravity(vector<Body*>& bodies);
    
    void setLoadDistributedBox(vector<Body*>& bodies, vector<Loads::LoadDistributedBox> loads);
 
    void setPrescribedPorePressureBox(vector<Body*>& bodies, vector<Loads::PressureBox> pressures);
 
    void setInitialPorePressureBox(vector<Body*>& bodies, vector<Loads::PressureBox> pressures);
 
    void setInitialPorePressureMaterial(vector<Body*>& bodies, vector<Loads::PressureMaterial> pressures);
 
    void updatePrescribedPorePressure(vector<Body*>* bodies);
    
    void setPrescribedPorePressureBoundaryForceBox(vector<Body*>& bodies, vector<Loads::PressureBoundaryForceBox> loads);
 
    void setInitialVelocity(vector<Body*>& bodies);
 
    void configureNodalPointLoads(Mesh* mesh);
 
    NodalPointLoadData& getNodalPointList();
};
 
#endif /* LOADS_H_ */