Influence of the shape and mass of a small robot when thrown to a dummy human head

Social robots have shown some efficacy in assisting children with autism and are now being considered as assistive tools for therapy. The physical proximity of a small companion social robot could become a source of harm to children with autism during aggressive physical interactions. A child exhibiting challenging behaviors could throw a small robot that could harm another child’s head upon impact. In this paper, we investigate the effects of the mass and shape of objects thrown on impact at different velocities on the linear acceleration of a developed dummy head. This dummy head could be the head of another child or a caregiver in the room. A total of 27 main experiments were conducted based on Taguchi’s orthogonal array design. The data were then analyzed using ANOVA and then optimized based on the signal-to-noise ratio. Our results revealed that the two design factors considered (i.e. mass and shape) and the noise factor (i.e. impact velocities) affected the response. Finally, confirmation runs at the optimal identified shape and mass (i.e. mass of 0.3 kg and shape of either cube or wedge) showed an overall reduction in the resultant peak linear acceleration of the dummy head as compared to the other conditions. These results have implications on the design and manufacturing of small social robots whereby minimizing the mass of the robots can aid in mitigating the potential harm to the head due to impacts.


Introduction
The global interest in robots is on the rise due to various advances in artificial intelligence and sensing technologies that made them more autonomous and intelligent ( Fig. 1) [17,42]. This opens up new possibilities in the application of robots in new areas [1]. These developments have implications and raise some safety concerns [2].
The interest in social robots and their clinical applications in assisting the therapy of children with Autism Spectrum Disorder (ASD) has grown considerably due to the reported evidences of their efficacy [15,21]. Unlike their neurotypical peers, children with ASD exhibit stereotypical and challenging behaviors. The reported prevalence rates of challenging behaviors among individuals with autism are high. Previous studies reported occurrence rates of greater than 49% and up to 69% [10,11,31].
Unlike typical toys, social robots have the ability to demonstrate emotions, establish social connections, display of personalities, use cues, and engage with partners at an interpersonal level [18,43,49]. The introduction of robots to children with ASD represents new challenge that must be taken into account. For example, children interacting with social robots might show some aggression [8,14,25]. Furthermore, robots are meant to elicit behaviors [19,21]. Hence, the introduction of such technology to children with ASD could represent a potential harm during the manifestation of unwanted behaviors (e.g. kicking, throwing, and banging [32]). For instance, a thrown small robot that hits the head might cause subconcussion or superficial injuries [2]. There is a need to evaluate social robots to achieve safe physical human-robot interaction during such scenarios.
When the head is subjected to bumps or blows and impacts with objects, a traumatic brain injury (TBI) might occur. In serious cases, TBI could lead to a disability or, in extreme cases, it could cause death. The occurrence of TBI among children was reported to cause challenges, such as disabilities and impairment in daily skills [9,39]. In 2013 alone, 2.8 million cases of TBI have been reported in the United States [47]. The cases of TBI among children were occurred due to getting struck by or against an object. A study analyzing mild brain injuries among children in Sweden for the years 1998 and 1999 has found that 47% of the cases occurred at home and during playtime or leisure activity and due to childcare products, which includes toys [23]. These figures have many implications on the design of robotic toys pertaining to the safety of the head. There is a need for further safety considerations and userfocused design to take into account the characteristics of special needs users, such as children with ASD [2,8]. The work in social robots safety is still limited [22,28], especially in relation to improving design aspects of small robots [13,16,29,48].
In this paper, we investigate two design factors (i.e. mass and shape) of a small robot that is subjected to different throwing velocities (i.e. noise factor) and understand their effects on the acceleration of a dummy's head. Furthermore, we identify the conditions of the design factors at which the response is minimized. This paper is structured as follows. Section 2 presents the design of experiment. Section 3 describes the materials and methods while Sect. 4 provides the results. Section 5 describes the discussion. Finally, Sect. 6 provides the conclusion.

Taguchi design
Experiments are usually conducted on processes and systems to understand the overall performance, to deduce the most influential parameters, and to determine the optimal settings to achieve the desired responses or overall goals [35]. Design of experiment (DOE) is an approach aimed at using the minimum amount of resources available while maximizing the amount of information obtained from a process through the selection of parameters to be investigated [35]. DOE helps in collecting different information by altering the optimal arrangement of parameters or factors to enhance product robustness [20,26]. There are many DOE techniques available. The choice depends on the investigated problem and the aim of the experiments conducted [20]. One of such techniques is the Taguchi method.
The method was developed in 1979 and was meant to be used as an off-line quality control tool to improve manufacturing products and goods [45]. Furthermore, it has been considered in different applications [7,27,34,46]. Taguchi DOE considers two types of variables, namely controllable and noise (or uncontrollable). Control factors can be controlled in production while noise  factors cannot be controlled, except experimentally. This method aims to improve the robustness of products against any variations in the noise factors by finding the optimal values of the controllable factors. Depending on the number of factors investigated, Taguchi DOE could take on different settings by considering different orthogonal arrays (OAs) [36]. The crossed array Taguchi design that was considered in this study provides robust solution by understanding the interaction between the control factors and the noise factors [20].

Experimental factors
The goal is to investigate whether the shape or the mass plays any role in affecting the response. Hence, two control factors are considered in this investigation for their possible influence in attenuating the peak acceleration of the head. Twenty seven (i.e. 9 × 3 ) experiments must be conducted to cover all the possible combination of all the factors based on the L 9 (3 2 ) OA ( Table 1). The considered control and noise factors are independent.
The levels of all the factors have been defined ( Table 2). The range for the mass is comparable to that of small robotic toys. For the sake of simplicity, the selection of the overall shapes of the form factor were limited to three basic 3D geometric shapes while ignoring any other features. For consistency, the noise levels were limited to low velocities divided into three levels (i.e. < 3 m/s) [37].

Impact setup
The impact setup that was used to conduct the experiments contained a 3D-printed head that was mounted on a frame (Fig. 2). The mass of the dummy head was made close to that of children's dummy heads [44]. To measure the linear acceleration of the head, an accelerometer was placed inside the head. The data was acquired at 20 kHz. The impact setup has been shown to give similar results to that of related studies. More indepth details about the experimental setup and validation can be found in our earlier studies [2,4,6,7]. Table 1 The standard L 9 (3 2 ) orthogonal array (OA) Run Control factors   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  . 2 The impact setup that was used to conduct the experiments showing the developed dummy head, data acquisition card, and the computer system that was used to acquire and store the data [2]

Impactors
We want to understand the influence of the mass and shape of a small robot on the resultant peak linear head acceleration due to an impact. Hence, 3D models of three basic shapes were considered (Fig. 3). The shapes were constructed using a 3D printer (Replicator 5th Generation, MakerBot Industries, USA). A clay material was used to adjust the mass of each shape according to the mass levels in Table 2. The center of mass was made sure to be balanced for all the impactors.

Procedures
All experiments were performed according to the L 9 (3 2 ) orthogonal array (OA) ( Table 1). For each noise level, 9 experiments were conducted that have covered all possible combinations of the control factors (See supplementary material and [3,5]). The objects were tied to the frame which provided a controlled condition in the execution of the experiments (Fig. 4). Furthermore, it has provided more consistent impact velocities by adjusting the drop height of an object. Three different drop locations generated three different noise levels. The impact velocities were based on the video analysis of the experiments ( Table 2). As for the confirmation runs, 18 experiments were conducted.
The impact velocities were estimated based on the slow-motion recordings of the experiments using a video analysis software (Tracker v5.0.7 [12]). The raw acceleration readings were obtained using a data acquisition card that then stored using a LabView script. Finally, the peak head acceleration for each experiment was calculated based on a MATLAB script. More detailed procedures can be found in our earlier studies [2,4,6,7].

Head severity index
Several head severity indices were used in the literature to study the potential harm to the head, such as Head Injury Criterion, 3 ms criterion, and peak head linear acceleration. Previous studies have considered the peak linear acceleration of the head to investigate concussive events due to impacts [41,51]. Among hockey players, an earlier study has reported the possibility of the occurrence of a concussion at 31.8 g [33]. Another study investigating football impacts identified the occurance of subconcussive events at 26 ± 20 g [40]. Dummies that report the head acceleration were also considered to evaluate the potential injuries to the head (e.g. in industrial robots and in sports [24,38,50]). Similarly, the impact setup in this study will use the peak linear head acceleration only as an indicator to sudy the influence of the considered factors. Fig. 3 The three basic 3D designs of the impactors that have been considered in the study. a Cube ( 10 × 10 × 10 cm 3 , length, width and height). b Cylinder ( 10 × 10 cm 2 , diameter and height).
c Wedge ( 10 × 10 cm 2 , length and height). For simplicity, the three basic featureless shapes were considered to isolate the contribution of the shape on the response Fig. 4 Samples of the experiments that were conducted in this study a For the cube. b For the cylinder. c For the wedge

ANOVA
To study the significance of each factor on the head's acceleration, a three-way analysis of variance (ANOVA) test was conduced on the inverse of the average responses of all the investigated factors. The level for statistical significance was set to p < 0.05 . All analyses were performed using Minitab (v19.1, Minitab Inc., USA).

Signal to noise ratio
The optimization technique considered is the signalto-noise (S/N) ratio. The aim of this study is reduce the head's acceleration. Hence, the appropriate S/N ratio was selected, and it is defined as [20]: where E is the expectation and y 2 i is the head's acceleration (i.e. the response).

Orthogonal array
The Taguchi L 9 orthogonal array was completed by finding the average response, the standard deviation, and the respective S/N ratio for each combination of the investigated factors ( Table 5). The range for the linear acceleration values was 2.72 to 13.03 g. The lowest response value was at a mass of 0.3 kg, wedge shape, and Level 1 impact velocity (i.e. A1-B3-X1). As for the highest response, it occurred at a mass of 0.5 kg, cylinder shape, and at an impact velocity Level 3 (i.e. A3-B2-X3). The lowest overall average response across all factors was 3.95 g with a S/N ratio of − 12.05 and was due to Level 1 velocity (Tables 3 and 4). On the other hand, Level 3 velocity scored the highest average response (ie. 9.70 g) with a S/N ratio of − 19.91.

ANOVA
The statistical significance of each factor on the average response due to varying their levels was investigated based on a three-way ANOVA test. The inverse transformation of the response was considered in this test. A post-hoc pairwise Tukey test was also conducted.

The object's mass
To understand the contribution of the mass on the head's acceleration, a three-way ANOVA was conducted ( Table 6). The test revealed a statistical significance due to altering the conditions of the mass on the response

The object's shape
To investigate the effects of the shape on the response, a three-way ANOVA test was conducted (Table 6). There was a significant difference due to the alteration of the shape on the response (F(2,8) = 7.03, p = 0.02 ) at p < 0.05 . A post-hoc Tukey test showed that Level 2 shape (M = 7.88, SD = 1.21) differed significantly from Level 3 shape (M = 6.69, SD = 1.49).

The object's impact velocity
To understand the contribution of the impact velocity, a three-way ANOVA was performed (Table 6). A significant difference was reported due to altering the con-   and Level 3 impact velocity (M = 9.70, SD = 2.07) were statistically different at p < 0.05.

Interactions
There was no significant interaction between the mass and the shape (i.e. A and B) on the response (F(4,8) = 1.22, p = 0.37 ) at p < 0.05 . There was no significant interaction between the mass and the impact velocity (i.e. A and X) on the response (F(4,8) = 2.70, at p = 0.11 ) at p < 0.05 . There was no significant interaction between the shape and the impact velocity (i.e. B and X) on the response (F(4,8) = 3.91, at p = 0.05 ) at p < 0.05.

Discussion
The responses of the 27 experiments conducted suggest that the control factors and noise factor have influenced the response. The velocity and mass of impactors have been reported to affect the magnitude of the peak linear head acceleration [24,30]. The average responses and the results of the ANOVA test were in accordance with these findings (Tables 3 and 6). As for the shape factor, it appears to affect the response the least ( Table 3). The average response for the cylinder shape was higher compared to the other shapes.

Optimization
The goal of Taguchi design in this study was to optimize the control factors in order to reduce the head linear acceleration. The optimization was accomplished by inspecting the average response and the corresponding average S/N ratio at every level of the control factors (Tables 3 and 4; Fig. 5). The optimization case in this study is to produce the lowest response and highest S/N ratio, hence, the focus is on the factor conditions that satisfy this criterion. The optimal conditions for the mass and shape were found (Fig. 5). The mass of 0.3 kg achieved the lowest response. As for control factor of the shape, cube (i.e. Level 1) and wedge (i.e. Level 3) scored comparably the best. For the noise factor, Level 1 (i.e. 1.14 ± 0.10 m/s) achieved the lowest response.
The identified conditions were 0.3 kg for the mass factor and cube or wedge for the shape (i.e. A1-B1 and A1-B3, respectively). Compared to other combinations, the identified optimal conditions generated relatively lower responses, for example, the individual responses at noise Level 2 and Level 3 (i.e. X2 and X3) for runs 1 and 3 (Table 5).

Confirmation runs
The confirmation runs are needed to validate the optimal conditions that were obtained in the previous section. To achieve that, 18 experiments were conducted (Table 6). For each of the identified conditions in the previous section, 3 runs were performed ( Table 7). The confirmation runs showed that the average response values were comparable to respective ones obtained in the complete OA  (Table 5), hence, confirming the optimal identified conditions for the control factors.

Limitations of the study
The investigation conducted in this study was limited to three factors conducted in a laboratory developed experimental setup with controlled conditions that might differ from an actual and more dynamic scenario (e.g. potential for secondary impacts). The values of the object's mass that were tested were limited to the range of 0.3-0.5 kg. In reality, lighter or heavier robotic toys could exist. Studying the effect of shape was limited to three basic featureless shapes. However, actual robotic forms could have other unusual shapes with a lot of detailed features. Due to the way the experimental setup was designed, the impact area for each shape was limited to achieve consistency between the experiments. The effects of different impact areas and their interactions with different shapes and surface's features were not investigated in the current work. The noise levels that were tested were limited to low velocity impacts. Hence, the generated head's accelerations from our experiments were limited to a small range (i.e. < 10 g), which is unlikely to cause any potential harm. However, the toy robot could be subjected to higher velocities in an actual scenario. Finally, the investigated severity index was limited to the head acceleration. Other severity indices (e.g. tissue injuries) could be considered to investigate different effects due to different design parameters.

Conclusion
During meltdowns, a child interacting with a social robotic toy could involuntary throw the robot to others and cause harm. In this study, the Taguchi parameter design method was used to identify the optimal design parameters for small form factor robotic toy in order to reduce the peak linear head acceleration due to impacts on the head. The investigated control factors were the mass and the shape. The impact velocity was the noise factor. Based on L 9 (3 2 ) orthogonal array, a total number of 27 experiments were conducted covering the possible combinations of the control factors and the noise factor. The optimal levels of the shape and mass that minimize the peak head linear acceleration were found based on the S/N ratio. The three-way ANOVA test revealed statistical significance for the control factors and the noise factor in influencing the head's acceleration. The optimal levels for the control factors were as 0.3 kg for the mass and cube or wedge for the shape. The confirmation runs at the optimal conditions for the control factors provided the best responses as compared to other conditions. For a given impact velocity, reducing the mass of an object appears to reduce the head acceleration considerably. This was more evident at higher impact velocities. This finding combined with our earlier findings of the application of soft materials could provide a better overall reduction in the potential harm [7]. Therefore, the designers of small social robots or robotic toys for special needs children, or even for neurotypical children, should try to minimize the mass of their products as much as possible while adding a layer of a suitable soft material. In addition to the optimization efforts at the design level, the manufacturers should include special warning labels for their products that are meant to be used by special needs users, especially those that exhibit challenging behaviors.
Optimization of a robot's design can benefit from the techniques provided by Taguchi methods as it gives an efficient and convenient way to assess and optimize the safe design of small social robots. Moreover, such benefits are magnified for two reasons. On one hand, the safety requirements of the target users, such as children with special needs who may have a tendency for meltdowns, are addressed. On the other hand, the manufacturers have done their due diligence in optimizing their design for minimizing the chances of harm to the users and avoid lawsuits or product recall later on.