倾斜铸造的动态建模过程和模具设计对铸件质量的影响大学毕业论文外文文献翻译及原文

1、毕业设计（论文） 外文文献翻译文献、资料中文题目:倾斜铸造的动态建模过程和模具设计对铸件质量的影响文献、资料英文题目：文献、资料来源： 文献、资料发表（出版）日期：院（部）：专 业：班 级:姓 名：学 号：指导教师: 翻译日期:2017. 02. 14课程题目先进铸造技术文献题目：modelling the dynamics of the tilt-casting process and the effect of the mould design on the casting qualityh. wang'g. djambazow, k.a. pericleous r.a. hard

2、ings m. wickinsbcentre for numerical modelling and process analysis, university of greenwich, london se10 9ls, uk birc in materials processing, university of birmingham, birmingham, b15 2tt, uabstractall titanium alloys are highly reactive in the molten condition and so are usually melted in a water

3、-cooled copper crucible to avoid contamination using processes such as induction skull melting (ism). these provide only limited superheat which，coupled with the surface turbulence inherent in most conventional mould filling processes，results in entrainment defects such as bubbles in the castings. t

4、o overcome these problems，a novel tilt-casting process has been developed in which the mould is attached directly to the ism crucible holding the melt and the two are then rotated together to achieve a tranquil transfer of the metal into the mould. from the modelling point of view，this process invol

5、ves complex three-phase flow, heat transfer and solidification. in this paper, the development of a numerical model of the tilt-casting process is presented featuring several novel algorithm developments introduced into a general cfd package (physica) to model the complex dynamic interaction of the

6、liquid metal and melting atmosphere. these developments relate to the front tracking and heat transfer representations and to a casting-specific adaptation of the turbulence model to account for an advancing solid front. calculations have been performed for a 0.4 m long turbine blade cast in a titan

7、ium aluminidc alloy using different mould designs. it is shown that the feeder/basin configuration has a crucial influence on the casting quality. the computational results are validated against actual castings and are used to support an experimental programme. although fluid flow and heat transfer

8、are inseparable in a casting，the emphasis in this paper will be on the fluid dynamics of mould filling and its influence on cast quality rather than heat transfer and solidification which has been reported elsewhere.keywordstilt-casting; mould design; 3-d computational model; casting process；1. intr

9、oductionthe casting process is already many centuries old，yet many researchers are still devoted to itschamber(b) ventceramic mouldpouringbasininduction - coilcopperqjfinger*rskullihtest barorblademoltenalloycopper basewater in, water out <=tzstudy. net shape casting is very attractive from the c

10、ost point of view compared to alternative component manufacturing methods such as forging or machining. however，reproducible quality is still an issue; the elimination of defects and control of microstructure drive research. casting involves first the filling of the mould and subsequently the solidi

11、fication of the melt. from the numerical modelling point of view，this simple sequence results in a very complex three-phase problem to simulate. a range of interactions of physical phenomena arc involved including free surface fluid flow as the mould fills，heterogeneous heat transfer from the metal

12、to the mould， solidification of the molten metal as it cools，and the development of residual stresses and deformation of the solidified component.in industry there are many variants of the casting process such as sand casting, investment casting， gravity，and low and high pressure die casting. in thi

13、s study，the investment casting process，also called lost-wax casting, has been investigated. one of the advantages of this process is that it is capable of producing (near) net shape parts，which is particularly important for geometrically complex and difficult-to-machine components. this process star

14、ts with making a ceramic mould which involves three main steps: injecting wax into a die to make a replica of the component and attaching this to a pouring basin and running system; building a ceramic shell by the application of several layers of a ceramic slurry and ceramic stucco to the wax assemb

15、ly; de-waxing and mould firing. the pouring of the casting is performed either simply under gravity (no control), or using a rapid centrifugal action 1 (danger of macro-segregation plus highly turbulent filling)，or by suction as in counter-gravity casting (e.g. the hitchiner process2)，or by tilt-cas

16、ting. in this study， tilt-casting was chosen in an attempt to achieve tranquil mould filling. tilt-casting was patented in 1919 by durville 3 and has been successfully used with sand castings4 and aluminium die castings5. in the impress project 6，a novel process has been proposed and successfully de

17、veloped to combine induction skull melting (ism) of reactive alloys with tilt-casting7，8，9 and 10，with a particular application to the production of turbine blades in titanium aluminidc alloys. as shown in fig. 1，this is carried out inside a vacuum chamber and the mould is pre-heated in situ to avoi

18、d misruns (incomplete mould filling due to premature solidification) and mould cracking due to thermal shock.tilt-casting process: (a) experimental equipment; (b) schematic view of the ism crucible and mould，showing the domed shape acquired by the molten metal; (c) different stages of mould filling

19、showing the progressive replacement of gas by the metal.the component(s) to be cast are attached to a pouring basin which also doubles as a source of metal to feed the solidification shrinkage. the components are angled on the basin to promote the progressive uni-directional flow of metal into the m

20、ould. as the metal enters the mould it displaces the gas and an escape route has to be included in the design so that the two counter-flowing streams arc not mixed leading to bubbles trapped in the metal. vents are also used to enable any trapped gas to escape. the feeder used to connect the mould t

21、o the crucible is normally in any casting the last portion of metal to solidify，so supplying metal to the mould to counter the effects of solidification shrinkage. in tilt-casting, the feeder is also the conduit for the tranquil flow of metal into the mould and also for the unhindered escape of gas.

22、 for this reason, the fluid dynamics of the mould feeder interface merit detailed study.as well as the mould/feeder design, the production of castings involves several other key parameters，such as the metal pouring temperature, initial mould temperature, selective mould insulation and the tilt cycle

23、 timing. all these parameters have an influence on the eventual quality of the casting leading to a very large matrix of experiments. modelling (once validated) is crucial in reducing the amount of physical experiments required. as mentioned above，the mathematical models are complex due to the fact

24、that this is a three-phase problem with two rapidly developing phase fronts (liquid/gas and solid/liquid). in this paper, a 3-d computational model is used to simulate the tilt-casting process and to investigate the effect of the design of the basin/feeder on the flow dynamics during mould filling a

25、nd eventually on casting quality.2. experimental descriptiondetails of the experimental setup have been published elsewhere 11，but for completeness a summary description is given here. fig. la shows an overall view of the equipment used to perform the casting. the induction skull melting (ism) coppe

26、r crucible is installed inside a vacuum chamber. to enable rotation, it is attached to a co-axial power feed，which also allows cooling water containing ethylene glycol to be supplied to the ism crucible and the induction coil. the coil supplies a maximum of 8 ka at a frequency of 6 khz. the crucible

27、 wall is segmented， so that the induction field penetrates through the slots (by inducing eddy currents into each finger segment) to melt the charge and at the same time repel the liquid metal away from the side wall to minimise the loss of superheat. a billet of tial alloy is loaded into the crucib

28、le before clamping on the ceramic shell mould. the mould is surrounded by a low thermal mass split-mould heater. after evacuating the vacuum chamber，the mould is heated to the required temperature (1200 °c maximum) and the vessel back-filled with argon to a partial pressure of 20 kpa prior to m

29、elting. this pressure significantly reduces the evaporative loss of the volatile aluminium contained in the alloy. the power applied to the induction coil is increased according to a pre-determined power vs. time schedule so that a reproducible final metal temperature is achieved. at the end of melt

30、ing (7-8 min), the mould heater is opened and moved away. the induction melting power is ramped down and，simultaneously，the ism crucible and mould are rotated by 180° using a programmable controller to transfer the metal into the mould. the mould containing the casting is held vertically as the

31、 metal solidifies and cools down.3. mathematical model3.1. fluid flow equationsthe modelling of the castingprocess has involved a number of complex computational techniques since there are a range of physical interactions to account for: free surface fluid flow, turbulence， heat transfer and solidif

32、ication, and so on. the fluid flow dynamics of the molten metal and the gas filling the rest of the space arc governed by the navicr-stokcs equations, and a 3d model is used to solve the incompressible time-dependent flow:(2)罢+.=0where u is the fluid velocity vector; p is the density; g is the fluid

33、 viscosity; su is a source term which contains body forces (such as gravitational force，a resistive force (darcy term) 12) and the influence of boundaries. there is a sharp, rapidly evolving，property interface separating metal and gas regions in these equations as explained below.3.2. free surface:

34、counter diffusion method (cdm)one of the difficulties of the simulation arises from the fact that two fluid media arc present during filling: liquid metal and resident gas and their density ratio is as high as 10,000:1. not only does the fluid flow problem need to be solved over the domain，but the m

35、odel also has to track the evolution of the interface of the two media with time. a scalar fluid marker o was introduced to represent the metal volume fraction in a control volume and used to track the interface of the two fluids，called the scalar equation algorithm (sea) by pericleous et al. 14. in

36、 a gas cell，0 = 0; in a metal cell，(i)= 1; for a partially filled cell (i) takes on an intermediate value which the interface of the two media crosses through. the dynamics of the interface are governed by the advection equation:(3)誓+ =0ixkpe/fek fermthe interface then represents a moving property d

37、iscontinuity in the domain, which has to be handled carefully to avoid numerical smearing. as in 14，an accurate explicit time stepping scheme such as that by van leer 15 may be used to prevent smearing. however, the scheme is then limited to extremely small time steps for stability, leading to very

38、lengthy computations. to overcome this problem, a new tracking method，the counter diffusion method (cdm) 11 and 16，was developed as a corrective mechanism to counter this numerical diffusion. this discretizes the free surface equation in a stable，fully implicit scheme which makes the computations an

39、 order of magnitude faster. the implementation assumes that an interface-normal counter diffusion flux can be defined for each internal face of the computational mesh and applied with opposite signs to elements straddling the interface as source terms for the marker variable. the equation for the fl

40、ux per unit area f can be written as:(4)f = cl»» h|( 1 中 少公where c is a scaling factor, a free parameter in cdm allowing the strength of the counter diffusion action to be adjusted，and n is the unit normal vector to the face in the mesh. of the two cells either side of the face，the one wit

41、h the lower value of the marker(i)d becomes the donor cell while the richer cell (i)a is the acceptor (in order to achieve the counter diffusion action). the proposed formula makes the counter diffusion action self-limiting as it is reduced to zero where the donor approaches zero (gas) and where the

42、 acceptor reaches unity (liquid). in this form，the adjustment remains conservative. quantitative validation of cdm against other vof type techniques is given in a later section of the paper for accuracy and efficiency.3.3. heat transfer and solidificationheat transfer takes place between the metal,

43、mould and gas，and between cold and hot metal regions as the mould filling is carried out. the heat flow is computed by a transient energy conservation equation:(5)where t is the temperature; k is the thermal conductivity; cp is the specific heat (properties can be functions of the local temperature

44、or other variables); st is the source term which represents viscous dissipation, boundary heat transfer and latent heat contributions when a phase change occurs. for the latter，a new marker variable fl is used to represent the liquid fraction of the metal with (1 - fl) being the volume fraction of s

45、olidified metal. voller et al. 13 used a non-linear temperature function to calculate the liquid fraction. in this study, the liquid fraction is assumed to be a linear function of the metal temperature:11 r>n(atl is the liquidus temperature and ts is the solidus temperature.3.4. lvel turbulence m

46、odel (applied to solid moving boundaries)even at low filling speeds，the reynolds number is such that the flow is turbulent. the lvel method of spalding 17 is chosen to compute the turbulence because of its mixing-length simplicity and robustness. lvel is an abbreviation of a distance from the neares

47、t wall (l) and thelocal velocity (vel). the approximate wall distance is solved by the eqs. (7) and (8):(7)where w is an auxiliary variable in the regions occupied by the moving fluid with boundary conditions w = 0 on all solid walls.£ =+|w1this distance and the local velocity are used in the c

48、alculation of the local reynolds number from which the local value of the turbulent viscosity vt is obtained using a universal non-dimensional velocity profile away from the wall. the effective turbulent viscosity is then computed from the following equation:(9) where k = 0.417 is the von karman con

49、stant，e = 8.6 is the logarithmic law constant 17 and u+ is determined implicitly from the local reynolds number reloc = ul/v with the magnitude of the local velocity u and the laminar kinematic viscosity v17. the lvel method was extended to moving solid boundaries and in particular to solidifying re

50、gions by setting w = 0 in every region that is no longer fluid and then solving eqs. (7) and (8) at each time step.in simulating the tilt-casting process, the geometry is kept stationary and the gravitational force vector is rotated to numerically model the tilt instead of varying the coordinates of

51、 the geometry. the rotating gravitational force vector appears in the source term of eq. (1) for the tilt-casting process. a mathematical expression relating the tilting speed to the tilting angle 0 has been used. since 0 is a function of time，the variable rotation speed is adjustable to achieve tra

52、nquil filling. this technique neglects rotational forces within the fluid (centrifugal，coriolis) since they are negligible at the slow rotation rates encountered in tilt-casting. finally, the numerical model of the tilt-casting process and the new algorithm developments were implemented in the gener

53、al cfd package (physica).4. description of simulations4.1. geometry, mould design and computational meshthe casting is a generic 0.4 m-long turbine blade typical of that used in an industrial gas turbine. fig. 2 shows three mould designs which comprise the blade, a feeder/basin and a cylindrical cru

54、cible. fig. 2a incorporates a separate cube-shaped feeder that partially links the root of the blade and the basin. fig. 2b is a variant in which the plane of the blade is rotated through 90°. in both cases, the computational mesh contains 31,535 elements and 38,718 points. six vents are locate

55、d on the platform and the shroud of the blade，as seen in fig. 2a and b. fig. 2c is an optimised design where the feeder and basin are combined to provide a smooth connection between the blade and the crucible. two vents are located in the last areas to be filled to help entrapped gas to escape from

56、the mould. mesh of the cruciblc-mould assembly for the three eases investigated.the mesh for the last case contains 30,185 elements and 37,680 vertices. as in all the cases presented, numerical accuracy depends on mesh fineness and also the degree of orthogonality. to ensure a mostly orthogonal mesh

57、 the various components of the assembly were created separately using a structured body-fitted mesh generator and then joined using a mixture of hcxahcdral and tetrahedral cells. the mesh was refined as necessary in thin sections (such as the blade itself or the shroud and base plates)，but not neces

58、sarily to be fine enough to resolve boundary layer details. for this reason the lvel turbulence model was used rather than a more usual two-equation model of turbulence that relies on accurate wall function representation. the practical necessity to run in parallel with the experimental programme al

59、so limited the size of the mesh used. as with all free surface tracking algorithms，the minimum cell size determines the time step size for the stable simulations. although the cdm method is implicit，allowing the time step to exceed the cell cfl limit, accuracy is then affected. with these restrictions，turnaround time for a complete til