meshUtil.h

Go to the documentation of this file.

/*
 * -------------------------------------------
 * Copyright (c) 2021 - 2024 Prashant K. Jha
 * -------------------------------------------
 * PeriDEM https://github.com/prashjha/PeriDEM
 *
 * Distributed under the Boost Software License, Version 1.0. (See accompanying
 * file LICENSE)
 */
 
#pragma once
 
#include "util/point.h"
#include "util/matrix.h"
#include "inp/decks/materialDeck.h"
#include <string>
#include <vector>
 
namespace fe {
 
// forward declare Mesh
class Mesh;
 
void createUniformMesh(fe::Mesh *mesh_p, size_t dim, std::pair<std::vector<double>, std::vector<double>> box, std::vector<size_t> nGrid);
 
void getCurrentQuadPoints(const fe::Mesh *mesh_p,
                          const std::vector<util::Point> &xRef,
                          const std::vector<util::Point> &u,
                          std::vector<util::Point> &xQuadCur,
                          size_t iNodeStart = 0,
                          size_t iQuadStart = 0,
                          size_t quadOrder = 1);
 
void getStrainStress(const fe::Mesh *mesh_p,
                     const std::vector<util::Point> & xRef,
                     const std::vector<util::Point> &u,
                     bool isPlaneStrain,
                     std::vector<util::SymMatrix3> &strain,
                     std::vector<util::SymMatrix3> &stress,
                     size_t iNodeStart = 0,
                     size_t iStrainStart = 0,
                     double nu = 0.,
                     double lambda = 0.,
                     double mu = 0.,
                     bool computeStress = false,
                     size_t quadOrder = 1);
 
void getMaxShearStressAndLoc(const fe::Mesh *mesh_p,
                             const std::vector<util::Point> & xRef,
                             const std::vector<util::Point> &u,
                             const std::vector<util::SymMatrix3> &stress,
                             double &maxShearStress,
                             util::Point &maxShearStressLocRef,
                             util::Point &maxShearStressLocCur,
                             size_t iNodeStart = 0,
                             size_t iStrainStart = 0,
                             size_t quadOrder = 1);
 
} // namespace fe