shells/fragmentation/CElement.h

Go to the documentation of this file.

// -*- C++ -*- 
//
// ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
//
//                                   Fehmi Cirak
//                        California Institute of Technology
//                           (C) 2003 All Rights Reserved
//
// <LicenseText>
//
// ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
//
#ifndef CELEMENT_H
#define CELEMENT_H

#include "CMaterial.h"

#include <iosfwd>
#include <functional>
#include <vector>


namespace shells {
    struct SElementS;
    struct SVertexS;
}


namespace fragment {
    class CElement;
    class ComputeInternalForceFunctor;

    // LE - left element, RI - right element
    enum ElementType {LE=0, RI=1};

    // parameters for thickness integration
    const int tgauss = 3;
    const double eta[3]  = {-0.5, 0, 0.5};
    const double weta[3] = {0.166666666666667, 0.666666666666667, 0.16666666666667};

    // parameters for edge integration
    const int    egauss = 3;
    const double wcsi[3] = {0.277777777777778, 0.444444444444444, 0.27777777777778};
}



class fragment::CElement {
public:
    CElement(shells::SElementS * const left, 
             shells::SVertexS * const leftVtxA, shells::SVertexS * const leftVtxB, 
             shells::SElementS * const right, 
             shells::SVertexS * const rightVtxA, shells::SVertexS * const rightVtxB,
             bool halfWeight): 
        _left(left), _leftVtxA(leftVtxA), _leftVtxB(leftVtxB),
        _right(right), _rightVtxA(rightVtxA), _rightVtxB(rightVtxB), 
        _halfWeight(halfWeight) 
        {
            _mat = NULL;
        } 
    
    void setMaterial(fragment::CMaterial * const mat);
    void computeResidual(const double thick);

    void ptrToHistory(double **data, int *size);
    void resetHistory(double *data);

    // for debugging
    void printCoordinates (std::ostream& os) const;
    void printActiveCoordinates (std::ostream& os) const;

private:
    CElement(const CElement &);
    const CElement & operator=(const CElement &);

private:
    shells::SElementS     * const _left;    
    shells::SVertexS      * const _leftVtxA;
    shells::SVertexS      * const _leftVtxB;
    
    shells::SElementS     * const _right;
    shells::SVertexS      * const _rightVtxA;
    shells::SVertexS      * const _rightVtxB;

    CMaterial             *_mat;    

    std::vector<double>    _internal;    
    
    bool                   _halfWeight;
};



class fragment::ComputeInternalForceFunctor : 
    public std::unary_function<fragment::CElement *, void > {
    // introduced in order to make the library usage similar to the shell element
public:
    ComputeInternalForceFunctor(const double thick):_thick(thick){}

    void operator()(CElement * const elem)  const
        {
            elem->computeResidual(_thick);
        }

private:
    const double     _thick;
};

#endif

---