a_fcriteria_in.cxx

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/>.
*/
#include "a_fcriteria_in.h"
#include "a_face_2d4.h"
 
#include "config.h"
 
//---------------------------------------------------------------------------
const std::string a_fcriteria_in::help()
{
        std::ostringstream o;
        o << a_fcriteria::help();
        o << "***************" << std::endl;
        o << "a_fcriteria_in:" << std::endl;
        o << "***************" << std::endl;
        o << "This class derives from a _fcriteria," << std::endl;
        o << "it implements a specific resistance criterion:" << std::endl;
        o << "The excentricity of the resultant of the force on a face is required to remain in between specific bounds" << std::endl;
        o << "max 'value': specify the bound (default: 1 =  pass in the section)" << std::endl;
        return o.str();
}
//---------------------------------------------------------------------------
bool a_fcriteria_in::ok(const a_face * f0) const
{
        const a_face_2d4 * f = cf_2d4(f0);
        double exrx = fabs(f->exrx());
        double exry = fabs(f->exry());
        // std::cerr << "ex: " << exrx << " " << exry << std::endl;
        if ((exrx<max_+verysmall_) and (exry <max_+verysmall_))
                return true;
        return false;
}
//---------------------------------------------------------------------------
double a_fcriteria_in::penalty(const a_face * f0) const
{
        const a_face_2d4 * f = cf_2d4(f0);
        double exrx = fabs(f->exrx());
        double exry = fabs(f->exry());
        if ((exrx<max_+verysmall_) and (exry <max_+verysmall_))
                return 0;
        if (exry>exrx)
                exrx = exry;
        return exrx-max_;
}
//---------------------------------------------------------------------------