b2property_contact_2d_tube.H

//------------------------------------------------------------------------
// b2property_contact_2d_tube.H --
//
// written by Thomas Ludwig
//
// Copyright (c) 2016
//     SMR Engineering & Development SA
//     2502 Bienne, Switzerland
//
// All Rights Reserved.  Proprietary source code.  The contents of
// this file may not be disclosed to third parties, copied or
// duplicated in any form, in whole or in part, without the prior
// written permission of SMR.
//------------------------------------------------------------------------
 
#ifndef _B2PROPERTY_CONTACT_2D_TUBE_H_
#define _B2PROPERTY_CONTACT_2D_TUBE_H_
 
#include "elements/properties/b2property_contact_rb.H"
#include "model/b2element.H"
 
namespace b2000 {
class PropertyContact2DTube : public PropertyContactRB {
public:
    bool check(const double coor[3], double coor_target[3]) const {
        double rcoor[2];
        sub_2(coor, centre, rcoor);
 
        const double r = norm_2(rcoor);
        const bool inside = (r_min > 0 && r <= r_min);
        const bool outside = (r_max > r_min && r >= r_max);
 
        if (!inside && !outside) { return false; }
 
        const double c = (inside ? r_min : r_max) / r;
        coor_target[0] = centre[0] + centre[0] + c * rcoor[0];
        coor_target[1] = centre[1] + centre[1] + c * rcoor[1];
        coor_target[2] = coor[2];
        return true;
    }
 
protected:
    double centre[2];
    double r_min;
    double r_max;
};
}  // namespace b2000
 
#endif /* _B2PROPERTY_CONTACT_2D_TUBE_H_ */