Geometric variation simulation and robust design for flexible cables and hoses

Abstract

This article presents how to do stability analysis, variation simulation and tolerance envelopes for deformable cables and hoses. In this way, the well-established methods for analyzing and minimizing the propagation of geometric variation are extended from rigid and sheet metal assemblies to large deformations of slender parts. This extension is highly motivated, not at least, by the increased amount of cables and hoses in electrified and hybrid solutions in the automotive industry. To secure the geometric correctness of the analysis, we use a nonlinear rod approach in the implementation of a simulation model of flexible cables that accounts for large deformations and supports a wide range of common cable clip types, cable joints and branches. A selection of measures and visualizations that reflect the spatial effects of geometric variation in deformable cables is incorporated in established variation analysis techniques. In particular, the concept of tolerance envelopes for rigid parts has been extended to also handle system of slender parts such as cables and hoses. This article extends previous research on robust design and variation simulation to cables subject to large deformations and has been successfully applied on two typical industrial cases: a wiring harness attached with various types of clips to a static surrounding and a cooler hose subjected to a high-frequency sampled engine motion.

Introduction

Geometric variation, originating from individual manufacturing and assembly processes, often propagates and accumulates during production, resulting in nonnominal products. Geometric quality problems are often discovered during preproduction or when the product is getting ready for market introduction. A change in the product or production concept at this stage often results in huge costs for product and/or production changes, market delays and bad publicity. Therefore, more and more efforts are made in early concept phases to virtually verify product and production concepts with respect to geometric variation.

Recently, we have seen a shift of focus of the automotive industry toward electrified and hybrid solutions; controlling and predicting the geometric variation in deformable parts such as electric harnesses and hoses become increasingly important as the competition for unoccupied space is intensified. For example, both conventional combustion engines and battery supplied electrical engines need to fit in an already densely packed vehicle. The placement of each component must be evaluated with respect to geometrical interference with other disciplines and their components considering component and assembly variations.

The geometry assurance process

The process aiming at controlling geometric variation throughout the whole product realization process is called the geometry assurance process. In the concept phase, the design concept and the production concept are developed. In the verification phase, the design and production concepts are physically verified by prototypes and test series. In the production phase, the product is in full production. The effect of geometric variation may be avoided, identified and controlled by a number of activities in these three main phases.

In the concept phase, the effect of variation is avoided by robust design concepts. ¹ Tolerances are simulated and verified with respect to fulfillment of product requirements; see, for example, studies by Gao et al. ² and Glancy and Chase. ³ Robust fixture design is considered by Cai et al. ⁴ and Wang. ⁵ The philosophy of robust design is combined with the philosophy of axiomatic design by Söderberg and Carlson ⁶ in order to achieve an assembly concept that is also easy to adjust and control.

In the verification phase, physical prototypes and test series are used to verify the concept geometrically. Prototypes verify the overall concept with respect to nominal fit and function, whereas test series verify the concept with respect to manufacturing and assembly variation under more realistic circumstances. These verification activities are quite expensive and time-consuming. Therefore, simulation techniques are being used more and more for nonnominal verification to reduce the needs for physical test series, see, for example, studies by Spensieri et al. ⁷ and Van den Berg et al. ⁸

In the production phase, inspection data are used for diagnosis and root cause analysis. Individual processes are monitored and kept under control using techniques for statistical process control. Overall trends are identified using univariate or multivariate techniques, see, for example, the study by Johnsson and Wichern. ⁹ Fault isolation in the auto-body assembly process using inline measurements has been addressed by, for example, Hu and Wu, ¹⁰ Ceglarek and Shi, ¹¹ Jin and Shi, ¹² Ding et al.,^13,14 Camelio et al. ¹⁵ and Carlson and Söderberg. ¹⁶

Variation simulation and robust design are an active research area, from which methods in order to control and reduce variation in rigid assemblies have evolved and are now also well established in the industry. When dealing with compliant parts, which is not that well established yet in industry, there are still many significant published results. Flexible assemblies have been treated by Liu et al. ¹⁷ Ceglarek and Shi ¹⁸ used a linear beam–based model for tolerance analysis in sheet metal assemblies. Sellem and Rivière, ¹⁹ Camelio et al. ²⁰ and Cai ²¹ used finite element analysis (FEA) for the same purpose. Hu and Wu ¹⁰ presented the “stream-of-variation” theory for automotive body assembly. Wärmefjord et al. ²² investigated the impact of spot welding sequences for variation in sheet metal assemblies. When the product involves complicated systems of deformable cables such as wiring harnesses, however, the linear FEA is not sufficient to properly capture the highly nonlinear geometric response in the material even when subject to small variations.

Scope of the article

The research of this article extends a set of well-established simulation methods for analyzing and controlling the propagation of geometric variation in assemblies to also handle deformable cables, hoses and wiring harnesses. In order to accomplish this, new components need to be in place. Crucial for the geometric correctness of the analysis are the choice and implementation of a simulation model for cables that accounts for large deformations and supports a wide range of common cable clip types, cable joints and branch breakouts. A computationally efficient implementation of the model allows then for carrying out a variation simulation by means of Monte Carlo iterations. In addition, a selection of critical measures and visualizations that reflect the spatial effects of geometric variation in deformable cables is incorporated in existing variation analysis techniques. In particular, a novel method to generate tolerance envelopes for compliant parts is introduced.

Structure of the article

The article is structured as follows. The next section describes Cosserat rods, the simulation model that is used in order to fast and accurately compute the static equilibrium of flexible cables. The proposed simulation methodology is presented in the section “Variation analysis”: different geometric similarity measures for flexible cables are discussed and the concepts of stability analysis, tolerance envelopes and full-scale variation simulation are introduced. The methodology is then applied to two different industrial cases in the “Case studies” section, and the results are finally concluded in the “Conclusion” section.

Cosserat rods

The cable simulation model must both account for large deformations and be computationally efficient. Cosserat rods are gaining widespread popularity for exactly these reasons and constitute a specialized topic geometrically nonlinear structural mechanics.^23,24 By using innovative mathematical techniques in combination with numerical procedures specifically tuned to the application problem, it is possible to retain the necessary physical accuracy even on very coarse grids. ²⁵

Rod kinematics

Cables are in this article regarded as deformable objects geometrically characterized as one-dimensional slender structures. This type of objects can express large overall deformations although locally the stresses and strains remain small. Under the assumption that the cable cross section always remains planar and rigid, the kinematics of a cable are captured in a framed curve.

A Cosserat rod of undeformed length L is a framed curve parameterized by arc length s

$$s:[0,L]\mapsto (\mathbf{\phi }\left(\mathbf{s}\right),\mathit{R}(\mathbf{s}))\in \mathrm{SE}\left(3\right)$$

(1)

$$\mathbf{\phi }$$ is the center curve of the cable piercing the centroids of the cross section along the rod and $$\mathit{R}\left(s\right)=({\mathit{d}}_1,{\mathit{d}}_2,{\mathit{d}}_3)$$ describes the cross-section orientation; $${\mathit{d}}_1$$ and $${\mathit{d}}_2$$ are the directors in the cross-section plane and $${\mathit{d}}_3={\mathit{d}}_1\times {\mathit{d}}_2$$ is orthogonal to the cable cross section (remark: for a Kirchoff rod, $${\mathit{d}}_3$$ is always tangent to the center curve; however, for the Cosserat rod, $${\mathit{d}}_3$$ can deviate from the center curve but with a shear force penalty induced by stiffness in the shear directions). The composition of position and orientation is an element in $$\mathrm{SE}\left(3\right)$$ , the Special Euclidean Group that describes rigid body transformations in $${ℝ}^3$$ . Any material point in the undeformed rod is then uniquely identified in the deformed rod by the deformation mapping $$\mathbf{\chi }:[0,\mathbf{L}]\times A\to {ℝ}^3$$

$$\mathbf{\chi }(\mathbf{s},{\mathbf{\xi }}_1,{\mathbf{\xi }}_2)=\mathbf{\phi }\left(\mathbf{s}\right)+{\mathbf{\xi }}_1·{\mathit{d}}_1\left(\mathbf{s}\right)+{\mathbf{\xi }}_2·{\mathit{d}}_2\left(\mathbf{s}\right)$$

(2)

where $$A\subset {ℝ}^2$$ is the cross-section domain (Figure 1).

Figure 1.

A Cosserat rod represented by a center curve and a frame field describing the cross-section orientation.

Mechanical equilibrium

To find the (static) mechanical equilibrium, one seeks a stationary point to the potential energy functional W of which the elastic part is composed of quadratic forms in terms of frame-invariant strain measures.

The stretching/shearing strain vector $$\Gamma =({\mathbf{\Gamma }}_1,{\mathbf{\Gamma }}_2,{\mathbf{\Gamma }}_3)$$ and curvature/torsion vector $$\mathbf{\Omega }=({\mathbf{\Omega }}_1,{\mathbf{\Omega }}_2,{\mathbf{\Omega }}_3)$$ in material coordinates read

$$\Gamma \left(\mathbf{s}\right)=\mathit{R}{\left(\mathbf{s}\right)}^{\mathbf{T}}\frac{\partial \mathbf{\phi }}{\partial \mathbf{s}}\left(\mathbf{s}\right)-{\mathit{e}}_3$$

(3a)

$$\widehat{\Omega }\left(s\right)=\mathit{R}{\left(s\right)}^T\frac{\partial \mathbf{\phi }}{\partial s}\mathit{R}\left(s\right)$$

(3b)

$${\Gamma }_{1,2}$$ are the shearing strain components in the 1- and 2-directions and $${\Gamma }_3$$ is the tension strain. $${\Omega }_{1,2}$$ are the bending curvature components and $${\Omega }_3$$ is the torsion corresponding to direction 3. Also, $$(\widehat{·})$$ is the “cross-product operator” identifying rotation axis vectors in $${ℝ}^3$$ with a skew 3 × 3 tensor and $${\mathit{e}}_3$$ is the third standard basis vector. Under the assumption of a hyperelastic constitutive relation (which can be justified since locally $$\Gamma$$ and $$\mathbf{\Omega }$$ are small), the stored shearing/stretching and bending/torsion energy densities in a deformed rod read

$$w^{\Gamma }\left(s\right)=\frac{1}{2}\Gamma {\left(\mathbf{s}\right)}^{\mathbf{T}}{\mathit{K}}^{\Gamma }\Gamma \left(\mathbf{s}\right)$$

(4a)

$$w^{\Omega }\left(s\right)=\frac{1}{2}\mathbf{\Omega }{\left(\mathbf{s}\right)}^{\mathbf{T}}{\mathit{K}}^{\mathbf{\Omega }}\mathbf{\Omega }\left(\mathbf{s}\right)$$

(4b)

for some stiffness tensors $${\mathit{K}}^{\Gamma }\text{and}{\mathit{K}}^{\mathbf{\Omega }}$$ . The (static) mechanic equilibrium for a single cable in the presence of a gravitational force field $$\mathit{g}$$ is found as the solution to the balance equations

$$\frac{\partial }{\partial s}\left(\mathit{R}{\mathit{K}}^{\Gamma }\Gamma \right)-K_{\rho }\mathit{g}=0$$

(5a)

$$\frac{\partial }{\partial s}\left(\mathit{R}{\mathit{K}}^{\Omega }\mathbf{\Omega }\right)+\frac{\partial \mathbf{\phi }}{\partial s}\times (\mathit{R}{\mathit{K}}^{\Gamma }\Gamma )=0$$

(5b)

or (equivalently and more efficiently) as the stationary states to the total potential energy

$$\partial W=0$$

(6a)

$$W=\underset{s=0}{\overset{L}{{\displaystyle \int }}}\left\{w^{\Gamma }\left(s\right)+w^{\Omega }\left(s\right)-K_{\rho }〈g,\phi (s)〉\right\}ds$$

(6b)

An adaptive finite difference discretization of a Cosserat rod coupled with an efficient quasi-Newton method to minimize the energy together with a set of boundary conditions forms the foundation of our cable simulation.

Boundary conditions

A wide range of common cable clips and different cable joints and branches are encountered in many industrial cases (Figure 2). It is geometric variation in the realization of these boundary conditions that we later wish to analyze. For this purpose, a general collection of boundary conditions needs to be supported.

Figure 2.

A wiring harness with clip and branch breakouts.

A cable clip can constrain the cable in different ways. For example, a clip can lock the position and the $$d_3$$ tangent direction of a cable but allow twists, mimicking the behavior of a steel ring. Table 1 lists a collection of different clips that have been implemented. A natural extension is that a set of clips could be connected; their relative locations are kept fixed, but they are free to move as a group of objects. In this way, any kind of cable branches and joints can be represented. All the boundary conditions are represented by elimination of degrees of freedom in the rod model.

Table 1.

A collection of implement cable clip constraints embedded in the simulation model.

Clip type	Position constraint	Frame constraint
Position and frame	$$\mathbf{\phi }\left(\mathbf{s}\right)={\mathbf{\phi }}_0$$	$$\mathbf{R}\left(s\right)={\mathbf{R}}_0$$
Position and direction	$$\mathbf{\phi }\left(\mathbf{s}\right)={\mathbf{\phi }}_0$$	$${\mathit{d}}_3\left(s\right)=v_0$$
Position	$$\mathbf{\phi }\left(\mathbf{s}\right)={\mathbf{\phi }}_0$$	Free
Free	Free	Free

An efficient implementation of a discrete Cosserat rod model together with the generalized set of boundary conditions allows us to gather large quantities of physically reliable cable simulation data to fully support the methodology presented in the next section.

Variation analysis

A manufacturing process is always afflicted by variation in various product design variables, causing deviation from nominal values in critical measures on the final product. Variation analysis aims to early in the geometry assurance process detect design and assembly solutions that are sensitive to variation and may cause problems later on during production.

Formally, we assume that the configuration of a product $$P$$ is determined by a set of stochastic design variables $$\mathit{x}=\{x_1,x_2,\dots ,x_n\}$$ , each $$x_i$$ belonging to a statistical distribution and having a dimensional tolerance $$[l_i,u_i]$$ . One or more critical measures $$\mathit{m}=\{m_1,m_2,\dots ,m_m\}$$ may be defined, measuring some property important for the functionality or the aesthetical appearance of $$P$$ . The relation between design variables and critical measures is called the design equation

$${\mathit{M}}_P\left(x\right)={\left(\begin{array}{c}\hfill m_1\left(x\right)m_2\left(x\right)\dots m_m\left(x\right)\hfill \end{array}\right)}^T$$

(7)

For deformable cables, design variables could mean the length of the cable, the orientation of the cable clips or material parameters, while critical measures of interest could be the maximum stress in the cable, normal force exerted on a cable clip or the shortest Euclidean distance between the cable and the surrounding.

This section describes a methodology for variation analysis of deformable cables with a focus on controlling the spatial effects of variation in geometric design variables. The methodology is tailored for early stages in the geometry assurance process when real prototypes and inspection data are not available and consists of the following:

• Stability analysis. Evaluation of the robustness of a cable design by linearization around the nominal design. The analysis indicates how design variables contribute to design sensitivity and how they may be updated in order to achieve a cable design less sensitive to fixture variation.

• Tolerance envelope analysis. Generation of the smallest volume enclosing all cable shapes satisfying a specified set of tolerances. A tolerance envelope can be used to validate that variation in a cable design does not interfere with components in the surrounding—or reversely—to serve as a forbidden region when considering the design and location of the surrounding components.

• Variation simulation. Geometric variation simulation by means of Monte Carlo iterations. The simulation captures nonlinearities and provides a good knowledge about the variation effects and serves as good feedback to the tolerancing of the design.

Geometric design variables

The geometric design variables considered in this article describe either the geometric part properties of the undeformed cable shape (cable length, radius, etc.) or fixturing variables that control the locations of the cable clips.

Typically, a cable clip is rigidly attached to a specific component, inheriting the location scheme of that component. Conventional assembly locating schemes such as the 3-2-1 system have been studied in the context of robust fixture design. ¹ We will additionally assume that there might be a variation in the assembly of the clip itself. For simplicity, we assume that the variation in clip orientation $$\mathit{R}$$ is locally expressed in terms of the roll, pitch and yaw (RPY) convention

$$\mathit{R}(\alpha ,\beta ,\gamma )={\mathit{R}}_z\left(\alpha \right)·{\mathit{R}}_y\left(\beta \right)·{\mathit{R}}_x\left(\gamma \right)$$

(8)

Here, $${\mathit{R}}_x$$ , $${\mathit{R}}_y$$ and $${\mathit{R}}_z$$ express rotations around the local $$x$$ , $$y$$ and $$z$$ axes with angles $$\alpha ,\beta$$ and $$\gamma$$ , respectively. This parameterization makes sense since small rotations are nearly commutative, and we are far away from the parametric singularity $$\beta =\pi /2$$ .

Variation measures

A variety of measures can be used to quantify the geometric variation in the final product. A measure can estimate the total variation for a body by, for example, statistical root sum square (RSS) methods, or it can relate to a dimension critical for the functionality or the aesthetical appearance of the product. In our study, two types of measures are of additional interest: shortest distance to a surrounding component and a geometrical similarity relation with respect to the nominal cable shape.

The cable response function $$\mathit{x}\to ({\mathbf{\phi }}_{\mathit{x}},{\mathbf{R}}_{\mathbf{x}})$$ maps the design variables $$\mathit{x}$$ to a rod in equilibrium. This relation could be one-to-many since there might exist multiple solutions satisfying the same set of boundary conditions. However, we may assume that locally around a cable shape in equilibrium, small perturbations in design parameters attract the cable to the “same” energy minima. We formally define the cable shape $$C\subset {ℝ}^3$$ as the trace of the deformation mapping $$\mathbf{\chi }:C=\mathbf{\chi }([0,\mathbf{L}]\times A)$$ and the cable curve $$c$$ as the trace of $$\mathbf{\phi }$$ , $$c=\mathbf{\phi }\left([0,\mathbf{L}]\right)$$ . The distance from a cable shape $$C$$ to a surrounding component $$S$$ is then defined as $$d^{*}={inf}_{u\in C,v\in S}\Vert \mathit{u}-v\Vert$$ . If the cable has a circular cross section of radius $$r$$ , then $$d^{*}$$ can be computed as $$d^{*}=max\{0,{inf}_{\mathit{u}\in c,\mathit{v}\in S}||\mathit{u}-v||-r\}$$ .

If there is little knowledge about the surrounding or a strong confidence that the nominal cable design is the best from a geometric perspective, then any geometric deviation from that is associated with a risk. Given a perturbation $$\tilde{\mathit{\chi }}$$ of the nominal cable shape function, we can compute the largest deviation of a material point as $$\underset{[0,L]\times A}{sup}\tilde{\chi }(s,{\xi }_1,{\xi }_2)-\mathbf{\chi }(\mathbf{s},{\mathbf{\xi }}_1,{\mathbf{\xi }}_2)$$ . This point wise similarity does indeed provide an upper bound on how much any material point in the cable has been displaced; however, it says little about the displacement of $$C$$ . Hence, it makes sense to also define a more general similarity measure. The Hausdorff distance between the two sets $$U$$ and $$V$$ is defined as

$$d^H(U,V)=max\{d_U\left(V\right),d_V\left(U\right)\},\text{where}{d}_U\left(V\right)=\underset{\mathit{u}\in U}{sup}{\inf }_{v\in \mathit{V}}||\mathit{u}-v||$$

(9)

Proximity queries between triangulated objects and/or space curves represented by polylines have been efficiently implemented, ²⁶ and Aspert et al. ²⁷ described how a discrete Hausdorff distance can be computed (Figure 3). The similarity measures presented in this section are true worst-case measures.

Figure 3.

(a) A distance measure between a cable shape and a static geometry. (b) The Hausdorff distance between a perturbed cable and its nominal shape.

Stability analysis

A stability analysis evaluates the robustness of the nominal design and indicates how design parameters contribute to variation and how a location scheme or design variable may be updated in order to achieve a cable design less sensitive to variation. The analysis is carried out by linearizing the design equation and studying the Jacobian, commonly referred to as the sensitivity matrix, at the nominal design

$$\left(\begin{array}{c}\Delta m_1\\ ⋮\\ \Delta m_m\end{array}\right)\approx \left(\begin{array}{ccc}\partial m_1/\partial x_1& \cdots & \partial m_1/\partial x_n\\ ⋮& \ddots & ⋮\\ \partial m_m/\partial x_1& \cdots & \partial m_m/\partial x_n\end{array}\right)\left(\begin{array}{c}\Delta x_1\\ ⋮\\ \Delta x_n\end{array}\right)$$

(10)

In the case of our Hausdorff similarity measure, we cannot evaluate the Jacobian since the measure is not even differentiable at the nominal design. For this measure, we will consider $\max \left\{\left|\partial m_i/\partial x_j^{-}\right|,\left|\partial m_i/\partial x_j^{+}\right|\right\}$ when evaluating the sensitivity matrix

Following the concept of Suh, ²⁸ a good coupled design is generally one with a distinct coupling between as few variables and with as little amplification as possible. A more robust design is achieved by an iterative procedure where the contribution to sensitivity is identified by the sensitivity matrix and resolved by updating the configuration of the product, for example, changing positions of the corresponding locators.

Tolerance envelopes

A tolerance envelope $$V^{*}$$ is the smallest volume in $${ℝ}^3$$ enclosing all cable shapes satisfying a specified set of design tolerances

$$V^{*}=\mathrm{argmin}|V|\text{such}\text{that}{C}_u\subset V\text{for}\text{all}u,\text{where}{u}_i\in [l_i,u_i]\text{for}i=1,\dots ,n$$

(11)

$$C_{\mathit{u}}$$ is the cable shape corresponding to design variable $$\mathit{x}$$ with values $$\mathit{u}$$ .

Lööf et al. ²⁹ generated a motion envelope for rigid bodies by computing the convex hull of a large set of variation simulation outcomes combined with convex partitioning. To generate a tolerance envelope of a triangulated cable shape, we instead deterministically sample the tolerance regions for all design variables and build up a collection of cable shape extrema, each with the same number of triangles generated via the deformation mapping $$\mathbf{\chi }$$ . For a fine enough sampling resolution, we assume that the variation of an individual triangle is always contained within the convex hull of the extreme triangle locations. The tolerance envelope is then the union of the convex hulls corresponding to the entire triangulation (Figure 4). This set can be efficiently generated with respect to memory consumption by implicitly tracing all triangle variations through a marching cubes algorithm. ³⁰

Figure 4.

(a) A set of tolerance extreme shapes and (b) the corresponding convex hull of one material triangle. (c) The entire tolerance envelope.

Variation simulation

A variation simulation analysis is carried out in order to gain knowledge about how the variation in the design variables propagates to the critical measures. It validates the robustness and serves as feedback to the tolerancing of the design. Typically, in a CAT software, the variation is color coded in order to get a good overview of areas with high and low degree of variation (see Figure 5).

Figure 5.

This figure illustrates two rear lamps of a car with different locating schemes. The left lamp has a more spread out set of locating points and as the color coding suggests is clearly the best locating scheme with respect to part variation.

To perform the variation simulation, the Monte Carlo simulation method (MC) is used. MC samples a large population of design variable values randomly (within their specified statistical distributions) and evaluates the critical measures in order to estimate their statistical properties. The method captures nonlinearities and allows any kind of statistical distribution as input parameter variation.

Case studies

Two industrial case studies have been conducted in order to test the methodology in real-life scenarios. The cases consist of highly detailed triangulated data and complex motions sampled with high frequency.

Wiring harness installation

This case considers a wiring harness of length 312 mm with a breakout at arc length position s = 258 mm. The three ends are mounted to parts of a car, see Figure 6. The task is to verify the robustness of the design with respect to the cable connector WH_Start as indicated in the figure. WH_Start has a rotational variation around the local y axis with tolerance ±5°.

Figure 6.

(a) Tolerance envelope and color-coded variation. (b) WH_Start and the sensitivity matrix of the initial design. (c) The final robust design at a safe distance from the car.

First, a tolerance envelope is generated to detect geometrical interference with other components within the estimated set of tolerances. It turns out that the tolerance envelope intersects the car geometry. A color coding indicates that the variation of the harness is indeed the largest in areas where the interference occurs. In order to improve the robustness of the design, the stability matrix is computed for the distance to the car D_CAR and for a similarity measure D_SIM. D_SIM is 31 mm/°, indicating that a perturbation of 1° in WH_Start affects the design of the cable with 31 mm. WH_Start is rotated by 3° around the local z axis toward the other harness ends, and we observe that both sensitivities are now reduced. The procedure is reiterated until we arrive at a design that is both robust and at a safe distance from geometric interference with other components.

Engine motion

In this scenario, a cooler hose of length 330 mm is attached to an engine, and the engine follows a recorded engine motion. The point of the hose connected to the cooler has a variation of size in the translatory degrees of freedom. Similar to the previous case, in order to analyze the motion, a tolerance envelope for the hose with respect to the entire motion is needed. To generate such a volume, we first divide the motion into a large number of motion time segments. For each such time segment $$[T_i,T_{i+1}]$$ , we collect the tolerance extreme shapes corresponding to $$T_i$$ and $$T_{i+1}$$ , respectively, and generate a tolerance envelope with the algorithm described in the section “Tolerance envelopes.” By incrementally adding tolerance envelopes of all motion segments and removing any enclosed triangles, a volume bounding all possible variations of the hose is finally generated. In addition, all critical measures are generalized to take their maxima over the entire motion. The iterative procedure of the previous case is applied in this generalized setting, achieving a robust and variation-insensitive design to the specific type of motion (see Figure 7).

Figure 7.

(a) A tolerance envelope for the entire motion indicates a potential interference with the surrounding (distance 4 mm). (b) Color coding of variation with respect to cable length (330 mm). (c) The tolerance envelope of the hose motion with respect to new cable length 300 mm.

Validation

The variation simulation methods described in this article have been carefully validated by the research community for rigid and sheet metal assemblies. Crucial to the correctness of the analysis when extending the methods to deformable cables is the underlying simulation model. The Cosserat rod model has been verified against higher order finite element method (FEM) simulation as well as against real measurements provided by Delphi. Figure 8 shows the deviation of computer tomography (CT) scans of three different specimens and their average centerline compared to the result of the simulation model. The maximum deviation from the average center curve of the specimens is approximately 1 mm and for the individual specimens at most 3 mm, which indicates a very good alignment between simulation and measurement data.

Figure 8.

Comparison of simulation model and CT scan data (courtesy of Delphi). The plot shows the deviation between the simulation and scans of three different material specimens, and the pictures show the setup from two different views (the bright shape is the simulated shape and the dark shape is the scans). The maximum deviation from the average centerline of the specimens is very small: approximately 1 mm.

Conclusion

This article has presented a set of simulation methods for analyzing and controlling the propagation of geometric variation in deformable cables and hoses. Efficient implementations of a cable simulation model, similarity measures for triangulated geometries and tolerance envelopes have been combined with state-of-the-art variation simulation methods in order to strengthen early stages of the geometry assurance process. The methodology was successfully applied on two industrial cases.

The main scientific contribution of this work is the extension of existing simulation methods and tolerance envelopes to handle systems of deformable slender objects where large deformations need to be accounted for.