vtf-logo

shells/applications/cylinderElastic/excentric/Excentricity.h

Go to the documentation of this file.
00001 // -*- C++ -*-
00002 //
00003 // ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00004 // 
00005 //                                   Fehmi Cirak
00006 //                        California Institute of Technology
00007 //                           (C) 2004 All Rights Reserved
00008 //
00009 // <LicenseText>
00010 //
00011 // ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00012 //
00013 #ifndef EXCENTRICITY_H
00014 #define EXCENTRICITY_H
00015 
00016 #include <functional>
00017 
00018 #include "shells/driverCC/SElementFunctors.h"
00019 
00020 
00021 namespace shells {
00022     struct SElementS;
00023 }
00024 
00025 
00026 namespace {    
00027     double excentricity(shells::SElementS* element)
00028         {
00029             // compute element center
00030             double center[3];
00031             shells::SElementTriangleAverage<double*> cfunc;
00032             cfunc(element, center);
00033             
00034             // apply some thickness perturbation in y-direction
00035             const double radius = 0.020195;
00036             const double fac = 1.0+center[1]/radius*0.05;
00037             
00038             return fac;
00039         }
00040 }
00041 
00042 #endif

Generated on Fri Aug 24 13:00:24 2007 for SFC Thin-Shell Finite Element Solver by  doxygen 1.4.7