Elements.h

Go to the documentation of this file.

/*=============================================================================
TexGen: Geometric textile modeller.
Copyright (C) 2006 Martin Sherburn
 
This program is free software; you can redistribute it and/or
modify it under the terms of the GNU General Public License
as published by the Free Software Foundation; either version 2
of the License, or (at your option) any later version.
 
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 General Public License for more details.
 
You should have received a copy of the GNU General Public License
along with this program; if not, write to the Free Software
Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA  02110-1301, USA.
=============================================================================*/
 
#pragma once
 
#include "Matrix.h"
 
namespace TexGen
{
 
    class CLASS_DECLSPEC CElement
    {
    public:
        CElement(const CMatrix &P = CMatrix(), int iOrder = 3);
        virtual ~CElement(void) {}
 
        void SetNodeCoordinates(const CMatrix &P);
        void SetIntegrationOrder(int iOrder);
 
        void GetKeMatrix(CMatrix& KeMatrix);
 
        void SetFibreDirection(XYZ FibreDirection) { m_FibreDirection = FibreDirection; }
        XYZ GetFibreDirection() { return m_FibreDirection; }
 
    protected:
        struct INTEGRATION_POINT
        {
            CMatrix Position;
            double dWeight;
        };
 
        void GetBMatrix(CMatrix& BMatrix, CMatrix& Position);
        virtual void GetAMatrix(CMatrix& AMatrix) = 0;
        virtual void GetXMatrix(CMatrix& XMatrix, CMatrix& Position) = 0;
        virtual void GetDMatrix(CMatrix& DMatrix) = 0;
 
        virtual void GetIntegrationPoints(vector<INTEGRATION_POINT> &IntegrationPoints) = 0;
 
        CMatrix m_P;    // These are the nodal coordinates
        int m_iIntegrationOrder;
        XYZ m_FibreDirection;
    };
 
 
    class CLASS_DECLSPEC CElementTri : public CElement
    {
    public:
        CElementTri(const CMatrix &P, int iOrder);
        double GetArea();
 
    protected:
        void GetIntegrationPoints(vector<INTEGRATION_POINT> &IntegrationPoints);
    };
 
 
    class CLASS_DECLSPEC CElementTriBending : public CElementTri
    {
    public:
        CElementTriBending(const CMatrix &P = CMatrix(), int iOrder = 3);
 
        void SetLongitudinalBendingModulus(double E1) { m_E1 = E1; }
        double GetLongitudinalBendingModulus() { return m_E1; }
 
        void SetTransverseBendingModulus(double E2) { m_E2 = E2; }
        double GetTransverseBendingModulus() { return m_E2; }
 
    protected:
        void GetAMatrix(CMatrix& AMatrix);
        void GetXMatrix(CMatrix& XMatrix, CMatrix& Position);
        void GetDMatrix(CMatrix& DMatrix);
 
        CMatrix GetvMatrix(CMatrix& Position);
        CMatrix GetTheta_xMatrix(CMatrix& Position);
        CMatrix GetTheta_yMatrix(CMatrix& Position);
        CMatrix GetEpsilon_xMatrix(CMatrix& Position);
        CMatrix GetEpsilon_yMatrix(CMatrix& Position);
        CMatrix GetEpsilon_xyMatrix(CMatrix& Position);
 
        double m_E1;
        double m_E2;
    };
 
 
    class CLASS_DECLSPEC CElementTriTension : public CElementTri
    {
    public:
        CElementTriTension(const CMatrix &P = CMatrix(), int iOrder = 3);
 
        void SetTensileStress(double E) { m_E = E; }
        double GetTensileStress() { return m_E; }
 
    protected:
        void GetAMatrix(CMatrix& AMatrix);
        void GetXMatrix(CMatrix& XMatrix, CMatrix& Position);
        void GetDMatrix(CMatrix& DMatrix);
 
        CMatrix GetvMatrix(CMatrix& Position);
        CMatrix GetEpsilon_xMatrix(CMatrix& Position);
        CMatrix GetEpsilon_yMatrix(CMatrix& Position);
 
        double m_E;
    };
 
};  // namespace TexGen