a_face.h

Go to the documentation of this file.

/*
Copyright 2010-2016 Pierre SMARS (smars@yuntech.edu.tw)
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, see <http://www.gnu.org/licenses/>.
*/
#ifndef _A_FACE_H_
#define _A_FACE_H_
 
#include <vector>
#include "a_wrench.h"
#include "a_point.h"
#include "a_segment.h"
#include "a_plane.h"
#include "a_contact.h"
#include "a_fcriteria.h"
 
class a_block;
 
class a_face
{
public:
        a_face();
        a_face(a_block * b, int ref_face);
        virtual ~a_face();
        static const std::string help();
        int nv() const {return r_vertices_.size();}
        void av(int ref) {r_vertices_.push_back(ref);}
        int rv(const int i) const;
        int rf() const {return ref_;}
        a_block * block() {return block_;}
        a_point v(const int i) const;
        a_point * vl(const int i);
        a_point * vl(const int i) const;
        a_point vc(const int i) const;
        a_point c() const;
        virtual a_point cl() const = 0;
        virtual void grow(const double& val) = 0;
        virtual a_point normale(const unsigned short i);
        virtual a_point point(const double& u, const double& v = 0.) const;
        virtual void slide(const double& u, const double& v);
        void deform(const unsigned short& i, const double& val);
        void exit(bool exit) {exit_face_ = exit;}
        bool exit() const {return exit_face_;}
        virtual double S() const = 0;
        double Sx() const;
        double Sy() const;
        double Sz() const;
        virtual a_point nx() const = 0;
        virtual a_point ny() const = 0;
        virtual a_point normal() const = 0;
        double dxM() const;
        double dxm() const;
        virtual double dyM() const;
        virtual double dym() const;
        //distance from max x to max y
        double dx() const;
        //distance from max y to max y
        double dy() const;
        virtual a_plane plane() const = 0; 
        virtual double N() const;
        virtual double N(const a_wrench&) const;
        virtual double Qx() const;
        virtual double Qy() const;
        virtual double Mt() const;
        virtual double Mx() const;
        virtual double Mx(const a_wrench&) const;
        virtual double My() const;
        virtual double My(const a_wrench&) const;
        double ax() const; 
        double ay() const; 
        virtual a_point x() const;
        virtual a_point x(const a_wrench&) const;
        double exx() const;
        double exx(const a_wrench&) const;
        double exy() const;
        double exy(const a_wrench&) const;
        double exn() const;
        void f(const a_wrench& f) {*f_=f;}
        void f(a_wrench * f);
        void f(const a_point& f);
        void f(const double fx, const double fy, const double fz);
        a_wrench f() const {return *f_;}
        a_wrench * fp() {return f_;}
        virtual void clear();
        void criteria(a_fcriteria * criteria) {criteria_=criteria;}
        bool ok() const;
        double penalty() const;
        void lface(a_face * f) {lface_ = f;}
        a_face * lface() const {return lface_;}
        a_block * lblock() const;
        bool in() const {return lface_ != 0;}
        bool out() const {return lface_ == 0;}
        a_point hinge() const;
        virtual int contacttype() const = 0;
        virtual a_point contactpoint() const = 0;
        virtual a_segment contactsegment() const = 0;
        virtual a_contact contactface() const = 0;
        virtual void read(std::istream &i);
        virtual void write(std::ostream &o) const;
        friend   std::istream&    operator>> (std::istream& i,       a_face& b);
        friend   std::ostream&    operator<< (std::ostream& o, const a_face& b);
 
        class bad_contact_type_error  {};
        class no_contact_error  {};
        class bad_indexing_error  {};
 
protected:
        void clear0();
        a_block * block_;
        int ref_;
        std::vector<int> r_vertices_;
        bool exit_face_;
        a_wrench * f_;
        bool f_original_;
        a_face * lface_;
        a_fcriteria * criteria_;
};
 
#endif