Home | ADCATS Info | Search | Site Map | Bulletin Board | Reports & Publications | Bibliography | Contact Us

AutoCAD 2D Analyzer:
CHAPTER 2: Analysis
Home : Manuals : AutoCAD - Analyzer - Analysis

2.0 Introduction

The tolerancing of machine parts is a necessary design function due to the variability of all manufacturing methods. Tolerancing is the assignment of upper and lower limits to basic dimensions to account for this variability.

The basis for the assignment of component tolerances should be the design requirements on critical assembly tolerances which assure proper performance. An assembly clearance or gap is an example of a performance parameter. If the clearance is too small, the parts may be difficult to assemble. If the clearance is too large, the parts may experience impact loads or misalign during operation. Such design parameters may be set by ANSI standards, trade group guidelines, government regulations, company design manuals, prototype testing or customer feedback.

In an assembly, the nominal dimensions of the component parts may be summed to determine the basic size of each clearance. However, the component tolerances may also accumulate, causing the clearance of some assemblies to exceed the design limits and thus fail to pass inspection. Hence, the individual component tolerances of the assembly must be chosen such that their sum does not exceed the critical assembly tolerance limits.

2.1 2-D and 3-D Tolerance Analysis

Tolerance accumulation or stack-up may be estimated from one of the expressions in Table 2.1. All three of these approximations are available in CATS. Which one you use depends upon customer requirements, process data available, and desired accuracy.

Table 2.1. Assembly Tolerance Accumulation Models

Worst Case
Assures 100% assembly acceptance if all parts meet component tolerances. Costly design model. Requires excessively tight component tolerances.
Root Sum Square
Assumes Normal distribution and +/-3Sigma tolerances. Some fraction of assemblies will not meet specification. May adjust ZASM to obtain desired acceptance fraction. Less costly. Permits looser component tolerances.
Six Sigma
Most realistic estimates. Accounts for process mean shifts and their long-term effects on assembly distribution.

In Table 2.1, dU is the predicted variation in the resultant assembly dimension; dxi is the variation in a component dimension; U/xi is the sensitivity that a variation in dxi has on U; and TASM is the assembly tolerance (or design limit) for variations in dU. ZASM is a target number of standard deviations at which the designer would like the assembly spec limits to lie. Cpi is the process capability ratio for a component dimension whose manufacturing process is assumed to be centered. Cpki is the process capability index for a component dimension, which includes a measure of long term drifts in the process mean.

CATS calculates the sensitivities from the AutoCATS 2-D vector model and predicts the tolerance accumulation of the assembly variables of interest. If specified limits have been set for an assembly variable, the computed distribution of the variable is used to calculate the number of assemblies which will be out of spec. This is shown graphically for a ZASM of 3.0 in Figure 2.1.

Figure 2.1. Determining the number of out-of-spec assemblies.

2.2 A Vector Model of an Assembly

Figure 2.2 shows a kinematic vector model overlaid on the cross section of a mechanical assembly. This device is a locking hub for holding a reel of computer tape on a tape drive. In this example, the plunger slides vertically downward while the wedge on its outer diameter pushes the arm outward until it locks against the inside diameter of a computer tape reel. In the extreme position shown, with the plunger resting on the base, we would like to assure an interference fit between the rubber pad and the reel (shown here as a gap or clearance for clarity).

Figure 2.2. Vector loop model of locking computer tape hub assembly.

Two loops describe the variation in this assembly, one closed and one open, as shown in Fig. 2.2. The open loop is one-dimensional and is set up to compute the variable Gap between the Arm and Reel. The closed loop is two-dimensional and is set up to compute RL, the resultant outer radius of the tape hub. The nominal value of the Gap and RL, along with dependent variables u and F, may be determined by an accurate CAD layout. However, the variations in the Gap, RL, u and F are desired. The resultant variation may be computed in terms of the independent variables and their sensitivities using a worst case or statistical model for tolerance accumulation.

Each 2-D vector loop equation may be resolved into three scalar equations, representing the sum of vector components in the X and Y directions and the sum of the relative rotations between adjacent vectors. The resulting system of equations may be expessed:

hx(xi,uj) = 0

hy(xi,uj) = 0

htheta(xi,uj) = 0

where xi are the manufactured component dimensions, uj are the unknown kinematic lengths and angles.

For the locking tape hub problem, the single closed loop shown in figure 2.2 gives the following set of nonlinear equations:

hx = b cos(90) + a cos(0) + u cos(theta) + r cos(theta - 90) + (e + i ) cos(theta + F - 90)

+ g cos(theta + F - 180) + RL cos(theta + F - 270) + h cos(theta + F - 180) = 0

hy = b sin(90) + a sin(0) + u sin(theta) + r sin(theta - 90) + (e + i ) sin(theta + F - 90)

+ g sin(theta + F - 180) + RL sin(theta + F - 270) + h sin(theta + F - 180) = 0

htheta = 90 - 90 + theta - 90 + F + 0 - 90 - 90 + 90 + 90 = 0

where a,b,e,g,h,i,r and theta are the independent manufactured dimensions and RL,u and F are the dependent kinematic dimensions.

The open loop equation is a one-dimensional equation for the resultant gap between the tape reel radius RT and the hub radius in locked position RL.

Gap = RT - RL

2.3 Linearized Solution for the Assembly Tolerances

There are two steps in solving the vector equations. First, all of the manufactured dimensions are set to their nominal values and the system of vector loop equations are solved for the nominal values of the kinematic variables and assembly resultants. Vector equations are generally nonlinear and must be solved by iterative methods. If the nominal values may be determined from a precise CAD layout, this step may be omitted.

Second, the equations are linearized for small variations about the nominal by Taylor's series expansion, retaining the first order derivatives. Each derivative is evaluated using the nominal dimensions of the dependent and independent variables. The component tolerances dxi are substituted into the equations and solve this system of linearized equations is solved for the corresponding variation in the kinematic variables and resultant assembly dimensions.

The linearized closed loop equations may be written as:

dhx(xi,uj,k) = Sigma(F(hx,xi)dxi) + Sigma(F(hx,uj)duj) + Sigma(F(hx,k)dk) = 0

dhy(xi,uj,k) = Sigma(F(hy,xi)dxi) + Sigma(F(hy,uj)duj) + Sigma(F(hy,k)dk) = 0

dhtheta(xi,uj,k) = Sigma(F(htheta,xi)dxi) + Sigma(F(htheta,uj)duj) + Sigma(F(htheta,k)dk) = 0

where dxi are the specified tolerances on the independent dimensions and duj are the resultant variations in the dependent assembly dimensions. The equations have also been expanded to include form and feature tolerances dk. The affects of form tolerances will be discussed in a later section. In matrix notation these equations become:

[A] {dx} + [B] {du} + [F] {d} = {0} Closed Loop Equations

where [A] is the matrix of sensitivities to variations in the independent variables xi

{dx} is the vector of specified variations of the independent variables dxi

[B] is the matrix of sensitivities to variations in the independent variables uj

{du} is the vector of unknown variations of the dependent variables duj

[F] is the matrix of sensitivities to variations in the form variables k

{d} is the vector of specified variations of the form variables dk

This system of equations may be solved by linear algebra for the unknowns duj:

{du} = -[B-1] [A] {dx} - [B-1] [F] {d}

from which expressions for the predicted accumulation of variations in the duj may be generated:

duj = Sigma(|F(uj,xi)|dxi) + Sigma(|F(uj,k)|dk) (Worst Case)

duj = (Sigma(F(uj,xi)dxi)2 + Sigma(F(uj,k)dk)½ (Statistical)

The linearized open loop equations in matrix form may be expressed:

{dg} = [C] {dx} + [E] {du} + [G] {d} Open Loop Equations

where C, E and G represent sensitivities to variations in dxi, duj and d[]k, respectively. dg is the non-zero vector of specified assembly gaps or rotational variations.

In general, the open loop equations do not require simultaneous solution. They are evaluated after eliminating the duj by substituting expressions from the closed loop equations.

{dg} = [C] {dx} - [E] ([B-1] [A] {dx} + [B-1] [F] {d[]}) + [G] {d[]}

= ([C] - [E] [B-1] [A] ) {dx} + ([G] - [E] [B-1] [F] ) {d[]}

from which expressions for the predicted accumulation of variations in the dgj may be generated:

dgj = Sigma(|F(gj,xi)|dxi) + Sigma(|F(gj,k)|dk) (Worst Case)

dgj = R(Sigma(F(gj,xi)dxi)2 + Sigma(F(gj,k)dk)2 (Statistical)

In the statistical expressions above, if the dxi and the dk represent three standard deviations of their respective statistical distributions, then the duj and dgj will also correspond to three standard deviations. These equations may be used to estimate the variance of the resultant assembly distributions. The resulting distribution may be used to compute the percent rejects as shown in Fig. 3.1.

Thus, the variational behavior of the whole population of assemblies may be determined by analyzing just one assembly statistically, rather than analyzing a large sample population one at a time, as is done in Monte Carlo Simulation.

2.4 Numerical Example

Nominal dimensions and tolerances of the tape hub lock :

Independent Dependent

a = 1.355 +/-0.0015 g = 0.493 +/-0.004 u = 0.319 +/- ?

b = 0.400 +/-0.006 h = 0.200 +/-0.008 RL = 1.864 +/- ?

r = 0.060 +/-0.002 theta = 75 +/-0.5 (deg) F = 15deg. +/- ?

e = 0.318 +/-0.003 = 0.2618 +/-0.0087 (rad) SPEC

i = 0.050 +/-0.002 RT = 1.856 +/-0.004 Gap = -0.010 +/- 0.006

Taking the partial derivatives of the three hub lock equations yields the following sensitivities (evaluated at the nominal dimensions):

The bold symbols above each column indicate with respect to which variable the derivative has been taken.

Solving for du in terms of the independent variables dx (ignoring form variations for now):

The solution for the dependent assembly variations in terms of the component tolerances:

{du dRL dF}T =

{db da dr de di dg dh dtheta}T

To estimate dRL, multiply the second row of the B-1A matrix times the vector of component tolerances, taking absolute values for worst case or root sum square for statistical.

dRL = |0.268 db| + |1.0 da| + |1.035 dr| + |1.0 (de+di)| +|0.268 (dg+dh)| + |0.3307 dtheta|

= 0.01548 (Worst Case)

dRL = (0.2682db2+1.0da2+1.0352dr+1.0(de2+di2)+0.2682(dg2+dh2)+0.33072dtheta2)½

= 0.00578 (Statistical)

Clearly, worst case is a much more conservative estimate of the variation.

By substituting trial values for the tolerances dxi, these expressions may be used for "what-if" studies and tolerance design. On the other hand, by substituting the design limits for the duj, tolerance allocation algorithms may be used to determine acceptable component tolerances for dxi.

To predict the variation in the resulting gap between the tape reel and hub, we must evaluate the open loop equation. Calculating the nominal gap:

Gap = RT - RL = 1.856 - 1.864 = -0.008 (interference)

dGap = dRT - dRL

The corresponding expressions for worst case and statistical variation:

dGap = |0.004| + |0.01548| = +/-0.0195 (Worst Case)

dGap = R(0.0042 + 0.005782 ) = +/-0.0070 (Statistical)

Calculated Gap limits: Min = -0.015 Max = -0.001 Mean = -0.008

Specified Gap limits: Min = -0.016 Max = -0.004 Mean = -0.010

% Rejects: Lower Limit: 0.03% Upper Limit: 4.36%

Thus, we see that the component tolerances as specified will yield gap limits which are too wide and the mean value will not be centered in the specified range. An unacceptable number of assemblies will have insufficient grip to lock a tape reel securely. Adjustments will be needed in the nominals to center the mean properly and some component tolerances will have to be tightened.

Examining the sensitivities and percent contribution of each component will help the designer to decide which tolerances to change. Chapters 5, 6, and 7 will use three separate example problems to show how to use the CATS tolerance design options to accomplish this automatically.

2.5 Approximating Form Variations

Form variations only introduce variation into an assembly at the points of contact between mating surfaces. Figure 2.3 shows a cylinder in contact with a plane. In one assembly, the cylinder might rest on a peak of the surface. In the next assembly, the cylinder could be down in a valley. Thus, the surface waviness in the plane results in a translational variation normal to the surface.

Similarly, the cylinder is probably not perfectly round, but exhibits waviness or lobing. If, during assembly, it is placed with a lobe at the contact point, the cylinder center will be higher than average. If a low point on the surface of the cylinder is in contact, the center will be lower. Thus, we see that both surfaces produce independent translational variations which are normal to the surface at the point of contact.

Figure 2.3
Translational variations

Figure 2.4 Rotational variationscaused by form variations. caused by form variations.

In contrast, a block on a plane produces rotational variation, as shown in Fig. 2.4, since one corner of the block may be higher than the other. The magnitude of the variation depends on the length of the block and the amplitude and period of the waviness.

2.6 Calculating the Contribution of Form Variations

Figure 2.5 shows the example assembly with feature tolerances defined at the mating surfaces. Each tolerance is specified by an ANSI Y14.5 feature control symbol which is added to the assembly drawing. The []i symbols identify the magnitude of feature tolerance as listed in Table 2.2. Corresponding algebraic terms are added to the kinematic assembly model and serve as independent sources of variation.

Figure 2.5. Feature tolerances added at mating surfaces of an assembly.

Table 2.2 Specified Feature Tolerances for Tape Hub Problem

Feature Feature Type Joint Type Variation Type Characteristic Length
1 perpendicularity rigid rotation 1.355
2 run out cylinder slider translation N/A
3 profile cylinder slider translation N/A
4 perpendicularity rectangular rotation 0.493
5 flatness planar rotation 0.359
6 flatness planar rotation 0.359
By introducing small translational or rotational variations at each point
of contact, the sensitivity matrix for form variations may be determined.

Solving for du:

{du} = -[B-1] [A] {dx} - [B-1] [F] {d}

Substituting for [A], [F] and [B-1], then taking the middle row of sensitivites to form the tolerance accumulation expression with form tolerances:

dRL =0.01548+|.694 d1|+|1.035 d2|+|1.035 d3|+|.397 d4|+|.397 d5|+|.397 d6|

= 0.02504 (Worst Case)

dRL =(.00578+0.2682d12+1.0d22+1.0352d32+1.0d42+0.2682d52+0.33072d62)½

= 0.00713 (Statistical)

The corresponding expressions for the gap:

dGap = |0.004| + |0.02504| = +/-0.02904 (Worst Case)

dGap = R(0.0042 + 0.007132 ) = +/-0.0082 (Statistical)

Calculated Gap limits: Min = -0.0162 Max = +0.0002 Mean = -0.008

Specified Gap limits: Min = -0.016 Max = -0.004 Mean = -0.010

% Rejects: Lower Limit: 0.13% Upper Limit: 7.08%

This represents nearly twice as many rejects as was predicted without form tolerances. Thus we see the necessity for including form variations in assembly tolerance analysis.

Please note that all of the calculations demonstrated here are performed automatically by the AutoCATS Analyzer. Once the engineering model has been created, the equations are derived, dependent variables separated, sensitivities calculated, and expressions for tolerance accumulation computed.

2.7 Advanced Six Sigma Analysis--Accounting for Process Mean Shifts

Six Sigma is the quality program instituted corporate-wide at Motorola, for which they received the first annual Baldridge National Quality award in 1988. Now, many of Motorola's suppliers and other major companies have adopted the same system.

The basic premise of the Six Sigma Program is: in order to achieve high quality in a complex product comprised of many components and processes, each component and process must be produced at significantly higher quality levels in order for the composite result to meet final quality standards.

Stated statistically, suppose there were 1000 dimensions or other characteristics of your product, any one of which could lower the quality of the finished product. If each characteristic were produced to +/-3Sigma quality (99.73% acceptable parts or 2700 defects per million), the resultant assemblies would be only (.9973)1000=.067 or 6.7% defect free. To have 99.73% defect free assemblies, you would need to produce each component to a quality of (.9973).001=.9999973 or 99.99973%, which is 2.7 defects per million.

To achieve the high quality levels required for world competition in the electronics industry, Motorola has mandated +/-6Sigma quality for all processes (.002 defects per million). However, they also recognize that shifts and drifts in the mean of the processes are expected, so they have introduced a modified process model which includes an allowance for accumulated mean shifts. The result is a net quality level of +/-4.5Sigma (3.4 defects per million).

Of course, tolerance analysis of assemblies is only one component of the complete Motorola quality management system, but the Six Sigma tolerance analysis model is so versatile and powerful that it can have a major impact on cutting production cost and improving quality.

Assigning tolerances based on statistical tolerance analysis alone is not sufficient to assure that reject rates will remain within quality limits. Naturally occurring shifts in the mean of a process can produce biased distributions which can result in increased assembly problems and a greater number of rejects than anticipated.

The process mean is difficult to control. Tool wear causes the mean to drift, creating a zone of uncertainty, as shown in Figure 2.6. Shifts can also occur due to setup error. Deliberate bias toward the least material condition may be introduced to account for tool wear. Bias toward the maximum material condition may be used to allow material for rework.

In an assembly, mean shifts can also accumulate just as tolerances do. The resulting assembly distribution can be pushed so far off center that significant rejects may occur even though all component parts are produced within tolerance specs.

Even when the mean shifts do not accumulate or bias the component distributions, they can still cause an increase in rejects. For large quantity production, the mean may drift around within the zone of uncertainty, causing the distribution to spread. The long-term distribution of the process will be greater than that determined by typical sampling procedures.

Figure 2.6. The location of the mean is not known precisely.


To account for all these effects, the designer needs to use a more complete model. CATS has incorporated the Motorola model for tolerance accumulation, which includes two mean shift factors, kstat and kdyn, to account for static and dynamic mean shifts, respectively. Static mean shifts describe non-centered processes. Dynamic mean shifts describe the long term increase in the spread of the distribution due to drifting about the center. Each shift factor expresses the mean shift as a fraction of the original component tolerance.

The new component means become:

xi' = xi + kstati Ti

The long term tolerance accumulation model may be expressed as:

dU = ZasmR(Sigma(f(Ti,3Cpki))2) <= Tasm

where, xi is the nominal component dimension, dU is the resultant assembly variation, Ti is the component tolerance. Zasm is the desired number of standard deviations of the assembly distribution to be included in the specified assembly tolerance Tasm. Set Zasm = 3.0 for the usual assumption of +/-3Sigma assembly tolerance limits. Cpki is the Process Capability Index used by manufacturing to quantify non-centered processes. Set Cpki = 1.0 for the usual assumption of +/-3Sigma component tolerance limits. Each component Cpk may be calculated from the expressions:

Cpk = Cp(1-kdyn)

Cp = f(USL - LSL,6Sigma)

Cp is the Process Capability Ratio. It is the ratio of the width of the tolerance zone (USL - LSL) to the 6Sigma width of the process. When the tolerance limits are +/-3Sigma, (USL - LSL) = 6Sigma and Cp = 1.0. When Cp = 1.0 and kdyn = 0, then Cpki = 1.0 and the component tolerances are +/-3Sigma and centered. If the process limits are set at +/-6Sigma and kdyn = 0.25, then the mean could be shifted as much as one fourth of the tolerance zone, or 1.5Sigma, as shown in Figure 2.7.

Figure 2.7. Effect of a mean shift on +/-6Sigma specification limits for all component kstat=0.25.

The Six Sigma model can provide a common ground for interaction between engineering and manufacturing. The model forces the engineer's attention on manufacturing considerations, and he must communicate with manufacturing to get the needed model data. He can convey the critical design parameters to manufacturing in a form that permits freedom to alter tolerances without violating design requirements. Manufacturing can communicate meaningfully to engineering in terms of the two quality assurance parameters most commonly used in statistical process control: mean and variance.

Summary of Six Sigma Advanced Analysis

1. The CATS Motorola Six Sigma model can perform statistical tolerance analysis of assemblies for any specified quality level. Individual quality levels for each component may also be specified in terms of the Process Capability Ratio (Cp) or the standard deviation of the process.

2. The mean shift for each component process may be specified and the resultant effect on the assembly yield analyzed. Two kinds of mean shift may be specified: static or dynamic. Static mean shift results from fixed bias in a process, such as setup or tooling errors, and is added to the mean of the component dimension. Dynamic mean shift is a random bias due to long-term variations in the process mean, such as tool wear or setup variations over numerous production runs. It causes an increase in the variance or the spread in the process distribution and a decrease in Cpk.

 PRO-E 2 D
Title | Overview | Modeling | Commands
Title | Overview | Analysis | Allocation | Interface
Verification: Overview

 AutoCAD 2 D
Title | Overview | Modeling | Commands
Analyzer: Title | Overview | Analysis | Allocation | Interface
Verification: Title | Overview

 Catia 3 D
Title | Overview | Modeling | Commands | Building a Tolerance Model

The ADCATS site: Home | ADCATS Info | Search | Site Map | Bulletin Board | Reports & Publications | Bibliography | Contact Us