2013 Vol.26(3)
Display Mode： |
2013, 27(3).
Abstract:
The flow characteristics of cavitation jets are essential issues among relevant studies. The physical properties of the jet are largely determined by the geometrical parameters of the nozzle. The structure and cavitation jets characteristics of the angular-nozzle and the self-resonating cavitation nozzle have been extensively studied, but little research is conducted in the central-body cavitation nozzle mainly because of its hard processing and the cavitation jet effect not satisfactory. In this paper, a novel central-body nozzle (a non-plunger central-body nozzle with square outlet) is studied to solve above problems. Submerged jets discharged from the novel central-body nozzle are simulated, employing the full cavitation model. The impact of nozzle configuration on jet properties is analyzed. The analysis results indicate that when central-body relative diameter keeps constant, there is an optimal contraction degree of nozzle’s outlet, which can induce intense cavitation in the jet. The central-body relative diameter also affects jet profiles. In the case of large central-body relative diameter, most of the bubbles settle in the jet core. On the contrary, a smaller relative diameter makes bubbles concentrate in the interface between the jet and its surrounding fluid. Moreover, the shorter outlet part allows the cavitation zone further extend in both the axial and racial directions. The research results further consummate the study on the central-body nozzles and the correlation between cavitation jet and the structure, and elementarily reveal the mechanism of cavitation jet produced in a non-plunger novel central-body nozzle and the effect of the structure parameters on the cavitation jet, moreover, provide the theoretical basis for the optimal design of the nozzle.
The flow characteristics of cavitation jets are essential issues among relevant studies. The physical properties of the jet are largely determined by the geometrical parameters of the nozzle. The structure and cavitation jets characteristics of the angular-nozzle and the self-resonating cavitation nozzle have been extensively studied, but little research is conducted in the central-body cavitation nozzle mainly because of its hard processing and the cavitation jet effect not satisfactory. In this paper, a novel central-body nozzle (a non-plunger central-body nozzle with square outlet) is studied to solve above problems. Submerged jets discharged from the novel central-body nozzle are simulated, employing the full cavitation model. The impact of nozzle configuration on jet properties is analyzed. The analysis results indicate that when central-body relative diameter keeps constant, there is an optimal contraction degree of nozzle’s outlet, which can induce intense cavitation in the jet. The central-body relative diameter also affects jet profiles. In the case of large central-body relative diameter, most of the bubbles settle in the jet core. On the contrary, a smaller relative diameter makes bubbles concentrate in the interface between the jet and its surrounding fluid. Moreover, the shorter outlet part allows the cavitation zone further extend in both the axial and racial directions. The research results further consummate the study on the central-body nozzles and the correlation between cavitation jet and the structure, and elementarily reveal the mechanism of cavitation jet produced in a non-plunger novel central-body nozzle and the effect of the structure parameters on the cavitation jet, moreover, provide the theoretical basis for the optimal design of the nozzle.
2013, 27(3).
Abstract:
Conventional analytical and numerical methods for the mechanical properties of helical threads are relied on many assumptions and approximations and thus hardly yield satisfied results. A parameterized 3D finite element model of bolted joints with real helical thread geometry is established and meshed with refined hexahedral elements. The Von Mises plasticity criterion, kinematic hardening rule of materials and interfacial contacts are employed to make it possible for the suggested model be able to approach real assembly conditions. Then, the mechanical properties of bolted joints with different thread pitches, thread numbers and modular ratios are investigated, including the contact pressure distribution at joint interfaces, the axial load distribution and stress concentration in screw threads during the loading and unloading process. Simulation results indicate that the load distribution in screw threads produced by the suggested model agrees well the results from CHEN’s photoelastic tests. In addition, an interesting phenomenon is found that tightening the bolt with a large preload first and then adjusting the clamping force by unloading can make the load distribution more uniform and reduce the maximum residual equivalent stress in thread roots by up to 40%. This research provides a simple and practical approach to constructing the 3D finite element model and predicting the mechanical properties of helical thread connection.
Conventional analytical and numerical methods for the mechanical properties of helical threads are relied on many assumptions and approximations and thus hardly yield satisfied results. A parameterized 3D finite element model of bolted joints with real helical thread geometry is established and meshed with refined hexahedral elements. The Von Mises plasticity criterion, kinematic hardening rule of materials and interfacial contacts are employed to make it possible for the suggested model be able to approach real assembly conditions. Then, the mechanical properties of bolted joints with different thread pitches, thread numbers and modular ratios are investigated, including the contact pressure distribution at joint interfaces, the axial load distribution and stress concentration in screw threads during the loading and unloading process. Simulation results indicate that the load distribution in screw threads produced by the suggested model agrees well the results from CHEN’s photoelastic tests. In addition, an interesting phenomenon is found that tightening the bolt with a large preload first and then adjusting the clamping force by unloading can make the load distribution more uniform and reduce the maximum residual equivalent stress in thread roots by up to 40%. This research provides a simple and practical approach to constructing the 3D finite element model and predicting the mechanical properties of helical thread connection.
2013, 27(3).
Abstract:
There are two separate traditional approaches to model contact problems: continuum and atomistic theory. Continuum theory is successfully used in many domains, but when the scale of the model comes to nanometer, continuum approximation meets challenges. Atomistic theory can catch the detailed behaviors of an individual atom by using molecular dynamics (MD) or quantum mechanics, although accurately, it is usually time-consuming. A multiscale method coupled MD and finite element (FE) is presented. To mesh the FE region automatically, an adaptive method based on the strain energy gradient is introduced to the multiscale method to constitute an adaptive multiscale method. Utilizing the proposed method, adhesive contacts between a rigid cylinder and an elastic substrate are studied, and the results are compared with full MD simulations. The process of FE meshes refinement shows that adaptive multiscale method can make FE mesh generation more flexible. Comparison of the displacements of boundary atoms in the overlap region with the results from full MD simulations indicates that adaptive multiscale method can transfer displacements effectively. Displacements of atoms and FE nodes on the center line of the multiscale model agree well with that of atoms in full MD simulations, which shows the continuity in the overlap region. Furthermore, the Von Mises stress contours and contact force distributions in the contact region are almost same as full MD simulations. The method presented combines multiscale method and adaptive technique, and can provide a more effective way to multiscale method and to the investigation on nanoscale contact problems.
There are two separate traditional approaches to model contact problems: continuum and atomistic theory. Continuum theory is successfully used in many domains, but when the scale of the model comes to nanometer, continuum approximation meets challenges. Atomistic theory can catch the detailed behaviors of an individual atom by using molecular dynamics (MD) or quantum mechanics, although accurately, it is usually time-consuming. A multiscale method coupled MD and finite element (FE) is presented. To mesh the FE region automatically, an adaptive method based on the strain energy gradient is introduced to the multiscale method to constitute an adaptive multiscale method. Utilizing the proposed method, adhesive contacts between a rigid cylinder and an elastic substrate are studied, and the results are compared with full MD simulations. The process of FE meshes refinement shows that adaptive multiscale method can make FE mesh generation more flexible. Comparison of the displacements of boundary atoms in the overlap region with the results from full MD simulations indicates that adaptive multiscale method can transfer displacements effectively. Displacements of atoms and FE nodes on the center line of the multiscale model agree well with that of atoms in full MD simulations, which shows the continuity in the overlap region. Furthermore, the Von Mises stress contours and contact force distributions in the contact region are almost same as full MD simulations. The method presented combines multiscale method and adaptive technique, and can provide a more effective way to multiscale method and to the investigation on nanoscale contact problems.
2013, 27(3).
Abstract:
There are many researches on the deformation behavior of wrought magnesium alloys, such as AZ31, AZ80, AZ91, and ZK60 magnesium alloys at different temperatures and strain rates, but few of them focuses on the deformation behavior of AZ41M and ZK60M alloys, especially under the twin-roll casting (TRC) state. Meanwhile, the existing researches only focus on the grain refinement law of the magnesium alloys under deformation conditions, the deformation mechanism has not been revealed yet. The hot compression behavior of AZ41M and ZK60M magnesium alloys under the temperature and strain rate ranges of 250–400 °C and 0.001–1 s–1 are studied by thermal simulation methods using Gleeble 1500 machine and virtual simulation using finite element analysis software. Simulation results show that sine hyperbolic law is the most suitable flow stress model for wider deformation conditions. The most reasonable selected deformation conditions of ZK60M alloy is 350 C/0.1 s–1 for TRC and 350 C/1 s–1 for conventional casting (CC), while AZ41M alloy is 300 C/0.01 s–1 for TRC and 350 C/0.1 s–1 for CC. Deformation behavior and dynamic recrystallization (DRX) mechanism of them are analyzed at the same deformation conditions. The microstructures of AZ41M and ZK60M alloys are observed at different deformed conditions by optical microscopy (OM) and electron back scatter diffraction (EBSD) and it reveals the flow behavior and deformation mechanism of them. Working harden and work soften contribute to the activation of basal, non-basal slip systems which promote DRX. The proposed research reveals the deformation behavior and mechanism of the AZ41M and ZK 60M magnesium alloys and concludes their optimized deformation parameters and processes and provides a theory basis for their manufacturing and application.
There are many researches on the deformation behavior of wrought magnesium alloys, such as AZ31, AZ80, AZ91, and ZK60 magnesium alloys at different temperatures and strain rates, but few of them focuses on the deformation behavior of AZ41M and ZK60M alloys, especially under the twin-roll casting (TRC) state. Meanwhile, the existing researches only focus on the grain refinement law of the magnesium alloys under deformation conditions, the deformation mechanism has not been revealed yet. The hot compression behavior of AZ41M and ZK60M magnesium alloys under the temperature and strain rate ranges of 250–400 °C and 0.001–1 s–1 are studied by thermal simulation methods using Gleeble 1500 machine and virtual simulation using finite element analysis software. Simulation results show that sine hyperbolic law is the most suitable flow stress model for wider deformation conditions. The most reasonable selected deformation conditions of ZK60M alloy is 350 C/0.1 s–1 for TRC and 350 C/1 s–1 for conventional casting (CC), while AZ41M alloy is 300 C/0.01 s–1 for TRC and 350 C/0.1 s–1 for CC. Deformation behavior and dynamic recrystallization (DRX) mechanism of them are analyzed at the same deformation conditions. The microstructures of AZ41M and ZK60M alloys are observed at different deformed conditions by optical microscopy (OM) and electron back scatter diffraction (EBSD) and it reveals the flow behavior and deformation mechanism of them. Working harden and work soften contribute to the activation of basal, non-basal slip systems which promote DRX. The proposed research reveals the deformation behavior and mechanism of the AZ41M and ZK 60M magnesium alloys and concludes their optimized deformation parameters and processes and provides a theory basis for their manufacturing and application.
Rule- and PSO Algorithm-based Dynamic Spatial Rescheduling Method for Hull Curved Block Construction
2013, 27(3).
Abstract:
Due to no effective rescheduling method in hull curved block construction planning, existing scheduling planning can’t be applied in practical production effectively. Two-dimensional layout and dynamic attributes of block construction planning are considered to develop a spatial rescheduling method, which is based on the spatial points searching rule and the particle swarm optimization(PSO) algorithm. A dynamic spatial rescheduling method is proposed to solve the manufacturing problem of rush-order blocks. Through spatial rescheduling, the rescheduling start time, the current processing information set and rescheduling blocks set can be obtained automatically. By using and updating the data of these sets, the rescheduling method combines the PSO algorithm with the spatial points searching rule to determine the rescheduling start time and layout of the blocks. Three types of dynamic events, including rush-order block delay, existing block delay and existing block position changes, are used to address problems with different function goals by setting different function weights. Finally, simulations based on three types of rush-order block events are performed to validate this method, including single rush-order block, multi rush-order blocks at the same time and multi rush-order blocks at different times. The simulation results demonstrate that this method can solve the rush-order block problems in hull block construction and reduce the interference to the existing manufacturing schedule. The proposed research provides a new rescheduling method and helps instruct scheduler to make production planning in hull block construction.
Due to no effective rescheduling method in hull curved block construction planning, existing scheduling planning can’t be applied in practical production effectively. Two-dimensional layout and dynamic attributes of block construction planning are considered to develop a spatial rescheduling method, which is based on the spatial points searching rule and the particle swarm optimization(PSO) algorithm. A dynamic spatial rescheduling method is proposed to solve the manufacturing problem of rush-order blocks. Through spatial rescheduling, the rescheduling start time, the current processing information set and rescheduling blocks set can be obtained automatically. By using and updating the data of these sets, the rescheduling method combines the PSO algorithm with the spatial points searching rule to determine the rescheduling start time and layout of the blocks. Three types of dynamic events, including rush-order block delay, existing block delay and existing block position changes, are used to address problems with different function goals by setting different function weights. Finally, simulations based on three types of rush-order block events are performed to validate this method, including single rush-order block, multi rush-order blocks at the same time and multi rush-order blocks at different times. The simulation results demonstrate that this method can solve the rush-order block problems in hull block construction and reduce the interference to the existing manufacturing schedule. The proposed research provides a new rescheduling method and helps instruct scheduler to make production planning in hull block construction.
2013, 27(3).
Abstract:
At present, the inner cutters of a full face rock tunnel boring machine (TBM) and transition cutter edge angles are designed on the basis of indentation test or linear grooving test. The inner and outer edge angles of disc cutters are characterized as symmetric to each other with respect to the cutter edge plane. This design has some practical defects, such as severe eccentric wear and tipping, etc. In this paper, the current design theory of disc cutter edge angle is analyzed, and the characteristics of the rock-breaking movement of disc cutters are studied. The researching results show that the rotational motion of disc cutters with the cutterhead gives rise to the difference between the interactions of inner rock and outer rock with the contact area of disc cutters, with shearing and extrusion on the inner rock and attrition on the outer rock. The wear of disc cutters at the contact area is unbalanced, among which the wear in the largest normal stress area is most apparent. Therefore, a three-dimensional model theory of rock breaking and an edge angle design theory of transition disc cutter are proposed to overcome the flaws of the currently used TBM cutter heads, such as short life span, camber wearing, tipping. And a corresponding equation is established. With reference to a specific construction case, the edge angle of the transition disc cutter has been designed based on the theory. The application of TBM in some practical project proves that the theory has obvious advantages in enhancing disc cutter life, decreasing replacement frequency, and making economic benefits. The proposed research provides a theoretical basis for the design of TBM three-dimensional disc cutters whose rock-breaking operation time can be effectively increased.
At present, the inner cutters of a full face rock tunnel boring machine (TBM) and transition cutter edge angles are designed on the basis of indentation test or linear grooving test. The inner and outer edge angles of disc cutters are characterized as symmetric to each other with respect to the cutter edge plane. This design has some practical defects, such as severe eccentric wear and tipping, etc. In this paper, the current design theory of disc cutter edge angle is analyzed, and the characteristics of the rock-breaking movement of disc cutters are studied. The researching results show that the rotational motion of disc cutters with the cutterhead gives rise to the difference between the interactions of inner rock and outer rock with the contact area of disc cutters, with shearing and extrusion on the inner rock and attrition on the outer rock. The wear of disc cutters at the contact area is unbalanced, among which the wear in the largest normal stress area is most apparent. Therefore, a three-dimensional model theory of rock breaking and an edge angle design theory of transition disc cutter are proposed to overcome the flaws of the currently used TBM cutter heads, such as short life span, camber wearing, tipping. And a corresponding equation is established. With reference to a specific construction case, the edge angle of the transition disc cutter has been designed based on the theory. The application of TBM in some practical project proves that the theory has obvious advantages in enhancing disc cutter life, decreasing replacement frequency, and making economic benefits. The proposed research provides a theoretical basis for the design of TBM three-dimensional disc cutters whose rock-breaking operation time can be effectively increased.
2013, 27(3).
Abstract:
Carrying capacity of the casing will reduce after the casing is worn, which seriously affects the subsequent well drilling, well completion, oil extraction and well repair. A lot of researches on calculation of casing wear collapse strength have been done, but few of them focus on collapsing failure mechanism, and influencing factors and law of collapse strength. So, significant difference between estimated value and actual value of collapse strength comes into being. By theoretical analysis, numerical simulation and actual test, the collapsing failure mechanism of casing wear as well as the influencing factors and laws of collapse strength are investigated, and the investigation results show that collapse of crescent casing wear belongs to “three hinged” instability. The severely-worn position on the casing is yielded into the plastic zone first then deformed greatly, which causes the plastic instability of the whole structure. The casing wear collapse strength presents changes of exponent, power function and linear trend with the residual casing wall thickness, wear radius and axial load, respectively. When the flexibility is less than 10°30 m, the borehole bending has less impact on casing collapse strength. Thus, the computation model for the casing wear collapsing strength is established by introducing wear radius coefficient and casing equivalent yield strength, at the same time, the model is tested. The test results show that the relative error for the computation model is less than 5%. The research results provide a basis for design of the casing string strength and evaluation of down-hole safety.
Carrying capacity of the casing will reduce after the casing is worn, which seriously affects the subsequent well drilling, well completion, oil extraction and well repair. A lot of researches on calculation of casing wear collapse strength have been done, but few of them focus on collapsing failure mechanism, and influencing factors and law of collapse strength. So, significant difference between estimated value and actual value of collapse strength comes into being. By theoretical analysis, numerical simulation and actual test, the collapsing failure mechanism of casing wear as well as the influencing factors and laws of collapse strength are investigated, and the investigation results show that collapse of crescent casing wear belongs to “three hinged” instability. The severely-worn position on the casing is yielded into the plastic zone first then deformed greatly, which causes the plastic instability of the whole structure. The casing wear collapse strength presents changes of exponent, power function and linear trend with the residual casing wall thickness, wear radius and axial load, respectively. When the flexibility is less than 10°30 m, the borehole bending has less impact on casing collapse strength. Thus, the computation model for the casing wear collapsing strength is established by introducing wear radius coefficient and casing equivalent yield strength, at the same time, the model is tested. The test results show that the relative error for the computation model is less than 5%. The research results provide a basis for design of the casing string strength and evaluation of down-hole safety.
2013, 27(3).
Abstract:
The failure of the key parts, such as gears, in cutter head driving system of tunneling boring machine has not been properly solved under the interaction of driving motors asynchronously and wave tunneling torque load. A dynamic model of multi-gear driving system is established considering the inertia effects of driving mechanism and cutter head as well as the bending-torsional coupling. By taking into account the nonlinear coupling factors between ring gear and multiple pinions, the influence for meshing angle by bending-torsional coupling and the dynamic load-sharing characteristic of multiple pinions driving are analyzed. Load-sharing coefficients at different rotating cutter head speeds and input torques are presented. Numerical results indicate that the load-sharing coefficients can reach up to 1.2–1.3. A simulated experimental platform of the multiple pinions driving is carried out and the torque distributions under the step load in driving shaft of pinions are measured. The imbalance of torque distribution of pinions is verified and the load-sharing coefficients in each pinion can reach 1.262. The results of simulation and test are similar, which shows the correctness of theoretical model. A loop coupling control method is put forward based on current torque master slave control method. The imbalance of the multiple pinions driving in cutter head driving system of tunneling boring machine can be greatly decreased and the load-sharing coefficients can be reduced to 1.051 by using the loop coupling control method. The proposed research provides an effective solution to the imbalance of torque distribution and synchronous control method for multiple pinions driving of TBM.
The failure of the key parts, such as gears, in cutter head driving system of tunneling boring machine has not been properly solved under the interaction of driving motors asynchronously and wave tunneling torque load. A dynamic model of multi-gear driving system is established considering the inertia effects of driving mechanism and cutter head as well as the bending-torsional coupling. By taking into account the nonlinear coupling factors between ring gear and multiple pinions, the influence for meshing angle by bending-torsional coupling and the dynamic load-sharing characteristic of multiple pinions driving are analyzed. Load-sharing coefficients at different rotating cutter head speeds and input torques are presented. Numerical results indicate that the load-sharing coefficients can reach up to 1.2–1.3. A simulated experimental platform of the multiple pinions driving is carried out and the torque distributions under the step load in driving shaft of pinions are measured. The imbalance of torque distribution of pinions is verified and the load-sharing coefficients in each pinion can reach 1.262. The results of simulation and test are similar, which shows the correctness of theoretical model. A loop coupling control method is put forward based on current torque master slave control method. The imbalance of the multiple pinions driving in cutter head driving system of tunneling boring machine can be greatly decreased and the load-sharing coefficients can be reduced to 1.051 by using the loop coupling control method. The proposed research provides an effective solution to the imbalance of torque distribution and synchronous control method for multiple pinions driving of TBM.
2013, 27(3).
Abstract:
The control study of vehicle semi-active suspension with magneto-rheological (MR) dampers has been attracted much attention internationally. However, a simple, real time and easy implementing semi-active controller has not been proposed for the MR full-vehicle suspension system, and a systematic analysis method has not been established for evaluating the multi-objective suspension performances of MR full-vehicle vertical, pitch and roll motions. For this purpose, according to the 7-degree of freedom (DOF) full-vehicle dynamic system, a generalized 7-DOF MR and passive full-vehicle dynamic model is set up by employing the modified Bouc-wen hysteretic force-velocity (F-v) model of the MR damper. A semi-active controller is synthesized to realize independent control of the four MR quarter-vehicle sub-suspension systems in the full-vehicle, which is on the basis of the proposed modified skyhook damping scheme of MR quarter-vehicle sub-suspension system. The proposed controller can greatly simplify the controller design complexity of MR full-vehicle suspension and has merits of easy implementation in real application, wherein only absolute velocities of sprung and unsprung masses with reference to the road surface are required to measure in real time when the vehicle is moving. Furthermore, a systematic analysis method is established for evaluating the vertical, pitch and roll motion properties of both MR and passive full-vehicle suspensions in a more realistic road excitation manner, in which the harmonic, rounded pulse and real road measured random signals with delay time are employed as different road excitations inserted on the front and rear two wheels, by considering the distance between front and rear wheels in full-vehicle. The above excitations with different amplitudes are further employed as the road excitations inserted on left and right two wheels for evaluating the roll motion property. The multi-objective suspension performances of ride comfort and handling safety of the proposed MR full-vehicle suspension are thus thoroughly evaluated by comparing with those of the passive full-vehicle suspension. The results show that the proposed controller can ideally improve multi-objective suspension performances of the ride comfort and handling safety. The proposed harmonic, rounded pulse and real road measured random signals with delay time and asymmetric amplitudes are suitable for accurately analyzing the vertical, pitch and roll motion properties of MR full-vehicle suspension system in a more realistic road excitation manner. This research has important theoretical significance for improving application study on the intelligent MR semi-active suspension.
The control study of vehicle semi-active suspension with magneto-rheological (MR) dampers has been attracted much attention internationally. However, a simple, real time and easy implementing semi-active controller has not been proposed for the MR full-vehicle suspension system, and a systematic analysis method has not been established for evaluating the multi-objective suspension performances of MR full-vehicle vertical, pitch and roll motions. For this purpose, according to the 7-degree of freedom (DOF) full-vehicle dynamic system, a generalized 7-DOF MR and passive full-vehicle dynamic model is set up by employing the modified Bouc-wen hysteretic force-velocity (F-v) model of the MR damper. A semi-active controller is synthesized to realize independent control of the four MR quarter-vehicle sub-suspension systems in the full-vehicle, which is on the basis of the proposed modified skyhook damping scheme of MR quarter-vehicle sub-suspension system. The proposed controller can greatly simplify the controller design complexity of MR full-vehicle suspension and has merits of easy implementation in real application, wherein only absolute velocities of sprung and unsprung masses with reference to the road surface are required to measure in real time when the vehicle is moving. Furthermore, a systematic analysis method is established for evaluating the vertical, pitch and roll motion properties of both MR and passive full-vehicle suspensions in a more realistic road excitation manner, in which the harmonic, rounded pulse and real road measured random signals with delay time are employed as different road excitations inserted on the front and rear two wheels, by considering the distance between front and rear wheels in full-vehicle. The above excitations with different amplitudes are further employed as the road excitations inserted on left and right two wheels for evaluating the roll motion property. The multi-objective suspension performances of ride comfort and handling safety of the proposed MR full-vehicle suspension are thus thoroughly evaluated by comparing with those of the passive full-vehicle suspension. The results show that the proposed controller can ideally improve multi-objective suspension performances of the ride comfort and handling safety. The proposed harmonic, rounded pulse and real road measured random signals with delay time and asymmetric amplitudes are suitable for accurately analyzing the vertical, pitch and roll motion properties of MR full-vehicle suspension system in a more realistic road excitation manner. This research has important theoretical significance for improving application study on the intelligent MR semi-active suspension.
Design and Optimization for the Occupant Restraint System of Vehicle Based on a Single Freedom Model
2013, 27(3).
Abstract:
Throughout the vehicle crash event, the interactions between vehicle, occupant, restraint system (VOR) are complicated and highly non-linear. CAE and physical tests are the most widely used in vehicle passive safety development, but they can only be done with the detailed 3D model or physical samples. Often some design errors and imperfections are difficult to correct at that time, and a large amount of time will be needed. A restraint system concept design approach which based on single-degree-of-freedom occupant-vehicle model (SDOF) is proposed in this paper. The interactions between the restraint system parameters and the occupant responses in a crash are studied from the view of mechanics and energy. The discrete input and the iterative algorithm method are applied to the SDOF model to get the occupant responses quickly for arbitrary excitations (impact pulse) by MATLAB. By studying the relationships between the ridedown efficiency, the restraint stiffness, and the occupant response, the design principle of the restraint stiffness aiming to reduce occupant injury level during conceptual design is represented. Higher ridedown efficiency means more occupant energy absorbed by the vehicle, but the research result shows that higher ridedown efficiency does not mean lower occupant injury level. A proper restraint system design principle depends on two aspects. On one hand，the restraint system should lead to as high ridedown efficiency as possible, and at the same time, the restraint system should maximize use of the survival space to reduce the occupant deceleration level. As an example, an optimization of a passenger vehicle restraint system is designed by the concept design method above, and the final results are validated by MADYMO, which is the most widely used software in restraint system design, and the sled test. Consequently, a guideline and method for the occupant restraint system concept design is established in this paper.
Throughout the vehicle crash event, the interactions between vehicle, occupant, restraint system (VOR) are complicated and highly non-linear. CAE and physical tests are the most widely used in vehicle passive safety development, but they can only be done with the detailed 3D model or physical samples. Often some design errors and imperfections are difficult to correct at that time, and a large amount of time will be needed. A restraint system concept design approach which based on single-degree-of-freedom occupant-vehicle model (SDOF) is proposed in this paper. The interactions between the restraint system parameters and the occupant responses in a crash are studied from the view of mechanics and energy. The discrete input and the iterative algorithm method are applied to the SDOF model to get the occupant responses quickly for arbitrary excitations (impact pulse) by MATLAB. By studying the relationships between the ridedown efficiency, the restraint stiffness, and the occupant response, the design principle of the restraint stiffness aiming to reduce occupant injury level during conceptual design is represented. Higher ridedown efficiency means more occupant energy absorbed by the vehicle, but the research result shows that higher ridedown efficiency does not mean lower occupant injury level. A proper restraint system design principle depends on two aspects. On one hand，the restraint system should lead to as high ridedown efficiency as possible, and at the same time, the restraint system should maximize use of the survival space to reduce the occupant deceleration level. As an example, an optimization of a passenger vehicle restraint system is designed by the concept design method above, and the final results are validated by MADYMO, which is the most widely used software in restraint system design, and the sled test. Consequently, a guideline and method for the occupant restraint system concept design is established in this paper.
2013, 27(3).
Abstract:
Although a great deal of research has been dedicated to the synthesis of spiral bevel gears, little related to reverse engineering can be found. An approach is proposed to reverse the machine-tool settings of the pinion of a spiral bevel gear drive on the basis of the blank and tooth surface data obtained by a coordinate measuring machine(CMM). Real tooth contact analysis(RTCA) is performed to preliminary ascertain the contact pattern, the motion curve, as well as the position of the mean contact point. And then the tangent to the contact path and the motion curve are interpolated in the sense of the least square method to extract the initial values of the bias angle and the higher order coefficients(HOC) in modified roll motion. A trial tooth surface is generated by machine-tool settings derived from the local synthesis relating to the initial meshing performances and modified roll motion. An optimization objective is formed which equals the tooth surface deviation between the real tooth surface and the trial tooth surface. The design variables are the parameters describing the meshing performances at the mean contact point in addition to the HOC. When the objective is optimized within an arbitrarily given convergence tolerance, the machine-tool settings together with the HOC are obtained. The proposed approach is verified by a spiral bevel pinion used in the accessory gear box of an aviation engine. The trial tooth surfaces approach to the real tooth surface on the whole in the example. The results show that the convergent tooth surface deviation for the concave side on the average is less than 0.5 μm, and is less than 1.3 μm for the convex side. The biggest tooth surface deviation is 6.7 μm which is located at the corner of the grid on the convex side. Those nodes with relative bigger tooth surface deviations are all located at the boundary of the grid. An approach is proposed to figure out the machine-tool settings of a spiral bevel pinion by way of reverse engineering without having known the theoretical tooth surfaces and the corresponding machine-tool settings.
Although a great deal of research has been dedicated to the synthesis of spiral bevel gears, little related to reverse engineering can be found. An approach is proposed to reverse the machine-tool settings of the pinion of a spiral bevel gear drive on the basis of the blank and tooth surface data obtained by a coordinate measuring machine(CMM). Real tooth contact analysis(RTCA) is performed to preliminary ascertain the contact pattern, the motion curve, as well as the position of the mean contact point. And then the tangent to the contact path and the motion curve are interpolated in the sense of the least square method to extract the initial values of the bias angle and the higher order coefficients(HOC) in modified roll motion. A trial tooth surface is generated by machine-tool settings derived from the local synthesis relating to the initial meshing performances and modified roll motion. An optimization objective is formed which equals the tooth surface deviation between the real tooth surface and the trial tooth surface. The design variables are the parameters describing the meshing performances at the mean contact point in addition to the HOC. When the objective is optimized within an arbitrarily given convergence tolerance, the machine-tool settings together with the HOC are obtained. The proposed approach is verified by a spiral bevel pinion used in the accessory gear box of an aviation engine. The trial tooth surfaces approach to the real tooth surface on the whole in the example. The results show that the convergent tooth surface deviation for the concave side on the average is less than 0.5 μm, and is less than 1.3 μm for the convex side. The biggest tooth surface deviation is 6.7 μm which is located at the corner of the grid on the convex side. Those nodes with relative bigger tooth surface deviations are all located at the boundary of the grid. An approach is proposed to figure out the machine-tool settings of a spiral bevel pinion by way of reverse engineering without having known the theoretical tooth surfaces and the corresponding machine-tool settings.
2013, 27(3).
Abstract:
Among most traditional piezo water cooling systems, piezoelectric valve pumps are adopted as their driving sources. The valves in these pumps induce problems of shock and vibration and also make their structure complicated, which is uneasy to minimize and reduce their reliability and applicability of the whole system. In order to avoid these problems caused by valve structure, a novel valveless piezoelectric pump is developed, which integrates both functions of transforming and cooling. The pump’s Y-shape tree-like construction not only increases the efficiency of cooling but also the system reliability and applicability. Firstly, a multistage Y-shape treelike bifurcate tube is proposed, then a valveless piezoelectric pump with multistage Y-shape treelike bifurcate tubes is designed and its working principle is analyzed. Then, the theoretical analysis of flow resistance characteristics and the flow rate of the valveless piezoelectric pump are performed. Meanwhile, commercial software CFX is employed to perform the numerical simulation for the pump. Finally, this valveless piezoelectric pump is fabricated, the relationship between the flow rates and driving frequency, as well as the relationship between the back pressure and the driving frequency are experimentally investigated. The experimental results show that the maximum flow rate is 35.6 mLmin under 100 V peak-to-peak voltage (10.3 Hz) power supply, and the maximum back pressure is 55 mm H2O under 100 V (9 Hz) power supply, which validates the feasibility of the valveless piezoelectric pump with multistage Y-shape treelike bifurcate tubes. The proposed research provides certain references for the design of valveless piezoelectric pump and improves the reliability of piezo water cooling systems.
Among most traditional piezo water cooling systems, piezoelectric valve pumps are adopted as their driving sources. The valves in these pumps induce problems of shock and vibration and also make their structure complicated, which is uneasy to minimize and reduce their reliability and applicability of the whole system. In order to avoid these problems caused by valve structure, a novel valveless piezoelectric pump is developed, which integrates both functions of transforming and cooling. The pump’s Y-shape tree-like construction not only increases the efficiency of cooling but also the system reliability and applicability. Firstly, a multistage Y-shape treelike bifurcate tube is proposed, then a valveless piezoelectric pump with multistage Y-shape treelike bifurcate tubes is designed and its working principle is analyzed. Then, the theoretical analysis of flow resistance characteristics and the flow rate of the valveless piezoelectric pump are performed. Meanwhile, commercial software CFX is employed to perform the numerical simulation for the pump. Finally, this valveless piezoelectric pump is fabricated, the relationship between the flow rates and driving frequency, as well as the relationship between the back pressure and the driving frequency are experimentally investigated. The experimental results show that the maximum flow rate is 35.6 mLmin under 100 V peak-to-peak voltage (10.3 Hz) power supply, and the maximum back pressure is 55 mm H2O under 100 V (9 Hz) power supply, which validates the feasibility of the valveless piezoelectric pump with multistage Y-shape treelike bifurcate tubes. The proposed research provides certain references for the design of valveless piezoelectric pump and improves the reliability of piezo water cooling systems.
2013, 27(3).
Abstract:
The existing investigations on piezoelectric materials containing an elliptic hole mainly focus on remote uniform tensile loads. In order to have a better understanding of the fracture behavior of piezoelectric materials under different loading conditions, theoretical and numerical solutions are presented for an elliptic hole in transversely isotropic piezoelectric materials subjected to uniform internal shearing forces based on the complex potential approach. By solving ten variable linear equations, the analytical solutions inside and outside the hole satisfying the permeable electric boundary conditions are obtained. Taking PZT-4 ceramic into consideration, numerical results of electro-elastic fields along the edge of the hole and axes, and the electric displacements in the hole are presented. Comparison with stresses in transverse isotropic elastic materials shows that the hoop stress at the ends of major axis in two kinds of material equals zero for the various ratios of major to minor axis lengths; If the ratio is greater than 1, the hoop stress in piezoelectric materials is smaller than that in elastic materials, and if the ratio is smaller than 1, the hoop stress in piezoelectric materials is greater than that in elastic materials; When it is a circle hole, the shearing stress in two materials along axes is the same. The distribution of electric displacement components shows that the vertical electric displacement in the hole and along axes in the material is always zero though under the permeable electric boundary condition; The horizontal and vertical electric displacement components along the edge of the hole are symmetrical and antisymmetrical about horizontal axis, respectively. The stress and electric displacement distribution tends to zero at distances far from the elliptical hole, which conforms to the conclusion usually made on the basis of Saint-Venant’s principle. Unlike the existing work, uniform shearing forces acting on the edge of the hole, and the distribution of electro-elastic fields inside and outside the elliptic hole are considered.
The existing investigations on piezoelectric materials containing an elliptic hole mainly focus on remote uniform tensile loads. In order to have a better understanding of the fracture behavior of piezoelectric materials under different loading conditions, theoretical and numerical solutions are presented for an elliptic hole in transversely isotropic piezoelectric materials subjected to uniform internal shearing forces based on the complex potential approach. By solving ten variable linear equations, the analytical solutions inside and outside the hole satisfying the permeable electric boundary conditions are obtained. Taking PZT-4 ceramic into consideration, numerical results of electro-elastic fields along the edge of the hole and axes, and the electric displacements in the hole are presented. Comparison with stresses in transverse isotropic elastic materials shows that the hoop stress at the ends of major axis in two kinds of material equals zero for the various ratios of major to minor axis lengths; If the ratio is greater than 1, the hoop stress in piezoelectric materials is smaller than that in elastic materials, and if the ratio is smaller than 1, the hoop stress in piezoelectric materials is greater than that in elastic materials; When it is a circle hole, the shearing stress in two materials along axes is the same. The distribution of electric displacement components shows that the vertical electric displacement in the hole and along axes in the material is always zero though under the permeable electric boundary condition; The horizontal and vertical electric displacement components along the edge of the hole are symmetrical and antisymmetrical about horizontal axis, respectively. The stress and electric displacement distribution tends to zero at distances far from the elliptical hole, which conforms to the conclusion usually made on the basis of Saint-Venant’s principle. Unlike the existing work, uniform shearing forces acting on the edge of the hole, and the distribution of electro-elastic fields inside and outside the elliptic hole are considered.
2013, 27(3).
Abstract:
Work on dynamic topology optimization of engineering structures for vibration suppression has mainly addressed the maximization of eigenfrequencies and gaps between consecutive eigenfrequencies of free vibration, minimization of the dynamic compliance subject to forced vibration, and minimization of the structural frequency response. A dynamic topology optimization method of bi-material plate structures is presented based on power flow analysis. Topology optimization problems formulated directly with the design objective of minimizing the power flow response are dealt with. In comparison to the displacement or velocity response, the power flow response takes not only the amplitude of force and velocity into account, but also the phase relationship of the two vector quantities. The complex expression of power flow response is derived based on time-harmonic external mechanical loading and Rayleigh damping. The mathematical formulation of topology optimization is established based on power flow response and bi-material solid isotropic material with penalization(SIMP) model. Computational optimization procedure is developed by using adjoint design sensitivity analysis and the method of moving asymptotes(MMA). Several numerical examples are presented for bi-material plate structures with different loading frequencies, which verify the feasibility and effectiveness of this method. Additionally, optimum results between topological design of minimum power flow response and minimum dynamic compliance are compared, showing that the present method has strong adaptability for structural dynamic topology optimization problems. The proposed research provides a more accurate and effective approach for dynamic topology optimization of vibrating structures.
Work on dynamic topology optimization of engineering structures for vibration suppression has mainly addressed the maximization of eigenfrequencies and gaps between consecutive eigenfrequencies of free vibration, minimization of the dynamic compliance subject to forced vibration, and minimization of the structural frequency response. A dynamic topology optimization method of bi-material plate structures is presented based on power flow analysis. Topology optimization problems formulated directly with the design objective of minimizing the power flow response are dealt with. In comparison to the displacement or velocity response, the power flow response takes not only the amplitude of force and velocity into account, but also the phase relationship of the two vector quantities. The complex expression of power flow response is derived based on time-harmonic external mechanical loading and Rayleigh damping. The mathematical formulation of topology optimization is established based on power flow response and bi-material solid isotropic material with penalization(SIMP) model. Computational optimization procedure is developed by using adjoint design sensitivity analysis and the method of moving asymptotes(MMA). Several numerical examples are presented for bi-material plate structures with different loading frequencies, which verify the feasibility and effectiveness of this method. Additionally, optimum results between topological design of minimum power flow response and minimum dynamic compliance are compared, showing that the present method has strong adaptability for structural dynamic topology optimization problems. The proposed research provides a more accurate and effective approach for dynamic topology optimization of vibrating structures.
2013, 27(3).
Abstract:
During the process of designing the mixed-flow pump impeller, the meridional flow passage shape directly affects the obtained meridional flow field, which then has an influence on the three-dimensional impeller shape. However, the meridional flow passage shape is too complicated to be described by a simple formula for now. Therefore, reasonable parameter selection for the meridional flow passage is essential to the investigation. In order to explore the effects of the meridional flow passage shape on the impeller design and the hydraulic performance of the mixed-flow pump, the hub and shroud radius ratio (HSRR) of impeller and the outlet diffusion angle (ODA) of outlet zone are selected as the meridional flow passage parameters. 25 mixed-flow pump impellers, with specific speed of 496 under the design condition, are designed with various parameter combinations. Among these impellers, one with HSRR of 1.94 and ODA of 90° is selected to carry out the model test and the obtained experimental results are used to verify accuracies of the head and the hydraulic efficiency predicted by numerical simulation. Based on SIMPLE algorithm and standard k-ε two-equation turbulence model, the three-dimensional steady incompressible Reynolds averaged Navier-Stokes equations are solved and the effects of different parameters on hydraulic performance of mixed-flow pump impellers are analyzed. The analysis results demonstrate that there are optimal values of HSRR and ODA available, so the hydraulic performance and the internal flow of mixed-flow pumps can be improved by selecting appropriate values for the meridional flow passage parameters. The research on these two parameters, HSRR and ODA, has further illustrated influences of the meridional flow passage shape on the hydraulic performance of the mixed-flow pump, and is beneficial to improving the design of the mixed-flow pump impeller.
During the process of designing the mixed-flow pump impeller, the meridional flow passage shape directly affects the obtained meridional flow field, which then has an influence on the three-dimensional impeller shape. However, the meridional flow passage shape is too complicated to be described by a simple formula for now. Therefore, reasonable parameter selection for the meridional flow passage is essential to the investigation. In order to explore the effects of the meridional flow passage shape on the impeller design and the hydraulic performance of the mixed-flow pump, the hub and shroud radius ratio (HSRR) of impeller and the outlet diffusion angle (ODA) of outlet zone are selected as the meridional flow passage parameters. 25 mixed-flow pump impellers, with specific speed of 496 under the design condition, are designed with various parameter combinations. Among these impellers, one with HSRR of 1.94 and ODA of 90° is selected to carry out the model test and the obtained experimental results are used to verify accuracies of the head and the hydraulic efficiency predicted by numerical simulation. Based on SIMPLE algorithm and standard k-ε two-equation turbulence model, the three-dimensional steady incompressible Reynolds averaged Navier-Stokes equations are solved and the effects of different parameters on hydraulic performance of mixed-flow pump impellers are analyzed. The analysis results demonstrate that there are optimal values of HSRR and ODA available, so the hydraulic performance and the internal flow of mixed-flow pumps can be improved by selecting appropriate values for the meridional flow passage parameters. The research on these two parameters, HSRR and ODA, has further illustrated influences of the meridional flow passage shape on the hydraulic performance of the mixed-flow pump, and is beneficial to improving the design of the mixed-flow pump impeller.
2013, 27(3).
Abstract:
Wear and mileage performance are the foremost performances for truck bus radial (TBR) tires. There are a lot of researches about the tire wear performance as well as the contact patch phenomenon by using finite element analysis (FEA) method or testing. But there is little published data on the correlations between the footprint geometry and the tread wear performance of tires. In this paper, an experiment on tire-ground performance of TBR tires is carried out by using Tekscan. The real-time changes of contact-area pressure distribution that occurred during the process of continuous load and unload are recorded. Three types of tires that act differently in behavior under normal usage are analyzed. A new method of researching in tire tread wear, which focuses on the geometrical characters of the footprint, is put forward. The experimental results of the three tires are described by using footprint geometrical characters. On the basis of studying the changing laws of footprint geometrical characters during the loading process and considering consumer survey and factory feedback information, the correlations between the geometrical character of footprints and tread destruction form are built. The analyzed results show that a greater contact area coefficient and a steady coefficient of contact result in a better wear performance for TBR tires. The footprint-shape coefficient changing laws in the process of loading are found to have a very good coincidence with the tread wear of the three types of tires. Tires with a smaller footprint-shape coefficient are likely to have an average tread wear while avoiding the shoulder wear first. The proposed research provides a new solution to predict tire-ground performance at the point of footprint and several useful references for improving tire design.
Wear and mileage performance are the foremost performances for truck bus radial (TBR) tires. There are a lot of researches about the tire wear performance as well as the contact patch phenomenon by using finite element analysis (FEA) method or testing. But there is little published data on the correlations between the footprint geometry and the tread wear performance of tires. In this paper, an experiment on tire-ground performance of TBR tires is carried out by using Tekscan. The real-time changes of contact-area pressure distribution that occurred during the process of continuous load and unload are recorded. Three types of tires that act differently in behavior under normal usage are analyzed. A new method of researching in tire tread wear, which focuses on the geometrical characters of the footprint, is put forward. The experimental results of the three tires are described by using footprint geometrical characters. On the basis of studying the changing laws of footprint geometrical characters during the loading process and considering consumer survey and factory feedback information, the correlations between the geometrical character of footprints and tread destruction form are built. The analyzed results show that a greater contact area coefficient and a steady coefficient of contact result in a better wear performance for TBR tires. The footprint-shape coefficient changing laws in the process of loading are found to have a very good coincidence with the tread wear of the three types of tires. Tires with a smaller footprint-shape coefficient are likely to have an average tread wear while avoiding the shoulder wear first. The proposed research provides a new solution to predict tire-ground performance at the point of footprint and several useful references for improving tire design.
2013, 27(3).
Abstract:
The suitable process parameters for a two-stage turbo air classifier are important for obtaining the ultrafine powder that has a narrow particle-size distribution, however little has been published internationally on the classification process for the two-stage turbo air classifier in series. The influence of the process parameters of a two-stage turbo air classifier in series on classification performance is empirically studied by using aluminum oxide powders as the experimental material. The experimental results show the following: 1) When the rotor cage rotary speed of the first-stage classifier is increased from 2 300 rmin to 2 500 rmin with a constant rotor cage rotary speed of the second-stage classifier, classification precision is increased from 0.64 to 0.67. However, in this case, the final ultrafine powder yield is decreased from 79% to 74%, which means the classification precision and the final ultrafine powder yield can be regulated through adjusting the rotor cage rotary speed of the first-stage classifier. 2) When the rotor cage rotary speed of the second-stage classifier is increased from 2 500 rmin to 3 100 rmin with a constant rotor cage rotary speed of the first-stage classifier, the cut size is decreased from 13.16 μm to 8.76 μm, which means the cut size of the ultrafine powder can be regulated through adjusting the rotor cage rotary speed of the second-stage classifier. 3) When the feeding speed is increased from 35 kgh to 50 kgh, the “fish-hook” effect is strengthened, which makes the ultrafine powder yield decrease. 4) To weaken the “fish-hook” effect, the equalization of the two-stage wind speeds or the combination of a high first-stage wind speed with a low second-stage wind speed should be selected. This empirical study provides a criterion of process parameter configurations for a two-stage or multi-stage classifier in series, which offers a theoretical basis for practical production.
The suitable process parameters for a two-stage turbo air classifier are important for obtaining the ultrafine powder that has a narrow particle-size distribution, however little has been published internationally on the classification process for the two-stage turbo air classifier in series. The influence of the process parameters of a two-stage turbo air classifier in series on classification performance is empirically studied by using aluminum oxide powders as the experimental material. The experimental results show the following: 1) When the rotor cage rotary speed of the first-stage classifier is increased from 2 300 rmin to 2 500 rmin with a constant rotor cage rotary speed of the second-stage classifier, classification precision is increased from 0.64 to 0.67. However, in this case, the final ultrafine powder yield is decreased from 79% to 74%, which means the classification precision and the final ultrafine powder yield can be regulated through adjusting the rotor cage rotary speed of the first-stage classifier. 2) When the rotor cage rotary speed of the second-stage classifier is increased from 2 500 rmin to 3 100 rmin with a constant rotor cage rotary speed of the first-stage classifier, the cut size is decreased from 13.16 μm to 8.76 μm, which means the cut size of the ultrafine powder can be regulated through adjusting the rotor cage rotary speed of the second-stage classifier. 3) When the feeding speed is increased from 35 kgh to 50 kgh, the “fish-hook” effect is strengthened, which makes the ultrafine powder yield decrease. 4) To weaken the “fish-hook” effect, the equalization of the two-stage wind speeds or the combination of a high first-stage wind speed with a low second-stage wind speed should be selected. This empirical study provides a criterion of process parameter configurations for a two-stage or multi-stage classifier in series, which offers a theoretical basis for practical production.
2013, 27(3).
Abstract:
Currently, the modeling of cutting process mainly focuses on two aspects: one is the setup of the universal cutting force model that can be adapted to a broader cutting condition; the other is the setup of the exact cutting force model that can accurately reflect a true cutting process. However, there is little research on the prediction of chatter stablity in milling. Based on the generalized mathematical model of inserted cutters introduced by ENGIN, an improved geometrical, mechanical and dynamic model for the vast variety of inserted cutters widely used in engineering applications is presented, in which the average directional cutting force coefficients are obtained by means of a numerical approach, thus leading to an analytical determination of stability lobes diagram (SLD) on the axial depth of cut. A new kind of SLD on the radial depth of cut is also created to satisfy the special requirement of inserted cutter milling. The corresponding algorithms used for predicting cutting forces, vibrations, dimensional surface finish and stability lobes in inserted cutter milling under different cutting conditions are put forward. Thereafter, a dynamic simulation module of inserted cutter milling is implemented by using hybrid program of Matlab with Visual Basic. Verification tests are conducted on a vertical machine center for Aluminum alloy LC4 by using two different types of inserted cutters, and the effectiveness of the model and the algorithm is verified by the good agreement of simulation result with that of cutting tests under different cutting conditions. The proposed model can predict the cutting process accurately under a variety of cutting conditions, and a high efficient and chatter-free milling operation can be achieved by a cutting condition optimization in industry applications.
Currently, the modeling of cutting process mainly focuses on two aspects: one is the setup of the universal cutting force model that can be adapted to a broader cutting condition; the other is the setup of the exact cutting force model that can accurately reflect a true cutting process. However, there is little research on the prediction of chatter stablity in milling. Based on the generalized mathematical model of inserted cutters introduced by ENGIN, an improved geometrical, mechanical and dynamic model for the vast variety of inserted cutters widely used in engineering applications is presented, in which the average directional cutting force coefficients are obtained by means of a numerical approach, thus leading to an analytical determination of stability lobes diagram (SLD) on the axial depth of cut. A new kind of SLD on the radial depth of cut is also created to satisfy the special requirement of inserted cutter milling. The corresponding algorithms used for predicting cutting forces, vibrations, dimensional surface finish and stability lobes in inserted cutter milling under different cutting conditions are put forward. Thereafter, a dynamic simulation module of inserted cutter milling is implemented by using hybrid program of Matlab with Visual Basic. Verification tests are conducted on a vertical machine center for Aluminum alloy LC4 by using two different types of inserted cutters, and the effectiveness of the model and the algorithm is verified by the good agreement of simulation result with that of cutting tests under different cutting conditions. The proposed model can predict the cutting process accurately under a variety of cutting conditions, and a high efficient and chatter-free milling operation can be achieved by a cutting condition optimization in industry applications.
2013, 27(3).
Abstract:
Hydrostatic spindles are increasingly used in precision machine tools. Thermal error is the key factor affecting the machining accuracy of the spindle, and research has focused on spindle thermal errors through examination of the influence of the temperature distribution, thermal deformation and spindle mode. However, seldom has any research investigated the thermal effects of the associated Couette flow. To study the heat transfer mechanism in spindle systems, the criterion of the heat transfer direction according to the temperature distribution of the Couette flow at different temperatures is deduced. The method is able to deal accurately with the significant phenomena occurring at every place where thermal energy flowed in such a spindle system. The variation of the motion error induced by thermal effects on a machine work-table during machining is predicated by establishing the thermo-mechanical error model of the hydrostatic spindle for a high precision machine tool. The flow state and thermal behavior of a hydrostatic spindle is analyzed with the evaluated heat power and the coefficients of the convective heat transfer over outer surface of the spindle are calculated, and the thermal influence on the oil film stiffness is evaluated. Thermal drift of the spindle nose is measured with an inductance micrometer, the thermal deformation data 1.35 μm after running for 4 h is consistent with the value predicted by the finite element analysis’s simulated result 1.28 μm, and this demonstrates that the simulation method is feasible. The thermal effects on the processing accuracy from the flow characteristics of the fluid inside the spindle are analyzed for the first time.
Hydrostatic spindles are increasingly used in precision machine tools. Thermal error is the key factor affecting the machining accuracy of the spindle, and research has focused on spindle thermal errors through examination of the influence of the temperature distribution, thermal deformation and spindle mode. However, seldom has any research investigated the thermal effects of the associated Couette flow. To study the heat transfer mechanism in spindle systems, the criterion of the heat transfer direction according to the temperature distribution of the Couette flow at different temperatures is deduced. The method is able to deal accurately with the significant phenomena occurring at every place where thermal energy flowed in such a spindle system. The variation of the motion error induced by thermal effects on a machine work-table during machining is predicated by establishing the thermo-mechanical error model of the hydrostatic spindle for a high precision machine tool. The flow state and thermal behavior of a hydrostatic spindle is analyzed with the evaluated heat power and the coefficients of the convective heat transfer over outer surface of the spindle are calculated, and the thermal influence on the oil film stiffness is evaluated. Thermal drift of the spindle nose is measured with an inductance micrometer, the thermal deformation data 1.35 μm after running for 4 h is consistent with the value predicted by the finite element analysis’s simulated result 1.28 μm, and this demonstrates that the simulation method is feasible. The thermal effects on the processing accuracy from the flow characteristics of the fluid inside the spindle are analyzed for the first time.
2013, 27(3).
Abstract:
The wheel-rail contact problems, such as the number, location and the track of contact patches, are very important for optimizing the spatial structure of the rails and lowering the vehicle-turnout system dynamics. However, the above problems are not well solved currently because of having the difficulties in how to determine the multi-contact, to preciously present the changeable profiles of the rails and to establish an accurate spatial turnout system dynamics model. Based on a high-speed vehicle-turnout coupled model in which the track is modeled as flexible with rails and sleepers represented by beams, the line tracing extreme point method is introduced to investigate the wheel-rail multiple contact conditions and the key sections of the blade rail, longer nose rail, shorter rail in the switch and nose rail area are discretized to represent the varying profiles of rails in the turnout. The dynamic interaction between the vehicle and turnout is simulated for cases of the vehicle divergently passing the turnout and the multi-point contact is obtained. The tracks of the contact patches on the top of the rails are presented and the wheel-rail impact forces are offered in comparison with the contact patches transference on the rails. The numerical simulation results indicate that the length of two-point contact occurrence of a worn wheel profile and rails is longer than that of the new wheel profile and rails; The two-point contact definitely occurs in the switch and crossing area. Generally, three-point contact doesn’t occur for the new rail profile, which is testified by the wheel-rails interpolation distance and the first order derivative function of the tracing line extreme points. The presented research is not only helpful to optimize the structure of the turnout, but also useful to lower the dynamics of the high speed vehicle-turnout system.
The wheel-rail contact problems, such as the number, location and the track of contact patches, are very important for optimizing the spatial structure of the rails and lowering the vehicle-turnout system dynamics. However, the above problems are not well solved currently because of having the difficulties in how to determine the multi-contact, to preciously present the changeable profiles of the rails and to establish an accurate spatial turnout system dynamics model. Based on a high-speed vehicle-turnout coupled model in which the track is modeled as flexible with rails and sleepers represented by beams, the line tracing extreme point method is introduced to investigate the wheel-rail multiple contact conditions and the key sections of the blade rail, longer nose rail, shorter rail in the switch and nose rail area are discretized to represent the varying profiles of rails in the turnout. The dynamic interaction between the vehicle and turnout is simulated for cases of the vehicle divergently passing the turnout and the multi-point contact is obtained. The tracks of the contact patches on the top of the rails are presented and the wheel-rail impact forces are offered in comparison with the contact patches transference on the rails. The numerical simulation results indicate that the length of two-point contact occurrence of a worn wheel profile and rails is longer than that of the new wheel profile and rails; The two-point contact definitely occurs in the switch and crossing area. Generally, three-point contact doesn’t occur for the new rail profile, which is testified by the wheel-rails interpolation distance and the first order derivative function of the tracing line extreme points. The presented research is not only helpful to optimize the structure of the turnout, but also useful to lower the dynamics of the high speed vehicle-turnout system.
2013, 27(3).
Abstract:
Setting engine emission targets to meet diesel car requirements is particularly important in engine performance development phase. Many researches are focused on associating vehicle performance with engine targets, but most work is done by testing, which is time and cost consuming, furthermore, the relationship of vehicle and engine will change when either engine or vehicle changes. A GT-Drive model to simulate New European Driving Cycle (NEDC) for passenger car is developed and calibrated by testing data, model precision is controlled within 5%. Time distribution of engine operating conditions when car running NEDC cycle has been analyzed, 10 critical major engine operating points are summarized according to running time proportion. Emission of NOx and smoke control regions containing these 10 points for target engine are set. Vehicle emissions are simulated and evaluated during engine development after engine performance test data are got, and engine combustion system layout and calibration are adjusted until vehicle targets are met. Vehicle is tested in chassis dynamometer finally, the testing results show a good agreement with the simulated results with an error of less than 5%, which proves that the emission value exchange of vehicle and engine is reliable. Performance target decomposition method for passenger car diesel presented can greatly shorten the development cycle and save costs.
Setting engine emission targets to meet diesel car requirements is particularly important in engine performance development phase. Many researches are focused on associating vehicle performance with engine targets, but most work is done by testing, which is time and cost consuming, furthermore, the relationship of vehicle and engine will change when either engine or vehicle changes. A GT-Drive model to simulate New European Driving Cycle (NEDC) for passenger car is developed and calibrated by testing data, model precision is controlled within 5%. Time distribution of engine operating conditions when car running NEDC cycle has been analyzed, 10 critical major engine operating points are summarized according to running time proportion. Emission of NOx and smoke control regions containing these 10 points for target engine are set. Vehicle emissions are simulated and evaluated during engine development after engine performance test data are got, and engine combustion system layout and calibration are adjusted until vehicle targets are met. Vehicle is tested in chassis dynamometer finally, the testing results show a good agreement with the simulated results with an error of less than 5%, which proves that the emission value exchange of vehicle and engine is reliable. Performance target decomposition method for passenger car diesel presented can greatly shorten the development cycle and save costs.
2013, 27(3).
Abstract:
The current match method of electric powertrain still makes use of longitudinal dynamics, which can’t realize maximum capacity for on-board energy storage unit and can’t reach lowest equivalent fuel consumption as well. Another match method focuses on improving available space considering reasonable layout of vehicle to enlarge rated energy capacity for on-board energy storage unit, which can keep the longitudinal dynamics performance almost unchanged but can’t reach lowest fuel consumption. Considering the characteristics of driving motor, method of electric powertrain matching utilizing conventional longitudinal dynamics for driving system and cut-and-try method for energy storage system is proposed for passenger cars converted from traditional ones. Through combining the utilization of vehicle space which contributes to the on-board energy amount, vehicle longitudinal performance requirements, vehicle equivalent fuel consumption level, passive safety requirements and maximum driving range requirement together, a comprehensive optimal match method of electric powertrain for battery-powered electric vehicle is raised. In simulation, the vehicle model and match method is built in Matlab/simulink, and the Environmental Protection Agency (EPA) Urban Dynamometer Driving Schedule (UDDS) is chosen as a test condition. The simulation results show that 2.62% of regenerative energy and 2% of energy storage efficiency are increased relative to the traditional method. The research conclusions provide theoretical and practical solutions for electric powertrain matching for modern battery-powered electric vehicles especially for those converted from traditional ones, and further enhance dynamics of electric vehicles.
The current match method of electric powertrain still makes use of longitudinal dynamics, which can’t realize maximum capacity for on-board energy storage unit and can’t reach lowest equivalent fuel consumption as well. Another match method focuses on improving available space considering reasonable layout of vehicle to enlarge rated energy capacity for on-board energy storage unit, which can keep the longitudinal dynamics performance almost unchanged but can’t reach lowest fuel consumption. Considering the characteristics of driving motor, method of electric powertrain matching utilizing conventional longitudinal dynamics for driving system and cut-and-try method for energy storage system is proposed for passenger cars converted from traditional ones. Through combining the utilization of vehicle space which contributes to the on-board energy amount, vehicle longitudinal performance requirements, vehicle equivalent fuel consumption level, passive safety requirements and maximum driving range requirement together, a comprehensive optimal match method of electric powertrain for battery-powered electric vehicle is raised. In simulation, the vehicle model and match method is built in Matlab/simulink, and the Environmental Protection Agency (EPA) Urban Dynamometer Driving Schedule (UDDS) is chosen as a test condition. The simulation results show that 2.62% of regenerative energy and 2% of energy storage efficiency are increased relative to the traditional method. The research conclusions provide theoretical and practical solutions for electric powertrain matching for modern battery-powered electric vehicles especially for those converted from traditional ones, and further enhance dynamics of electric vehicles.
2013, 27(3).
Abstract:
A lot of research has been focused on the necking process during the plastic deformation of sheet metals, but the localized necking is rarely distinguished form diffused necking by experiments, due to the limit of measurement equipment and method. Quenching and Partitioning (Q&P) steel is a 3rd generation advanced high strength steel (AHSS). Its good combination of high strength and ductility ensures potential application in automobile industry. Uniaxial tensile tests of QP980 steel sheet at five strain rates are performed to investigate the necking process and the effect of strain rate on necking behavior of Q&P steel. Digital image correlation (DIC) method is applied during tensile tests, and evolutions of major strain, minor strain and normal strain distributions along gauge section of the tensile specimens are obtained. The diffused and localized necking strains are determined according to SWIFT necking theory and HILL necking theory respectively. The test results indicate that with the increasing of strain rate in the investigated range, the diffused necking strain decreases from 0.152 to 0.120 and localized necking strain decreases from 0.245 to 0.137. Meanwhile, the difference of the two strains decreases form 0.096 to 0.017. Thus it can be concluded that strain rate has an influence on both necking strains during the deformation of QP980 steel sheet. Diffused and localized necking strains are determined by uniaxial tensile tests with the aid of DIC technique and the effect of strain rate on necking strains is evaluated.
A lot of research has been focused on the necking process during the plastic deformation of sheet metals, but the localized necking is rarely distinguished form diffused necking by experiments, due to the limit of measurement equipment and method. Quenching and Partitioning (Q&P) steel is a 3rd generation advanced high strength steel (AHSS). Its good combination of high strength and ductility ensures potential application in automobile industry. Uniaxial tensile tests of QP980 steel sheet at five strain rates are performed to investigate the necking process and the effect of strain rate on necking behavior of Q&P steel. Digital image correlation (DIC) method is applied during tensile tests, and evolutions of major strain, minor strain and normal strain distributions along gauge section of the tensile specimens are obtained. The diffused and localized necking strains are determined according to SWIFT necking theory and HILL necking theory respectively. The test results indicate that with the increasing of strain rate in the investigated range, the diffused necking strain decreases from 0.152 to 0.120 and localized necking strain decreases from 0.245 to 0.137. Meanwhile, the difference of the two strains decreases form 0.096 to 0.017. Thus it can be concluded that strain rate has an influence on both necking strains during the deformation of QP980 steel sheet. Diffused and localized necking strains are determined by uniaxial tensile tests with the aid of DIC technique and the effect of strain rate on necking strains is evaluated.
2013, 27(3).
Abstract:
Few function about 3D tool radius compensation is applied to generating executable motion control commands in the existing computer numerical control (CNC) systems. Once the tool radius is changed, especially in the case of tool size changing with tool wear in machining, a new NC program has to be recreated. A generic 3D tool radius compensation method for multi-axis peripheral milling in CNC systems is presented. The offset path is calculated by offsetting the tool path along the direction of the offset vector with a given distance. The offset vector is perpendicular to both the tangent vector of the tool path and the orientation vector of the tool axis relative to the workpiece. The orientation vector equations of the tool axis relative to the workpiece are obtained through homogeneous coordinate transformation matrix and forward kinematics of generalized kinematics model of multi-axis machine tools. To avoid cutting into the corner formed by the two adjacent tool paths, the coordinates of offset path at the intersection point have been calculated according to the transition type that is determined by the angle between the two tool path tangent vectors at the corner. Through the verification by the solid cutting simulation software VERICUT with different tool radiuses on a table-tilting type five-axis machine tool, and by the real machining experiment of machining a soup spoon on a five-axis machine tool with the developed CNC system, the effectiveness of the proposed 3D tool radius compensation method is confirmed. The proposed compensation method can be suitable for all kinds of three- to five-axis machine tools as a general form.
Few function about 3D tool radius compensation is applied to generating executable motion control commands in the existing computer numerical control (CNC) systems. Once the tool radius is changed, especially in the case of tool size changing with tool wear in machining, a new NC program has to be recreated. A generic 3D tool radius compensation method for multi-axis peripheral milling in CNC systems is presented. The offset path is calculated by offsetting the tool path along the direction of the offset vector with a given distance. The offset vector is perpendicular to both the tangent vector of the tool path and the orientation vector of the tool axis relative to the workpiece. The orientation vector equations of the tool axis relative to the workpiece are obtained through homogeneous coordinate transformation matrix and forward kinematics of generalized kinematics model of multi-axis machine tools. To avoid cutting into the corner formed by the two adjacent tool paths, the coordinates of offset path at the intersection point have been calculated according to the transition type that is determined by the angle between the two tool path tangent vectors at the corner. Through the verification by the solid cutting simulation software VERICUT with different tool radiuses on a table-tilting type five-axis machine tool, and by the real machining experiment of machining a soup spoon on a five-axis machine tool with the developed CNC system, the effectiveness of the proposed 3D tool radius compensation method is confirmed. The proposed compensation method can be suitable for all kinds of three- to five-axis machine tools as a general form.
2013, 27(3).
Abstract:
The unified analysis of omni-directional wheeled mobile manipulator(OWMM) through the dimensionally nonhomogeneous Jacobian matrix may lead to unreliable results. The existing researches focus on the integrated performance evaluation of OWMM, without taking the influence of velocity difference into consideration. This paper presents a new approach to formulate the dimensionally homogeneous Jacobian matrix and a new index for OWMM. First, the universal transformational matrix of the link coordinate frame is derived based on double quaternion. The degree of locomotion of a mobile platform and the degrees of manipulation of a manipulator are treated equally as the joints of a redundant composite robot. Then the integrated modeling of OWMM is established, and the non-dimensional Jacobian matrix is obtained from the above matrix. Next, by examining the concept of directional manipulability, a new index is proposed to evaluate the directional manipulability of OWMM along the specified task direction. Furthermore, for the same task and the time of simulation, the kinematic performance of the fixed base operation of the manipulator and the operation of OWMM are analyzed by numerical simulation. The results suggest that the proposed approach is equivalent to the method of traditional kinematic modeling, with the error of numerical solution being 10–7, and the varying frequency of DMTR is higher than that of the condition number. This indicates that DMTR is more effective to reflect the variation of the task direction. The proposed method can be used to analyze the manipulating capability of OWMM, and has a simple structure and generality in the unified analysis of OWMM.
The unified analysis of omni-directional wheeled mobile manipulator(OWMM) through the dimensionally nonhomogeneous Jacobian matrix may lead to unreliable results. The existing researches focus on the integrated performance evaluation of OWMM, without taking the influence of velocity difference into consideration. This paper presents a new approach to formulate the dimensionally homogeneous Jacobian matrix and a new index for OWMM. First, the universal transformational matrix of the link coordinate frame is derived based on double quaternion. The degree of locomotion of a mobile platform and the degrees of manipulation of a manipulator are treated equally as the joints of a redundant composite robot. Then the integrated modeling of OWMM is established, and the non-dimensional Jacobian matrix is obtained from the above matrix. Next, by examining the concept of directional manipulability, a new index is proposed to evaluate the directional manipulability of OWMM along the specified task direction. Furthermore, for the same task and the time of simulation, the kinematic performance of the fixed base operation of the manipulator and the operation of OWMM are analyzed by numerical simulation. The results suggest that the proposed approach is equivalent to the method of traditional kinematic modeling, with the error of numerical solution being 10–7, and the varying frequency of DMTR is higher than that of the condition number. This indicates that DMTR is more effective to reflect the variation of the task direction. The proposed method can be used to analyze the manipulating capability of OWMM, and has a simple structure and generality in the unified analysis of OWMM.
E-mail Alerts
