Safety experiments for small robots investigating the potential of soft materials in mitigating the harm to the head due to impacts

There is a growing interest in social robots to be considered in the therapy of children with autism due to their effectiveness in improving the outcomes. However, children on the spectrum exhibit challenging behaviors that need to be considered when designing robots for them. A child could involuntarily throw a small social robot during meltdown and that could hit another person's head and cause harm (e.g. concussion). In this paper, the application of soft materials is investigated for its potential in attenuating head's linear acceleration upon impact. The thickness and storage modulus of three different soft materials were considered as the control factors while the noise factor was the impact velocity. The design of experiments was based on Taguchi method. A total of 27 experiments were conducted on a developed dummy head setup that reports the linear acceleration of the head. ANOVA tests were performed to analyze the data. The findings showed that the control factors are not statistically significant in attenuating the response. The optimal values of the control factors were identified using the signal-to-noise (S/N) ratio optimization technique. Confirmation runs at the optimal parameters (i.e. thickness of 3 mm and 5 mm) showed a better response as compared to other conditions. Designers of social robots should consider the application of soft materials to their designs as it help in reducing the potential harm to the head.


Introduction
The interest in robots is increasing globally as estimated by the International Federation of Robotics (IFR) [1]. The application of robots is extending to new areas, such as that in healthcare. Most notably is the application of social robots in therapy sessions with children with autism, which has been reported to improve the overall outcomes [2]. However, such children exhibit a multitude of challenging behaviors that could raise some safety concerns when a robot is present in their vicinity [3]. The occurrence rates of challenging behaviors are high (e.g. 49% up to 69% [4][5] [6]), and that have many consequences on the services and treatments provided to them [7].
The ability to convey emotions, exhibiting different personalities, employing communication cues, and forming social relationships are some of the traits that make social robots different from typical toys [8][9] [10]. Introducing a new stimuli (e.g. a social robot) that is meant to elicit behaviors could provoke some unwanted behaviors among children with autism. For example, they might show some aggression toward the robots [11][12] [13]. Kicking, throwing, hitting, and banging are some of the challenging behaviors that could potentially cause harm during the interactions with social robots (Fig. 1). For example, the throwing of a small social robot to the head could cause superficial injuries, subconcussion or even concussion in extreme cases [3]. The occurrence of such behaviors implies the need for safer robotic designs [11] [14].
The establishment of safety standards in different fields of robotics is making a notably advances. However, the progress in establishing safety standards in relation to social robots and robotic toys is still lacking [15] [16]. Some of the existing safety standards in toys can be readily imported to cover some fundamental design and safety aspects. For example, the ISO 8124 standard [17]. Safety aspects of the mechanical and physical properties of toys are covered in part one of this standard while part two and three covers flammability and migration of certain elements, respectively. More rigorous Safety experiments for small robots investigating the potential of soft materials in mitigating the harm to the head due to impacts A PREPRINT Table 1: The standard L 9 (3 2 ) orthogonal array (OA) that was used in this study. A  B   1  1  1  2  1  2  3  1  3  4  2  1  5  2  2  6  2  3  7  3  1  8  3  2  9 3 3

Control factors
design considerations are needed that consider the unwanted behaviors exhibited by children with autism. For example, considering a scenario where the robotic toy is thrown to the head.
In this paper, we use Taguchi design of experiment (DOE) method to investigate a way to reduce the harm to the head by studying the influence of two control factors (i.e. soft material storage modulus and its thickness) and one noise factor (i.e. throwing velocity) of a small form factor toy on the resultant head's acceleration. Furthermore, the optimal levels of the investigated control factors that help in reducing the response are identified. This paper is organized as follows. Section 2 describes the Taguchi method while Section 3 provides the materials and methods. Section 4 presents the results of the study. Section 5 provides discussion while Section 6 concludes the work.

Background
The understanding of a process is usually done through conducting a series of experiments that provide information about the processes, such as the most influential parameters, the optimal settings, and the overall performance [18]. To maximize the information obtained while minimizing the resources needed, design of experiment (DOE) approach is usually used to study processes. Taguchi's method is one of the DOE techniques that was developed in 1979 to improve the quality of goods [19]. This technique has been considered in the optimization of different applications, such as strain measurements, wear studies, and turning parameters [20][21] [22].
There are two types of process factors that are considered in Taguchi method. The first factor is controllable at the product design level and referred to as the control factor. The other factor is known as the noise factor that is uncontrollable during production and at product design level. However, noise factor can be simulated experimentally.
Taguchi method aims to make products more robust by optimizing the control factors to reduce the sensitivity to noise factors. Furthermore, Taguchi method provides the quality loss function to measure losses in relation to product's variation from the target performance. In Taguchi DOE, different Orthogonal Arrays (OAs) can be considered depending on the number of factors being investigated [23].

Experimental design
Taguchi design method has been considered in this study because it considers the noise variables and because it greatly reduces the number of experiments. Furthermore, the method provides a robust parameter design by finding the optimal values of the control variables to reduce the investigated problem sensitivity to noise variables [24]. The Orthogonal Arrays (OAs) of this method provides more information to understand the relationship between the control variables and the noise variables, which is essential for a robust design. Additionally, Taguchi method provides performance measures that are needed in this study to optimize the design, such as the signal-to-noise (S/N) ratio.
In this study, two control factors are investigated for their potential in reducing the linear acceleration of the head. Experiments are conducted based on an L 9 (3 2 ) Taguchi OA (Table 1). A total of 27 (i.e. 9 × 3) experiments have to be conducted that consider the three levels of the control and noise factors. The considered control factors can be adjusted at the product design level while the noise factor is dependent on the real life scenario (e.g. throwing). Finally, the selected factors are independent while the measured output (i.e. head's acceleration) is dependent. Safety experiments for small robots investigating the potential of soft materials in mitigating the harm to the head due to impacts A PREPRINT The levels of the two control factors (i.e. material thickness and storage modulus) and noise factor (i.e. impact velocity) have been defined (Table 2). To achieve consistency, the mass and the shape of the impactor were kept the same throughout the experiments. The mass of the impactor was kept at 0.4 kg, which is within the expected range of the targeted applications (i.e. small robotic toys). The shape of the impactor was cylindrical without any features on the surface. Finally, the impact velocities used were limited to low velocities to achieve more consistency in terms of the noise levels (i.e. less than 3 m/s) [25].

Experimental Setup
The experimental setup was based on a 3D printed dummy head situated with nylon coated wire ropes in a dedicated frame (Fig. 2). The head was embedded with an accelerometer that measures the linear acceleration of the head. The accelerometer sensor was interfaced to the computer through a data acquisition card. The experimental setup has been validated in giving comparable results to those of previous studies under similar testing conditions. More detailed description of the experimental setup can be viewed in our earlier studies [3][26] [27].

Impactor
A 3D printed cylindrical object was used as an impactor in this study. The dimensions of the impactor were (10 × 10 cm 2 , height and diameter). A 3D printer (Replicator 5 th Generation, MakerBot Industries, USA) was used to build the object. Clay was used to fill the impactor to reach 0.4 kg. The soft materials (Ecoflex OO-30 & Dragon skin FX-Pro, Smooth-On, USA) were prepared according to manufacturer's instructions. The soft materials were prepared in molds of different thicknesses and then rectangular (5 × 8 cm 2 ) samples of each were attached to the impactor covering the area of impacts (Fig. 3).

Dynamic Mechanical Analysis
The soft materials properties were studied using a dynamic mechanical analyzer (RSA-G2, TA instruments, USA ; Fig. 4). The dynamic mechanical analysis (DMA) is a common test to measure the properties (i.e. elastic and viscous) of a material. The properties were studied by applying a stress (e.g. sinusoidal) and measuring the resultant strain and the phase difference between the input and output. A frequency sweep tests were conducted to study the storage modulus. In these tests, the frequency was varied from 0.1 Hz to 100 Hz while the strain and temperature kept constant. The storage modulus readings for each material were generated (Fig. 5). The values of storage modulus at 1 Hz for each material were considered in the analysis. This frequency is believed to be at which the high rate of challenging behaviors might occur [28]. Safety experiments for small robots investigating the potential of soft materials in mitigating the harm to the head due to impacts A PREPRINT Figure 2: The experimental setup used in this study [3]. Safety experiments for small robots investigating the potential of soft materials in mitigating the harm to the head due to impacts A PREPRINT Figure 3: Sample of the experiments conducted.

Impact Experiments
A total of 27 main experiments (i.e. 9 experiments at each noise level) were conducted (See supplementary material and data [29][30]).The experiments covered different combinations of the control factors according to the L 9 (3 2 ) orthogonal array (OA) ( Table 1). To achieve consistent impact velocities (i.e. noise levels), the impactor was tied to the frame allowing it to swing freely (Fig. 3). By altering the drop height, three different impact velocities levels were generated in the range of 1 -3 m/s ( Table 2). As for the confirmation runs, a total of 24 runs at each noise level were conducted after identifying the optimal control factors levels from the main experiments.

Peak Linear Head Acceleration
The peak linear acceleration of the head during impacts has been used as one of the biomechanical measures to investigate concussive events [31] [32]. In a study on young football players (age 7 -8 years) [33], the reported peak linear accelerations during impacts could fall anywhere in the range from 10 g to 111 g, and that is believed to be the range at which concussions could potentially occur based on the football-related incidents [34]. Furthermore, concussion could occur at a low peak linear acceleration of 31.8 g as reported by another study [35]. Another study analyzed 62,974 head acceleration data points and identified subconcussive impacts with an average head acceleration Safety experiments for small robots investigating the potential of soft materials in mitigating the harm to the head due to impacts A PREPRINT Figure 4: The Dynamic Mechanical Analyzer device that was used to analyze the soft materials.
of 26 ± 20 g with a median of 19 g while the average head acceleration for concussive impacts was 105 ± 27 g with a median of 103 g [36].
Laboratory settings using dummies are commonly used to evaluate the potential harm to a human by simulating a certain scenario, such as that used to evaluate harm levels due to impact in industrial robots [37]. Similar setups have been adopted in other fields to quantify harm due to impacts, such as that in sports to evaluate potential harm and the efficacy of protective gears [38][39] [40]. Similarly, a laboratory setup generating the peak linear head acceleration severity index will be used in this study as an indicator to understand the effects on the head linear acceleration due to altering some of the parameters of the impacting objects.

Data Analysis
Two main analyses were conducted on the generated results. The first analysis was a one-way analysis of variance (ANOVA). This test was conducted on the average responses of the peak head acceleration at each level of each factor to understand the contribution of each factor on the resultant response. The level for statistical significance was set to p < 0.05. The second analysis was based on the signal-to-noise (S/N) ratio and it was used for optimization. The goal of this analysis is to minimize the response variable, hence, the smaller-the-better case was considered and the corresponding S/N ratio was used; it is defined as follows [24]: where E is the expected value and y 2 i is the response. All analyses were performed using Minitab (v17.0, Minitab Inc., USA).

Orthogonal array
A total of 27 responses of the peak linear head acceleration were recorded and the corresponding average value, standard deviation, and S/N ratio for each combination were calculated to complete the Taguchi L 9 orthogonal array ( Table 3). The obtained linear acceleration values were in the range from 2.42 to 10.75 g due to different levels of control and noise factors. The lowest linear acceleration value obtained corresponds to a thickness of 3 mm, Ecoflex, and Level 1 impact velocity (i.e. A2-B1-X1) while the highest linear acceleration value corresponds to a thickness of 1 mm, Ecoflex, and Level 3 impact velocity (i.e. A1-B1-X3). The average response and the S/N ratio due to varying the level of each factor were tabulated (Table 4 and Table 5). The lowest average linear head acceleration was 3.18 g with S/N ratio of -10.10 occurred at Level 1 impact velocity while the highest average linear head acceleration was 10.10 g with S/N ratio of -20.09 due to Level 3 impact velocity.

ANOVA Tests
One-way ANOVA tests for each factor were conducted to understand if there is a significant difference due to varying the conditions on the resultant average peak linear head acceleration. In case the ANOVA test reported a significant difference, a post-hoc Tukey test was conducted. Safety experiments for small robots investigating the potential of soft materials in mitigating the harm to the head due to impacts A PREPRINT

Effect of the material's thickness
A one-way ANOVA was conducted to compare the effect of varying the three conditions of the thickness on the response ( Table 6). The test revealed that there was no significant difference due to varying the thickness on the resultant average linear head acceleration for the three conditions (F(2,24) = 0.04, p = 0.96) at the p< 0.05.

Effect of the material's storage modulus
A one-way ANOVA was conducted to compare the effect of the three different conditions of the storage modulus on the head peak acceleration ( Table 7). The test revealed that there was no significant difference due to changing the storage modulus on the resultant average linear head acceleration for the three conditions (F(2,24) = 0.02, p = 0.98) at p< 0.05. Safety experiments for small robots investigating the potential of soft materials in mitigating the harm to the head due to impacts A PREPRINT

Effect of the object's impacting velocity
A one-way ANOVA was conducted to compare the effect of the three different levels of the impact velocity on the head peak acceleration ( Table 8). The test revealed that there was a significant difference due to varying the impact velocity on the resultant average linear head acceleration for the three conditions (F(2,24) = 515.63, p = 0.00) at the p< 0.05. A post-hoc Tukey test showed that velocity Level 1 (M = 3.18, SD = 0.42), velocity Level 2 (M = 6.16, SD = 0.28), and velocity Level 3 (M = 10.10, SD = 0.61) were different significantly at p< 0.05.

Analysis
The alteration of the control factors and noise factor levels have an effect on the resultant head accelerations (Table  3). No definite trend can be observed between the thickness (i.e. Factor A) and the resultant response by visually investigating the orthogonal array. For example, looking at the response values in the noise factor columns X1 -X3 from row 1 to 9, the registered head acceleration value appear to remain consistent. Similar observation can be made for Factor B. One the other hand, the noise factor levels seem to affect the response significantly. For example, examining columns X1, X2, and X3 reveals that the response increases proportionally with the applied impact velocity (i.e. noise factor) as supported by the relatively large standard deviations across each row.
The peak linear head acceleration has been reported to be influenced by the impact velocity of an impactor [41] [37]. The increasing average values of the response at each level of the impact velocity supports these findings (Table 4). Furthermore, ANOVA test results and post-hoc Tukey findings showed that the impact velocity has a significant effect on the resultant peak head acceleration. On the other hand, no similar conclusion can be made for the control factors. For example, the reported ANOVA results for the effect of the material's thickness revealed that there was no statistical significance. This could be attributed to the relatively small selected thickness range and to the small difference between each level.

Optimization
The second goal of this study is to find the optimal values that reduces the response. This is achieved by investigating the mean value of the resultant head acceleration and the corresponding mean S/N ratio for each factor (Table 4 and Table 5). The plots were generated for the mean response at each factor level and the corresponding mean S/N ratios for a better visual comparison (Fig. 6). The criterion for selecting the optimal conditions is based on finding the levels that produce the lowest response and highest S/N ratio.
The best conditions for the thickness and storage modulus were identified based on the lowest generated head linear acceleration and highest S/N ratio (Fig. 6). For the control factor of material's thickness (i.e. Factor A), 3 mm (i.e. Level 2) and 5 mm (i.e. Level 3) achieved closely the best results. As for control factor of material's storage modulus (i.e. Factor B), ecoflex (i.e. Level 1) and clay (i.e. Level 3) scored closely the best results. Even though ANOVA tests on both of the control factors reported no significance in affecting the peak head acceleration, the identified levels for each factor provided the lowest responses and highest S/N ratios as compared to other conditions. Hence, these conditions  were selected as the optimal values. As for the noise factor, 1 -1.2 m/s impact velocities (i.e Level 1) scored the best results. The optimized conditions for the control factors investigated were 3 mm or 5 mm for the material's thickness factor and ecoflex or clay for the material's storage modulus (i.e. A2-B1, A2-B3, A3-B1, and A3-B3). Interestingly, the selected optimal control factors' levels produced relatively lower average peak head accelerations compared to other conditions at even higher noise levels ( Safety experiments for small robots investigating the potential of soft materials in mitigating the harm to the head due to impacts A PREPRINT

Confirmation Tests
After finding the optimal levels for each control factor, the last stage of Taguchi design is to perform the confirmation runs. The goal of this study is to minimize the peak head acceleration due to an object being thrown at the head by investigating the effect of two control factors, namely the material's thickness and its storage modulus. Hence, the optimal levels obtained in the previous section should produce relatively smaller head accelerations as compared to other conditions. Confirmation runs are needed to confirm these findings. To ensure that the optimal levels are robust and applicable to different noise scenarios, two confirmation runs will be conducted at every noise level.
A total of 24 confirmation runs were conducted at the optimal control factors' levels. For each control and noise factors combination, 2 runs were conducted and the corresponding mean values for each were calculated (Table 9). Comparing the results of the confirmation runs to that obtained from the main experiments, the average values were very close to respective ones obtained in the complete Taguchi orthogonal array (Table 3). Hence, the confirmation runs confirmed that the selected optimal levels produced the lowest peak head accelerations.

Limitations of the Study
This study considered only the application of three soft materials while there are many other candidates that could be considered, for example, PDMS [42]. The effects of the added mass of the soft materials were ignored (e.g Less than 0.05 kg). However, this added mass might influence the results significantly, especially when larger area is covered (e.g. covering the whole object with a soft material) or larger thickness is considered (i.e. greater than 5 mm). For consistency, the shape of the object was limited to one shape while the velocity of impacts was limited to low range. However, different shapes of robotic toys exist and higher impact velocities might occur in realistic scenarios. Other severity indices could have been considered to measure different potential harm. For example, measuring the soft tissue injuries and quantify the potential of soft materials in mitigating it.

Conclusion
In this study, the influence of an added soft material to an object on the linear acceleration of the head upon impact has been investigated. The Taguchi L 9 (3 2 ) orthogonal array design has been used to plan the 27 main experiments that were conducted. The control factors were the thickness and the storage modulus of three different soft materials. The noise factor was the impact velocity. The significance of each factor has been identified based on ANOVA tests while the optimal levels for the control factors were identified based on the analysis of S/N ratio. ANOVA tests showed that the control factors were not statistically significant in influencing the linear acceleration of the head. On the other hand, ANOVA test of the noise factor revealed that it was statistically significant. Material thickness of 3 mm and 5 mm achieved the best results. This implies that the application of a higher thickness of soft material will attenuate the head's acceleration better. Ecoflex and clay have achieved better response as compared to dragon skin. Confirmation runs at the optimal identified conditions achieved better responses as compared to other conditions.
The control factors, especially for the thickness of a soft material, may provide reduction in the overall head's acceleration. Hence, the manufacturers and designers of small robotic toys for special needs children should consider adding a layer of safe and soft material to their products to improve safety and to reduce any potential harm (e.g. subconcussions and superficial injuries). Additionally, they need to investigate different soft materials to find suitable materials that provide robustness and ease of application to their products while improving the safety aspects. Considering soft materials will also open for the possibility of embedding different sensors [43][44] [45]. The investigation conducted in this study used Taguchi method to design the experiments in which it provided a convenient, cost effective, and efficient way in the assessment and optimization of product designs concerning safety aspects. The emerging field of companion social robots can benefit from Taguchi methods in its approach to optimize a robot s design.