Experimental and Numerical Investigation on High-Pressure Centrifugal Pumps:Ultimate Pressure Formulation,Fatigue Life Assessment and Topological Optimization of Discharge Section

A high percentage of failure in pump elements originates from fatigue.This study focuses on the discharge section behavior,made of ductile iron,under dynamic load.An experimental protocol is established to collect the strain under pressurization and depressurization tests at specific locations.These experimental results are used to formulate the ultimate pressure expression function of the strain and the lateral surface of the discharge section and to validate finite element modeling.Fe-Safe is then used to assess the fatigue life cycle using different types of fatigue criteria(Coffin-Manson,Morrow,Goodman,and Soderberg).When the pressure is under 3000 PSI,pumps have an unlimited service life of 107 cycles,regardless of the criterion.However,for a pressure of 3555 PSI,only the Morrow criterion denotes a significant decrease in fatigue life cycles,as it considers the average stress.The topological optimization is then applied to the most critical pump model(with the lowest fatigue life cycle)to increase its fatigue life.Using the solid isotropic material with a penalization approach,the Abaqus Topology OptimizationModule is employed.The goal is to reduce the strain energy density while keeping the volume within bounds.According to the findings,a 5%volume reduction causes the strain energy density to decrease from 1.06 to 0.66106 J/m^(3).According to Morrow,the fatigue life cycle at 3,555 PSI is 782,425 longer than the initial 309,742 cycles.

Large Eddy Simulations of Rayleigh-Bénard Convection in a Cubic Cavity

This paper conducts a Large Eddy Simulation (LES) of Rayleigh Bénard convection in a cubic cavity based on the WMLES S-Omega subgrid-scale model. For a cubic cavity with a vertical temperature difference of 6.7°C and 20°C, the velocity pulsation profiles and the mean velocity profiles of the vertical section in the middle of the cubic cavity were simulated, respectively. And they are consistent with the experiment results. Furthermore, the mean velocity field of the vertical cross-section in the middle of the cavity was calculated. Structures of the mean velocity field in the two cases are similar. A counterclockwise large vortex is found to occupy the cavity, and there are two small clockwise vortices in the lower left and upper right corners, and the mean velocity fields at two different temperature differences are consistent with the experimental results. The two-dimensional instantaneous temperature field and mean temperature field with different cross-sections in the z-direction, as well as the three-dimensional instantaneous isothermal surface structure, indicate that the large-scale circulation motion within the cubic cavity is moving diagonally. In addition, the structure of the mean streamline also illustrates this viewpoint. For the reverse vortex formed at two corners in the mean streamline structure, we used the Q criterion to identify and obtain two vortex structures similar to boomerangs. The basic turbulent structure in RB thermal convection includes the rising and falling plumes generated by buoyancy effects.

Preliminary Model Study for Forecasting a Hot Weather Process in Guangdong Province Using CMA-TRAMS

In this paper, the CMA-TRAMS tropical high-resolution system was used to forecast a typical hot weather process in Guangdong, China with different horizontal resolutions and surface coverage. The results of resolutions of 0.02° and 0.06° were presented with the same surface coverage of the GlobeLand30 V2020, companies with the results of resolution 0.02° with the USGS global surface coverage. The results showed that, on the overall assessment the 2 km model performed better in forecasting 2 m temperature, while the 6 km model was more accurate in predicting 10 m wind speed. In the evaluation of representative stations, the 2 km model performed better in forecasting 2 m temperature and 2 m relative humidity at the coastal stations, and the 2 km model was also better in forecasting 2 m pressure at the representative stations. However, the 6 km model performed better in forecasting 10 m wind speed at the representative stations. Furthermore, the 2 km model, owing to its higher horizontal resolution, presented a more detailed stratification of various meteorological field maps, allowing for a more pronounced simulation of local meteorological element variations. And the use of the surface coverage data of the GlobeLand30 V2020 improved the forecasting of 2 m temperature, and 10 m wind speed compared to the USGS surface coverage data.

Hopf-flip bifurcation of high dimensional maps and application to vibro-impact systems

This paper addresses the problem of Hopf-flip bifurcation of high dimensional maps. Using the center manifold theorem, we obtain a three dimensional reduced map through the projection technique. The reduced map is further transformed into its normal form whose coefficients are determined by that of the original system. The dynamics of the map near the Hopf-flip bifurcation point is approximated by a so called “time-2τ^2 map” of a planar autonomous differential equation. It is shown that high dimensional maps may result in cycles of period two, tori T^1 （Hopf invariant circles）, tori 2T^1 and tori 2T^2 depending both on how the critical eigenvalues pass the unit circle and on the signs of resonant terms＇ coefficients. A two-degree-of-freedom vibro-impact system is given as an example to show how the procedure of this paper works. It reveals that through Hopf-flip bifurcations, periodic motions may lead directly to different types of motion, such as subharmonic motions, quasi-periodic motions, motions on high dimensional tori and even to chaotic motions depending both on change in direction of the parameter vector and on the nonlinear terms of the first three orders.

A hyperbolic Lindstedt-Poincare method for homoclinic motion of a kind of strongly nonlinear autonomous oscillators

A hyperbolic Lindstedt-Poincare method is presented to determine the homoclinic solutions of a kind of nonlinear oscillators, in which critical value of the homoclinic bifurcation parameter can be determined. The generalized Lienard oscillator is studied in detail, and the present method＇s predictions are compared with those of Runge-Kutta method to illustrate its accuracy.

CHAOTIC MOTIONS AND LIMIT CYCLE FLUTTER OF TWO-DIMENSIONAL WING IN SUPERSONIC FLOW

Based on the piston theory of supersonic flow and the energy method, the flutter motion equations of a two-dimensional wing with cubic stiffness in the pitching direction are established. The aeroelastic system contains both structural and aerodynamic nonlinearities. Hopf bifurcation theory is used to analyze the flutter speed of the system. The effects of system parameters on the flutter speed are studied. The 4th order Runge-Kutta method is used to calculate the stable limit cycle responses and chaotic motions of the aeroelastic system. Results show that the number and the stability of equilibrium points of the system vary with the increase of flow speed. Besides the simple limit cycle response of period 1, there are also period-doubling responses and chaotic motions in the flutter system. The route leading to chaos in the aeroelastic model used here is the period-doubling bifurcation. The chaotic motions in the system occur only when the flow speed is higher than the linear divergent speed and the initial condition is very small. Moreover, the flow speed regions in which the system behaves chaos axe very narrow.

Uniaxial Time-Dependent Ratcheting of SS304 Stainless Steel at High Temperatures

The uniaxial time-dependent strain cyclic behaviors and ratcheting of SS304 stainless steel were studied at high temperatures （350 ℃ and 700 ℃）. The effects of straining and stressing rates, holding time at the peak and/or valley of each cycle in addition to ambient temperature on the cyclic softening/hardening behavior and ratcheting of the material were discussed. It can be seen from experimental results that the material presents remarkable time dependence at 700 ℃, and the ratcheting strain depends greatly on the stressing rate, holding time and ambient temperature. Some significant conclusions are obtained, which are useful to build a constitutive model describiog the time-dependent cyclic deformation of the material.

Simulation of landslide run-out by considering frictional heating and thermal pressurization

Some of the remarkable characteristics of natural landslides, such as surprisingly long travel distances and high velocities, have been attributed to the mechanisms of frictional heating and thermal pressurization. In this work, this mechanism is combined with a depth-averaged model to simulate the long runout of landslides in the condition of deformation. Some important factors that influence frictional heating and thermal pressurization within the shear zone are further considered, including velocity profile and pressurization coefficient. In order to solve the coupled equations, a combined computational method based on the finite volume method and quadratic upwind interpolation for convective kinematics scheme is proposed. Several numerical tests are performed to demonstrate the feasibility of the computational scheme, the influence of thermal pressurization on landslide run-out, and the potential of the model to simulate an actual landslide.

NONLINEAR VIBRATION ANALYSIS OF VISCOELASTIC CABLE WITH SMALL SAG

Both the inplane and out-of-plane transverse vibrations of aviscoelastic cable subjected to an initial stress distributinguniform on the cross section are studied. The constitution of thecable material is as- sumed to be of the hereditary integral type.The partial differential -integral equations of motion are derivedfirst. Then by applying Galerkin's method, the governing equationsare reduced to a set of second-order non- linear differential-integral equations which are solved by finite difference numericalintegration procedures. Finally, the effects of the viscosityparameter and the elastic parameter on the transient amplitudes ofthe first Mode are investigated by numerical simulation.

Numerical Simulation and Experimental Investigation on Performance of the Wells Turbine in Irregular Oscillating Flow

The Wells turbine is an axial-flow air-turbine designed to extract energy from ocean waves. An important consideration is the self-starting capability of the Wells turbine, a phenomenon encountered where the turbine accelerate by itself up to a certain speed for the best turbine performance. In order to clarify the self-starting characteristic and running performance of the Wells turbine in an irregular oscillating flow, a numerical simulation process is established in this paper on the rational assumption of quasi-steady flow conditions, Both self-starting characteristics and running performance are obtained through the numerical simulation and subsequently compared with the experimental data achieved on a computer-controlled oscillating flow test rig which could realize some irregular oscillating flow according to the specified spectrum. Results show that the self-starting time decreases with the increase of the significant wave height and the mean frequency of the irregular oscillating flow, Therefore, it is possible to predict accurately the performance of the Wells turbine by computer simulation.

Comparison of Drags on Fish Cages of Different Shapes

The effects of the direction of current on the drag on fish cages are studied in the present paper. The drags on cages of different shapes, including cylindrical, tnmcated conical, cuboidal and hexagonal, are compared. The drag on the tnmcated conical net is smaller than that on other shapes of the same area. This net shape with a small apex angle is suggested for the design of fish cages.

Nonlinear flutter of a two-dimension thin plate subjected to aerodynamic heating by differential quadrature method

The problem of nonlinear aerothermoelasticity of a two-dimension thin plate in supersonic airflow is examined. The strain-displacement relation of the von Karman＇s large deflection theory is employed to describe the geometric non-linearity and the aerodynamic piston theory is employed to account for the effects of the aerodynamic force. A new method, the differential quadrature method （DQM）, is used to obtain the discrete form of the motion equations. Then the Runge-Kutta numerical method is applied to solve the nonlinear equations and the nonlinear response of the plate is obtained numerically. The results indicate that due to the aerodynamic heating, the plate stability is degenerated, and in a specific region of system parameters the chaos motion occurs, and the route to chaos motion is via doubling-period bifurcations.

CALCULATION OF THE DAMPING OF THE Zn-27Al ALLOY BASED ON THE MICRO INTERFACE SLIDING MODEL

The microstructures of the Zn-27Al alloy after modification, solid-solution treatment, and natural aging were studied. It was clarified why the damping properties of Zn-27Al alloys, after treatment, had advanced most on the basis of analyzing the microstructures. Approximate expressions have been educed, which can be used to quantificationally work out the damping of the Zn-27Al alloy on the basis of the micro interface sliding model. By comparing the testing damping properties of the foundry Zn-27Al alloys and the Zn-27Al alloys after modification, solid solution, and natural aging, it was shown that the expressions were rational.

In-line response of vertical cylinders in regular and random waves

The in-line response of a vertical flexibly mounted cylinder in regular and random waves is reported. Both theoretical analyses and experimental measurements have been performed. The theoretical predictions are based on the Morison equation which is solved by the incremental harmonic balance method. Experiments are then performed in a wave flume to determine the accuracy of the Morison equation in predicting the in-line response of the cylinder in regular and random waves. The interaction between waves and vibrating cylinders are investigated.

Experimental Research on Behavior of Composite Material Projectile Penetrating Concrete Target

Projectile made of carbon fiber composite material shell and metal warhead penetrates concrete target at speeds of 336,m/s,447,m/s and 517,m/s.The angles between the perpendicu-lar of target surface and projectile axis are 0°and 30°.The thickness of concrete target is 200,mm and the compression strength is 30 MPa.The experimental results indicate that the strength of composite material structure is high.Composite projectile can go through concrete tar-get without fiber segregation and breakage.The percent fill is 18.5% in the composite material projectile.It is about twice as that of metal projectile,if the density of metal is taken as 7.8,g/cm3.Comparing with metal projectile,low-density,high-strength composite material can lessen projec-tile weight,improve charge-weight ratio of detonator and enhance destructive powder.

UNCONVENTIONAL HAMILTON-TYPE VARIATIONAL PRINCIPLES FOR DYNAMICS OF REISSNER SANDWICH PLATE

According to the basic idea of classical yin-yang complementarity and modern dual-complementarity, in a simple and unified way proposed by Luo（1987）, some uncon ventional Hamilton-type variational principles for dynamics of Reissner sandwich plate can be established systematically. The unconventional Hamilton-type variation principle can fully characterize the initial boundary value problem of this dynamics. In this paper, an important integral relation is given, which can be considered as the generalized principle of virtual work in mechanics. Based on this relation, it is possible not only to obtain the principle of virtual work in dynamics of Reissner sandwich plate, but also to derive systematically the complementary functionals for fivefield, two-field and one-field unconventional Hamilton-type variational principles by the generalized Legender transformations. Furthermore, with this approach, the intrinsic relationship among the various principles can be explained clearly.

Wave Dissipation Characteristics of A Mountain-Type Breakwater

One mountain-type breakwater consisting of two inclined plates and one vertical plate is proposed based on several types of traditional free surface breakwaters, including the horizontal plate, curtain wall, and trapezoidal barriers. The interaction between the regular waves and the fixed free surface mountain-type breakwater is measured in one wave flume(15.0 m×0.6 m×0.7 m). The wave propagation, reflection, and transmission process are simulated using the VOF method and the hybrid SAS/laminar method. The simulated wave profiles are consistent with the experimental observations. For waves with a length smaller than four times width of the mountain-type breakwater, the reflected wave amplitudes are slightly larger than those of the vertical-plate breakwater, while the wave transmission coefficients are all smaller than 0.5, and the wave loss coefficients are larger than 0.7. The wave energy is dissipated by wave breaking on the windward inclined plate, and turbulent flow around the vertical plate and the leeward inclined plate.

3-D distribution of tensile stress in rock specimens for the Brazilian test

It is claimed that the formula used for calculating the tensile strength of a disk-shaped rock specimen in the Brazilian test is not accurate, because the formula is based on the 2-dimensional elastic theory and only suitable for very long or very short cylinders. The Matlab software was used to obtain the 2-dimensional distribution of stress in the rock specimen for Brazilian test. Then the 2-dimensional stress distribution in Brazilian disk was analyzed by the Marc FEM software. It can be found that the results obtained by the two software packages can verify each other. Finally, the 3-dimensional elastic stress in the specimen was calculated. The results demonstrate that the distribution of stress on the cross section of the specimen is similar to that in 2-dimension. However, the value of the stress on the cross section varies along the thickness of the specimen and the stress is bigger when getting closer to the end of the specimen. For the specimen with a height-to-diameter ratio of 1 and a Poisson＇s ratio of 0.25, the tensile strength calculated with the classical 2-D formula is 23.3% smaller than the real strength. Therefore, the classical 2-D formula is too conservative.

Phase Behaviors in Bi-phase Simulation of Powder Segregation in Metal Injection Molding

Powder segregation induced by mold filling is an important phenomenon that affects the final quality of metal injection molding （MIM）. The prediction of segregation in MIM requires a bi-phase flow model to describe distinctly the flows of metallic powder and polymer binder. Viscous behaviors for the flows of each phase should hence be determined. The coefficient of interaction between the flows of two phases should also be evaluated. However, only viscosity of the mixed feedstock is measurable by capillary tests. Wall sticking is supposed in the traditional model for capillary tests, while the wall slip is important to be taken into account in MIM injection. Objective of the present paper is to introduce the slip effect in bi-phase simulation, and search the suitable way to determine the viscous behaviors for each phase with the consideration of wall slip in capillary tests. Analytical and numerical methods were proposed to realize such a specific purpose. The proposed method is based on the mass conservation between the capillary flows in mono-phase model for the mixed feedstock and in bi-phase model for the flows of two phases. Examples of the bi-phase simulation in MIM were realized with the software developed by research team. The results show evident segregation, which is valuable for improving the mould designs.

Application of Reproducing Kernel Particle Method in an Analysis of Elasto-plastic Deformation Under Taylor Impact

The Reproducing Kernel Particle Method (RKPM) is one of several new meshless numerical methods de- veloped internationally in recent years. The ideal elasto-plastic constitutive model of material under a Taylor impact is characterized by the Jaumann stress- and strain-rates. An updated Lagrangian format is used for the calculation in a nu- merical analysis. With the RKPM, this paper deals with the calculation model for the Taylor impact and deduces the control equation for the impact process. A program was developed to simulate numerically the Taylor impact of projec- tiles composed of several kinds of material. The simulation result is in good accordance with both the test results and the Taylor analysis outcome. Since the meshless method is not limited by meshes, it is believed to be widely applicable to such complicated processes as the Taylor impact, including large deformation and strain and to the study of the dy- namic qualities of materials.