Benchmark solution ===================== The reference solution, which is analytical, is calculated as follows. Either .. image:: images/Object_17.svg :width: 30 :height: 19 .. _RefImage_Object_17.svg: the x-axis of the mass in the fixed coordinate system and .. image:: images/Object_18.svg :width: 30 :height: 19 .. _RefImage_Object_18.svg: the abscissa of the vibrating support in this same frame of reference. Initially, it is assumed that the mass is adherent to its support. It then remains so for some time after the initial moment .. image:: images/Object_19.svg :width: 30 :height: 19 .. _RefImage_Object_19.svg: . It therefore undergoes the acceleration imposed by the rigid support, i.e. .. image:: images/Object_20.svg :width: 30 :height: 19 .. _RefImage_Object_20.svg: . The tangential force exerted by the mass on the support is then .. image:: images/Object_21.svg :width: 30 :height: 19 .. _RefImage_Object_21.svg: (zero at the initial moment, which justifies the initial hypothesis that initially, the mass is adherent to its support). The mass remains adherent as long as .. image:: images/Object_22.svg :width: 30 :height: 19 .. _RefImage_Object_22.svg: . Si .. image:: images/Object_23.svg :width: 30 :height: 19 .. _RefImage_Object_23.svg: , the mass therefore remains adherent to its support indefinitely, and its movement is exactly the same as this one. By introducing the dimensionless coefficient .. image:: images/Object_24.svg :width: 30 :height: 19 .. _RefImage_Object_24.svg: , the condition of permanent adherence is written .. image:: images/Object_25.svg :width: 30 :height: 19 .. _RefImage_Object_25.svg: . The acceleration curve of the mass, as of the support, then has the following appearance as a function of time: .. image:: images/10003ECC00001DB60000123C1446B75E5FDD2E5D.svg :width: 30 :height: 19 .. _RefImage_10003ECC00001DB60000123C1446B75E5FDD2E5D.svg: As for speed, it looks as follows (the only primitive with zero mean): .. image:: images/10003BF800001DB60000123C62D873916373B464.svg :width: 30 :height: 19 .. _RefImage_10003BF800001DB60000123C62D873916373B464.svg: Si .. image:: images/Object_26.svg :width: 30 :height: 19 .. _RefImage_Object_26.svg: , there is a smaller time .. image:: images/Object_27.svg :width: 30 :height: 19 .. _RefImage_Object_27.svg: Such as .. image:: images/Object_28.svg :width: 30 :height: 19 .. _RefImage_Object_28.svg: . This smallest time is necessarily such that .. image:: images/Object_29.svg :width: 30 :height: 19 .. _RefImage_Object_29.svg: , which makes it possible to remove the absolute value in the previous expression, and to obtain an explicit expression .. image:: images/Object_30.svg :width: 30 :height: 19 .. _RefImage_Object_30.svg: . In particular, .. image:: images/Object_31.svg :width: 30 :height: 19 .. _RefImage_Object_31.svg: . After this moment, the mass slides to the left with respect to the support, so it verifies the dynamic equation .. image:: images/Object_32.svg :width: 30 :height: 19 .. _RefImage_Object_32.svg: , or .. image:: images/Object_33.svg :width: 30 :height: 19 .. _RefImage_Object_33.svg: . Its speed therefore increases linearly with time, starting at .. image:: images/Object_34.svg :width: 30 :height: 19 .. _RefImage_Object_34.svg: Of the negative value .. image:: images/Object_35.svg :width: 30 :height: 19 .. _RefImage_Object_35.svg: (in fact, .. image:: images/Object_36.svg :width: 30 :height: 19 .. _RefImage_Object_36.svg: ). .. image:: images/Object_37.svg :width: 30 :height: 19 .. _RefImage_Object_37.svg: **Movement for** .. image:: images/Object_38.svg :width: 30 :height: 19 .. _RefImage_Object_38.svg: **, "stick-slip" regime, succession of grip and slipping** Necessarily, for a certain amount of time .. image:: images/Object_39.svg :width: 30 :height: 19 .. _RefImage_Object_39.svg: satisfying .. image:: images/Object_40.svg :width: 30 :height: 19 .. _RefImage_Object_40.svg: , the speed of the mass becomes equal to the speed of the support again. At this moment, the movement becomes adherent again *if and only if* the acceleration experienced by the mass at the beginning of the adhesion is less than in absolute value than .. image:: images/Object_41.svg :width: 30 :height: 19 .. _RefImage_Object_41.svg: . The translation of this condition is examined below. To begin with, we express the value of .. image:: images/Object_42.svg :width: 30 :height: 19 .. _RefImage_Object_42.svg: . The weather .. image:: images/Object_43.svg :width: 30 :height: 19 .. _RefImage_Object_43.svg: Satisfy the equation .. image:: images/Object_44.svg :width: 30 :height: 19 .. _RefImage_Object_44.svg: , or .. image:: images/Object_45.svg :width: 30 :height: 19 .. _RefImage_Object_45.svg: , or .. image:: images/Object_46.svg :width: 30 :height: 19 .. _RefImage_Object_46.svg: . This equation, transcendent, allows the determination of .. image:: images/Object_47.svg :width: 30 :height: 19 .. _RefImage_Object_47.svg: depending on .. image:: images/Object_48.svg :width: 30 :height: 19 .. _RefImage_Object_48.svg: and .. image:: images/Object_49.svg :width: 30 :height: 19 .. _RefImage_Object_49.svg: , or finally, taking into account the expression of .. image:: images/Object_50.svg :width: 30 :height: 19 .. _RefImage_Object_50.svg: , the determination of .. image:: images/Object_51.svg :width: 30 :height: 19 .. _RefImage_Object_51.svg: depending on the physical parameters of the system .. image:: images/Object_52.svg :width: 30 :height: 19 .. _RefImage_Object_52.svg: and .. image:: images/Object_53.svg :width: 30 :height: 19 .. _RefImage_Object_53.svg: . If the acceleration of support in .. image:: images/Object_54.svg :width: 30 :height: 19 .. _RefImage_Object_54.svg: is less in absolute value than .. image:: images/Object_55.svg :width: 30 :height: 19 .. _RefImage_Object_55.svg: , the movement then remains adherent until a moment .. image:: images/Object_56.svg :width: 30 :height: 19 .. _RefImage_Object_56.svg: for which the acceleration of the support and the mass reach the value .. image:: images/Object_57.svg :width: 30 :height: 19 .. _RefImage_Object_57.svg: , an instant that, for reasons of clear symmetries on the graphs above, satisfies exactly .. image:: images/Object_58.svg :width: 30 :height: 19 .. _RefImage_Object_58.svg: . The mass then begins a sliding phase until an instant. .. image:: images/Object_59.svg :width: 30 :height: 19 .. _RefImage_Object_59.svg: , after which the movement is repeated periodically. It is understood that for sufficiently small values of .. image:: images/Object_60.svg :width: 30 :height: 19 .. _RefImage_Object_60.svg: , the movement will not be able to become a member over time .. image:: images/Object_61.svg :width: 30 :height: 19 .. _RefImage_Object_61.svg: , because the acceleration of the mass would exceed the threshold .. image:: images/Object_62.svg :width: 30 :height: 19 .. _RefImage_Object_62.svg: . So there is a critical value. .. image:: images/Object_63.svg :width: 30 :height: 19 .. _RefImage_Object_63.svg: Such as for .. image:: images/Object_64.svg :width: 30 :height: 19 .. _RefImage_Object_64.svg: , the movement of the mass passes without an adherence phase from a sliding to a sliding in the opposite direction. A reflection on the continuity of the mass velocity response function with respect to the parameter .. image:: images/Object_65.svg :width: 30 :height: 19 .. _RefImage_Object_65.svg: Show that for .. image:: images/Object_66.svg :width: 30 :height: 19 .. _RefImage_Object_66.svg: , the subsequent movement is always slippery ("slip-slip" regime, with alternating meanings). For .. image:: images/Object_67.svg :width: 30 :height: 19 .. _RefImage_Object_67.svg: , the movement periodically alternates phases of grip and sliding. The critical value .. image:: images/Object_68.svg :width: 30 :height: 19 .. _RefImage_Object_68.svg: Admits a simple analytic expression. In fact, for .. image:: images/Object_69.svg :width: 30 :height: 19 .. _RefImage_Object_69.svg: , the moments .. image:: images/Object_70.svg :width: 30 :height: 19 .. _RefImage_Object_70.svg: and .. image:: images/Object_71.svg :width: 30 :height: 19 .. _RefImage_Object_71.svg: are confused. So .. image:: images/Object_72.svg :width: 30 :height: 19 .. _RefImage_Object_72.svg: And the equation .. image:: images/Object_73.svg :width: 30 :height: 19 .. _RefImage_Object_73.svg: becomes .. image:: images/Object_74.svg :width: 30 :height: 19 .. _RefImage_Object_74.svg: . By squaring, we get .. image:: images/Object_75.svg :width: 30 :height: 19 .. _RefImage_Object_75.svg: , or .. image:: images/Object_76.svg :width: 30 :height: 19 .. _RefImage_Object_76.svg: . For .. image:: images/Object_77.svg :width: 30 :height: 19 .. _RefImage_Object_77.svg: , the movement is **only asymptotically** periodic. The suite .. image:: images/Object_78.svg :width: 30 :height: 19 .. _RefImage_Object_78.svg: Moments of change in the direction of sliding check .. image:: images/Object_79.svg :width: 30 :height: 19 .. _RefImage_Object_79.svg: when .. image:: images/Object_80.svg :width: 30 :height: 19 .. _RefImage_Object_80.svg: tends to infinity. The figure below shows the typical shape (broken line) of the mass's speed in the slip-slip situation. .. image:: images/Object_81.png :width: 36 :height: 19 .. _RefImage_Object_81.png: **Movement for** .. image:: images/Object_82.svg :width: 30 :height: 19 .. _RefImage_Object_82.svg: **: "slip-slip" regime, no adherence** Let's summarize the findings: We have the dimensionless coefficient .. image:: images/Object_83.svg :width: 30 :height: 19 .. _RefImage_Object_83.svg: and its critical value* such as .. image:: images/Object_84.svg :width: 30 :height: 19 .. _RefImage_Object_84.svg: . +------------------------------------------------------------------------------------------------------------------------------+-----------------------------------------------------------------------------------------------+ |If |the established regime is of the "stick-slip" type: alternation of adhesion and sliding phases;| + + + | | | + .. image:: images/Object_85.svg + + | :width: 30 | | + :height: 19 + + | | | + + + | | | + + + | | | +------------------------------------------------------------------------------------------------------------------------------+-----------------------------------------------------------------------------------------------+ |If |the established regime is of the "slip-slip" type: alternating permanent sliding; | + + + | | | + .. image:: images/Object_86.svg + + | :width: 30 | | + :height: 19 + + | | | + + + | | | + + + | | | + + + | | | +, + + | | | +------------------------------------------------------------------------------------------------------------------------------+-----------------------------------------------------------------------------------------------+ |If |the established regime is of the "stick" type: permanent adhesion with the base. | + + + | | | + .. image:: images/Object_87.svg + + | :width: 30 | | + :height: 19 + + | | | + + + | | | + + + | | | + + + | | | +, + + | | | +------------------------------------------------------------------------------------------------------------------------------+-----------------------------------------------------------------------------------------------+ In the following analytical calculation/*Aster* comparison results, the amplitude choices .. image:: images/Object_88.svg :width: 30 :height: 19 .. _RefImage_Object_88.svg: are such that these three situations are visited. In fact, we take :math:`m\mathrm{=}1\mathit{kg}`, :math:`g\mathrm{=}10m\mathrm{/}{s}^{2}`, :math:`\mu \mathrm{=}0.1`, :math:`{a}_{0}\mathrm{=}15m\mathrm{/}{s}^{2}`, :math:`{a}_{0}\mathrm{=}1.5m\mathrm{/}{s}^{2}`, :math:`{a}_{0}\mathrm{=}1.01m\mathrm{/}{s}^{2}`, and :math:`{a}_{0}\mathrm{=}0.99m\mathrm{/}{s}^{2}`. The wear power is physically zero during the adhesion phases. In *Code_Aster*, with the dyna_ VIBRA operator used here, adhesion is not detected because the integration of the movement is done by regularizing the law of friction. Compliance with the zero result of the wear power during adhesion phases required the introduction of a criterion on the sliding speed, so that below a certain value, it must be considered zero, and the adherent movement. You can consult the reference documentation Wear Calculation Operator/Archard Model [:ref:`R7.04.10 `]. During the sliding phases, the wear power follows the law .. image:: images/Object_101.svg :width: 30 :height: 19 .. _RefImage_Object_101.svg: , where .. image:: images/Object_102.svg :width: 30 :height: 19 .. _RefImage_Object_102.svg: is the relative sliding speed of the mass on the support. In the stick-slip regime situation, in which the movement becomes strictly periodic after a finite amount of time, the wear energy during a half-period is exactly .. image:: images/Object_103.svg :width: 30 :height: 19 .. _RefImage_Object_103.svg: .. image:: images/Object_104.svg :width: 30 :height: 19 .. _RefImage_Object_104.svg: . The transcendent formulation of .. image:: images/Object_105.svg :width: 30 :height: 19 .. _RefImage_Object_105.svg: apparently does not make it possible to simplify the expression of this energy of wear and tear. The average wear power .. image:: images/Object_106.svg :width: 30 :height: 19 .. _RefImage_Object_106.svg: Is simply wear energy .. image:: images/Object_107.svg :width: 30 :height: 19 .. _RefImage_Object_107.svg: above divided by the half-period of the response .. image:: images/Object_108.svg :width: 30 :height: 19 .. _RefImage_Object_108.svg: . In the case of a movement that is always slippery ( .. image:: images/Object_109.svg :width: 30 :height: 19 .. _RefImage_Object_109.svg: ), the integration interval to take is of the form .. image:: images/Object_110.svg :width: 30 :height: 19 .. _RefImage_Object_110.svg: with .. image:: images/Object_111.svg :width: 30 :height: 19 .. _RefImage_Object_111.svg: big enough, so that .. image:: images/Object_112.svg :width: 30 :height: 19 .. _RefImage_Object_112.svg: is close enough to the limit value .. image:: images/Object_113.svg :width: 30 :height: 19 .. _RefImage_Object_113.svg: . We can avoid the numerical calculation by recurrence of this sequence, knowing that the mean of the asymptotic speed is zero. In fact, the following .. image:: images/Object_114.svg :width: 30 :height: 19 .. _RefImage_Object_114.svg: Has a finite limit .. image:: images/Object_115.svg :width: 30 :height: 19 .. _RefImage_Object_115.svg: . The properties satisfied by .. image:: images/Object_116.svg :width: 30 :height: 19 .. _RefImage_Object_116.svg: are illustrated in the following figure: .. image:: images/10000000000002090000015B49E0265E1983B019.png :width: 36 :height: 19 .. _RefImage_10000000000002090000015B49E0265E1983B019.png: The line segment has the equation .. image:: images/Object_117.svg :width: 30 :height: 19 .. _RefImage_Object_117.svg: , And for .. image:: images/Object_118.svg :width: 30 :height: 19 .. _RefImage_Object_118.svg: , the speed .. image:: images/Object_119.svg :width: 30 :height: 19 .. _RefImage_Object_119.svg: Or take the opposite value .. image:: images/Object_120.svg :width: 30 :height: 19 .. _RefImage_Object_120.svg: , which gives the equation .. image:: images/Object_121.svg :width: 30 :height: 19 .. _RefImage_Object_121.svg: , either .. image:: images/Object_122.svg :width: 30 :height: 19 .. _RefImage_Object_122.svg: ; whose solution is .. image:: images/Object_123.svg :width: 30 :height: 19 .. _RefImage_Object_123.svg: . Note that we find although for .. image:: images/Object_124.svg :width: 30 :height: 19 .. _RefImage_Object_124.svg: , the acceleration of the support calculated at time .. image:: images/Object_125.svg :width: 30 :height: 19 .. _RefImage_Object_125.svg: Give the limit value .. image:: images/Object_126.svg :width: 30 :height: 19 .. _RefImage_Object_126.svg: . Indeed .. image:: images/Object_127.svg :width: 30 :height: 19 .. _RefImage_Object_127.svg: . In the case of the always sliding movement, the wear energy during an asymptotic period is exactly given by the formula .. image:: images/Object_128.svg :width: 30 :height: 19 .. _RefImage_Object_128.svg: which can be explained according to the previous calculation, by taking .. image:: images/Object_129.svg :width: 30 :height: 19 .. _RefImage_Object_129.svg: and .. image:: images/Object_130.svg :width: 30 :height: 19 .. _RefImage_Object_130.svg: , which gives .. image:: images/Object_131.svg :width: 30 :height: 19 .. _RefImage_Object_131.svg: , either .. image:: images/Object_132.svg :width: 30 :height: 19 .. _RefImage_Object_132.svg: . The average wear power (over a period) asymptotic is then .. image:: images/Object_133.svg :width: 30 :height: 19 .. _RefImage_Object_133.svg: . The following Maple program allows the calculation of the exact wear power in a specified time interval, as well as the plot of the graph showing the convergence of the mass speed function to a periodic limit function, for any value of physical and excitation parameters such as the regime is of the slip-slip type (:math:`\eta \mathrm{\le }{\eta }^{\text{*}}`), and the exact value of the average wear power over a period (the only useful value for what we are interested in) in the case of stick-slip. # This program does the calculation, on the transitory part # of the start of the signal, of the exact wear power, # up to a time specified at the start of the program. Digits: = 20: pi: = evalf (Pi): T: = 1: # period of support movement omega: = 2*pi/T: rpm: = 4: tmax: = 12: # duration of the transient in question ncycle: = floor (tmax/T) +2: # calculation iteration number of ti [i] and tf [i] Nmax: = 100*ncycle: # to replace the sin function with a broken line m: = 1: g: = 10: Mu: = 0.1: a0: = 1.5: eta: =mu*g/a0: omega: = 2*pi/T: etarstar: = 2/sqrt (pi^2+4): ti [1] := 1/omega*arcsin (eta): dX: = t -> -a0/omega*cos (omega*t): dxminus [0] := dX (t): lines x: = [:ref:`ti[1 `], dX (ti [1])]: Wear: = 0: # there is no wear during the adhesion phase [:ref:`0, ti[1 <0, ti[1>`]] # # Note that ti [i+1] is necessarily in the interval [:ref:`i*T-T/4,i*T+T/2 `] # and that tf [i] is necessarily in the interval [:ref:`i*T-3*T/4, i*T `]. # These two intervals overlap, but we still have tf [i] `], dX (ti [i+1])]: tinf: = max (tf [i], tmin): tsup: = min (ti [i+1], tmax): If tinf mu*g* (t-theta) +dX (theta): Vginfini: = xinfinity - xX: Eumoyana: = - int (M*G*vginInfinity (t), t=theta.. (theta+pi/omega)): Eumoyanaana: = m*g*a0/omega^2*sqrt (4-eta^2*pi^2): Pumoyana: = 2*Eumoyana/T: Pumoyanaana: = 2*Eumoyanaana/T: Powder: = Wear/ (tmax-tmin); elif (eta>etaetoile and eta<1) then # stick-slip regime lines x: = [:ref:`ti[1 `], dX (ti [1])]: dxplus [1] := mu*g* (t-ti [1]) + subs (t=ti [1], dxminus [0])): f [1] := fsolve (dX (t) =dxplus [1], t= (T-3*T/4).. T): dxplus: = unapply (dxplus [1], t): Vg: = dxplus - dX: Eu: = -int (M*g*vg (t), t=ti [1].. tf [1]): Pusuremoy: = 2*EU/T; else # permanent adherence regime Eu: = 0; End: The *Aster* solution considered is the calculation of the average wear power during a transitory phase ranging from :math:`4` to :math:`\mathrm{11,99}\mathit{secondes}` (from :math:`8\pi \mathrm{/}\omega` to :math:`24\pi \mathrm{/}\omega`). The wear energy during this transitory period differs somewhat from the average wear energy (asymptotic) over this period (both in a stick-slip and slip-slip situation). It is therefore appropriate, in order to compare it precisely with the *Aster* results, to make an exact calculation of this energy in the time interval :math:`\mathrm{[}\mathrm{4s},\mathrm{11,99}s\mathrm{]}`. For :math:`{a}_{0}\mathrm{=}\mathrm{15m}\mathrm{/}{s}^{2}`, the asymptotic mean wear power is :math:`\mathrm{15,1146144886}\mathit{Watt}` while the average wear power over the time interval :math:`\mathrm{[}\mathrm{4s},\mathrm{11,99}s\mathrm{]}` is :math:`\mathrm{15,257521794}\mathit{Watt}`. It is this last value that constitutes the reference result. **Note:** As an average power calculation, the wear power calculated over an interval does not have to increase with the duration of the interval. If you add to the interval a period of time over which there is adhesion, the average wear power will be lower. Benchmark results ---------------------- .. csv-table:: "**Acceleration value** :math:`\mathit{max.}\mathit{a0}` **(** :math:`{\mathit{ms}}^{\mathrm{-}2}` **)**", "**Average wear power value** **On the range** :math:`\mathrm{[}\mathrm{4s},\mathrm{11,99}s\mathrm{]}` **, in Watt**" "15 (slip-slip)", "15,26709959" "1.5 (stick-slip)", "0.40906245" "1.01 (stick-slip)", "2.261641E-4" "0.99 (stick)", "0" Uncertainty about the solution --------------------------- Quasi-analytical solution (presence of transcendent equations solved numerically with arbitrary precision). Bibliographical references --------------------------- 1. B. WESTERMO, F. UDWADIA: Periodic Response of a sliding oscillator system to harmonic excitation. Earthquake Engineering and Structural Dynamics Vol 14 135-146 (1983) 2. *Code_Aster* documentation [:ref:`R7.04.10 `]