mech_output_writer.h

Go to the documentation of this file.

/*
 * Copyright (c) 2014-2015:  G-CSC, Goethe University Frankfurt
 * Author: Raphael Prohl
 * 
 * This file is part of UG4.
 * 
 * UG4 is free software: you can redistribute it and/or modify it under the
 * terms of the GNU Lesser General Public License version 3 (as published by the
 * Free Software Foundation) with the following additional attribution
 * requirements (according to LGPL/GPL v3 §7):
 * 
 * (1) The following notice must be displayed in the Appropriate Legal Notices
 * of covered and combined works: "Based on UG4 (www.ug4.org/license)".
 * 
 * (2) The following notice must be displayed at a prominent place in the
 * terminal output of covered works: "Based on UG4 (www.ug4.org/license)".
 * 
 * (3) The following bibliography is recommended for citation and must be
 * preserved in all covered files:
 * "Reiter, S., Vogel, A., Heppner, I., Rupp, M., and Wittum, G. A massively
 *   parallel geometric multigrid solver on hierarchically distributed grids.
 *   Computing and visualization in science 16, 4 (2013), 151-164"
 * "Vogel, A., Reiter, S., Rupp, M., Nägel, A., and Wittum, G. UG4 -- a novel
 *   flexible software system for simulating pde based models on high performance
 *   computers. Computing and visualization in science 16, 4 (2013), 165-179"
 * 
 * This program is distributed in the hope that it will be useful,
 * but WITHOUT ANY WARRANTY; without even the implied warranty of
 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
 * GNU Lesser General Public License for more details.
 */
 
#ifndef MECH_OUTPUT_WRITER_H_
#define MECH_OUTPUT_WRITER_H_
 
using namespace std;
 
namespace ug{
namespace SmallStrainMechanics{
 
 
template<typename TDomain>
class MechOutputWriter
{
  private:
    typedef typename domain_traits<TDomain::dim>::grid_base_object TBaseElem;
 
  public:
    static const int dim = TDomain::dim;
 
  public:
    MechOutputWriter():
      m_quadOrder(0), m_bIP_values_written(false), m_stressEV(false),
      m_normalStress(false){};
 
    void stress_eigenvalues_at(const number CoordX,
        const number CoordY, const number CoordZ)
    {
      m_stressEV = true;
      m_evalPointEV[0] = CoordX;
      m_evalPointEV[1] = CoordY;
      if (dim == 3)
        m_evalPointEV[2] = CoordZ;
    }
 
    void normal_stresses_at(const number CoordX,
        const number CoordY, const number CoordZ)
    {
      //  TODO:check, if point is within the domain!
      //  if not, set m_normalStress to false!
      m_normalStress = true;
      m_evalPointNormStress[0] = CoordX;
      m_evalPointNormStress[1] = CoordY;
      if (dim == 3)
        m_evalPointNormStress[2] = CoordZ;
    }
 
    ~MechOutputWriter(){};
 
  public:
    void preprocess();
    void pre_timestep()
    { m_bIP_values_written = false;}
 
    template<typename TFEGeom>
    void post_timestep_elem(const number time, SmartPtr<TDomain> dom,
        TFEGeom& geo, TBaseElem* elem, const LocalVector& u);
 
    void post_timestep(const number time);
    void postprocess();
 
    void material_law(SmartPtr<IMaterialLaw<TDomain> > spMatLaw)
      {m_spMatLaw = spMatLaw;}
    void quad_order(const int quadOrder)
      {m_quadOrder = quadOrder;}
 
    void plastIP_elem(LocalVector& locPlastIP, TBaseElem* elem,
        const LocalVector& locU, SmartPtr<TDomain> dom);
 
    void eqPlastStrain_elem(LocalVector& locEqPlastStrain, TBaseElem* elem, const LocalVector& locU,
        SmartPtr<TDomain> dom);
 
    void normal_stress_strain_elem(LocalVector& locSigma, LocalVector& locEps, LocalVector& locStressFunc,
        LocalVector& locStrainFunc, TBaseElem* elem, const LocalVector& locU, SmartPtr<TDomain> dom);
 
    void invariants_kirchhoff_stress_elem(LocalVector& locInvarKirchhoffStress,
        TBaseElem* elem, const LocalVector& u, SmartPtr<TDomain> dom);
 
  private:
    template<typename TFEGeom>
    void next_ips_to_point(vector<size_t>& vNextIP, const MathVector<dim>& point,
        const TFEGeom& geo);
 
    template<typename TFEGeom>
    void stress_eigenvalues_near_point(const number time,
        TFEGeom& geo, const LocalVector& u);
 
    template<typename TFEGeom>
    void normal_stress_near_point(const number time,
        TFEGeom& geo, const LocalVector& u);
 
  private:
    SmartPtr<IMaterialLaw<TDomain> > m_spMatLaw;
 
    int m_quadOrder;
 
    bool m_bIP_values_written;
 
    bool m_stressEV;
    bool m_normalStress;
 
    MathVector<dim> m_evalPointEV;
    MathVector<dim> m_evalPointNormStress;
 
    FILE* m_fileStressEV;
 
};
 
//  some global output functions
 
template<typename TGridFunction>
void normal_stresses_strains(MechOutputWriter<typename TGridFunction::domain_type>& mechOut,
    TGridFunction& sigma, TGridFunction& epsilon, TGridFunction& stressFunc,
    TGridFunction& strainFunc, TGridFunction& u);
 
template <typename TGridFunction>
void plast_ip(TGridFunction& plastIP,
    MechOutputWriter<typename TGridFunction::domain_type>& mechOut,
    TGridFunction& u);
 
template <typename TGridFunction>
void equiv_plast_strain(TGridFunction& eqPlastStrain,
    MechOutputWriter<typename TGridFunction::domain_type>& mechOut,
    TGridFunction& u);
 
template <typename TGridFunction>
void invariants_kirchhoff_stress(TGridFunction& invarKirchhoffStress,
    MechOutputWriter<typename TGridFunction::domain_type>& mechOut,
    TGridFunction& u);
 
// end group small_strain_mechanics
 
} //end of namespace SmallStrainMechanics
} //end of namespace ug
 
#include "mech_output_writer_impl.h"
 
#endif /* MECH_OUTPUT_WRITER_H_ */