The present invention relates to a method for producing the drilled holes with accurate geometrical and positional parameters of an air-cooled components, in particular for the production of film cooling holes with accurate geometrical and positional values in a blade of a gas or steam turbine. Preferably in an aero-engine hollow turbine blade with multiple drilled holes introduced to ensure proper cooling air flow through the cooling holes.
To achieve high levels of efficiency with a high output at the same time, gas turbines are operated with an extremely inlet temperature, which often exceeds the permissible material temperature of turbine blading. It is essential to strengthen the bearing working temperature ability of the turbine blade, and to use certain cooling strategies to make it possible nevertheless for the blades to be operated reliably and with long life-time. Development in cooling technique resulted in complicated cooling scheme, multiple (104-106) fine-shaped holes (0.3-1.2 mm) with divergent oblique angles (10°-75°) exist on such air-cooled components in which cooling air flows through the airfoil and discharged through cooling holes. Consequently, the flow of the cooling fluid as a film is determined considerably by the shape and location of the hole, which means the cooling holes need to be straight, accurate, and exactly positioned. For example, diametrical tolerances for air-cooled holes are generally on the order of approximately ±0.002 inch (±0.05 mm) or less.
A wide variety of methods in the state-of-art for producing drilled cooling holes for gas turbine blades, for instance, or in other embodiments, precision drilling techniques such as the laser beam machining, the electrical discharge machining (EDM), and the electrical discharge machining (ECM), etc., are typically used to producing the holes. Specifically, ultrashort pulsed lasers for developing in cooling holes drilling applications with minimized recast layer and negligible heat-affected zone are widely used.
Although the methods for cooling holes drilling known from the state of art have a significant performance on drilling the holes with high efficiency, trial-and-error methods are still implemented in drilling operations. Firstly it is because the actual wall thickness at the designed location of drilled holes deviates from the ideal wall thickness due to the that the components such as gas turbine blades with intricate geometry and tight tolerances are often fabricated by investment casting, which is a multiplex nonlinear physical processes coupled with geometry, material properties and boundary conditions. The component, in the embodiment of gas turbine blades, the dimensions produced exhibit non-uniform deformation due to the deformation of wax pattern, expansion of shell, and alloy shrinkage. In addition, holding a thin-walled components such as blades with sculptured surfaces firmly depends on the fixture. Once the turbine blade is localized thoroughly, and restrained by the locators and clamp of the fixture, the component is processed to generate drilled holes with geometric features according to the cooling requirements, In general, the hole feature may have deviation in terms of the localization error. Furthermore, the deviation may occur when clamping the component in the machining process, including the movement and deformation of the component, etc.
As is apparent from the deviations of the ideal drilled holes features discussed above regarding the geometries and positions of drilled hole, the defects such as the occurrence of blind holes, laser burns on the opposite wall, as well as the lack of consistency accuracy of holes. This causes a change in the through-flow of the cooling fluid in comparison with the intended through-flow.
Document U.S. Pat. No. 0,229,759A1 describes a method for the parametric production of a drilled hole in a component, particularly in a gas turbine blade, comprising the steps of: Measuring an actual wall thickness of the component at the location of drill hole; Adjusting the wall thickness on the basis of adjusted parametric drilled hole geometrical value according to measurement. As such, the actual wall thickness of drilled holes can be determined based on the corrected value.
Document U.S. Pat. No. 0,183,325 A1 discloses a method for producing holes in a component, in particular of turbo-machines. The method utilizes 3D model of the actual geometry of the component, adopting each hole on the basis of actual geometry, then the production program is generated for each individual hole. As a result, the process quality and the quality of the holes can be improved.
Document U.S. Pat. No. 6,339,879 B1 teaches a method of sizing and forming a cooling hole in a gas turbine engine component. The method generally entails drilling a hole in the surface of a component, then measuring the thickness of the recast layer surrounding the hole. The thickness of the additive layer of the diffusion coating that will deposit on a corresponding recast surface of the component can be predicted based on an inverse relationship determined by the measured thickness of the recast layer during drilling process. As a result, an appropriately-oversized hole can then be formed in the component.
Although these documents teach various methods for producing drilled hole or drilled holes in a gas turbine blade, or the like, there is no intent to establish, or knowledge about the method for producing drilled holes parametrically in a gas turbine engine component with accurate geometrical and positional values with the consideration of the deviations generated before and after machining.
In view of the above, it would be desirable if an improved process were available for producing the cooling holes of gas turbine engine components with the accurate geometrical and positional parameters.
Accordingly, It is therefore an object of the invention to solve the problem mentioned above. To provide a method for the production of a drilled hole in a component parametrically, it being intended that the drilled hole produced should correspond to predetermined geometries and positions to reach the achieve pre-established cooling airflow conditions. In particular, this method is suitable for producing an air-cooled drilled cooling hole in a component of a gas or steam turbine, preferably in a blade of a gas or steam turbine in an accurate manner.
The present invention provides a method entails producing the drilled holes in a surface region of a substrate. In accordance with this method, a component is formed by investment casting process, determining the surface error of the component generated during the process, modeling the deformation characteristics of the surface error, including shrinkage, bending, and torsion; then fixing the component by fixture system, calculating the deviation resulted from the transformation of the reference frame of the component, and calculating the deviation caused by clamping deformation. During the component drilling process, predicting the deformation due to the drilling force. As such, the surface deviation, comprising the component manufacturing error, locating positional error, clamping deformation, drilling error can be determined. As a result, the geometrical and positional parameters describing actual values of drilled holes can be determined parametrically. Knowing the actual parameters of drilled holes, an appropriately cooling performance can be achieved in the component accurately.
The method according to the invention is suitable in particular for producing drilled holes in a component of gas turbine parametrically, preferably in a blade of a gas or steam turbine. Furthermore, the process of this invention is useful for newly manufactured gas turbine blade and presumably for repaired ones, wherein the drilled holes required accurate geometries and positions. As a result, a notable advantage of the present invention is that it succeeds in maintaining required cooling requirement while significantly reducing the necessary labor and time to producing cooling holes due to trial-and-error method widely adopted nowadays.
These and other features, aspects, and advantages of the present invention will become better understood when the following detailed description is read with reference to the accompanying drawings in which like characters represent like parts throughout the drawings, wherein:
The present invention will be described in reference to a high pressure gas turbine blade [1] illustrated in
According to this invention, the airfoil section [2] with cooling holes [5] through a wall [6] of the high pressure gas turbine blade [1] section is shown in
In the case of a drilled cooling holes [5] that is to be produced, the cooling holes [5] comprises in the direction of flow [7] of the cooling fluid an supply section [8], and a diffuser section (not shown) adjoining the supply section [8]. For noncircular cooling hole, the supply section [8] and diffuser section have different cross-sections to obtain cooling efficiency improvement. Since the circular drilled hole cross-sections can be produced more easily from a production engineering viewpoint than noncircular drilled hole cross-sections, moreover, the noncircular cooling hole is the circular cooling hole evolution of the through-flow characteristics at the inlet region, which make the circular holes the most representative types of the drilled holes. In one embodiment, the supply section [8] and the diffuser section may also, have a constant cross-section. In an effort to provide a concise description of the embodiment, the embodiments for producing noncircular cooling holes may not be described in the specification. It is contemplated that a variety of holes with different shapes and sizes can be accurately formed using the processes disclosed herein.
Here in
However, the geometrical and positional values of the cooling holes [5] deviate from the ideal values, deviations of the actual cooling holes [5] from the ideal values hole also occur and, as a consequence, in some cases considerable changes of the through-flow characteristics of the drilled cooling hole occur, even in the case of deviations are within the production tolerances. For example, the length of the supply section [8] which represents the depth of the hole L given here as a specified value may be changed. The occurrence of the deviations are owing to two obvious seasons: Firstly, due to the high pressure gas turbine blade [1] with intricate geometry and tight tolerances are often fabricated by investment casting, which is a multiplex nonlinear physical processes coupled with geometry, material properties and boundary conditions. Therefore, the turbine blade dimensions produced exhibit non-uniform deformation due to the deformation of wax pattern, expansion of shell, and shrinkage of super-alloy. Secondly, holding a thin-walled high pressure gas turbine blade [1] with sculptured surfaces firmly depends on the manufacturing fixture. Fixture system is used to establish the datum reference frame of the surface of turbine blade to be processed with respect to the reference frame of a drilling tool. Once the high pressure gas turbine blade [1] is localized thoroughly, and restrained by the locators and clamp of the fixture, the high pressure gas turbine blade [1] is processed to generate cooling holes [5] with geometric features according to the cooling requirements. In general, the cooling holes [5] feature may have geometric errors in terms of its form and position related to the component datum reference frame. If there exists a misalignment error between the component datum reference frame and drilling tool reference frame, also known as localization error. A localization error is essentially caused by a deviation in the position of the contact point between a locator and the component surface from its nominal specification, and it plays an important role in the deviation of the airfoil section [2] besides the geometrical and positional values of the cooling holes [5].
In general, the deviation of geometrical and positional values of the cooling holes [5] generated before and during the machining process.
E=E
b
+E
a=(ei+ec+el)+(em+ed) (1)
To make allowance during the production of drilled holes for deviations of the geometry and positions actually produced from hole values ideally specified, the invention, need to determine each deviation as it is depicted in
For the determination of manufacturing error [15] ei, the manufacturing error [15] often exhibits non-uniform distribution, this is primarily due to the nonhomogeneous stress induced by different cooling rate generated by complex structures. Hence, to establish the geometrical and positional parameters of drilled hole, including the open angle of the diffuser [12] to surface tangent, the nominal thickness [9], and the position of the hole center [11], the determination of manufacturing error [15] during investment casting becomes essential.
The manufacturing error [15] depends on the integrated operation of many factors. It is not only closely related to the properties of materials having unstable compositions but also concerned with the casting structure and forming conditions. The non-uniform manufacturing error [15], also known as the displacement field, includes different deformation characteristics, consisting of shrinkage, torsion, and bending. For any point on the surface of the casting, the displacement filed Di can be described as:
D
i
=D
si
+D
bi
+D
ti (2)
where Dsi, Dbi, and Dti represent the deformation characteristics of shrinkage, bending, and torsion, respectively. Since the accurate geometrical and positional parameters of cooling holes [5] significantly depend on the displacement between actual and ideal models of the high pressure gas turbine blade [1], it is essential to modeling each deformation characteristic coupled in the displacement field.
The method for determining the torsional deformation Dti is described as: A series of the 2D section [21] of the high pressure gas turbine blade [1] are obtained from the ideal and actual models, illustrated in
The determination of geometrical centroid is based on the discrete points of each 2D section [21]. Given the set of discrete points of the 2D section [21] {p1, p2, . . . , pn}, let pi(xi, yi, zi) be the coordinates of the discrete points, i=1, 2, . . . , n, and let M be the geometrical centroid of pi, which can be expressed as:
To determine the deformation characteristics of shrinkage and bending Dsi, Dbi in the manufacturing error [15] em, Dti must be eliminated and compensated. The method for compensating torsional deformation characteristic is described as:
P
i
′=P
i
+D
ti (4)
where Pi and Pi′ represent any point on the surface of 2D section of actual model [26] and the 2D section of ideal model [25] without torsional deformation, respectively. Hence, Pc′=Pc+Dti, where Pc and Pc′ represent the geometrical centroid of actual section [22] and the geometrical centroid of ideal section [23] without torsional deformation, respectively. Theoretically, the geometrical centroid of the 2D section of actual model [26] and the 2D section of ideal model [25] coincide, Pc′=Mc, after compensation.
The bending deformation characteristic Dbi is represented by the angular displacement of the actual model relative to the ideal model, which is described by torsional angle, as illustrated in
P
i
″=P
i
′·R (5)
where R represents the rotation matrix, and θ is the torsional angle [27].
With the determination of torsional and bending deformation Dti and Dbi, the manufacturing error [15] em of high pressure gas turbine blade [1] only contains shrinkage deformation Dsi. The shrinkage deformation has important influence on the determination the nominal thickness [9] of the cooling holes [5], which can be discussed in the step of the correction of geometrical and positional parameters of cooling holes [5].
In the step, locating error el is determined. In the fixture model, a locator is used for providing unique and accurate location in the both aspects of position and rotation with respect to a fixture reference frame. For the localization of high pressure gas turbine blade [1] prepared for cooling holes [5] drilling, a 3-2-1 locating principle is applied with 6 locators. Define CSw as the component-embedded coordinate system, a locator i (i=1, 2, . . . , 6) contacted the component surface at position ri along the surface normal ni, if a small displacement of locating position δr generates the perturbations in position and rotation of the reference frame of component by δq, which CSw is transformed to CSw′. According to the kinematics principle of rigid body, the relationship between δr and δq is given as:
δr=Gδq (7)
where δr=[δr1, δr2, . . . δri, . . . δrn]T is the locating positional error, specifically, δri is the locating error of the locator i, which represents the distance from the nominal of the locator ni to the component surface, δq represents the coordinate transformation from CSw to CSwt.
Thus, δq can be expressed as a matrix consisting of displacement derivatives of translation and rotation:
δq=[δxw,δyw,δzw,δαw,δβw,δγw]T (8)
where [δxw, δyw, δzw] are the translation displacements of the CSw in the X, Y, and Z directions; [δαw, δβw, δγw] denote three independent angles representing the orientation of CSw relative to CSw′, G is the Jacobian matrix in three-dimensional space, which is calculated as:
Therefore, regarding the transformation δq, the displacement of the point on the Component surface, denoted as eli can be written as:
e
li
=T(δaw)·T(δβw)·(δγw)·R(ri)−Pi (10)
where R(ri)=[xi+δxw, yi+δyw, zi+δzw]T, Pi=[xi, yi, zi]T, T(δαw), T(δβw), and T(δγw) represent the 3×3 rotation matrix in terms of three angles.
Clamping deformation, ec, which is the deviation of component arising from localized elastic deformation at fixture-component contacts owing to clamping is known to affect significantly the component location accuracy. The finite-element method is employed to modeling the component-fixture contact mechanism and the surface error under the static clamping circumstances, the method for determining ec is demonstrated in
Although the following example refers to nickel-based high pressure gas turbine blade [1] and stainless steel fixture system [32], the same process would apply for other materials, such as metals and alloys of metals, as would be understood by those of skill in the art. Examples of suitable materials include, but are not limited to, aluminum, magnesium, nickel, and their alloys.
In the exemplary embodiment, the high pressure gas turbine blade [1] and the fixture system [32] are meshed by the available commercially software, such as Hypermesh from Altair Engineering, ANSYS Meshing from ANSYS Corp, etc., can be used for modeling the meshing model [31]. For the fixture system, the contact model [33] between high pressure gas turbine blade [1] and fixture system [32] needs to be determined. Locators and clamps of the contact model [33] are assumed to have the spherical tipped shape, the contact area between component and clamps is relatively small compared to the dimensions of component-fixture system, hence, the contact region can be presumed to have the lumped contact stiffness, the deformation at the contact point can be solved by the parameters of the component-fixture system, including the radius of the locator at the contact region, the Young's modulus and the Poisson's ratio of the component and fixture materials, in this embodiment, are the Young's modulus and the Poisson's ratio of the nickel-based super-alloy and stainless steel. The meshing model [31] of high pressure gas turbine blade [1], the contact model [33], combined with the applied clamping force [34], the static model [35] can be obtained, which is used for modeling the clamping deformation [16] ec.
When machining the cooling holes [5] on the sculptured surface of high pressure gas turbine blade [1] using mechanical drilling operation, the structural deformation under dynamic drilling forces could be a damageable factor to the final performance. The deviation Ea caused by machining is composed of the error caused by component movement em and component deformation ed, which can be regarded as the sum of the responses both of geometric error and structural deformation.
To determine the deviation Ea induced by machining process, finite-element analysis can be used. During this analysis, the high pressure gas turbine blade [1] is held firmly by the fixture system [32], the locators and clamps are modeled as rigid body with area contact, the relative deformation of the component-fixture is regard as zero due to the clamping deformation is analyzed independently.
Predicting the deformation of high pressure gas turbine blade [1] during the machining, the machining forces in a mechanical drilling operation include a torque force for tool rotation, a shear force at the cutting edge of the hole, and an axial force due to feeding. Since the forces in machining operation are time and position-dependent, the simulation is realized by a time-discrete modeling method. The machining process is divided into adequate time steps Δt, due to the short machining interval, each time step (ti, ti+1) is assumed to follow the static equilibrium condition as:
[K]r(ti)=f(ti)i=1,2, . . . ,T (11)
where K denotes the structural stiffness matrix, r(ti) is the displacement vector, f (ti) is the vector of instantaneous force, i is the ith time step.
The instantaneous force [F] can be written as:
where n represents the number of the total calculation points. In this embodiment, n is chosen as the node of the meshing model [31]. Therefore, the deviation during machining Ea represented by the displacement vectors can be determined using Eq. (11) by static solution based on the available commercial finite element software package, such as ABAQUS from Dassault Group, ANSYS from ANSYS Corp, etc.
To further address the method for determining the deviations during machining process, different machining approaches need to be considered. Different machining approach has obvious different performance of the caused deviation error, for the non-traditional machining (NTM) techniques, such as electrical discharge machining (EDM), electrical discharge machining (EDM) and laser drilling processing, one of the advantages of NTM techniques is that the mechanical energy is not necessarily used to provide material removal, the generated residual stress can be minimized, hence, the deviations induced by machining process based on NTM techniques is smaller compared with the deviation caused by mechanical machining. However, the method for determining the deviations described is not intended to be limiting of the machining approach used. It can be used for the determination of deviation during different machining approaches, including, but not limit to mechanical drilling, electrical discharge machining, electrical discharge machining, laser drilling, etc.
For the purpose of production of the cooling holes [5] for the deviations E=Eb+Ea determined actually produced from the ideal model of high pressure gas turbine blade [1], the geometrical and positional values of the cooling holes [5] can expediently be determined from the parametrical correlations, and transferred to the control system of the drilling apparatus.
To develop a parametric model for describing the cooling holes [5] with corrected values, three parameters need to determined: the hole center [11], Pc, the nominal thickness [9], T, and the open angle of the diffuser [12], α, which can be represented by the nominal of the direction of flow [7], n.
For the determination of hole center [11] Pc, let CSw be the global coordinate system fixed on high pressure gas turbine blade [1] without deformation. Pc=[xc, yc, zc]T denotes the ideal coordinate of the hole center [11], Pe1=[xe1, ye1, ze1]T is the coordinate of the hole center [11] after high pressure gas turbine blade [1] investment casting process finishes, on the basis of discussion on modeling of manufacturing error em of high pressure gas turbine blade [1], the displacement from the hole center [11] Pc to Pe1 is consist of the deformation characteristics of torsion and bending. Therefore, the coordinate transformation from Pc to Pe1 can be written as:
P
c
e1
=R
e1
·P
i
+T
e1(i=1,2, . . . ,n) (13)
where Re1 and Te1 represent the rotation matrix and translation matrix of the component, respectively, which can be determined by decoupling the torsional and bending deformation. i denotes the ith discrete point on the surface of component.
Deviations occur when the high pressure gas turbine blade [1] is fixed and clamped, CS′w is the component fixed coordinate system, due to locating and clamping, hole center [11] Pe1 changes to Pe2, which can be described as:
P
c
e2
=R
e2
·P
c
e1
+T
e2(i=1,2, . . . ,n) (14)
where Re2 and Te2 denote the rotation matrix and translation matrix of the component generated by locating and clamping, which can be determined by Eq. (10) and Eq. (11), respectively.
For the deviation caused by machining process, hole center [11] Pie2 will be changing to Pe3 which can be determined using finite-element method as illustrated in
The method for the determination of the direction of flow [7] n of the cooling holes [5] is described as: let the neighborhood points [36]: P1, P2, P3, and P4 delineate four discrete points located in the neighborhood of the hole center [11], M1, M2, M3, and M4 are the middle points [37] of the lines [38]
The actual nominal thickness [9] of the component is essential for determining the length of supply section [8], which is also the depth of the cooling holes [5]. To calculate the actual nominal thickness [9], the actual and ideal models of high pressure gas turbine blade [1] are transformed into a series of 2D cross-sectional airfoil section [2] by cutting planes perpendicular to the blade length direction, hence, the calculated nominal thickness [9] can be regarded as the distance of the intersection points between the intersecting line and curve. The method of calculating the nominal thickness [9] is described as: firstly, the outer contour of the airfoil section [2] Cω (μ) is dispersed into n points, for example, n=200, denoted by Pi (i=1, 2, . . . , n), the unit normal vector along the inner normal direction at each point is Ni (i=1, 2, . . . , n), where point Pi is the starting point and Ni is the direction vector. Line L can be constructed as
L
i
=P
i
+tN
i(i,1,2, . . . ,n) (16)
The parametric equation for constructing the inner profile of the airfoil section [2] is given by
r=r(μ) (17)
The intersection points of line Li and inner profiles Q1i and Q2i are determined by Eq. (16) and Eq. (17). By definition, Ti=|PiQi1| is the nominal thickness [9].
To modeling the hole depth, which is represented as supply section [8], on the basis of the establishment of the actual nominal thickness [9], the deviation of nominal thickness [9] can be determined by the shrinkage deformation Dsi. The non-uniform shrinkage can be represented by a scale coefficient K called shrinkage ratio, represented as:
where K is the shrinkage ratio, D is the dimension of the 2D section of ideal model [25] of the high pressure gas turbine blade [1], and M is the dimension of the 2D section of actual model [26] of the high pressure gas turbine blade [1]. Therefore, the actual depth of the drilled cooling holes [5] can be determined the actual nominal thickness [9] with the consideration of the direction of flow [7].
To make allowance during the production of the cooling holes [5] for deviations determined from the ideally specified, the invention, as it is depicted in the flow diagram schematically represented in
Having described the embodiment of the present invention in detail, it will be apparent that modifications and variations are possible without departing from the scope of the embodiment defined in the appended claims. More specifically, although some aspects of the described embodiment of the present invention are identified herein as preferred or particularly advantageous, it is contemplated that the embodiment of the present invention is not necessarily limited to these preferred aspects.
This non-provisional application claims the benefit under 35 U.S.C. § 119(e) to U.S. Provisional Application No. 62/404,877 filed in the United States on Oct. 6, 2016, the entire contents of which are hereby incorporated by reference.
Number | Date | Country | |
---|---|---|---|
62404877 | Oct 2016 | US |