# Physical pendulum experiment re-investigated with an accelerometer sensor

Cyril Dauphin, Frédéric Bouquet

► To cite this version:

Cyril Dauphin, Frédéric Bouquet. Physical pendulum experiment re-investigated with an accelerometer sensor. Papers in physics, 2018, 10, pp.100008. 10.4279/PIP.100008. hal-04128659

HAL Id: hal-04128659

https://hal.science/hal-04128659

Submitted on 2 Jun 2024

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_1.png" /></figure>

# Physical pendulum experiment re-investigated with an accelerometer sensor

C. Dauphin,<fnref n="1" />,<fnref n="2" />* F. Bouquet<fnref n="3" />

We have conducted a compound pendulum experiment using Arduino and an associated two-axis accelerometer sensor as measuring device. We have shown that the use of an accelerometer to measure both radial and orbital accelerations of the pendulum at different positions along its axis offers the possibility of performing a more complex analysis compared to the usual analysis of the pendulum experiment. In this way, we have shown that this classical experiment can lead to an interesting and low-cost experiment in mechanics.

## I. Introduction

The physical pendulum experiment is the typical one to introduce the physics of oscillating systems. The usual aim of the analysis of this experiment is to determine the pendulum period and damping factor by using an angular position sensor [1, 2].

We have conducted the pendulum experiment by using a two axis accelerometer sensor. Such sensor has already been used by Fernandes et al. (2017) [3] but their study focused on the analysis of the time variation of the radial acceleration to investigate large-angle anharmonic oscillations.

Here, we have used the accelerometer to measure both radial and orbital accelerations of the pendulum at different positions along its axis, which offers the possibility of performing a more complex analysis compared to the usual single measurement of the pendulum period.

Furthermore, we use a microcontroller and an associated two axis accelerometer sensor to acquire the data. Thus, we have used this simple and low cost experiment compared to the ready to use commercial one to introduce a richer theory and data analysis tools that can lead to an interesting experiment in mechanics.

We have described the theoretical analysis of this experiment and present an example of a possible experimental setup, the analysis of the measured radial and orbital acceleration in order to acquire the moment of inertia, the center of mass, the damping factor and the period of the pendulum.

## II. Example of experimental setup

We use an Arduino [4] and a two-axis accelerometer sensor as measuring devices. The Arduino is an interesting choice for an experiment, as it is an easy-to-use and low-cost microcontroller, with a large user community. Even if Arduino was not initially developed as a physicist tool, it can be used in various contexts of experimental physics activities (e.g., see references [5–10]).

<footnote>*E-mail: cyril.dauphin@villebon-charpak.fr</footnote>

<footnote n="1">¹ Institut Villebon-Georges Charpak, Université Paris-Sud (Bât. 490) - Rue Hector Berlioz - 91400 Orsay, France.</footnote>

<footnote n="2">² Département de Physique, University Paris-Sud, Université Paris-Saclay, 91405 Orsay cedex, France.</footnote>

<footnote n="3">³ Laboratoire de Physique des Solides, CNRS, University Paris-Sud, Université Paris-Saclay, 91405 Orsay cedex, France.</footnote>
Received: 24 July 2018, Accepted: 9 October 2018

Edited by: A. Martí, M. Monteiro

Licence: Creative Commons Attribution 4.0

DOI: http://dx.doi.org/10.4279/PIP.100008

<figure><img src="image_2.png" /></figure>

# Physical pendulum experiment re-investigated with an accelerometer sensor

C. Dauphin,<fnref n="1" />,<fnref n="2" />* F. Bouquet<fnref n="3" />

We have conducted a compound pendulum experiment using Arduino and an associated two-axis accelerometer sensor as measuring device. We have shown that the use of an accelerometer to measure both radial and orbital accelerations of the pendulum at different positions along its axis offers the possibility of performing a more complex analysis compared to the usual analysis of the pendulum experiment. In this way, we have shown that this classical experiment can lead to an interesting and low-cost experiment in mechanics.

## I. Introduction

The physical pendulum experiment is the typical one to introduce the physics of oscillating systems. The usual aim of the analysis of this experiment is to determine the pendulum period and damping factor by using an angular position sensor [1, 2].

We have conducted the pendulum experiment by using a two axis accelerometer sensor. Such sensor has already been used by Fernandes et al. (2017) [3] but their study focused on the analysis of the time variation of the radial acceleration to investigate large-angle anharmonic oscillations.

Here, we have used the accelerometer to measure both radial and orbital accelerations of the pendulum at different positions along its axis, which offers the possibility of performing a more complex analysis compared to the usual single measurement of the pendulum period.

Furthermore, we use a microcontroller and an associated two axis accelerometer sensor to acquire the data. Thus, we have used this simple and low cost experiment compared to the ready to use commercial one to introduce a richer theory and data analysis tools that can lead to an interesting experiment in mechanics.

We have described the theoretical analysis of this experiment and present an example of a possible experimental setup, the analysis of the measured radial and orbital acceleration in order to acquire the moment of inertia, the center of mass, the damping factor and the period of the pendulum.

## II. Example of experimental setup

We use an Arduino [4] and a two-axis accelerometer sensor as measuring devices. The Arduino is an interesting choice for an experiment, as it is an easy-to-use and low-cost microcontroller, with a large user community. Even if Arduino was not initially developed as a physicist tool, it can be used in various contexts of experimental physics activities (e.g., see references [5–10]).

<footnote>*E-mail: cyril.dauphin@villebon-charpak.fr</footnote>

<footnote n="1">¹ Institut Villebon-Georges Charpak, Université Paris-Sud (Bât. 490) - Rue Hector Berlioz - 91400 Orsay, France.</footnote>

<footnote n="2">² Département de Physique, University Paris-Sud, Université Paris-Saclay, 91405 Orsay cedex, France.</footnote>

<footnote n="3">³ Laboratoire de Physique des Solides, CNRS, University Paris-Sud, Université Paris-Saclay, 91405 Orsay cedex, France.</footnote>
<figure><img src="image_2.png" /><figcaption>Figure 1: (a) An example of an experimental setup. (b) An accelerometer sensor attached to the centerline of the pendulum, with one of its axes parallel to the centerline.</figcaption></figure>

<figure><img src="image_3.png" /><figcaption>Figure 2: Radial (red) and orbital (black) accelerations measured by the accelerometer obtained with the experimental setup displayed in Fig. 1. Inset shows the temporal evolution of both accelerations between 40 s and 45 s. In order to have both curves on the same plot, orbital acceleration is shifted by a constant offset of 9 in the inset.</figcaption></figure>

The experimental setup that we used here is shown in Fig. 1(a). The accelerometer sensor is a microelectromechanical inertial sensor which is precisely calibrated by using the values +g, 0g and −g for each axis with g = 9.8 m s⁻².

The pendulum used in this experiment is composed of a bar on which masses can be attached to different positions. The accelerometer is attached to the bar and positioned in such a way that one of its measurement axes lies parallel to the bar (Fig. 1(b)). Special care should be taken so that the wires connecting the accelerometer to the board are flexible enough in order not to damp the pendulum.

Figure 2 shows an example of data acquired by the accelerometer. The main features of the graph are:

* the radial acceleration measurement decreases and goes to g as t approaches the infinity.

* when the radial acceleration reaches its maximum values, the orbital acceleration is essentially equal to zero (inset of Fig. 2).

* the radial acceleration is asymmetric about the straight line a = g contrary to the orbital acceleration that is symmetric about the line a = 0.

* the period of the radial acceleration oscillations is twice that of the orbital acceleration oscillations (inset of Fig. 2).

In the next section, we will present the theory that explains these main features.

### III. Theory

#### i. Expression of the acceleration components measured by the accelerometer sensor

Figure 3 shows a pendulum sketch with notations that will be used throughout this paper. O is the pivot point, G is the pendulum mass center and A the accelerometer position. L and r stands for the distance between O and G and O and A, respectively. We note $I_λ$ the moment of inertia of the pendulum about the Oz axis.

Applying the angular momentum theorem to the pendulum and considering viscous damping leads to:

$$I_{\lambda}\ddot{\theta} = -MgL \sin \theta - \gamma\dot{\theta} \quad (1)$$

where θ is the angle between the pendulum axis OG and the vertical axis, M is the mass of the system and γ is the coefficient of friction. We note
<figure><img src="image_3.png" /><figcaption>Figure 1: (a) An example of an experimental setup. (b) An accelerometer sensor attached to the centerline of the pendulum, with one of its axes parallel to the centerline.</figcaption></figure>

<figure><img src="image_4.png" /><figcaption>Figure 2: Radial (red) and orbital (black) accelerations measured by the accelerometer obtained with the experimental setup displayed in Fig. 1. Inset shows the temporal evolution of both accelerations between 40 s and 45 s. In order to have both curves on the same plot, orbital acceleration is shifted by a constant offset of 9 in the inset.</figcaption></figure>

The experimental setup that we used here is shown in Fig. 1(a). The accelerometer sensor is a microelectromechanical inertial sensor which is precisely calibrated by using the values +g, 0g and −g for each axis with g = 9.8 m s⁻².

The pendulum used in this experiment is composed of a bar on which masses can be attached to different positions. The accelerometer is attached to the bar and positioned in such a way that one of its measurement axes lies parallel to the bar (Fig. 1(b)). Special care should be taken so that the wires connecting the accelerometer to the board are flexible enough in order not to damp the pendulum.

Figure 2 shows an example of data acquired by the accelerometer. The main features of the graph are:

* the radial acceleration measurement decreases and goes to g as t approaches the infinity.

* when the radial acceleration reaches its maximum values, the orbital acceleration is essentially equal to zero (inset of Fig. 2).

* the radial acceleration is asymmetric about the straight line a = g contrary to the orbital acceleration that is symmetric about the line a = 0.

* the period of the radial acceleration oscillations is twice that of the orbital acceleration oscillations (inset of Fig. 2).

In the next section, we will present the theory that explains these main features.

### III. Theory

#### i. Expression of the acceleration components measured by the accelerometer sensor

Figure 3 shows a pendulum sketch with notations that will be used throughout this paper. O is the pivot point, G is the pendulum mass center and A the accelerometer position. L and r stands for the distance between O and G and O and A, respectively. We note $I_λ$ the moment of inertia of the pendulum about the Oz axis.

Applying the angular momentum theorem to the pendulum and considering viscous damping leads to:

$$I_{\lambda} \ddot{\theta} = -MgL \sin \theta - \gamma \dot{\theta} \quad (1)$$

where θ is the angle between the pendulum axis OG and the vertical axis, M is the mass of the system and γ is the coefficient of friction. We note
<figure><img src="image_4.png" /><figcaption>Figure 3: Sketch of the pendulum and notations used in the text. The pivot is at O, the center of mass at G and the accelerometer sensor is at A.</figcaption></figure>

$I_λ = αML^2$ where the numerical factor $α$ depends on the type of pendulum ($α = 1$ for a simple pendulum and $α ≠ 1$ for a physical pendulum). We introduce the damping factor $κ = \frac{γ}{2I_λ}$ to obtain:

$$ \ddot{θ} = -\frac{g}{αL} \sin θ - 2κ\dot{θ} \quad (2) $$

Forces acting on the proof mass ($m_{sensor}$) inside the accelerometer are its weight and the inertial force due to its movement. Thus, the radial and orbital components of these forces in the non-inertial reference frame of the pendulum are given by:

$$ F_r = m_{sensor}r\dot{θ}^2 + m_{sensor}g \cos θ \quad (3) $$

$$ F_θ = -m_{sensor}r\ddot{θ} - m_{sensor}g \sin θ \quad (4) $$

The acceleration components, as measured by the accelerometer at a distance $r$ from the pivot, are then given by:

$$ a_r = r\dot{θ}^2 + g \cos θ \quad (5) $$

$$ a_θ = -r\ddot{θ} - g \sin θ \quad (6) $$

Including the expression of $\ddot{θ}$ (Eq. (2)) into the expression of $a_θ$ gives:

$$ a_θ = \frac{g}{α} \left( \frac{r}{L} - α \right) \sin θ + 2κr\dot{θ} \quad (7) $$

We choose $θ = θ_0$ and $\dot{θ} = 0$ as initial conditions of the pendulum movement. We only consider here

small damping ($κ << ω$) and the small angle approximation ($θ_0 << \frac{π}{2}$), which leads from Eq. (2) to:

$$ θ = θ_0 e^{-κt} \cos(ωt) \quad (8) $$

$$ \dot{θ} = -ωθ_0 e^{-κt} \sin(ωt) \quad (9) $$

$$ \ddot{θ} = -ω^2θ_0 e^{-κt} \cos(ωt) \quad (10) $$

with the pendulum angular frequency $ω = \sqrt{ω_0^2 - κ^2}$ and $ω_0 = \sqrt{g/αL}$. Thus, the acceleration components measured by the accelerometer for $κ << ω$ and $θ_0 << \frac{π}{2}$ are given by:

$$ a_r = rω^2θ_0^2 e^{-2κt} \sin^2(ωt) + g \cos(θ_0e^{-κt} \cos(ωt)) \quad (11) $$

$$ a_θ = -2κrωθ_0e^{-κt} \sin(ωt) + \frac{g}{α} \left( \frac{r}{L} - α \right) \sin(θ_0e^{-κt} \cos(ωt)) \quad (12) $$

With $θ_0 << \frac{π}{2}$, we can further approximate these two expressions by:

$$ a_r = rω^2θ_0^2 e^{-2κt} \sin^2(ωt) + g(1 - θ_0^2e^{-2κt} \cos^2(ωt)) \quad (13) $$

$$ a_θ = -2κrωθ_0e^{-κt} \sin(ωt) + \frac{g}{α} \left( \frac{r}{L} - α \right) θ_0 e^{-κt} \cos(ωt) \quad (14) $$

Some comments can be made about Eqs. (11) and (12). We first focus on the radial acceleration $a_r$.

* $\cos(\cos(x))$ and $\sin^2(x)$ are both $π$-periodic functions. So $a_r$ is a $\frac{T}{2}$ periodic function with $T = \frac{2π}{ω}$ being the pendulum period. Indeed, the pendulum reaches its maximum velocity and so its maximum radial acceleration each time $θ$ is equal to 0.

* $g \cos(θ_0e^{-κt} \cos(ωt))$ varies between between $g \cos(θ_0)$ and $g$ (blue curve in Fig. 4(b)) and goes to $g$ as $t$ approaches infinity. Physically, this function represents the projection of $\vec{g}$ onto the pendulum axis OG.

* $rω^2θ_0^2e^{-2κt} \sin^2(ωt)$ varies between $rω^2θ_0^2$ and 0 and the upper envelope of this function (red
<figure><img src="image_5.png" /><figcaption>Figure 3: Sketch of the pendulum and notations used in the text. The pivot is at O, the center of mass at G and the accelerometer sensor is at A.</figcaption></figure>

$I_λ = αML^2$ where the numerical factor $α$ depends on the type of pendulum ($α = 1$ for a simple pendulum and $α ≠ 1$ for a physical pendulum). We introduce the damping factor $κ = \frac{γ}{2I_λ}$ to obtain:

$$ \ddot{θ} = -\frac{g}{αL} \sin θ - 2κ\dot{θ} \quad (2) $$

Forces acting on the proof mass ($m_{sensor}$) inside the accelerometer are its weight and the inertial force due to its movement. Thus, the radial and orbital components of these forces in the non-inertial reference frame of the pendulum are given by:

$$ F_r = m_{sensor}r\dot{θ}^2 + m_{sensor}g \cos θ \quad (3) $$

$$ F_θ = -m_{sensor}r\ddot{θ} - m_{sensor}g \sin θ \quad (4) $$

The acceleration components, as measured by the accelerometer at a distance $r$ from the pivot, are then given by:

$$ a_r = r\dot{θ}^2 + g \cos θ \quad (5) $$

$$ a_θ = -r\ddot{θ} - g \sin θ \quad (6) $$

Including the expression of $\ddot{θ}$ (Eq. (2)) into the expression of $a_θ$ gives:

$$ a_θ = \frac{g}{α} \left( \frac{r}{L} - α \right) \sin θ + 2κr\dot{θ} \quad (7) $$

We choose $θ = θ_0$ and $\dot{θ} = 0$ as initial conditions of the pendulum movement. We only consider here

small damping ($κ << ω$) and the small angle approximation ($θ_0 << \frac{π}{2}$), which leads from Eq. (2) to:

$$ θ = θ_0 e^{-κt} \cos(ωt) \quad (8) $$

$$ \dot{θ} = -ωθ_0 e^{-κt} \sin(ωt) \quad (9) $$

$$ \ddot{θ} = -ω^2θ_0 e^{-κt} \cos(ωt) \quad (10) $$

with the pendulum angular frequency $ω = \sqrt{ω_0^2 - κ^2}$ and $ω_0 = \sqrt{g/αL}$. Thus, the acceleration components measured by the accelerometer for $κ << ω$ and $θ_0 << \frac{π}{2}$ are given by:

$$ a_r = rω^2θ_0^2 e^{-2κt} \sin^2(ωt) + g \cos(θ_0 e^{-κt} \cos(ωt)) \quad (11) $$

$$ a_θ = -2κrωθ_0 e^{-κt} \sin(ωt) + \frac{g}{α} \left( \frac{r}{L} - α \right) \sin(θ_0 e^{-κt} \cos(ωt)) \quad (12) $$

With $θ_0 << \frac{π}{2}$, we can further approximate these two expressions by:

$$ a_r = rω^2θ_0^2 e^{-2κt} \sin^2(ωt) + g (1 - θ_0^2 e^{-2κt} \cos^2(ωt)) \quad (13) $$

$$ a_θ = -2κrωθ_0 e^{-κt} \sin(ωt) + \frac{g}{α} \left( \frac{r}{L} - α \right) θ_0 e^{-κt} \cos(ωt) \quad (14) $$

Some comments can be made about Eqs. (11) and (12). We first focus on the radial acceleration $a_r$.

* $\cos(\cos(x))$ and $\sin^2(x)$ are both $π$-periodic functions. So $a_r$ is a $\frac{T}{2}$ periodic function with $T = \frac{2π}{ω}$ being the pendulum period. Indeed, the pendulum reaches its maximum velocity and so its maximum radial acceleration each time $θ$ is equal to 0.

* $g \cos(θ_0 e^{-κt} \cos(ωt))$ varies between between $g \cos(θ_0)$ and $g$ (blue curve in Fig. 4(b)) and goes to $g$ as $t$ approaches infinity. Physically, this function represents the projection of $\vec{g}$ onto the pendulum axis OG.

* $rω^2θ_0^2 e^{-2κt} \sin^2(ωt)$ varies between $rω^2θ_0^2$ and 0 and the upper envelope of this function (red
<figure><img src="image_5.png" /><figcaption>Figure 4: (a) Calculated radial acceleration versus time. (b) In red, graph of the $r\omega^2\theta_0^2e^{-2\kappa t}\sin^2(\omega t) + g$ contribution; in blue, graph of the $g\cos(\theta_0e^{-\kappa t}\cos(\omega t))$ contribution. L = 30 cm, $\theta_0 = 15^\circ$, r = 45 cm, $\kappa = 0.1 s^{-1}$ and $\alpha = 1$ for both panels.</figcaption></figure>

<figure><img src="image_6.png" /><figcaption>Figure 5: (a) Calculated orbital acceleration versus time. (b) In red, graph of the $\frac{g}{\alpha}(\frac{r}{L}-\alpha)\sin(\theta_0 e^{-\kappa t}\cos(\omega t))$ contribution; in blue, graph of the part $-2\kappa r\omega\theta_0 e^{-\kappa t}\sin(\omega t)$. L = 30 cm, $\theta_0 = 15^\circ$, r = 45 cm, $\kappa = 0.1 s^{-1}$ and $\alpha = 1$ for both panels.</figcaption></figure>

curve in Fig. 4(b), note that it has been displaced by $g$ decreases exponentially as $e^{-2\kappa t}$. Physically, this function represents the acceleration due to the radial centrifugal force felt by the sensor (as expected, the acceleration due to the radial centrifugal force is maximum when the pendulum is vertical). This point and the previous one explains the asymmetry of the radial acceleration about the straight line $a_r = g$, as observed in Fig. 2.

We can now study the orbital acceleration $a_\theta$:

* $\sin(\cos(x))$ and $\sin(x)$ are both $2\pi$-periodic functions. Therefore, $a_\theta$ is a T-periodic func-

tion with $T = \frac{2\pi}{\omega}$. Indeed, the angular acceleration $\ddot{\theta}$ reaches its highest and lowest values each time the pendulum passes through its highest points. This point explains that the period of the radial acceleration is twice that of the orbital one.

* $\frac{g}{\alpha}(\frac{r}{L}-\alpha)\sin(\theta_0 e^{-\kappa t}\cos(\omega t))$ goes to 0 as $t$ approaches infinity and the $\sin(\cos(x))$ function implies that the upper and lower envelope of this part of the $a_\theta$ function are symmetric about the straight line of equation $a = 0$ (red curve in Fig. 5(b)).

* $-2\kappa r\omega\theta_0 e^{-\kappa t}\sin(\omega t)$ goes to zero as $t$ ap-
<figure><img src="image_6.png" /><figcaption>Figure 4: (a) Calculated radial acceleration versus time. (b) In red, graph of the $r\omega^2\theta_0^2e^{-2\kappa t}\sin^2(\omega t) + g$ contribution; in blue, graph of the $g\cos(\theta_0e^{-\kappa t}\cos(\omega t))$ contribution. L = 30 cm, $\theta_0 = 15^\circ$, r = 45 cm, $\kappa = 0.1 s^{-1}$ and $\alpha = 1$ for both panels.</figcaption></figure>

<figure><img src="image_7.png" /><figcaption>Figure 5: (a) Calculated orbital acceleration versus time. (b) In red, graph of the $\frac{g}{\alpha}(\frac{r}{L}-\alpha)\sin(\theta_0 e^{-\kappa t}\cos(\omega t))$ contribution; in blue, graph of the part $-2\kappa r\omega\theta_0 e^{-\kappa t}\sin(\omega t)$. L = 30 cm, $\theta_0 = 15^\circ$, r = 45 cm, $\kappa = 0.1 s^{-1}$ and $\alpha = 1$ for both panels.</figcaption></figure>

curve in Fig. 4(b), note that it has been displaced by $g$ decreases exponentially as $e^{-2\kappa t}$. Physically, this function represents the acceleration due to the radial centrifugal force felt by the sensor (as expected, the acceleration due to the radial centrifugal force is maximum when the pendulum is vertical). This point and the previous one explains the asymmetry of the radial acceleration about the straight line $a_r = g$, as observed in Fig. 2.

We can now study the orbital acceleration $a_\theta$:

* $\sin(\cos(x))$ and $\sin(x)$ are both $2\pi$-periodic functions. Therefore, $a_\theta$ is a T-periodic func-

tion with $T = \frac{2\pi}{\omega}$. Indeed, the angular acceleration $\ddot{\theta}$ reaches its highest and lowest values each time the pendulum passes through its highest points. This point explains that the period of the radial acceleration is twice that of the orbital one.

* $\frac{g}{\alpha}(\frac{r}{L}-\alpha)\sin(\theta_0 e^{-\kappa t}\cos(\omega t))$ goes to 0 as $t$ approaches infinity and the $\sin(\cos(x))$ function implies that the upper and lower envelope of this part of the $a_\theta$ function are symmetric about the straight line of equation $a = 0$ (red curve in Fig. 5(b)).

* $-2\kappa r\omega\theta_0 e^{-\kappa t}\sin(\omega t)$ goes to zero as $t$ ap-
<figure><img src="image_7.png" /><figcaption>Figure 6: Radial (red) and orbital (black) accelerations obtained from Eqs. (11) and (12) with $\theta_0 = 20^\circ$, $L = 30$ cm, $r = 34$ cm, $\alpha = 1.2$ and $\kappa = 0.01 \text{ s}^{-1}$. Inset shows the temporal evolution of radial (red) and orbital (black) accelerations between 40 s and 45 s. In order to have both curves on the same plot, orbital acceleration is shifted by a constant offset of 9 in the inset.</figcaption></figure>

proaches infinity. The upper and lower envelopes of this part of the $a_\theta$ function are symmetric about the straight line $a_\theta = 0$ (blue curve in Fig. 5(b)) and decrease exponentially as $e^{-\kappa t}$.

Figure 6 displays the radial (red) and orbital (black) accelerations obtained from Eqs. (11) and (12). We can see that the main features inferred from the data (Fig. 2) are well reproduced.

## ii. Experimentally accessible quantities

The parameters describing the pendulum and its motion can be derived from the measurements of $a_r$ and $a_\theta$.

* $a_r(t = 0) = a_{r,min} = g \cos \theta_0$. Thus, the measured value of $a_r$ at $t = 0$ leads to the value of $\theta_0$.

* A fit to the measured $a_r$ upper envelope with an exponential function allows us to determine the damping factor $\kappa$ of the pendulum from Eqs. (11) and (12). While, for small deflection angles $\theta$, see Eqs. (13) and (14), exponential fit of any envelope of measured $a_r$ or $a_\theta$ allows us to determine the damping factor.

* The sensor position OA can be measured with great accuracy. Thus, the position of the pendulum center of mass and moment of inertia are the two quantities which are difficult to determine experimentally. Here, we use the fit of the temporal evolution of $a_r$ and $a_\theta$ to determine the product $\alpha L$.

## iii. Impact of $\alpha$ and $r$ on the measured radial and orbital accelerations

Figure 7 displays the evolution of the radial and orbital accelerations with time for different values of the moment of inertia (from $\alpha = 1$ (Panel (a)) to $\alpha = 2.5$ (Panel (d))). At the difference of the radial acceleration, we can see that the orbital acceleration depends strongly on $\alpha$. Indeed, $a_\theta$ expression at $t = 0$ leads to $a_\theta(t = 0) = (\frac{r}{\alpha L} - 1) g \sin \theta_0$, which is an inverse function of $\alpha$. We can also note that $a_\theta(t = 0) < 0$ if $\alpha > \frac{r}{L}$ and $a_\theta(t = 0) > 0$ if $\alpha < \frac{r}{L}$. Thus, value of $a_\theta$ at $t = 0$ gives information on the $\alpha$ value.

Figure 8 shows the evolution of the radial and orbital accelerations with time for different values of $r$. The distance OA increases from Panel (a) to Panel (d). As expected, the amplitude of the radial acceleration increases with larger OA values as the centrifugal force acting on the proof mass increases and the amplitude of the orbital acceleration decreases with larger OA values as the rate of variation of $\dot{\theta}$ decreases with this distance.

In particular, acceleration components measured by the accelerometer attached to the position of the point O are given by:

$$a_r = g \cos (\theta_0 e^{-\kappa t} \cos(\omega t)) \quad (15)$$

$$a_\theta = -g \sin (\theta_0 e^{-\kappa t} \cos(\omega t)) \quad (16)$$

In this case, the only force acting on the mass inside the accelerometer is its weight and the expressions of the acceleration do not depend on $\alpha$. Accelerometer sensor is used in this case as an angular position sensor.

We also note that the acceleration components measured by the accelerometer attached to the position $r = L\alpha$ are given by:
<figure><img src="image_8.png" /><figcaption>Figure 6: Radial (red) and orbital (black) accelerations obtained from Eqs. (11) and (12) with $\theta_0 = 20^\circ$, $L = 30$ cm, $r = 34$ cm, $\alpha = 1.2$ and $\kappa = 0.01 s^{-1}$. Inset shows the temporal evolution of radial (red) and orbital (black) accelerations between 40 s and 45 s. In order to have both curves on the same plot, orbital acceleration is shifted by a constant offset of 9 in the inset.</figcaption></figure>

proaches infinity. The upper and lower envelopes of this part of the $a_\theta$ function are symmetric about the straight line $a_\theta = 0$ (blue curve in Fig. 5(b)) and decrease exponentially as $e^{-\kappa t}$.

Figure 6 displays the radial (red) and orbital (black) accelerations obtained from Eqs. (11) and (12). We can see that the main features inferred from the data (Fig. 2) are well reproduced.

## ii. Experimentally accessible quantities

The parameters describing the pendulum and its motion can be derived from the measurements of $a_r$ and $a_\theta$.

* $a_r(t = 0) = a_{r,min} = g \cos \theta_0$. Thus, the measured value of $a_r$ at $t = 0$ leads to the value of $\theta_0$.

* A fit to the measured $a_r$ upper envelope with an exponential function allows us to determine the damping factor $\kappa$ of the pendulum from Eqs. (11) and (12). While, for small deflection angles $\theta$, see Eqs. (13) and (14), exponential fit of any envelope of measured $a_r$ or $a_\theta$ allows us to determine the damping factor.

* The sensor position OA can be measured with great accuracy. Thus, the position of the pendulum center of mass and moment of inertia are the two quantities which are difficult to determine experimentally. Here, we use the fit of the temporal evolution of $a_r$ and $a_\theta$ to determine the product $\alpha L$.

## iii. Impact of $\alpha$ and $r$ on the measured radial and orbital accelerations

Figure 7 displays the evolution of the radial and orbital accelerations with time for different values of the moment of inertia (from $\alpha = 1$ (Panel (a)) to $\alpha = 2.5$ (Panel (d))). At the difference of the radial acceleration, we can see that the orbital acceleration depends strongly on $\alpha$. Indeed, $a_\theta$ expression at $t = 0$ leads to $a_\theta(t = 0) = (\frac{r}{\alpha L} - 1) g \sin \theta_0$, which is an inverse function of $\alpha$. We can also note that $a_\theta(t = 0) < 0$ if $\alpha > \frac{r}{L}$ and $a_\theta(t = 0) > 0$ if $\alpha < \frac{r}{L}$. Thus, value of $a_\theta$ at $t = 0$ gives information on the $\alpha$ value.

Figure 8 shows the evolution of the radial and orbital accelerations with time for different values of $r$. The distance OA increases from Panel (a) to Panel (d). As expected, the amplitude of the radial acceleration increases with larger OA values as the centrifugal force acting on the proof mass increases and the amplitude of the orbital acceleration decreases with larger OA values as the rate of variation of $\dot{\theta}$ decreases with this distance.

In particular, acceleration components measured by the accelerometer attached to the position of the point O are given by:

$$a_r = g \cos (\theta_0 e^{-\kappa t} \cos(\omega t)) \quad (15)$$

$$a_\theta = -g \sin (\theta_0 e^{-\kappa t} \cos(\omega t)) \quad (16)$$

In this case, the only force acting on the mass inside the accelerometer is its weight and the expressions of the acceleration do not depend on $\alpha$. Accelerometer sensor is used in this case as an angular position sensor.

We also note that the acceleration components measured by the accelerometer attached to the position $r = L\alpha$ are given by:
<figure><img src="image_8.png" /><figcaption>Figure 7: The calculated radial (red) and orbital (black) acceleration with time for $α = 1$ (a), $α = 1.5$ (b), $α = 2.0$ (c) and $α = 2.5$ (d). In order to have both curves on the same plot, orbital acceleration is shifted by a constant offset of 9.2. L = 30 cm, $θ_0 = 15°$, r = 45 cm and $κ = 0.1$ for all panels.</figcaption></figure>

<figure><img src="image_9.png" /><figcaption>Figure 8: The calculated radial (red) and orbital (black) accelerations with time for different positions of the accelerometer relative to the center of mass G. In order to have both curves on the same plot, the orbital acceleration is shifted by a constant offset of 9.2. The accelerometer positions are (a) r = 0 cm, (b) r = 30 cm, (c) r = 45 cm, (d) r = 60 cm. L = 30 cm, $θ_0 = 15°$, $α = 2$ and $κ = 0.1$ for all panels.</figcaption></figure>
<figure><img src="image_9.png" /><figcaption>Figure 7: The calculated radial (red) and orbital (black) acceleration with time for $α = 1$ (a), $α = 1.5$ (b), $α = 2.0$ (c) and $α = 2.5$ (d). In order to have both curves on the same plot, orbital acceleration is shifted by a constant offset of 9.2. L = 30 cm, $θ_0 = 15°$, r = 45 cm and $κ = 0.1$ for all panels.</figcaption></figure>

<figure><img src="image_10.png" /><figcaption>Figure 8: The calculated radial (red) and orbital (black) accelerations with time for different positions of the accelerometer relative to the center of mass G. In order to have both curves on the same plot, the orbital acceleration is shifted by a constant offset of 9.2. The accelerometer positions are (a) r = 0 cm, (b) r = 30 cm, (c) r = 45 cm, (d) r = 60 cm. L = 30 cm, $θ_0 = 15°$, $α = 2$ and $κ = 0.1$ for all panels.</figcaption></figure>
<figure><img src="image_10.png" /><figcaption>Figure 9: Experimental setup used to study a bar pendulum.</figcaption></figure>

<figure><img src="image_11.png" /><figcaption>Figure 10: Radial (red) and orbital (green) accelerations measured by the accelerometer in the configuration shown in Fig. 9. Solid lines are the best fit of the data with $\kappa$ and $\alpha$ as free parameters. The best fit is given by $1/\kappa = 12.2$ s and $\alpha = 1.56$.</figcaption></figure>

$$a_r = \alpha\omega^2\theta_0^2e^{-2\kappa t}\sin^2(\omega t) + g\cos(\theta_0e^{-\kappa t}\cos(\omega t)) \quad (17)$$

$$a_{\theta} = -2\kappa\alpha L\omega\theta_0 e^{-\kappa t}\sin(\omega t) \quad (18)$$

In this case, component of $\ddot{\theta}$ due to the gravity force is counterbalanced by the orbital component of the force of gravity acting on the accelerometer sensor. Thus, $a_{\theta}$ is then directly proportional to the pendulum angular velocity.

After having shown that Eqs. (11) and (12) explain the features observed experimentally, we will now use them to retrieve the pendulum parameters $\kappa$, $\alpha$ and $L$ for different experimental setups.

# IV. Example of data analysis

## i. Example with a pendulum bar

We first focus on the pendulum shown in Fig. 9 to derive its physical parameters. We use a bar of a 45 cm length and mass 45 g with the pivot at 16.5 cm from the center of mass and the accelerometer at 33 cm.

Figure 10 shows the radial and orbital accelerations measured by the accelerometer after the pendulum has been displaced from the equilibrium position to an initial angle of 22°.

We have analyzed the data using Eqs. (11) and (12) with $\kappa$ and $\alpha$ as free parameters. The results of the fit are shown as black curves in Fig. 10. We have derived the best fit for $1/\kappa = 12.2$ s and $\alpha = 1.56$.

Assuming the pendulum to be a simple homogeneous slab, we calculate that the inertia moment of the bar about the rotation axis $\lambda$ is equal to $1.93 \times 10^{-3}$ kg m$^2$, which leads to $\alpha = 1.62$. The holes in the bar, which are used for attaching the masses, together with the mass of the accelerometer, explain the difference between this and with the $\alpha$ value obtained from the data fit. This experiment allows us to determine with a good accuracy the moment of inertia of the pendulum.

## ii. Retrieval of the mass center position

The pendulum of the previous subsection is a symmetric bar, thus, its mass center position can be determined precisely. In the general case, the mass center position can be difficult to determine and we can fit the data by using Eqs. (11) and (12) with $\kappa$, $\alpha$ and $L$ as free parameters and infer the position of the center of mass.

As an example of such analysis, we use a pendulum.
<figure><img src="image_11.png" /><figcaption>Figure 9: Experimental setup used to study a bar pendulum.</figcaption></figure>

<figure><img src="image_12.png" /><figcaption>Figure 10: Radial (red) and orbital (green) accelerations measured by the accelerometer in the configuration shown in Fig. 9. Solid lines are the best fit of the data with $\kappa$ and $\alpha$ as free parameters. The best fit is given by $1/\kappa = 12.2$ s and $\alpha = 1.56$.</figcaption></figure>

$$a_r = \alpha\omega^2\theta_0^2e^{-2\kappa t}\sin^2(\omega t) + g\cos(\theta_0e^{-\kappa t}\cos(\omega t)) \quad (17)$$

$$a_{\theta} = -2\kappa\alpha L\omega\theta_0 e^{-\kappa t}\sin(\omega t) \quad (18)$$

In this case, component of $\ddot{\theta}$ due to the gravity force is counterbalanced by the orbital component of the force of gravity acting on the accelerometer sensor. Thus, $a_{\theta}$ is then directly proportional to the pendulum angular velocity.

After having shown that Eqs. (11) and (12) explain the features observed experimentally, we will now use them to retrieve the pendulum parameters $\kappa$, $\alpha$ and $L$ for different experimental setups.

# IV. Example of data analysis

## i. Example with a pendulum bar

We first focus on the pendulum shown in Fig. 9 to derive its physical parameters. We use a bar of a 45 cm length and mass 45 g with the pivot at 16.5 cm from the center of mass and the accelerometer at 33 cm.

Figure 10 shows the radial and orbital accelerations measured by the accelerometer after the pendulum has been displaced from the equilibrium position to an initial angle of 22°.

We have analyzed the data using Eqs. (11) and (12) with $\kappa$ and $\alpha$ as free parameters. The results of the fit are shown as black curves in Fig. 10. We have derived the best fit for $1/\kappa = 12.2$ s and $\alpha = 1.56$.

Assuming the pendulum to be a simple homogeneous slab, we calculate that the inertia moment of the bar about the rotation axis $\lambda$ is equal to $1.93 \times 10^{-3}$ kg m$^2$, which leads to $\alpha = 1.62$. The holes in the bar, which are used for attaching the masses, together with the mass of the accelerometer, explain the difference between this and with the $\alpha$ value obtained from the data fit. This experiment allows us to determine with a good accuracy the moment of inertia of the pendulum.

## ii. Retrieval of the mass center position

The pendulum of the previous subsection is a symmetric bar, thus, its mass center position can be determined precisely. In the general case, the mass center position can be difficult to determine and we can fit the data by using Eqs. (11) and (12) with $\kappa$, $\alpha$ and $L$ as free parameters and infer the position of the center of mass.

As an example of such analysis, we use a pendulum.
<figure><img src="image_12.png" /><figcaption>Figure 11: Experimental setup used to study compound pendula, with added masses.</figcaption></figure>

dulum composed of a bar on which several masses can be attached to at different positions to acquired data that corresponds to pendulums of different values of $\alpha$ and $L$ (Fig. 11 (a) and (b)).

Figure 12 (a) and (b) displays the radial and orbital accelerations measured by the accelerometer after each pendulum shown in Fig. 11 has been displaced by an initial angle of 20.5°.

We fit the data by using Eqs. (11) and (12) with $\kappa$, $\alpha$ and $L$ as free parameters. Results of the fit are shown as the black curves in the insets of Fig. 12. Best fit are obtained with $1/\kappa = 111$ s, $\alpha = 1.162$ and $L = 0.293$ m for the pendulum of Fig. 11(a) and $1/\kappa = 101$ s, $\alpha = 6.55$ and $L = 0.074$ m for the pendulum of Fig. 11(b).

For the examples displayed in Fig. 12, the evolution of the orbital acceleration from Panels (a) to (b) shows an increase of $\alpha$, which is consistent with the fact that pendulum configuration goes from a configuration close to a simple pendulum ($\alpha \approx 1$ Fig. 11(a)) to a compound pendulum ($\alpha > 1$ Fig. 11(b)).

Figure 12 shows that the angular acceleration (in green) is much more sensitive to the value of $\alpha$ than the radial acceleration (in red). Therefore, orbital acceleration is a good quantity to measure and to

fit in order to determine the moment of inertia of a pendulum.

### iii. Impact of $r$ on the acquired data

We have shown in section III that the orbital acceleration is very sensitive to the accelerometer sensor position with respect to the pendulum rotation axis. As this position is precisely known, we can perform several measurements with different positions of the accelerometer to improve the accuracy and/or the precision of the derived pendulum parameters. As an example of such analysis, we use the experimental setups shown in Fig. 13.

Figure 14 (a) and (b) display the radial and orbital accelerations measured by the accelerometer after each pendulum shown in Fig. 13 has been displaced by an initial angle of 20.5°.

Fits of the data using Eqs. (11) and (12) with $\kappa$, $\alpha$ and $L$ as free parameters are shown as black curves in Fig. 14. Best fit are obtained with $1/\kappa = 111$ s, $\alpha = 1.162$ and $L = 29.3$ cm for the pendulum of Fig. 13(a) and with $1/\kappa = 100$ s, $\alpha = 1.162$ and $L = 29.2$ cm for the pendulum of Fig. 13(b).

The position of the accelerometer does not affect the calculated values of the moment of inertia of the pendulum and only slightly affects the center of mass position. The new configuration of the wires connecting the accelerometer in Fig. 13(b) changes slightly the $\kappa$ value. Thus, performing a second measurement with a different position of the accelerometer allow us to be more confident in the results retrieved from the first one.

We can also note that the orbital acceleration increases with lower OA values, therefore the orbital acceleration fit precision is better when the accelerometer is in the position of Fig. 13(b), while the precision of the radial acceleration fit is better when the accelerometer is in the position of Fig. 13(a).

## V. Conclusions

We have shown that the pendulum experiment analyzed with an accelerometer sensor leads to a theoretical study richer than the classical one. We have derived theoretical expressions for the radial and orbital acceleration data recorded by an accelerometer and separated the contributions from the pen-
<figure><img src="image_13.png" /><figcaption>Figure 11: Experimental setup used to study compound pendula, with added masses.</figcaption></figure>

dulum composed of a bar on which several masses can be attached to at different positions to acquired data that corresponds to pendulums of different values of $\alpha$ and $L$ (Fig. 11 (a) and (b)).

Figure 12 (a) and (b) displays the radial and orbital accelerations measured by the accelerometer after each pendulum shown in Fig. 11 has been displaced by an initial angle of 20.5°.

We fit the data by using Eqs. (11) and (12) with $\kappa$, $\alpha$ and $L$ as free parameters. Results of the fit are shown as the black curves in the insets of Fig. 12. Best fit are obtained with $1/\kappa = 111$ s, $\alpha = 1.162$ and $L = 0.293$ m for the pendulum of Fig. 11(a) and $1/\kappa = 101$ s, $\alpha = 6.55$ and $L = 0.074$ m for the pendulum of Fig. 11(b).

For the examples displayed in Fig. 12, the evolution of the orbital acceleration from Panels (a) to (b) shows an increase of $\alpha$, which is consistent with the fact that pendulum configuration goes from a configuration close to a simple pendulum ($\alpha \approx 1$ Fig. 11(a)) to a compound pendulum ($\alpha > 1$ Fig. 11(b)).

Figure 12 shows that the angular acceleration (in green) is much more sensitive to the value of $\alpha$ than the radial acceleration (in red). Therefore, orbital acceleration is a good quantity to measure and to

fit in order to determine the moment of inertia of a pendulum.

### iii. Impact of $r$ on the acquired data

We have shown in section III that the orbital acceleration is very sensitive to the accelerometer sensor position with respect to the pendulum rotation axis. As this position is precisely known, we can perform several measurements with different positions of the accelerometer to improve the accuracy and/or the precision of the derived pendulum parameters. As an example of such analysis, we use the experimental setups shown in Fig. 13.

Figure 14 (a) and (b) display the radial and orbital accelerations measured by the accelerometer after each pendulum shown in Fig. 13 has been displaced by an initial angle of 20.5°.

Fits of the data using Eqs. (11) and (12) with $\kappa$, $\alpha$ and $L$ as free parameters are shown as black curves in Fig. 14. Best fit are obtained with $1/\kappa = 111$ s, $\alpha = 1.162$ and $L = 29.3$ cm for the pendulum of Fig. 13(a) and with $1/\kappa = 100$ s, $\alpha = 1.162$ and $L = 29.2$ cm for the pendulum of Fig. 13(b).

The position of the accelerometer does not affect the calculated values of the moment of inertia of the pendulum and only slightly affects the center of mass position. The new configuration of the wires connecting the accelerometer in Fig. 13(b) changes slightly the $\kappa$ value. Thus, performing a second measurement with a different position of the accelerometer allow us to be more confident in the results retrieved from the first one.

We can also note that the orbital acceleration increases with lower OA values, therefore the orbital acceleration fit precision is better when the accelerometer is in the position of Fig. 13(b), while the precision of the radial acceleration fit is better when the accelerometer is in the position of Fig. 13(a).

## V. Conclusions

We have shown that the pendulum experiment analyzed with an accelerometer sensor leads to a theoretical study richer than the classical one. We have derived theoretical expressions for the radial and orbital acceleration data recorded by an accelerometer and separated the contributions from the pen-
<figure><img src="image_13.png" /><figcaption>Figure 12: Radial (red) and orbital (green) accelerations measured by the accelerometer. (a) and (b) correspond to the setups displayed in Fig. 11 (a) and (b), respectively. Insets in (a) and (b) show the temporal evolution of the radial (red) and orbital (green) accelerations and the best fit (black) between 35 s and 45 s. In order to have both curves on the same plot, orbital acceleration is shifted by a constant offset of 9 in the inset of (a) and by 8 in (b).</figcaption></figure>

<figure><img src="image_14.png" /><figcaption>Figure 13: Experimental setup used to study the impact of the acceleration sensor position on the acquired data.</figcaption></figure>

dulum angular motion and the gravitational force on the proof mass.

We have also shown that the orbital acceleration is an interesting data to retrieve the moment of inertia of a pendulum.

The possibility to have different positions of the sensor allows us to perform several measurements with the same pendulum to improve the accuracy and/or the precision of the derived pendulum parameters.

In this paper, we have only focused on the classical pendulum but the device used here could also be applied to more complex systems such as chaotic pendulums.

*Acknowledgements* - This work benefited from a grant “Pédagogie innovante” from IDEX Paris-Saclay. The Institut Villebon - Georges Charpak was awarded the Initiative of Excellence in Innovative Training in March 2012 (IDEFI IVICA: 11-IDFI-002), and is supported by the Paris-Saclay Initiative of Excellence (IDEX Paris-Saclay: 11-IDEX-0003). The authors would like to thank the referee for his/her helpful comments.
<figure><img src="image_14.png" /><figcaption>Figure 12: Radial (red) and orbital (green) accelerations measured by the accelerometer. (a) and (b) correspond to the setups displayed in Fig. 11 (a) and (b), respectively. Insets in (a) and (b) show the temporal evolution of the radial (red) and orbital (green) accelerations and the best fit (black) between 35 s and 45 s. In order to have both curves on the same plot, orbital acceleration is shifted by a constant offset of 9 in the inset of (a) and by 8 in (b).</figcaption></figure>

<figure><img src="image_15.png" /><figcaption>Figure 13: Experimental setup used to study the impact of the acceleration sensor position on the acquired data.</figcaption></figure>

dulum angular motion and the gravitational force on the proof mass.

We have also shown that the orbital acceleration is an interesting data to retrieve the moment of inertia of a pendulum.

The possibility to have different positions of the sensor allows us to perform several measurements with the same pendulum to improve the accuracy and/or the precision of the derived pendulum parameters.

In this paper, we have only focused on the classical pendulum but the device used here could also be applied to more complex systems such as chaotic pendulums.

*Acknowledgements* - This work benefited from a grant “Pédagogie innovante” from IDEX Paris-Saclay. The Institut Villebon - Georges Charpak was awarded the Initiative of Excellence in Innovative Training in March 2012 (IDEFI IVICA: 11-IDFI-002), and is supported by the Paris-Saclay Initiative of Excellence (IDEX Paris-Saclay: 11-IDEX-0003). The authors would like to thank the referee for his/her helpful comments.
<figure><img src="image_15.png" /><figcaption>Figure 14: Radial (red) and orbital (green) accelerations measured by the accelerometer. (a) and (b) correspond to the setups displayed in Fig. 13(a) and (b), respectively. Insets in (a) and (b) show the temporal evolution of the radial (red), and orbital (green) accelerations and the best fit (black) between 35 s and 45 s. In order to have both curves on the same plot, orbital acceleration is shifted by a constant offset of 9 in the insets of (a) and (b).</figcaption></figure>

[1] http://www2.vernier.com/sample_labs/PHYS-AM-18-physical_pendulum.pdf

[2] https://physics fullerton edu/files/Labs/225/Physical Pendulum V4 1 1.pdf

[3] J C Fernandes, P J Sebastião, L N Gonçalves, A Ferraz, Study of large-angle anharmonic oscillations of a physical pendulum using an acceleration sensor, Eur. J. Phys. 38, 2017.

[4] https://www arduino cc

[5] Š Kubínová, J Šlégr, *Physics demonstrations with the Arduino board*, Phys. Educ. 50, 472 (2015).

[6] C Galeriu, S Edwards, G Esper, *An Arduino investigation of simple harmonic motion*, Phys. Teach. 52, 157 (2014).

[7] K Atkin, *Construction of a simple low-cost teslameter and its use with Arduino and Maker-Plot software*, Phys. Educ. 51, 024001 (2016).

[8] C Petry et al., *Project teaching beyond Physics: Integrating Arduino to the laboratory*, In Technologies Applied to Electronics Teaching (TAEE), pp. 1-6, IEEE (2016).

[9] R Henaff et al., *A study of kinetic friction: The Timoshenko oscillator*, Am. J. Phys. 86, 174 (2018).

[10] F Bouquet, J Bobroff, M Fuchs-Gallezot, L Maurines, *Project-based physics labs using low-cost open-source hardware*, Am. J. Phys. 85, 216 (2017).
<figure><img src="image_16.png" /><figcaption>Figure 14: Radial (red) and orbital (green) accelerations measured by the accelerometer. (a) and (b) correspond to the setups displayed in Fig. 13(a) and (b), respectively. Insets in (a) and (b) show the temporal evolution of the radial (red), and orbital (green) accelerations and the best fit (black) between 35 s and 45 s. In order to have both curves on the same plot, orbital acceleration is shifted by a constant offset of 9 in the insets of (a) and (b).</figcaption></figure>

[1] http://www2.vernier.com/sample_labs/PHYS-AM-18-physical_pendulum.pdf

[2] https://physics.fullerton.edu/files/Labs/225/Physical_Pendulum_V4_1(1).pdf

[3] J C Fernandes, P J Sebastião, L N Gonçalves, A Ferraz, *Study of large-angle anharmonic oscillations of a physical pendulum using an acceleration sensor*, Eur. J. Phys. **38**, 2017.

[4] https://www.Arduino.cc

[5] Š Kubínová, J Šlégr, *Physics demonstrations with the Arduino board*, Phys. Educ. **50**, 472 (2015).

[6] C Galeriu, S Edwards, G Esper, *An Arduino investigation of simple harmonic motion*, Phys. Teach. **52**, 157 (2014).

[7] K Atkin, *Construction of a simple low-cost teslameter and its use with Arduino and Maker-Plot software*, Phys. Educ. **51**, 024001 (2016).

[8] C Petry et al., *Project teaching beyond Physics: Integrating Arduino to the laboratory*, In Technologies Applied to Electronics Teaching (TAEE), pp. 1-6, IEEE (2016).

[9] R Henaff et al., *A study of kinetic friction: The Timoshenko oscillator*, Am. J. Phys. **86**, 174 (2018).

[10] F Bouquet, J Bobroff, M Fuchs-Gallezot, L Maurines, *Project-based physics labs using low-cost open-source hardware*, Am. J. Phys. **85**, 216 (2017).