以上為資料預(yù)覽概圖,下載文件后為壓縮包資料文件?!厩逦?,無水印,可編輯】dwg后綴為cad圖,doc后綴為word格式,需要請自助下載即可,有疑問可以咨詢QQ 197216396 或 11970985
Jou rnal of Zhejiang University SCIENCE A ISSN 1 673 565X Print ISSN 1 862 1 775 Online www zju edu cn jzus www springerlink com E mail jzus zju edu ca Lf ef al J Zhejiang Univ SciA 2007 8 7 1077 1083 Single gate optimization for plastic inj ection mold LI Ji quan十 LI De qun GUO Zhi ying LV Hai yuan Department ofPlasticity Technology Shanghai Jiao Tong University Shanghai 200030 China E mail hutli 163 corn Received Nov 22 2006 revision accepted Mar 1 9 2007 1077 Abstract This paper deals with a methodology for single gate location optimization for plastic injection mold The objective of the gate optimization is to minimize the warpage of injection molded parts because warpage is a crucial quality issue for most injection molded parts while it is influenced greatly by the gate location Feature warpage is defined as the ratio of maximum displacement on the feature surface to the projected length of the feature surface to describe part warpage The optimization is combined with the numerical simulation technology to find the optimal gate location in which the simulated annealing algorithm is used to search for the optimum Finally an example is discussed in the paper and it can be concluded that the proposed method is effective Key words Injection mold Gate location Optimization Feature warpage doi 10 163 1 jzus 2007 A1077 Document code A CLC number TQ320 66 INTRODUCTION Plastic injection molding is a widely used com plex but highly efficient technique for producing a large variety of plastic products particularly those with high production requirement tight tolerance and complex shapes The quality ofinjection molded parts is a function of plastic material part geometry mold structure and process conditions The most important part of an injection mold basically is the following three sets of components cavities gates and runners and cooling system Lam and Seow 2000 and Jin and Lain 2002 achieved cavity balancing by varying the wall thick ness of the part A balance filling process within the cavity gives an evenly distributed pressure and tem perature which can drastically reduce the warpage of the part But the cavity balancing is only one of the important influencing factors of part qualities Espe cially the part has its functional requirements and its thicknesses should not be varied usually Project No 50675080 supported by the National Natural Science Foundation of China From the pointview ofthe injection mold design a gate is characterized by its size and location and the runner system by the size and layout The gate size and ruIlller layout are usually determined as constants Relatively gate locations and runner sizes are more flexible which can be varied to influence the quality of the part As a result they are often the design pa rameters for optimization Lee and Kim r l 996a optimized the sizes of runners and gates to balance runner system for mul tiple inj ection cavities The runner balancing was described as the differences of entrance pressures for a multi cavity mold with identical cavities and as differences of pressures at the end of the melt flow path in each cavity for a family mold wim different cavity volumes and geometries The methodology has shown uniform pressure distributions among the cavities during the entire molding cycle of multiple cavities mold Zhai et a1 2005a1 presented the two gate loca tion optimization of one molding cavity by an ef前一 cient search method based on pressure gradient PGSS and subsequently positioned weld lines to the desired locations by varying runner sizes for 維普資訊 1078 L e a1 d Zhejiang Univ SciA 2007 8 7 1077 1083 multi gate parts Zhai et a1 2006 As large volume part multiple gates are needed to shorten the maxi mum flow path with a corresponding decrease in injection pressure The method is promising for de sign of gates and runners for a single cavity wich multiple gates Many of injection molded parts are produced with one gate whether in single cavity mold or in multiple cavities mold Therefore the gate location of a single gate is the most common design parameter for optimization A shape analysis approach was pre sented by COHrbebaisse and Garcia 2002 by which the best gate location of iniection molding was esti mated Subsequently they developed this methodo1 ogy further and applied it to single gate location op timization of an L shape example Courbebaisse 2005 It is easy to use and not time consuming while it only serves the turning of simple flat parts wich uniforii1 thickness Pandelidis and Zou r 1 990 presented the opti mization ofgate location by indirect quality measures relevant to warpage and material degradation which is represented as weighted sum of a temperature dif ferential term an over pack term and a frictional overheating term Warpage is influenced by the above factors but the relationship between them is not clear Therefore the optimization effect is restricted by the determination of the weighting factors Lee and Kim r l 996b developed an automated selection method of gate location in which a set of initial gate locations were proposed by a designer and then the optimal gate was located by the a acent node evaluation method The conclusion to a great extent depends much on the human designer s intuition because the first step of the method is based on the designer s proposition So the result is to a large ex tent limited to the designer s experience Lam and Jin 200 l developed a gate location optimization method based on the minimization of the Standard Deviation of Flow Path Length SD L and Standard Deviation of Filling Time SD T during the molding filling process Subsequently Shen et 口 004a 2004b optimized the gate location design by minimizing the weighted sum of filling pressure filling time difference between different flow paths temperature difference and over pack percentage Zhai et a1 2005b investigated optimal gate location with evaluation criteria of iniection pressure at the end of filling These researchers presented the obi ec tive functions as performances of ini ection molding filling operation which are correlated with product qualities But the correlation between the perforii1 ances and qualities is very complicated and no clear relationship has been observed between them yet It is also difficult to select appropriate weighting factors for eachterm A new objective function is presented here to evaluate the warpage of inj ection molded parts to optimize gate location To measure part quality di rectly this investigation defines feature warpage to evaluate part warpage which is evaluated from the flow plus warpage simulation outputs of Moldflow Plastics Insight MPI software The obj ective func tion is minimized to achieve minimum deformation in gate location optimization Simulated annealing al gorithm is employed to search for the optimal gate location An example is given to illustrate the effec tivity of the proposed optimization procedure QUALITY MEASURES FEATURE WARPGE Definition of feature warpage To apply optimization theory to the gate design quality measures of the part must be specified in the first instance The term quality may be referred to many product properties such as mechanical thermal electrical optical ergonomical or geometrical prop erties There are two types of part quality measures direct and indirect A model that predicts the proper ties from numerical simulation results would be characterized as a direct quality measure In contrast an indirect measure of part quality is correlated with target quality but it cannot provide a direct estimate Ofchat quality For warpage the indirect quality measures in related works are one of performances of iniection molding flowing behavior or weighted sum of those The performances are presented as filling time dif ferential along different flow paths temperature dif ferential over pack percentage and so on It is ob vious that warpage is influenced by these perforii1 ances but the relationship between warpage and these performances is not clear and the determination of these weighting factors is rather difficult Therefore the optimization with the above Obiective function 維普資訊 Lf ef al J Zhejiang Univ SciA 2007 8 1077 1083 probably will not minimize part warpage even with perfect optimization technique Sometimes improper weighting factors will result in absolutely wrong re sults Some statistical quantities calculated from the nodal displacements were characterized as direct quality measures to achieve minimum deformation in related optimization studies The statistical quantities are usually a maximum nodal displacement an av erage oftop 1 0 percentile nodal displacements and an overall average nodal displacement Lee and Kim l 995 I 996b These nodaI displacements are easy to obtain from the simulation results the statistical val ues to some extents representing the deformation But the statistical displacement cannot effectively describe the deformation of the in i ection molded parts In industry designers and manufacturers usually pay more attention to the degree of part warpage on some specific features than the whole deformation of the in iection molded parts In this study feature warpage is defined to describe the deformation of the injection parts The feature warpage is the ratio of the maximum displacement of the feature surface to the proj ected length of the feature surface Fig 1 h 100 L where y is the feature warpage h is the maximum displacement on the feature surface deviating from the reference platform and L is the proj ected length of the feature surface on a reference direction paralleling the reference platform Surface ceplane Fig 1 The definition of feature warpage For complicated features only plane feature discussed here the feature warpage is usually sepa rated into two constituents on the reference plane which are represented on a 2D coordinate system h l0 l0 0o o 1079 where yy are the constituent feature warpages in the X Y direction and Lx L are the projected lengths of the feature surface on Y component Evaluation of feature warpage Atier the determination of target feature com bined with corresponding reference plane and pro iection direction the value of L can be calculated immediately from the part with the calculating method of analytic geometry Fig 2 L is a constant for any part on the specified feature surface and pro iected direction But the evaluation of h is more com plicated than that ofL y Feature Fig 2 The projected length evaluation Simulation of iniection molding process is a common technique to forecast the quality of part de sign mold design and process settings The results of warpage simulation are expressed as the nodal de flections on Z component and the nodal displacement Wis the vector length ofvector sum of and Wz k where k are the unit vectors on Z component The h is the maximum displacement of the nodes on the feature surface which is correlated with the normal orientation of the reference plane and can be derived from the results of warpage simulation To calculate h the deflection of fth node is evaluated firstly as follows COSO cosfl W z cosy一 3 where is the deflection in the normal direction of the reference plane of ith node W WI1 Wiz are the deflections on X Y Z component of ith node 伐 B are the angles ofnormal vector ofthe reference A and B are the terminal nodes of the feature to proiecting direction Fig 2 and are the deflections of nodesA and 維普資訊 1080 Lf ef al dZhefiangUnivSciA 2007 8 1077 1083 f COSa cosfl cosy l COSa cosfl Gz cosy where WA WAy are the deflections on Z component ofnode Wsx and Wez are the de flections on X Y Z component ofnode COiA and lB are the weighting factors of the terminal node deflec tions calculated as follows 2 iA 1一LiA L co 1一 J iA where LiA is the proj ector distance between fth node and node A Ultimately h is the maximum of the absolute value of h max l I I 1 In industry the inspection of the warpage is carried out with the help ofa feeler gauge while the measured part should be placed on a reference plat form The value of h is the maximum numerical reading of the space between the measured part sur face and the reference platform GATE L0CAT10N 0PTIMIZAT10N PR0BLEM FORMATION The quality term warpage means the perma nent deformation of the part which is not caused by an applied load It is caused by differential shrinkage throughout the part due to the imbalance of polymer flow packing cooling and crystallization The placement of a gate in an iniection mold is one of the most important variables of the total mold design The quality of the molded part is greatly af fected by the gate location because it influences the manner that the plastic f1ows into the mold cavity Therefore different gate locations introduce inho mogeneity in orientation density pressure and temperature distribution accordingly introducing different value and distribution of warpage Therefore gate location is a valuable design variable to minimize the injection molded part warpage Because the cor relation between gate location and warpage distribu tion is to a large extent independent of the melt and mold temperature it is assumed that the molding conditions are kept constant in this investigation The inj ection molded part warpage is quantified by the feature warpage which was discussed in the previous section The single gate location optimization can thus be formulated as fo11ows Minimize mini X Subjectto g X P Po一1 0 X X i 1 2 N where is the feature warpage P is the injection pressure at the gate position P0 is the allowable in iection pressure of injection molding machine or the allowable injection pressure specified by the designer or manufacturer X is the coordinate vector of the candidate gate locations is the node on the finite element mesh model of the part for injection molding process simulation N is the total number of nodes In the finite element mesh model of the part every node is a possible candidate for a gate There fore the total number of the possible gate location n is a function of the total number of nodes N and the total number of gate locations to be optimized N N一1 一 一1 In this study only the single gate location problem is investigated SIMULATED ANNEALING ALG0RITHM The simulated annealing algorithm is one of the most powerful and popular meta heuristics to solve optimization problems because of the provision of good global solutions to real world problems The algorithm is based upon that of Metropolis et a1 1 953 which was originally proposed as a means to find an equilibrium configuration of a collection of atoms at a given temperature The connection be tween this algorithm and mathematical minimization was first noted by Pincus r1970 but it was Kirkpatrick et a1 1983 who proposed that it formed the basis of an optimization technique for combina tional and other problems To apply the simulated annealing method to op 維普資訊 Lf ef al d Zhejiang Univ SciA 2007 8 1077 1083 timization problems the obiective function f is used as an energY function E Instead of finding a low energY configuration the problem becomes to seek an approximate global optimal solution The configura tions Of me values of design variables are substituted f or the energy cOnfiguratiOns Of the body and the control parameter for the process is substituted for temperature A random number generator is used as a way ofgenerating new values for the design variables It is obvious that this algorithm iust takes the mini mization problems into account Hence while per forming a maximization problem the obj ective func tion is multiplied by f l 1 to obtain a capable form The maj or advantage of simulated annealing algorithm over omer methods is the ability to avoid being trapped at local minima This algorithm em ploys a random search which not only accepts changes that decrease oNective function but also accepts some changes that increase it The Iatter are accepted with a probability P P e J where is the increase off k is Boltzman s constant and T is a control parameter which by analogy with the original application is known as the system temperature irrespective of the objective function involved In the case of gate location optimization the implementation ofthis algorithm is illustrated in Fig 3 and this algorithm js detailed as follows 1 SA algorithm starts from an initial gate loca tion Xold with an assigned value of the tempera ture parameter T fthe temperature counter k is initially set to zero Proper control parameter 0 c 1 in annealing process and Markov chain enera e are given 2 SA algorithm generates a new gate location Xnew in the neighborhood of Xou and the value of the objective functionflX is calculated 3 The new gate location will be accepted with probability determined bv me acceptance function 唧 min 1 exp 一k f Xn 一f Xo A uniform random variable Punif is generated in 0 1 If Punif 尸accept Xnew is accepted otherwise it is rejected 1081 f4 This process is repeated for a large enough number of iterations enerate for The sequence of trial gate locations generated in this way is known as Markov chain 5 A new Markov chain is then generated starting from the last accepted gate location in the previous Markov chain for a reduced temperature Tk l CTk and the same process continues for de creasing values of temperature until the algorithm stops System initialization f O 朋 O J generate X ld J 1 t 0 1 一 i一 Generate Xnew if一 Fig 3 The flow chart of the simulated annealing algorithm APPLICATION AND DISCUSSION The application to a presented in this section complex industrial part is to illustrate the proposed quality measure and optimization methodology The part is provided by a manufacturer as shown in Fig 4 In mis part the flatness of basal surface is the most important profile precision requirement Therefore the feature warpage is discussed on basal surface in which reference platform is specified as a horizontal plane attached to the basal surface and the longitu dinal direction is specified as pr0iected reference 維普資訊 1082 Lj ef al d Zhejiang Univ SciA 2007 8 1077 1083 direction The parameter h is the maximum basal surface deflection on the normal direction namely the vertical direction and the parameter L is the projected length of the basal surface to the longitudinal direc tion Fig 4 Industrial part provided by the manufacturer The material of the part is Nylon Zytel 101L 30 EGF DuPont Engineering Polymer The molding conditions in the simulation are listed in Table 1 Fig 5 shows the finite element mesh model of the part employed in the numerical simulation It has 1469 nodes and 2492 elements The objective func tion namely feature warpage is evaluated by Eqs 1 3 6 The h is evaluated from the results of Flow Warp Analysis Sequence in MPI by Eq 1 and the is measured on the industrial part immediately L 20 50 mm Table 1 The molding conditions in the simulation Conditions Values Fill time s Melt temperature C Mold temperature C Packing time s Packing pressure offilling pressure 2 5 295 70 10 80 MPI is the most extensive software for the in iection molding simulation which can recommend the best gate location based on balanced flow Gate loca tion analysis is an effective tool for gate location de sign besides empirical method For this part the gate location analysis of MPI recommends that the best gate location is near node N7459 as shown in Fig 5 The part warpage is simulated based on this recom mended gate and thus the feature warpage is evaluated 7 5 15 which is a great value In trial manufactur ing part warpage is visible on the sample work piece This is unacceptable for the manufacturer The great warpage on basal surface is caused by the uneven orientation distribution of the glass fiber as shown in Fig 6a Fig 6a shows that the glass fiber orientation changes from negative direction to posi tive direction because of the location of the gate par ticularly the greatest change of the fiber orientation appears near the gate The great diversifi