ContactModel< T >::Impl Class Reference

#include <ContactModel_impl.h>

Collaboration diagram for ContactModel< T >::Impl:

Public Types
typedef Array< T >	TArray
typedef Array< int >	IntArray
typedef Vector< T, 3 >	VectorT3
typedef Vector< double, 3 >	VectorD3
typedef Array< VectorT3 >	VectorT3Array

Public Member Functions
	Impl (const GeomFields &geomFields, ContactFields &contactFields, const MeshList &meshes)
void	init ()
ContactModelConstants< T > &	getConstants ()
void	computeSolidSurfaceForce (const StorageSite &solidFaces, bool perUnitArea)

Private Attributes
const MeshList	_meshes
const GeomFields &	_geomFields
ContactFields &	_contactFields
ContactModelConstants< T >	_constants

Detailed Description

template<class T>
class ContactModel< T >::Impl

Definition at line 17 of file ContactModel_impl.h.

Member Typedef Documentation

template<class T >

typedef Array<int> ContactModel< T >::Impl::IntArray

Definition at line 22 of file ContactModel_impl.h.

template<class T >

typedef Array<T> ContactModel< T >::Impl::TArray

Definition at line 21 of file ContactModel_impl.h.

template<class T >

typedef Vector<double, 3> ContactModel< T >::Impl::VectorD3

Definition at line 24 of file ContactModel_impl.h.

template<class T >

typedef Vector<T,3> ContactModel< T >::Impl::VectorT3

Definition at line 23 of file ContactModel_impl.h.

template<class T >

typedef Array<VectorT3> ContactModel< T >::Impl::VectorT3Array

Definition at line 25 of file ContactModel_impl.h.

Constructor & Destructor Documentation

template<class T >

ContactModel< T >::Impl::Impl ( const GeomFields &  geomFields, ContactFields &  contactFields, const MeshList &  meshes  )

inline

Definition at line 28 of file ContactModel_impl.h.

                              :
    _meshes(meshes),
    _geomFields(geomFields),
    _contactFields(contactFields)
    { }

Member Function Documentation

template<class T >

void ContactModel< T >::Impl::computeSolidSurfaceForce ( const StorageSite &  solidFaces, bool  perUnitArea  )

inline

Definition at line 41 of file ContactModel_impl.h.

References StorageSite::getCount(), and Array< T >::zero().

   {
     const int nSolidFaces = solidFaces.getCount();

     boost::shared_ptr<VectorT3Array>
      forcePtr( new VectorT3Array(nSolidFaces));
     
     VectorT3Array& force = *forcePtr;

     force.zero();

     _contactFields.force.addArray(solidFaces,forcePtr);

     const VectorT3Array& solidFaceCentroid =
       dynamic_cast<const VectorT3Array&>(_geomFields.coordinate[solidFaces]);

     //const VectorT3Array& solidFaceArea =
     // dynamic_cast<const VectorT3Array&>(_geomFields.area[solidFaces]);

     const TArray& solidFaceAreaMag =
      dynamic_cast<const TArray&>(_geomFields.areaMag[solidFaces]);

     
     vector<NearestCell> solidFacesNearestCell(nSolidFaces);
   

     //for each face in solidFaces, find out the nearest neighbor on substrate
     /* 
     const int numMeshes = _meshes.size();
     for (int n=0; n<numMeshes; n++)
     {
         const Mesh& mesh = *_meshes[n];
         foreach(const FaceGroupPtr fgPtr, mesh.getBoundaryFaceGroups())
         {
             const FaceGroup& fg = *fgPtr;
             if (fg.id == 5){
             const StorageSite& faces = fg.site;
             const VectorT3Array& fluidFaceCoords =
               dynamic_cast<const VectorT3Array&>(_geomFields.coordinate[faces]);

             FILE *fp = fopen("./fluidface.dat","w");
             for (int i=0; i<fg.site.getCount();i++){
               fprintf(fp, "%i\t%e\t%e\t%e\n", i, fluidFaceCoords[i][0],  fluidFaceCoords[i][1], fluidFaceCoords[i][2]);
             }
             fclose(fp);

             KSearchTree searchTree(fluidFaceCoords);

             Array<int> fluidNeighbors(1);
     
             for(int f=0; f<nSolidFaces; f++)
             {
               const VectorT3& xf = solidFaceCentroid[f];
               searchTree.findNeighbors(xf, 1, fluidNeighbors);
               const int c = fluidNeighbors[0];
               const VectorT3& xc = fluidFaceCoords[c];
               const double distanceSquared = mag2(xf-xc);
               NearestCell& nc = solidFacesNearestCell[f];
               if ((nc.mesh == 0) || (nc.distanceSquared > distanceSquared))
                 {
                   nc.mesh = &mesh;
                   nc.cell = c;
                   nc.distanceSquared = distanceSquared;
                 }
             }
             }
         }
     }
     */
     //for each face in solidFaces, calculate the contact force based on the nearest distance
     const double H = 0.23e-20;
     const double B = 3529e3;
     const double alpha = 0.1127;
     const double gamma = 22.69e9;
     const double alpha01 = 1.6e-9;
     const double alpha02 = 1.99e-9; 

     const double thickness = _constants["thickness"];
     const double gap = _constants["gap"];
     
     double cloestDistance = 1.;
     cout << "gap " << gap << " thickness "<< thickness << endl;
     for(int f=0; f<nSolidFaces; f++)       {
       //const double distance = sqrt(solidFacesNearestCell[f].distanceSquared);
       const VectorT3& xf = solidFaceCentroid[f];
       const double distance = gap + thickness*0.5 + xf[2];
       if (distance < cloestDistance)
         cloestDistance = distance;
       //force[f] = -H/(6*PI) * ((1-alpha)/pow(distance,3) + alpha/pow(distance-alpha01,3)) 
       //                            + B*exp(-(distance-alpha02)*gamma);
       force[f][2] = B*exp(-(distance-alpha02)*gamma);

       if (!perUnitArea){
         force[f] *= solidFaceAreaMag[f];
       }       
     }  

     cout << "cloest distance between beam and substrate  " << cloestDistance << endl;
     /*
     FILE *fp2 = fopen("./force.dat","w");
     for(int f=0; f<nSolidFaces; f++)       {
       fprintf(fp2, "%i\t%e\t%e\t%e\n", f, force[f][0],force[f][1],force[f][2]);
     }
     fclose(fp2); 
     FILE *fp1 = fopen("./beamCoords.dat","w");
     for(int f=0; f<nSolidFaces; f++)       {
       fprintf(fp1, "%i\t%e\t%e\t%e\n", f, solidFaceCentroid[f][0], solidFaceCentroid[f][1],solidFaceCentroid[f][2]); 
     }
     fclose(fp1);
     */
   }

template<class T >

ContactModelConstants<T>& ContactModel< T >::Impl::getConstants ( )

inline

Definition at line 39 of file ContactModel_impl.h.

{return _constants;}

template<class T >

void ContactModel< T >::Impl::init ( )

inline

Definition at line 36 of file ContactModel_impl.h.

   {}

Member Data Documentation

template<class T >

ContactModelConstants<T> ContactModel< T >::Impl::_constants

private

Definition at line 156 of file ContactModel_impl.h.

template<class T >

ContactFields& ContactModel< T >::Impl::_contactFields

private

Definition at line 155 of file ContactModel_impl.h.

template<class T >

const GeomFields& ContactModel< T >::Impl::_geomFields

private

Definition at line 154 of file ContactModel_impl.h.

template<class T >

const MeshList ContactModel< T >::Impl::_meshes

private

Definition at line 153 of file ContactModel_impl.h.

---

The documentation for this class was generated from the following file:

• ContactModel_impl.h