shells/applications/cylinderElastic/excentric/ImperfectThickness.h

Go to the documentation of this file.

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

#include <cmath>
#include <map>
#include <iterator>
#include <cassert>

#include "shells/driverCC/SElementFunctors.h"


namespace shells {
    struct SElementS;
}


class ImperfectThickness {
private:
    typedef std::map<double, double>     DDMap;
    typedef std::pair<double, double>    DDPair;
    DDMap     _angleThickness;

    static const double _pi = 3.14159265;

public: 
    ImperfectThickness() {
        _angleThickness.insert(std::make_pair(-1.e-8 , 0.0008204));
        _angleThickness.insert(std::make_pair(0.25*_pi, 0.0008230));
        _angleThickness.insert(std::make_pair(0.50*_pi, 0.0008560));
        _angleThickness.insert(std::make_pair(0.75*_pi, 0.0009093));
        _angleThickness.insert(std::make_pair(1.00*_pi, 0.0009246));
        _angleThickness.insert(std::make_pair(1.25*_pi, 0.0009169));
        _angleThickness.insert(std::make_pair(1.50*_pi, 0.0008814));
        _angleThickness.insert(std::make_pair(1.75*_pi, 0.0008433));
        _angleThickness.insert(std::make_pair(2.00*_pi+1.e-8 , 0.0008204));
    }
        
    // compute the thickness at the angle phi
    double thicknessPhi(double phi)
        {               
            DDMap::iterator itl = _angleThickness.lower_bound(phi);
            assert(itl != _angleThickness.begin());
            DDMap::iterator itb = itl;
            --itb;
        
            const double slope = (itl->second-itb->second)/(itl->first-itb->first);             
            const double thickness =  itb->second +(phi-itb->first)*slope;
                
            return thickness;
        }
        
    double thickness(shells::SElementS* element)
        {
            // compute element center
            double center[3];
            shells::SElementTriangleAverage<double*> cfunc;
            cfunc(element, center);
                
            double phi = std::atan2(center[1], center[0]);

            if (center[1] < 0) phi = 2.0*_pi+phi;
                
            return thicknessPhi(phi);
        }
};


#endif

---