linear_strain_energy_response_utils.h

Go to the documentation of this file.

//    |  /           |
//    ' /   __| _` | __|  _ \   __|
//    . \  |   (   | |   (   |\__ `
//   _|\_\_|  \__,_|\__|\___/ ____/
//                   Multi-Physics
//
//  License:         BSD License
//                   license: OptimizationApplication/license.txt
//
//  Main author:     Reza Najian Asl,
//                   Suneth Warnakulasuriya
//
 
#pragma once
 
// System includes
#include <variant>
#include <vector>
 
// Project includes
#include "includes/define.h"
#include "includes/model_part.h"
#include "expression/container_expression.h"
 
// Application includes
 
namespace Kratos
{
 
 
class KRATOS_API(OPTIMIZATION_APPLICATION) LinearStrainEnergyResponseUtils
{
public:
 
    using IndexType = std::size_t;
 
    using GeometryType = ModelPart::ElementType::GeometryType;
 
    using PhysicalFieldVariableTypes = std::variant<const Variable<double>*, const Variable<array_1d<double, 3>>*>;
 
    using ContainerExpressionType = std::variant<ContainerExpression<ModelPart::NodesContainerType>::Pointer, ContainerExpression<ModelPart::ConditionsContainerType>::Pointer, ContainerExpression<ModelPart::ElementsContainerType>::Pointer>;
 
 
    static double CalculateValue(ModelPart& rEvaluatedModelPart);
 
    static void CalculateGradient(
        const PhysicalFieldVariableTypes& rPhysicalVariable,
        ModelPart& rGradientRequiredModelPart,
        ModelPart& rGradientComputedModelPart,
        std::vector<ContainerExpressionType>& rListOfContainerExpressions,
        const double PerturbationSize);
 
private:
 
    template<class TEntityType>
    static double CalculateEntityStrainEnergy(
        TEntityType& rEntity,
        Matrix& rLHS,
        Vector& rRHS,
        Vector& rX,
        const ProcessInfo& rProcessInfo);
 
    template<class TEntityType>
    static void CalculateStrainEnergyEntitySemiAnalyticShapeGradient(
        TEntityType& rEntity,
        Vector& rX,
        Vector& rRefRHS,
        Vector& rPerturbedRHS,
        typename TEntityType::Pointer& pThreadLocalEntity,
        ModelPart& rModelPart,
        const double Delta,
        const Variable<array_1d<double, 3>>& rOutputGradientVariable);
 
    static void CalculateStrainEnergySemiAnalyticShapeGradient(
        ModelPart& rModelPart,
        const double Delta,
        const Variable<array_1d<double, 3>>& rOutputGradientVariable);
 
    static void CalculateStrainEnergyLinearlyDependentPropertyGradient(
        ModelPart& rModelPart,
        const Variable<double>& rPrimalVariable,
        const Variable<double>& rOutputGradientVariable);
 
    static void CalculateStrainEnergySemiAnalyticPropertyGradient(
        ModelPart& rModelPart,
        const double Delta,
        const Variable<double>& rPrimalVariable,
        const Variable<double>& rOutputGradientVariable);
 
};
 
}