Computational methods for the investigation of ski boots ergonomics

Ski boots are known to cause vasoconstriction in the wearer’s lower limbs and, thus, cause a “cold leg” phenomenon. To address this problem, this work provides a computational framework for analysing interactions between the ski boot and the lower limb. The geometry of the lower limb was derived from magnetic resonance imaging and computed tomography techniques and anthropometric data. The geometry of the ski boot shell was obtained by means of three-dimensional computer aided design models from a manufacturer. Concerning the ski boot liner, laser scanning techniques were implemented to capture the geometry of each layer. The mechanical models of the ski boot and the lower limb were identified and validated by means of coupled experimental investigations and computational analyses. The computational models were exploited to simulate the buckling process and to investigate interaction phenomena between the boot and the lower limb. Similarly, experimental activities were performed to further analyse the buckling phenomena. The obtained computational and experimental results were compared regarding both interaction pressure and displacements between the buckle and the corresponding buckle hooks. These comparisons provided reasonable agreement (mean value of discrepancy between the model and mean experimental results in the tibial region: 20%), underlining the model’s capability to correctly interpret results from experimental measurements. Results identified the critical areas of the leg, such as the tibial region, the calcaneal region of the foot and the anterior sole, which may suffer the most due to the hydrostatic pressure and compressive strain exerted on them. The results highlight that computational methods allow investigation of the interaction phenomena between the lower leg and ski boot, potentially providing an effective framework for a more comfortable and ergonomic design of ski boots.


Introduction
Alpine skiing is one of the most popular winter sports [1].From a biomechanical viewpoint, ski boots are the most important piece of ski equipment [2].Many factors determine ski boot comfort and ergonomics, particularly liner materials and structures.Previous studies demonstrated the influence of ski boot structure on heat and moisture transfer phenomena [3,4].Furthermore, ski boots are known to generate vasoconstrictive phenomena on the wearer's lower limbs, thus causing a "cold leg" sensation [5][6][7] and compromising ergonomics and comfort.
Ski boot buckling determines strain and stress fields within lower limb soft tissues, which may cause the vasoconstriction phenomena.Tools that enable the analysis of such stress and strain fields can be highly relevant.Experimental activities may provide preliminary information about the interaction phenomena between the ski boot and lower limb.For example, pressure sensors can be applied to measure the interaction pressure between the lower limb and ski boot at specific points [8,9].On the other hand, computational methods are extensively used for structural investigation of different parts of the ski boot leading to more refined and detailed results in a more efficient way.Moreover, computational approach enables a wider and deeper investigation, providing the continuous distribution of interaction pressure and effective data about stress and strain within biological tissues [10].
The aim of this work was to develop a finite element model of the interaction between the ski boot and the lower limb, to asses and to analyse pressures, strains and stresses characterising this interaction.Mechanical stimulation of soft tissue is mainly influenced by the following two factors: (a) the ski boot's geometrical configuration and materials, and (b) the wearers' lower limb anatomy and morphometry.A reliable computational framework must account for such parameters, with the goal to improve ski boot's comfort and user's ergonomics.Experimental tests were developed to validate the computational model.

Materials and methods
The in silico investigation of ski boot ergonomics and comfort has been developed by coupling experimental and computational activities.A general overview of the activities is reported in Fig. 1.The preliminary actions pertained to the development of lower limb and ski boot computational models.The coupling of such models entailed the computational analysis of ski boot buckling.Contemporarily, experimental measurements of the interface pressure between the lower limb and ski boot from buckling have been developed.The comparison of experimental and computational results made it possible to assess the dependability of the models.Finally, post-processing of computational results provided information about the mechanical stimulation of lower limb biological tissues, which allow evaluating ski boot comfort and ergonomics.

Finite element model of the lower limb
To obtain the three-dimensional computer aided design (3D CAD) model of the lower limb, computed tomography (CT) (NewTom/NTVGiMK4, 3 mA, 110 kV, 0.3 mm cubic voxels) and magnetic resonance imaging (MRI) (3 T Siemens Magnetom Trio, fat-suppressed 3D dual echo steady state (DESS) sequence with 0.7 mm isotropic resolution, size matrix 320 × 576 × 104, TE/TR = 5.5/13 ms, flip angle Fig. 1 Flowchart of the activities: computational modelling of lower limb (blue boxes; Sect.2.1) and ski boot (yellow boxes; Sect.2.2); models coupling and computational analysis of ski boot buckling (green boxes; Sect.2.3); experimental investigation of interface pres-sure (pink box; Sect.2.4); model validation and exploitation for the evaluation of ski boot comfort and ergonomics (grey boxes; Sects.3,4) 25 degrees, TA = 5 min) images of a Caucasian male (age: 30 years; height: 175 cm; weight: 77 kg; foot EU size: 42) were acquired in a no-load set-up.The specific person was informed about the purpose and procedures of the study, and signed an informed consent for processing the resulting models from anonymised biomedical images.In details, biomedical images were processed by means of segmentation techniques to obtain stereolithography (STL) triangulated models.With the process of segmentation, models obtained may be characterised by irregular triangulated surfaces, with gaps and scattering due to the stratification of the biomedical images.For this purpose, the first step involved the regularisation of the surface (Simpleware, Synopsys Corporate, Mountain View, California, USA), starting from the operation of smoothing that solves the minor defects, while maintaining the overall shape.On the other side, the surfaces that presented relevant discontinuity were selected, deleted and replaced with a regular shape that continuously followed the trend of the surrounding region.Subsequently, triangulated surfaces were interpolated by means of free-form surfaces (UGS NX, Siemens PLM Software, Plano, Texas, USA).The process led to free-form 3D CAD models of the main biological components of the lower limb, such as bones, adipose tissue, calf muscle, Achilles' tendon and skin.
The availability of average morphometric data [11][12][13][14][15][16][17] and the use of morphing techniques [18,19] made it possible to obtain an accurate and reliable geometry of a typical adult male (lower leg length = 530 mm [16]; calf circumference = 395 mm [13]; long heel length = 126 mm [14]; foot length = 260 mm [17]; ball width = 102 mm [17]) (Fig. 2a).Furthermore, morphometric data permitted to add biological structures that are not easily and accurately identifiable through MRI or CT images, such as cartilages [20] and ligaments [21].The main cartilages and ligaments of the ankle joints were 3D modelled, while truss elements were used to provide simplified links among the other bony elements, with particular regard to the metatarsal and phalangeal regions [22].The model was defined accounting for the possibility to morph the geometry on the basis of subject specific characteristics, with specific regard to the anatomical properties reported in Fig. 2a.On the other side, the specific dimension of some biological structures, such as ligaments and cartilages, must be defined based on anatomical considerations and data.The initial angle between the floor and tibia was defined by considering a relaxed position that is usually assumed during the ski boot wearing and buckling, such as 77° that was experimentally measured, as fully reported in Sect.2.4.
The resulting 3D CAD model was discretised in the pre-processing environment Abaqus CAE 2018 (Dassault Systèmes, France) using four-node linear tetrahedral finite elements.With regard to the different regions of the model, different mesh discretization levels and techniques were assumed.Specific sensitivity analyses have been performed Fig. 2 Virtual solid and finite element modelling of the lower limb: a anthropometric reference data (lower leg length = 530 mm [16]; calf circumference = 395 mm [13]; long heel length = 126 mm [14]; foot length = 260 mm [17]; ball width = 102 mm [17]); b complete geometry (blue: skin and subcutaneous loose tissue; red: bony elements; light brown: calf muscle), c internal view of bony elements (red), calf muscle (light brown) and Achilles tendon (azure), d details of some of the foot structures (azure: Achilles tendon; blue: ligaments; green: cartilages; red: plantar fascia; violet: plantar tissues; yellow: skin), details of e cartilages (dark grey) and f ligaments (dark grey) of the tibiotalar joint aiming to define the optimal element size, by considering both the quality of the computational results and the computational effort of the model.In this sense, the average size of the elements varied to discretize in the best way the geometry and the level of details of each model subcomponent (ranging between 0.05 and 1.00 mm) (Fig. 2b-f).
Concerning the mechanical behaviour of the biological tissues, the soft tissues were differentiated from the hard ones, giving the latter a rigid formulation [23].This choice did not alter the model's dependability because the forces generated due to the lower limb and the ski boot interaction do not strain the bone structure substantially.Moreover, the rigid body assumptions reduced the computational effort of the model, permitting to analyse the strains of soft tissues, which are mainly involved in vasoconstriction and whose magnitude is typically greater than that of bony elements.
Concerning the soft tissues, their mechanical response must account for large strains and displacements.Furthermore, such tissues are known to exhibit non-linear elastic behaviour as well as time dependent effects.For these reasons, these tissues were assigned hyperelastic or viscohyperelastic behaviour using both isotropic and anisotropic formulations depending on the specific tissue [24]: where S is the second Piola-Kirchhoff stress tensor, C is the right Cauchy Green strain tensor and q i are viscous variables.W 0 is a strain energy function that specifies the instantaneous hyperelastic behaviour of the tissue, and can assume isotropic or anisotropic conformations.W 0m is the isotropic contribution, while W 0f provides anisotropic features from ( 1) fibres families that are initially oriented along directions a 0 j [25]: Standard differential equations define the evolution of viscous variables, which specify relaxed stress components, depending on viscous parameters, such as relative stiffens γ i and relaxation timeτ i [26]: With regard to the different tissues, the specific model formulation and the identification of the associated constitutive parameters required relevant efforts.They mainly concerned to the specific mathematical definition of the strain energy function of Eq. ( 2) on the basis of tissue histology and preliminary data about tissue mechanics, the identification of the associated constitutive parameters by means of the inverse analysis of mechanical tests [10], and the computational implementation of the constitutive formulation by means of "ad hoc" computational subroutines.Specific references to constitutive analyses and experimental investigations are reported in Table 1.
Finally, frictionless contact conditions and tie constraints defined the interactions among soft and bony elements, depending on the specific situation.Specifically, tie constraint ensured continuity among skin and adipose tissues, between bone and cartilages or ligaments/tendons insertions, and between skin, calf muscle and Achilles' tendon.The final finite element discretization of the lower limb structures led to about 3 million four-node tetrahedral elements, 700,000 nodes and 2.5 million variables.

Muscular tissue of the calf Anisotropic hyperelastic -passive behaviour Tensile tests along different loading directions [52] Structural compression tests [53] Skin
Anisotropic hyperelastic [54] Tensile tests along different loading directions [55] Compression tests [56,57] Before exploiting the computational model of the lower limb for predictive investigations, its reliability must be assessed.The evaluation of such dependability usually requires analysing further experimentations, which must have been performed at the structure level.The comparison between model and experimental results finally claims the capability of the model to actually predict the mechanical behaviour of the biological structure.Because of the complexity of the overall model and the lack of experimentations on the overall lower limb, specific experimental activities on different anatomical structures were investigated.To validate the mechanical model of the calcaneal region, computational analyses simulated indentation tests, with regard to both in vitro [59] and in vivo [61,62] conditions.Fig 3a, b report the experimental and the computational set-up, respectively, while in Fig. 3c the comparison between model and experimental results is provided.The dependability of the ankle joint model required simulating typical clinical tests, like the dorsiflexion and the plantarflexion tests [39], and the anterior drawer and the talar tilt tests [63].Again, Fig. 3d-f provides an overview of the experimental and computational experimentations, while experimental and model results are compared in Fig. 3g, h.Aiming to further assess the effectiveness of the foot model, the different steps of the gait cycle were analysed [22].

Finite element model of the ski boot
The development of the 3D CAD model of the ski boot involved processing of aesthetic virtual solid models (Fig. 4a) that were provided by a boot manufacturer, with particular regard to shell components.The manufacturer provided an STL model that they mainly use for marketing purposes, which was interpolated by means of free-form surfaces.Laser scanner techniques were also exploited to develop a model of the liner and the cuff.Even though the procedure was performed on a specific product, it can be applied to model other ski boots.The size of the ski boot was 26 cm (Mondopoint sizing system), the internal width in the wider section (LAST) was 100 mm and the nominal flex index was 130, as stated by the manufacturer.
As previously reported, the provided geometric models of the shells were rendered using STL techniques and interpolated by means of free-form surfaces (Fig. 4b).The  3D models were discretised using four-node tetrahedral elements with an average size that ranged between 0.4 and 2.0 mm of mesh seed, depending on the region meshed (Fig. 4c).Furthermore, the use of tie constraints, multi-point constraints and the introduction of frictional contact conditions (using 0.2 static friction coefficient [64][65][66]) allowed for the correct description of the interactions among shell components (Fig. 4d,e).
Concerning the liner and the cuff, the model must envisage the different building layers and materials to correctly interpret their different mechanical behaviour.In the first place, it was decided to acquire geometry and inner structure by means of a CT scan but, due to the similar interaction of the different materials with X-ray, the approach failed in discriminating the different layers.For this reason, laser scanning techniques were used on a dummy lower limb on which each layer was placed sequentially.Manufacturing of the dummy model required 3D printing (Materialise NV, Leuven, Belgium; PC-ABS) of the previously developed 3D CAD model of the lower limb, accounting for the external surface, such as the skin surface, of the model only.After positioning a layer, the geometry was acquired.Subsequently, Boolean operations led to preliminary STL models of the layers (Fig. 4f, h), whose subsequent processing provided 3D CAD (Fig. 4g, i) and finite element (Fig. 4j) models.Tie constraints ensured the continuity among finite element models of the different layers.Finite element discretization was performed by means of four-node tetrahedral elements, whose size ranged between 0.2 and 1.0 mm.
The mechanical behaviour of the rigid thermoplastic components of the shells and of the liner (thermoplastic polyurethane Covestro Desmopan 3860 and 3856, respectively) was identified using technical data provided by the manufacturer.These components followed the assumption of linear elastic isotropic behaviour (tensile modulus 500 MPa for Desmopan 3860 and 310 MPa for Desmopan 3856 at 20 °C temperature, while 0.45 Poisson ratio was assumed for both materials), and the properties agreed with those reported elsewhere [67].Mechanical properties of rigid thermoplastic components were assumed at 20 °C temperature, because the experimentations on the liner soft foam and fabric materials were performed at room temperature, as described in the following.
Concerning the foam and the coupled foam/fabric materials of the liner and the cuff, the materials were tested distinctively by means of quasi-static tensile (EN ISO 1798:2008; dog bone specimen shape; 40 and 10 mm length and width, respectively, of the central region; 10 mm minimum thickness; 500 mm/min loading rate) and quasi-static compression (EN ISO 3386-1:1997/A1 2010; cylindrical or parallelepiped shape; 2500 mm 2 minimum specimen section area; 10 mm minimum specimen thickness; 100 mm/min loading rate) tests.With specific regard to thermoformed elements, experimentations were performed on specimens that were cut out from the manufactured components.The analysis of the experimental results suggested the use of a first-order Ogden hyperelastic formulation [25].The inverse analyses of such compression and tensile tests led to the materials' parameters (Fig. 4k-m).More details on experimental activities on liner materials are reported in Online Resource 1.
The final finite element discretization of the ski boot components led to about 2 million elements (four-node tetrahedral brick elements and three-node triangular shell elements), 600,000 nodes and 2.3 million variables.

Computational analysis of ski boot buckling process
CAD virtual translations and rotations made it possible to position the ski boot liner and the cuff within the shells (Fig. 5a).Friction contact conditions defined the interaction among the liner, cuff and shells (0.3 static friction coefficient for interaction between shells and liner/cuff, 0.4 static friction coefficient for self-interaction between liner and cuff surfaces).Again, CAD virtual position procedures, together with minimal Boolean operations, allowed for positioning of the lower limb model within the ski boot (Fig. 5b).Considering the presence of a sock, a friction contact condition specified the interactions between lower limb and ski boot assuming a 0.4 static friction coefficient.The specific static friction coefficients were assumed depending on the specific materials that build up the analysed ski boot [64][65][66].Nevertheless, the computational approach makes it possible to analyse the influence of different friction properties on the interaction phenomena between the ski boot and lower limb.Indeed, many different factors (materials, surface characteristics, temperature, etc.) may influence friction conditions.The final model of the system composed of the ski boot and lower limb resulted in about 5 million elements, 1.3 million nodes and more than 5 million variables.
To simulate the buckling of the ski boot, the buckles and buckle hooks were simplified to include only their centre points, connected to the shell nodes by means of multi-point constraints.The connection of centre points by means of connector elements (Fig. 5c) made it possible to describe the buckling process by means of reducing the length of the connectors (Fig. 5d) [68].The toe and heel lugs were fully constrained, as for the top region of the tibial head.The complexity of interactions phenomena among model components and the complexity of materials behaviour suggested assuming the explicit dynamic approach [69] of the general-purpose finite element code Abaqus Explicit 2018 (Dassault Systèmes, France).The mechanical behaviour of biological tissues was specified by the previously described visco-hyperelastic formulations, which were implemented by means of specifically developed VUMAT Fortran routines [68].
With regard to the step duration, the buckling process requires about 30 s to develop.On the other hand, the specific configuration of the model led to small stable time increments within the explicit algorithm [69], as low as 10 -7 s, causing a high computational effort in solving the quasi-static process.To reduce the computational effort, mass scaling techniques were adopted.Monitoring nodal accelerations, overall model kinetic and residual energies made it possible to ensure the reliability of the analysis.Computational calculations were performed by means of a high-performance computing server equipped with two Intel Xeon E7 8890 v4 (Intel Corporation, Santa Clara, CA), 256 GB RAM and solid-state hard disk drives.The analysis of the buckling process engaged 40 cores (80 threads) and lasted about 160 h.An animation of the buckling process from the model is included in Online Resource 2.

Experimental analysis of ski boot buckling process
Aiming at acquiring the necessary data regarding the pressure exerted on the lower limb during the buckling process, a set-up consisting of a portable acquisition system and an inflatable cuff was specifically developed.The inflatable cuff (Kinder Armumfang, Riester, Germany) was placed on the tibia, with the distal part at the level of the malleoli.In this way, mean average pressures were taken over a large anatomical area (~ 28 cm 2 ) to have an indication of the general compression generated by the progressive buckling process.The inflatable cuff was connected to a silicon pressure sensor (MPX5050DP pressure sensor, Freescale Semiconductor, USA) utilising temperature compensation and operating with a pressure range of 50 kPa.The signal was acquired by an Arduino Nano microcontroller (Arduino.cc)using the internal AD 10 bit converter whose minimal measurable pressure change was 0.049 kPa.The system was calibrated by comparing the measures with those from a manometer (LEX 1, Italmanometri, Reggio Emilia, Italy; pressure measurement range: -1 ÷ 2 bar, resolution: 0.1 mbar).All the experiments were performed at room temperature.Before experimentally measuring the interaction pressure, the angle between the tibia and floor was evaluated by means of a goniometer (RVS Goniometer, Baseline Evaluation Instruments, New York, USA), providing results ranging between 74° and 79°.Angle measurements were performed by assuming a relax position of the tester.
To perform the experimental analysis of the ski boot's buckling, a coding procedure was assumed.The buckles were numbered starting from the foot tip toward the top of the cuff (Fig. 6a) and each hook was numbered according to their proximity to the buckle hand (Fig. 6b).Thus, it was possible to codify each possible buckling configuration with four numbers representing the buckle (by the position of the number) and the buckle hook chosen (with the value of the number).This coding system allowed for the acquisition of data regarding the different buckling configurations, in terms of both pressure and displacements between the buckle and the corresponding buckle hooks (Fig. 6c).
Buckles displacements were measured with reference to the initial configuration of the unbuckled ski boot.In this sense, the relative displacement exactly defines the shortening of connector elements (Fig. 4c) that was applied to computationally simulate the buckling process.In detail, buckling was done by lacing up the first and the second buckles on the second hook, as configuration 2200.Afterward, the process was carried on the third and fourth buckles, leading to configuration 2222.Following the previous step, the first and the second buckles were adjusted to reach configuration 3422, and finally the third and the fourth buckles were tightened to the fourth hook, reaching configuration 3444 (Fig. 6c).
Experimental information was necessary to perform computational analyses whose results can be compared with experimental data.Aiming to provide general data, measurements were performed on six participants, whose anthropometric parameters, such as foot length, ball width, calf circumference, long heel length and lower leg length, were similar (± 10%) to the developed computational model of the lower limb (Fig. 2a).The study was approved by the Ethical Committee of the University of Verona (2019-UNVR-CLE-0434163, Department of Neurosciences, Biomedicine and Movement, 28th November 2019) and performed in accordance with the Declaration of Helsinki.All the participants were informed about the test purpose and procedures of the study, and signed an informed consent for the scientific implementation, application and processing of experimental results and data.

Results
The computational analyses of the buckling process provided data about the forces that were required to shorten the connector elements (Fig. 6c), as the buckling forces, which are summarised in Table 2.
Concerning the interaction pressure between the lower limb and ski boot at the tibial region, experimental data and computational results were compared (Fig. 7).
The distributions of the interaction pressure between the lower limb and ski boot and compressive strain within the lower limb tissues are reported in Fig. 8 for the different buckling configurations.The maximum values of interaction pressure and compressive strain were reached for configurations 3422 and 3444 in correspondence of the soft tissues of the metatarsal and calcaneal regions and the anterior lower leg region.Configuration 3422 showed an interaction pressure of 25 kPa and a compressive strain of 18%, while configuration 3444 showed a high increment of interaction pressures (up to 50 kPa) and compressive strains (over 36%).
Figure 9 provides the distribution of hydrostatic pressure in the most severe buckling configuration (3444), together with a representation of blood pressure within the vascular system [70].Analysing the hydrostatic pressure exerted on the lower limb during configuration 3444, it is possible to observe that the local hydrostatic pressure reached values exceeding 75 mmHg, in correspondence of the calcaneal and anterior regions (Fig. 9a, b).To better evaluate the influence of hydrostatic pressure on soft tissues functionality, typical values of blood pressure in the human organism are reported (Fig. 9c).

Discussion
The reported activities highlight the potentiality and the effectiveness of computational biomechanics for the investigation of the interaction phenomena between a ski boot and lower limb.Coupled experimental and computational activities demonstrated the dependability of the lower limb model Fig. 8 Results from the computational analysis of ski boot buckling at different configurations: a contours of interaction pressure between lower limb and ski boot on the skin; b contours of compressive strain on a middle section of the lower limb Fig. 9 Investigation of hydrostatic pressure: a distribution of the hydrostatic pressure within the lower limb and b details on a section; c blood pressure within the vascular system [70] (Fig. 3) and methodologies (Fig. 7).With specific regard to the ski boot, the computational investigation of the buckling process provided values of buckling forces (Table 2) that were similar to data from experimental activities reported elsewhere [71].Computational methods can also provide more information and data than experimental activities, like buckling forces, interaction pressure between the lower limb and ski boot (Fig. 7a), and stress and strain fields within biological tissues (Fig. 8b).
Mechanical stimuli determine many different biological effects.For example, high compressive strain values can induce tissue damage.Prolonged compression of soft tissues adjacent to bony prominences can lead to tissues degeneration, resulting in a condition called pressure-related deep tissue injury, which can develop into deep pressure ulcers and wounds.Previous studies have demonstrated that compression strains over 30% can induce significant tissue damage [72].High hydrostatic pressure values (Fig. 9) highlight possible occlusion of blood vessels, providing a source of vasoconstriction and explaining the "cold leg" sensation experienced by some ski boot users [5,6].Such mechanobiological factors can affect the ski boot performance, comfort and ergonomics [7].In this sense, the proposed approach defines a novel methodology for the design and the reliability assessment of ski boot geometry and materials.
Great efforts pertained to the definition of a reliable and accurate model of the lower limb, by considering the geometrical modelling of the different biological structures and the mechanical characterisation of the different biological tissues.With particular regard to parameters reported in Fig. 2a, the model is suitable for morphing procedures, aiming at interpreting the anatomical configuration of different skiers.Furthermore, tissue parameters can be modified to interpret specific situations, with particular regard to pathologic ones.On the other side, ski boot modelling was performed by considering a specific product.The action made it possible to illustrate the potentialities of the computational approach.The same methodology can be extended to analyse the interaction phenomena between the lower limb and different ski boots.
The presented results are promising, but it is worth noticing some critical points in the approach presented here.To reach an appreciable level of reliability, the model was discretised with many finite elements, thus leading to computationally heavy models and simulations.The process for the acquisition of the liner's geometry, due to the complex geometries and number of layers, also remains a complex and long process.A further limitation pertains to the Velcro strap that usually characterises ski boots, and was not included within the model.Nevertheless, the Velcro strap mainly acts on the upper region of the cuff and it was not tied during experimental measurements of interaction pressure.
The proposed model has been exploited to analyse a static condition, such as the interaction phenomena that occur between the lower limb and ski boot in a steady state situation.The activities aimed at analysing the influence of incremental ski boot buckling on mechanical stimulation of biological tissues.Both the proposed experimental and computational methods can be further exploited to analyse dynamic conditions.Indeed, plantar-and dorsiflexion that occur during skiing affect interactions between the biological structures and the ski boot.Further activities are under development to analyse these topics.Furthermore, aiming to a more accurate reliability assessment of interaction pressure measurements, pressure sensors could be located also in other positions, such as the metatarsal and the calcaneal region of the foot, and the anterior sole, as reported in [73].
Another limitation of the proposed investigations pertains to temperature.Ski boots ergonomics was here analysed at room temperature, but they are usually worn outdoors in sub-zero temperature (Celsius) conditions.Unfortunately, the experimental facilities available made it possible to characterise the mechanical behaviour of liner materials at room temperature only.Furthermore, the experimental evaluation of interaction pressure during the buckling process was performed at room temperature, and such data were mandatory for the reliability assessment of the overall computational model.Anyway, the paper aims at providing a novel computational approach and framework for the investigation of ski boot functionality and interaction phenomena with lower limb tissues.Given materials parameters at a specific temperature, the model should allow for analyses at such temperature, although its accuracy would need to be reassessed at this temperature.

Conclusions
The developed activities highlight the potentialities of computational biomechanics for both the investigation of sport devices reliability and the design of better equipment.Computational methods provide data and information that barely can be achieved by means of experimental activities.Furthermore, computational models, once they have been developed, can be exploited to analyse a broad scenario of situations, considering many conformations of the device, many anatomical configurations of the user and many loading conditions.With specific regard to the here developed investigations, aside from being a reliable base for the modelling of the lower limb and ski boot interaction, the reported approach and models are also suitable to be coupled with morphing techniques based on acquiring the skier's morphometric data by means of laser scanning.This non-invasive technique should allow for the acquisition of the most relevant data for the designing of a tailored ski boot, or for the evaluation of ski boot suitability for the specific user.

Fig. 3
Fig. 3 Examples of coupled experimental and computational activities for lower limb model validation.Indentation of calcaneal region: a experimental set-up by Erdemir et al. (2009) [58] b contour of compressive strain from computational analysis; c comparison of computational results (continuous lines) and median experimental data (empty circles): computational simulations have been performed assuming the principal set of constitutive parameters (middle grey

Fig. 4
Fig. 4 3D CAD, finite element modelling of ski boot shell and liner components and constitutive analysis: a aesthetic model of the ski boot; b 3D CAD free-form models of shell components; c finite element models of shell components; d coupling of lower and upper shell components: friction contact (0.2 friction coefficient) characterises the typical interaction condition; tie constraint ensured the coupling by screws or nails in the bottom region (blue rectangle); the red rectangle highlights the region of sub-plot e, where multipoint constraints limit the relative motion of upper and lower shells to a rotation around a specific local axis (violet arrow); f STL model

Fig. 5
Fig. 5 Coupling of computational models and analysis of the buckling process: a coupling of 3D CAD models of shells, liner and cuff; b detail of the finite element model of the coupling of ski boot and

Fig. 6 2 Fig. 7
Fig. 6 Experimental coding procedure of ski boot buckling: a definition of buckles numbers; b definition of buckle hooks numbers; c relative displacement between buckle and buckle hooks, which have been measured with reference to the initial unbuckled configuration This article is a part of Topical Collection in Sports Engineering on Winter Sports, edited by Dr. Aimee Mears, Dr. David Pearsall, Dr. Irving Scher and Dr. Carolyn Steele.