Dimensional control strategy and products distortions identification

More documents

Recommendations

Info

properties [3].2.1.b Inductive approachIt consists in proposing models to identifydistortions coming from experiment, mainly frommetrology. So, we need to define the number anddistribution of points so that significant surfaces forthe optimization criterion would be measured withthe smallest possible loss of information [4]. Thenumber of points is subject to part size, tolerancespecifications and measurement uncertainty. Thetheoretical smallest sampling size can be used forideal geometric shape identification. But, if the formdeviation is unknown (that is our case), we have toincrease the sampling size so that estimated value ofthe measurement converges to the “true” value ofform error [5]. As for points distribution, for notcreating local weighting disturbing the optimizationmethod, it is important to dispose them uniformly.Next, we need an optimization criterion to resolvethis “data-overabundant” system. In our case ofwarped surfaces identification, a point-per-pointcomparison between measured and theoreticalgeometry is needed [6]. We chose the least squarescriterion because it gives a very stable result and isless sensitive to asperities effects [7].2.2 Developed method2.2.a Introducing C-ring test partThe C-ring type sample is currently used because itsgeometry allows us to amplify distortion which issuitable for the metrological analysis [8] [9]. We usea 100 mm long C-ring with a 16 mm wide opening.Outer and inner cylinders (respectively Ø70 and Ø45mm) are 11 mm off-centered. These dimensionsauthorize thermocouples without disturbing thermalflow; they minimize side effects and allow obtainingrequired cooling velocity [10].First, a metrological analysis occurs at threemanufacturing states: after machining, after stressrelieving and after high pressure gas quench. Then,we present a method to correlate C-ring measuredmesh with f.e.m. one in order to compare simulatedistortions field with measurement’s one.2.2.b Experiment approachMeasurement strategyIn order to characterise the distortions significanttypes, we define a fine mesh for the C-ring. We keepin mind that many measurement points may increasethe appearance of abnormal points but we coulddetect them via a graphic interface (see section 3).Points are evenly distributed on all surfaces (figure1), but not on the putting plane (bottom one).Outer cylinder1722 points1722 pointsInner cylinderUpper planeBottom plane110 pointsRight line38 pointsLeft lineFig. 1.Fine mesh for points probing21 points21 pointsData processing strategyOur mathematic approach is to optimize, for eachpoint i, the measure error ε i between theoreticalpoint T i and measured point M i , orthogonallyprojected onto theoretical normal N i (equation (1)).⎛ ε1⎞⎜ ⎟ with :r M⎜ ⎟r ε = εiεi= ( ΟΤi− ΟΜi).Νi(1)⎜ ⎟⎜ M ⎟ Νi= 1⎝εn ⎠In referenceframer( Ο,i, j)We also define the matrix M ph of optimizationphenomena (equation (2)). Its scalars correspond tothe phenomena elementary amplitude effect (e Phnj )on the theoretical points written in column. Thedescription unit of phenomena would be close totheir estimated value, so the millimetre.ΜphΤxΤyK PhΟΤ1eTxeTyK ePh1 1=ΟΤMjeMTxeMej TyK K Mj PhΟΤM eMeMKKeMnTxnTynPhnn1n jnnwith Τx=eeTTx jy j= Τx.Ν= Τy.ΝjjΤy= 1mm= Ν . x= Ν . yTo dissociate phenomena, those must bemathematically independent and linearlysuperimposable. In metrology, the size order of thedefects is often small compared with nominaldimensions of the measured part and so, we will usethe small displacements assumption.Finally, equation (3) gives the residual r to minimizeusing the least square criteria thanks to any solver.jj(2)
= ε −a. Ph1− b.Ph2− c.Phi−...− n.Phnwith elementary quantities (a, b, c,..., n) of significantdisplacements and distortions phenomena ( Phi,..., Phn)(3)Experimental deviceOur coordinate measuring machine (CMM) has anindication error for length measurement defined byE=(3.5+L/350)µm, L in mm.Extension barTP20 probeStar stylusC-ringFig. 2. Experimental control devicePositioning deviceWe use a dynamic probe (TP20) with a star stylusand an extension bar (figure 2). Combined withprobe head rotations, we can access all surfaces. But,even with rotations, we do not obtain a uniformdistribution of points for the inner cylinder. We thuschecked the accuracy of our device compared with asingle TP20 probe (without extension bar and starstylus) in measuring a 44.99mm standard ring(table1). To succeed, we optimize the 108 points ofCMM with Tx-Ty and dilatation vectors.Table1. CMM measurements of a standard ring diameterValues in mm Single TP20 Our TP20 deviceDiameter 44.989 44.989Range 0.006 0.0052.2.c Processing of finite element simulationWe use Forge 2005 software to make 3D-quenchsimulations for which the input data are partlyexperimentally determined (the phase transformationkinetics, the heat transfer coefficient) and partlyfrom the literature (the mechanical data) [12].Data processing strategyFor a direct comparison between experimental andsimulation results, we use the method previouslypresented. For instance, this method is applied foronly the cylinders. By using “2D-cutting planes”, weexport sections’ boundary geometry at the sameheights as for CMM measured C-ring. However, wemust keep in mind that CMM measurement pointsare obtained by probing the surface following thetheoretical normals. But during simulation, nodesmove on all directions. So, we have to find theintersection point S 12 between theoretical normals T 1and finite element frontier (S 1 S 2 ).Fig. 3. Intersection point S 12 between (S 1 S 2 ) and D linesUncertainty due to the methodThe method generates a quantifiable error ε s due tothe linear approximation between two consecutivenodes of the finite element frontier (figure 3). Tominimise this error, we look for the 2 simulationpoints (S 1 and S 2 ) whose angular coordinates are asclose as possible to theoretical point (T 1 ) ones. Infigure 4, we show that higher is the sampling size;higher is the accuracy. This last also increase whenthe diameter size is smaller for the same samplingsize.Error in mm betweentheoretical andidentifyed diameters0,0400,0300,0200,0100,000Evolution of diameter identification error Evolution of diameter identification e0,040Outer section0,030Inner section0,0200,0100,00082 91 121 181 362 723Number of points from a cylinder section82 91 121 181 362 723Number of points from a cylinder sectionCase 1: Simulation nodes are theoretical Case 2: Simulation nodes are theoreticalpoints of outer and inner diameters points of dilated diameters (0.5mm)Fig. 4. Uncertainty due to data processing of simulation pointsInitial geometry for quench simulationWe must take into account C-ring dilation due toheating (at 930 °C). We consider for the studiedsteel grade a linear dilatation value of 22.3 µm/m/°Cin austenitic phase. Then, we create the same meshas for experimental analysis, consisting of 98sections with 400 nodes for each cylinder’sperimeter. This fine mesh aims to minimize theuncertainty of our method.3 FIRST RESULTSA software based on the described method wasdeveloped. Thus, numerical and 3D-identification ofdistortions phenomena can be calculated. For abetter visual understanding, a scale factor is appliedon the errors onto the theoretical normals (figure 5).Theses measure errors include displacementsphenomena, like translations, rotations and offcentringof the inner diameter (due to machining).By deleting them, we did what is often called best
Page 1: Dimensional control strategy and pr

Dimensional control strategy and products distortions identification

You also want an ePaper? Increase the reach of your titles

Delete template?

Save as template?