Welcome to Adams/Solver Subroutines - Kxcad.net
Welcome to Adams/Solver Subroutines - Kxcad.net
Welcome to Adams/Solver Subroutines - Kxcad.net
You also want an ePaper? Increase the reach of your titles
YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.
<strong>Welcome</strong> <strong>to</strong> <strong>Adams</strong>/<strong>Solver</strong> <strong>Subroutines</strong><br />
The force has two components: a spring or stiffness component and a damping or viscous component.<br />
The stiffness component is a function of the pe<strong>net</strong>ration of the floating part in<strong>to</strong> the restricting part. The<br />
stiffness opposes the pe<strong>net</strong>ration.<br />
The damping component of the force is a function of the speed of pe<strong>net</strong>ration multiplied by a damping<br />
coefficient. The damping opposes the direction of relative motion. To prevent a discontinuity in the<br />
damping force at contact, the damping coefficient is, by definition, a cubic step function of the<br />
pe<strong>net</strong>ration. Therefore, at zero pe<strong>net</strong>ration, the damping coefficient is always zero. The damping<br />
coefficient achieves a maximum, cmax, at a user-defined pe<strong>net</strong>ration, d. Even though the points of<br />
contact between the floating part and the restricting part may change as the system moves, <strong>Adams</strong>/<strong>Solver</strong><br />
always exerts the force between the I and the J markers. Examples of systems you can model with the<br />
BISTOP function include a ball rebounding between two walls and a slider moving in a slot. The slider<br />
is the floating body, and the part containing the slot is the restricting body. As long as the slider remains<br />
within the confines of the slot, there is no force acting on the slider. But if the slider tries <strong>to</strong> move beyond<br />
the slot, a force turns on, effectively preventing the slider's escape.<br />
The following summarizes the BISTOP function:<br />
• When x1 x x2, force = 0.<br />
• When x < x1, p = x1 - x and the force is positive.<br />
• When x > x2, p = x - x2 and the force is negative.<br />
• When p < d, the damping coefficient is a cubic step function of the pe<strong>net</strong>ration.<br />
• When p d, the damping coefficient is cmax.<br />
The values of k, e, cmax, and d depend on the materials used in the two parts and on the shapes<br />
of the parts.<br />
The following equation defines BISTOP:<br />
STOP<br />
Min k ( x1 – x)<br />
e ( ⋅ – STEP ( x, x1 – d, cmax, x10) ⋅x′0<br />
, )<br />
0<br />
Max – k ( x – x2) e ⎧<br />
x < x1 ⎪<br />
= ⎨<br />
x1 ≤x ≤x2<br />
⎪<br />
⎩ ( ⋅ – STEP ( x, x2, 0, x2 + d, cmax ) ⋅x′0<br />
, ) x > x2 See an explanation of the STEP function (C++ or FORTRAN).<br />
Caution: • BISTOP is only used <strong>to</strong> determine forces or <strong>to</strong>rques.<br />
• When e is less than or equal <strong>to</strong> one, the rate of change of the force is discontinuous<br />
at contact. This may cause convergence problems.<br />
BODY_MASS_PROPERTY<br />
The BODY_MASS_PROPERTY utility subroutine accepts specifiers for either a single part or the entire<br />
ADAMS model and returns <strong>to</strong> the caller the mass properties (mass, center of mass and inertia tensor) for<br />
that part or model.<br />
19