shells/applications/cylinderElastic/excentric/Excentricity.h

Go to the documentation of this file.

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

#include <functional>

#include "shells/driverCC/SElementFunctors.h"


namespace shells {
    struct SElementS;
}


namespace {    
    double excentricity(shells::SElementS* element)
        {
            // compute element center
            double center[3];
            shells::SElementTriangleAverage<double*> cfunc;
            cfunc(element, center);
            
            // apply some thickness perturbation in y-direction
            const double radius = 0.020195;
            const double fac = 1.0+center[1]/radius*0.05;
            
            return fac;
        }
}

#endif

---