modelBearingFriction

Coulomb friction in bearings

Extends from Modelica.Mechanics.Rotational.Interfaces.PartialElementaryTwoFlangesAndSupport2 (Partial model for a component with two rotational 1-dim. shaft flanges and a support used for textual modeling, i.e., for elementary models), Rotational.Interfaces.PartialFriction (Partial model of Coulomb friction elements), Modelica.Thermal.HeatTransfer.Interfaces.PartialElementaryConditionalHeatPortWithoutT (Partial model to include a conditional HeatPort in order to dissipate losses, used for textual modeling, i.e., for elementary models).

Information

This element describes Coulomb friction in bearings, i.e., a frictional torque acting between a flange and the housing. The positive sliding friction torque "tau" has to be defined by table "tau_pos" as function of the absolute angular velocity "w". E.g.

 w | tau
---+-----
 0 |   0
 1 |   2
 2 |   5
 3 |   8

gives the following table:

tau_pos = [0, 0; 1, 2; 2, 5; 3, 8];

Outside of the table, extrapolation through the last two table entries is used. It is assumed that the negative sliding friction force has the same characteristic with negative values. Friction is modelled in the following way:

When the absolute angular velocity "w" is not zero, the friction torque is a function of w and of a constant normal force. This dependency is defined via table tau_pos and can be determined by measurements, e.g., by driving the gear with constant velocity and measuring the needed motor torque (= friction torque).

When the absolute angular velocity becomes zero, the elements connected by the friction element become stuck, i.e., the absolute angle remains constant. In this phase the friction torque is calculated from a torque balance due to the requirement, that the absolute acceleration shall be zero. The elements begin to slide when the friction torque exceeds a threshold value, called the maximum static friction torque, computed via:

maximum_static_friction = peak * sliding_friction(w=0)  (peak >= 1)

This procedure is implemented in a "clean" way by state events and leads to continuous/discrete systems of equations if friction elements are dynamically coupled which have to be solved by appropriate numerical methods. The method is described in (see also a short sketch in UsersGuide.ModelingOfFriction):

Otter M., Elmqvist H., and Mattsson S.E. (1999):
Hybrid Modeling in Modelica based on the Synchronous Data Flow Principle. CACSD'99, Aug. 22.-26, Hawaii.

More precise friction models take into account the elasticity of the material when the two elements are "stuck", as well as other effects, like hysteresis. This has the advantage that the friction element can be completely described by a differential equation without events. The drawback is that the system becomes stiff (about 10-20 times slower simulation) and that more material constants have to be supplied which requires more sophisticated identification. For more details, see the following references, especially (Armstrong and Canudas de Wit 1996):

Armstrong B. (1991):
Control of Machines with Friction. Kluwer Academic Press, Boston MA.
Armstrong B., and Canudas de Wit C. (1996):
Friction Modeling and Compensation. The Control Handbook, edited by W.S.Levine, CRC Press, pp. 1369-1382.
Canudas de Wit C., Olsson H., Åström K.J., and Lischinsky P. (1995):
A new model for control of systems with friction. IEEE Transactions on Automatic Control, Vol. 40, No. 3, pp. 419-425.

Parameters

TypeNameDefaultDescription
BooleanuseSupport (from PartialElementaryTwoFlangesAndSupport2)false= true, if support flange enabled, otherwise implicitly grounded
Real[:,2]tau_pos[0, 1]Positive sliding friction characteristic [N.m] as function of w [rad/s] (w>=0)
SmoothnesssmoothnessSmoothness.LinearSegmentsSmoothness of interpolation in tau_pos table
Realpeak1Peak for maximum friction torque at w==0 (tau0_max = peak*tau_pos[1,2])
IntegerUnknown (from PartialFriction)3Value of mode is not known
IntegerFree (from PartialFriction)2Element is not active
IntegerForward (from PartialFriction)1w_relfric > 0 (forward sliding)
IntegerStuck (from PartialFriction)0w_relfric = 0 (forward sliding, locked or backward sliding)
IntegerBackward (from PartialFriction)-1w_relfric < 0 (backward sliding)
BooleanuseHeatPort (from PartialElementaryConditionalHeatPortWithoutT)false= true, if heatPort is enabled
Advanced
SI.AngularVelocityw_small (from PartialFriction)1.0e10Relative angular velocity near to zero if jumps due to a reinit(..) of the velocity can occur (set to low value only if such impulses can occur)

Connectors

TypeNameDefaultDescription
Flange_aflange_a (from PartialElementaryTwoFlangesAndSupport2)Flange of left shaft
Flange_bflange_b (from PartialElementaryTwoFlangesAndSupport2)Flange of right shaft
Supportsupport (from PartialElementaryTwoFlangesAndSupport2)Support/housing of component
HeatTransfer.Interfaces.HeatPort_aheatPort (from PartialElementaryConditionalHeatPortWithoutT)Optional port to which dissipated losses are transported in form of heat

Components

TypeNameDefaultDescription
SI.AngularVelocityw_relfric (from PartialFriction)Relative angular velocity between frictional surfaces
SI.AngularAccelerationa_relfric (from PartialFriction)Relative angular acceleration between frictional surfaces
SI.Torquetau0 (from PartialFriction)Friction torque for w_relfric=0 and forward sliding
SI.Torquetau0_max (from PartialFriction)Maximum friction torque for w_relfric=0 and locked
Booleanfree (from PartialFriction)= true, if frictional element is not active
Realsa (from PartialFriction)Path parameter of friction characteristic tau = f(a_relfric)
BooleanstartForward (from PartialFriction)= true, if w_relfric=0 and start of forward sliding
BooleanstartBackward (from PartialFriction)= true, if w_relfric=0 and start of backward sliding
Booleanlocked (from PartialFriction)= true, if w_rel=0 and not sliding
Integermode (from PartialFriction)Mode of friction (-1: backward sliding, 0: stuck, 1: forward sliding, 2: inactive, 3: unknown)
SI.PowerlossPower (from PartialElementaryConditionalHeatPortWithoutT)Loss power leaving component via heatPort (> 0, if heat is flowing out of component)
SI.AnglephiAngle between shaft flanges (flange_a, flange_b) and support
SI.TorquetauFriction torque
SI.AngularVelocitywAbsolute angular velocity of flange_a and flange_b
SI.AngularAccelerationaAbsolute angular acceleration of flange_a and flange_b