After graduating, I felt that I still needed to learn a lot, so I decided to start my PhD in biomechanics at the same faculty. In addition, I am currently studying at the University School of Physical Education in Wrocław, at athlete trainer speciality. I also train sprinting.
I am a person who is into new challenges and is really happy to cooperate and learn from others. That is why during my studies I participated in a few foreign mathematics programs. One of those was 33rd edition of the ECMI Modelling Week which took place in Darmstadt, Germany. During this event, we were working in small, international groups (under a supervision of the instructor) on projects based on the real-life problems. Our task was to predict water consumption by households, based on earlier meter readings. We received real data collected recurrently from big villas with gardens in Portugal. To solve this problems the knowledge gained during mathematical studies on such courses as time series was needed. We also learned a lot from our supervisor and from ourselves.
Moreover, at the turn of August and September this year, I participated in Applied Mathematics Summer School organized by King Abdullah University of Science and Technology (KAUST) in the Kingdom of Saudi Arabia. Here, together with other students, we proposed a price-formation model, with a population consisting of a finite number of agents storing and trading a commodity. I believe that participating in such projects and seeing how one can use (sometimes quite abstract) university knowledge to overcome industrial problems can be highly motivating for the future studies.
My master’s thesis, written under the supervision of Professor Łukasz Płociniczak, concerned modelling human movement. We consider the classic spring-mass model of running which is built upon an inverted elastic pendulum (cf. [1-3]). In a natural way, during the contact phase, we get the governing system of two nonlinear ordinary differential equations:
![\left[\begin{array}{c} L^2\ddot{\theta} + 2L\dot{L}\dot{\theta}\\ \ddot{L} \end{array}\right] = \left[\begin{array}{c} L\sin\theta\\ L\dot{\theta}^2+K(1-L)-\cos\theta \end{array}\right],](https://s0.wp.com/latex.php?latex=%5Cleft%5B%5Cbegin%7Barray%7D%7Bc%7D+L%5E2%5Cddot%7B%5Ctheta%7D+%2B+2L%5Cdot%7BL%7D%5Cdot%7B%5Ctheta%7D%5C%5C+%5Cddot%7BL%7D+%5Cend%7Barray%7D%5Cright%5D+%3D+%5Cleft%5B%5Cbegin%7Barray%7D%7Bc%7D+L%5Csin%5Ctheta%5C%5C+L%5Cdot%7B%5Ctheta%7D%5E2%2BK%281-L%29-%5Ccos%5Ctheta+%5Cend%7Barray%7D%5Cright%5D%2C+&bg=ffffff&fg=000&s=0&c=20201002)
where
Figure 1. Schematic of the main model
The initial conditions for the polar coordinate system have the form:
with
where
Problem. Let
This means that the spring stiffness
Next, we derive an approximate solution assuming a small angle of attack and a small spring compression:
![\left[\begin{array}{c} \widetilde\theta(T)\\ \widetilde{L}(T) \end{array}\right] = \left[\begin{array}{c} \theta_d\sinh{T} - \alpha\cosh{T}\\ 1-{\omega}^{-1}L_d\sin{\omega T}- {\omega^{-2}}(1 - \cos{\omega T}) \end{array}\right],](https://s0.wp.com/latex.php?latex=%5Cleft%5B%5Cbegin%7Barray%7D%7Bc%7D+%5Cwidetilde%5Ctheta%28T%29%5C%5C+%5Cwidetilde%7BL%7D%28T%29+%5Cend%7Barray%7D%5Cright%5D+%3D+%5Cleft%5B%5Cbegin%7Barray%7D%7Bc%7D+%5Ctheta_d%5Csinh%7BT%7D+-+%5Calpha%5Ccosh%7BT%7D%5C%5C+1-%7B%5Comega%7D%5E%7B-1%7DL_d%5Csin%7B%5Comega+T%7D-+%7B%5Comega%5E%7B-2%7D%7D%281+-+%5Ccos%7B%5Comega+T%7D%29+%5Cend%7Barray%7D%5Cright%5D%2C+&bg=ffffff&fg=000&s=0&c=20201002)
where
where
Figure 2. Comparison between the numerically calculated value of
Numerical calculations showed very good compatibility of the approximation with the exact solution (see Figure 2). In addition, measurements of the running parameters of professional athletes allowed to calibrate the model.
The studied problem turned out to be interesting enough that together with Professor Płociniczak we wrote a paper on the analysis of equations related to the system with the use of Poincaré – Lindstedt series, see [4]. A perturbative expansion with respect to the large spring stiffness is relevant for the solution of the Problem. Moreover, we used the obtained asymptotic estimates to prove that there exists a unique solution to the aforementioned boundary value problem and provided the same approximation to the sought stiffness.
[1] Blickhan, R. (1989). The spring-mass model for running and hopping. Journal of Biomechanics, 22, 1217-1227.
[2] Geyer, H., Seyfarth, A. & Blickhan, R. (2005). Spring-mass running: simple approximate solution and application to gait stability. Journal of Theoretical Biology, 232, 315-328.
[3] McMahon, T. A. & Cheng, G. C. (1990). The mechanics of running: How does stiffness couple with speed? Journal of Biomechanics, 23, 65-78.
[4] Płociniczak, Ł. & Wróblewska, Z. (2019). Solution and asymptotic analysis of a boundary value problem in the spring-mass model of running. Online available at https://arxiv.org/abs/1903.02247(pdf).
By Zofia Wróblewska, student, Faculty of Pure and Applied Mathematics, Wrocław University of Science and Technology