SolidMechanicsFields.hpp

Go to the documentation of this file.

/*
 * ------------------------------------------------------------------------------------------------------------
 * SPDX-License-Identifier: LGPL-2.1-only
 *
 * Copyright (c) 2016-2024 Lawrence Livermore National Security LLC
 * Copyright (c) 2018-2024 TotalEnergies
 * Copyright (c) 2018-2024 The Board of Trustees of the Leland Stanford Junior University
 * Copyright (c) 2023-2024 Chevron
 * Copyright (c) 2019-     GEOS/GEOSX Contributors
 * All rights reserved
 *
 * See top level LICENSE, COPYRIGHT, CONTRIBUTORS, NOTICE, and ACKNOWLEDGEMENTS files for details.
 * ------------------------------------------------------------------------------------------------------------
 */
 
#ifndef GEOS_PHYSICSSOLVERS_SOLIDMECHANICS_SOLIDMECHANICSFIELDS_HPP_
#define GEOS_PHYSICSSOLVERS_SOLIDMECHANICS_SOLIDMECHANICSFIELDS_HPP_
 
#include "common/DataLayouts.hpp"
#include "mesh/MeshFields.hpp"
 
namespace geos
{
namespace fields
{
 
namespace solidMechanics
{
 
using array2dLayoutTotalDisplacement = array2d< real64, nodes::TOTAL_DISPLACEMENT_PERM >;
using arrayView2dLayoutTotalDisplacement = arrayView2d< real64, nodes::TOTAL_DISPLACEMENT_USD >;
using arrayViewConst2dLayoutTotalDisplacement = arrayView2d< real64 const, nodes::TOTAL_DISPLACEMENT_USD >;
 
using array2dLayoutIncrDisplacement = array2d< real64, nodes::INCR_DISPLACEMENT_PERM >;
using arrayView2dLayoutIncrDisplacement = arrayView2d< real64, nodes::INCR_DISPLACEMENT_USD >;
using arrayViewConst2dLayoutIncrDisplacement = arrayView2d< real64 const, nodes::INCR_DISPLACEMENT_USD >;
 
using array2dLayoutStrain = array2d< real64, cells::RANK2_TENSOR_PERM >;
using arrayView2dLayoutStrain = arrayView2d< real64, cells::RANK2_TENSOR_USD >;
using arrayViewConst2dLayoutStrain = arrayView2d< real64 const, cells::RANK2_TENSOR_USD >;
 
using array2dLayoutAvgStress = array2d< real64, cells::RANK2_TENSOR_PERM >;
using arrayView2dLayoutAvgStress = arrayView2d< real64, cells::RANK2_TENSOR_USD >;
using arrayViewConst2dLayoutAvgStress = arrayView2d< real64 const, cells::RANK2_TENSOR_USD >;
 
using array2dLayoutVelocity = array2d< real64, nodes::VELOCITY_PERM >;
using arrayView2dLayoutVelocity = arrayView2d< real64, nodes::VELOCITY_USD >;
using arrayViewConst2dLayoutVelocity = arrayView2d< real64 const, nodes::VELOCITY_USD >;
 
using array2dLayoutAcceleration = array2d< real64, nodes::ACCELERATION_PERM >;
using arrayView2dLayoutAcceleration = arrayView2d< real64, nodes::ACCELERATION_USD >;
using arrayViewConst2dLayoutAcceleration = arrayView2d< real64 const, nodes::ACCELERATION_USD >;
 
 
DECLARE_FIELD( totalDisplacement,
               "totalDisplacement",
               array2dLayoutTotalDisplacement,
               0,
               LEVEL_0,
               WRITE_AND_READ,
               "Total displacements at the nodes" );
 
DECLARE_FIELD( incrementalDisplacement,
               "incrementalDisplacement",
               array2dLayoutIncrDisplacement,
               0,
               LEVEL_3,
               WRITE_AND_READ,
               "Incremental displacements for the current time step on the nodes" );
 
DECLARE_FIELD( averageStrain,
               "averageStrain",
               array2dLayoutStrain,
               0,
               LEVEL_0,
               WRITE_AND_READ,
               "Quadrature averaged strain components in cell" );
 
DECLARE_FIELD( averageStress,
               "averageStress",
               array2dLayoutAvgStress,
               0,
               LEVEL_0,
               WRITE_AND_READ,
               "Quadrature averaged stress components in cell" );
 
DECLARE_FIELD( averagePlasticStrain,
               "averagePlasticStrain",
               array2dLayoutStrain,
               0,
               LEVEL_0,
               WRITE_AND_READ,
               "Average plastic strain in cell" );
 
DECLARE_FIELD( velocity,
               "velocity",
               array2dLayoutVelocity,
               0,
               LEVEL_0,
               WRITE_AND_READ,
               "Current velocity on the nodes" );
 
DECLARE_FIELD( acceleration,
               "acceleration",
               array2dLayoutAcceleration,
               0,
               LEVEL_1,
               WRITE_AND_READ,
               "Current acceleration on the nodes. This array also is used "
               "to hold the summation of nodal forces resulting from the governing equations" );
 
DECLARE_FIELD( externalForce,
               "externalForce",
               array2d< real64 >,
               0,
               LEVEL_0,
               WRITE_AND_READ,
               "External forces on the nodes. This includes any boundary"
               " conditions as well as coupling forces such as hydraulic forces" );
 
DECLARE_FIELD( mass,
               "mass",
               array1d< real64 >,
               0,
               LEVEL_0,
               WRITE_AND_READ,
               "Mass on the nodes" );
 
DECLARE_FIELD( velocityTilde,
               "velocityTilde",
               array2d< real64 >,
               0,
               NOPLOT,
               WRITE_AND_READ,
               "Velocity predictors on the nodes" );
 
DECLARE_FIELD( uhatTilde,
               "uhatTilde",
               array2d< real64 >,
               0,
               NOPLOT,
               WRITE_AND_READ,
               "Incremental displacement predictors on the nodes" );
 
DECLARE_FIELD( contactForce,
               "contactForce",
               array2d< real64 >,
               0,
               NOPLOT,
               WRITE_AND_READ,
               "Contact force" );
 
}
 
}
 
}
 
#endif // GEOS_PHYSICSSOLVERS_SOLIDMECHANICS_SOLIDMECHANICSFIELDS_HPP_