Molecular simulation study of the stabilization process of NEPE propellant

In this reported study, the density functional theory(DFT) was used at the(U)B3LYP/6-311G(d,p) level to investigate the stabilization process of the nitrate ester plasticized polyether propellant(NEPE). Molecular simulations were conducted of the reaction that generates NO_(2), the autocatalytic and aging reaction triggered by the NO_(2), and the nitrogen dioxide absorption reaction of the stabilizers during the propellent stabilization process. These simulations were derived using the transition-state theory(TST)and variational transition-state theory(VTST). The simulation results suggested that the stabilization of the NEPE propellant consisted of three stages. First, heat and NO_(2) were generated during the denitrification reaction of nitroglycerine(NG) and 1,2,4-butanetriol trinitrate(BTTN) in the NEPE propellant.Second, nitroso products were generated by the reactions of N-Methyl-4-nitroaniline(MNA) and 2-nitrodiphenylamine(2NDPA) with NO_(2). Third, the stabilizers were exhausted and the autocatalytic reaction of NG and BTTN and the aging reaction of polyethylene glycol(PEG) were triggered by the heat and NO_(2)generated in the first stage. By comparing the energy barriers of the various reactions, it was found that the NO_(2)generated from the denitrification reaction significantly reduced the reaction energy barrier to 105.56-126.32 kJ/mol, also increased the reaction rate constant, and decreased the thermal stability and energetic properties of the NEPE propellant. In addition, the NO_(2)also weakened the mechanical properties of the NEPE propellant by attacking the-CH2groups and the O atoms in the PEG molecular chain. The energy barriers of the reactions of MNA and 2NDPA with NO_(2)(94.61-133.61 k J/mol) were lower than those of the autocatalytic and decomposition reactions of NG, BTTN, and the aging reactions of PEG(160.30-279.46 kJ/mol). This indicated that, by eliminating NO_(2), the stabilizer in the NEPE propellant can effectively prevent NO_(2)from reacting with the NG, BTTN, and PEG in the NEPE propellant. Consequently, this would help maintain the energy and mechanical properties of the NEPE propellant, thereby improving its thermal stability.

FEM Analysis of Spring-backs in Age Forming of Aluminum Alloy Plates

The age forming technology, characterized by huge spring-backs, has been developed to manufacture large integral wing-skin panel parts, which necessitates devising a method of predicting spring-backs. A 7B04-T7451 aluminum alloy creep test in tension is accomplished at 155 ℃, and the creep curves are obtained. The material constants of the mechanism-based creep constitutive equations are determined through experiments. The age forming process and the spring-backs of 7B04 aluminum alloy plates are analyzed using the commercial finite element software ABAQUS. The effects of plate thickness and forming time on spring-backs are researched. The spring-backs decrease with the increase of plate thickness and forming time. The test results verify the reliability of the finite element method （FEM） analysis.

A CUMULATIVE FATIGUE DAMAGE RULE UNDER THE ALTERNATIVE OF CORROSION OR CYCLIC LOADING

Fatigue damage increases with the applied loading cycles in a cumulative manner and the material deteriorates with the corrosion time. A cumulative fatigue damage rule under the alternative of corrosion or cyclic loading was proposed. The specimens of aluminum alloy LY12-CZ soaked in corrosive liquid for different times were tested under the constant amplitude cyclic loading to obtain S-N curves. The test was carried out to verify the proposed cumulative fatigue damage rule under the different combinations among corrosion time, loading level, and the cycle numbers. It was shown that the predicted residual fatigue lives showed a good agreement with the experimental results and the proposed rule was simple and can be easily adopted.

MECHANISM ANALYSIS OF THICKNESS EFFECT ON MIXED MODE Ⅰ/Ⅱ FRACTURE OF LC4-CS ALUMINUM ALLOY

Mixed mode Ⅰ/Ⅱ fracture erperiments of LC4-CS aluminum alloy were conductedby using tension--shear specimens with thicknesses of 2, 4, 8 and 14mm. Fracturemechanisms of thickness effect on mixed mode Ⅰ/Ⅱ fracture were first examined fromfracture surface morphology to correlate with the macroscopic fracture behavior andstress state. It is found that specimen thickness has a strong influence on mixed modefracture. As thickness varies from thin to thick the macroscopic fracture surfacesappear the characteristics of plane stress state (2mm, 4mm--thick specimen), three--dimensional stress state (8mm--thick specimens), and plane strain state (14mm--thickspecimens), respectively. The specimens of all kinds of thicknesses are typical of ten-sile type failure under mode Ⅰ loading condition and shear type failure under mode Ⅱloading condition. Two distinct features coexist on the fracture surfaces under mixedmode loading conditions, and the corresponding proportion varies with loading mix-ity. Void--growth processes are the failure mechanism in both predominately tensile-and shears--type fractures. The size and depth of dimples on the fracture surface varygreatly with thickness. Therefore, it is extraordinary necessary to take into accountthe thickness effect when a mixed mode fracture criterion is being established.

CALCULATION OF ELASTIC DAMPING CHARACTERISTICS OF ROTOR SUPPORT MADE OF METAL RUBBER MATERIAL UNDER VARIABLE LOADS

A metal rubber（MR） dry friction damper was designed based on the load supported by the rotor. An experimental apparatus for obtaining hysteresis loops of support under the precession load was designed. The elastic-damping characteristics of the ring-shaped MR damper used as a rotor support under variable loads were presented by studying the hysteresis loops of the damper. The vibration rigidity and the energy dissipation coefficient were calculated from the hysteresis loops, based on the description of the deformation process of the MR element with simple structure in a dimensionless coordinating system. The calculation results showed that the energy dissipation coefficient in the inner of MR element and on the boundary between the damper and the frame of the rotor support were approximately equal. The comparison of the hysteresis loops for a precession load and a one-axial load indicated a large difference when the coefficient of the energy dissipation and the stiffness of the MR damper were concerned.

Mechanical properties and microstructure of composites produced by recycling metal chips

In this study, the processing and mechanical properties of porous metal matrix composites(MMCs) composed of spheroidal cast iron chips(GGG40) and bronze chips(CuSn10) and formed by hot isostatic pressing were investigated. Bronze chips(CuSn10) were used as a matrix component, and spheroidal cast iron(GGG40) chips were used as a reinforcement component. The MMCs were produced with different CuSn10 contents(90 wt%, 80 wt%, 70 wt%, and 60 wt%). The hot isostatic pressing process was performed under three different pressures and temperatures. The produced MMCs were characterized using density tests, Brinell hardness tests, and compression tests. In addition, the consolidation mechanism was investigated by X-ray diffraction(XRD) analysis and scanning electron microscopy. The test results were compared with those for bulk CuSn10 and bulk GGG40. Mechanical tests results revealed that the metallic chips can be recycled by using hot pressing and that the mechanical properties of the produced MMCs were similar to those of bulk CuSn10. XRD and microscopy studies showed that no intermetallic compounds formed between the metallic chips. The results showed that the CuSn10 and GGG40 chips were consolidated by mechanical interlocking.

FUZZY ARITHMETIC AND SOLVING OF THE STATIC GOVERNING EQUATIONS OF FUZZY FINITE ELEMENT METHOD

The key component of finite element analysis of structures with fuzzy parameters, which is associated with handling of some fuzzy information and arithmetic relation of fuzzy variables, was the solving of the governing equations of fuzzy finite element method. Based on a given interval representation of fuzzy numbers, some arithmetic rules of fuzzy numbers and fuzzy variables were developed in terms of the properties of interval arithmetic. According to the rules and by the theory of interval finite element method, procedures for solving the static governing equations of fuzzy finite element method of structures were presented. By the proposed procedure, the possibility distributions of responses of fuzzy structures can be generated in terms of the membership functions of the input fuzzy numbers. It is shown by a numerical example that the computational burden of the presented procedures is low and easy to implement. The effectiveness and usefulness of the presented procedures are also illustrated.

Dynamic Simulation and Hemolysis Evaluation of the Regurgitant Flow over a Tilting-Disc Mechanical Heart Valve in Pulsatile Flow

Regurgitation in the heart diastolic phase represents a critical flow condition associated with many heart valve design considerations. The finite volume method, the Low-Reynolds-Number k-ω turbulent model and sliding mesh model are employed to solve and compare the complex flow field and the torque in each case. The end results expected from a cardiovascular CFD analysis are not limited only to the flowfield investigations. More importantly, it needs an evaluation criterion to judge if the design is acceptable as considered from a broader blood cell damage or activation perspective. In this study, blood cell damage index developed based on stress-time empirical rule and Lagrangian particle tracking is introduced to assess the viscous and turbulence-induced stresses effect to the blood cells.

Hot Cracking Susceptibility of 800H and 825 Nickel-Base Superalloys during Welding via Spot Varestraint Test

Hot cracking susceptibility of fillers 52 and 82 in 800H and 825 nickel-base superalloys was discussed using the Spot Varestraint test.The fillers of 52 and 82 were added into nickel-base superalloys via a gas tungsten arc welding(GTAW).Experimental results showed that the hot cracking sensitivity of the nickel-base superalloys with filler at high temperature was lower than that without filler.The hot cracking sensitivity had a slight effect when the filler 82 was added.The total length of crack was increased,the liquid-solid(L-S)two-phase range is higher so that the hot cracking susceptibility will be raised.The morphologies of cracks included the intergranular crack in the molten pool,molten pool of solidification cracking,heat-affected zone of intergranular cracks,and transgranular crack in the heat-affected zone.

NUMERICAL SOLUTION OF THE UNSTEADY EULER EQUATIONS FOR AIRFOILS USING APPROXIMATE BOUNDARY CONDITIONS

This paper presents an efficient numerical method for solving the unsteady Euler equations on stationary rectilinear grids. Boundary conditions on the surface of an airfoil are implemented by using their first-order expansions on the mean chord line. The method is not restricted to flows with small disturbances since there are no restrictions on the mean angle of attack of the airfoil. The mathematical formulation and the numerical implementation of the wall boundary conditions in a fully implicit time-accurate finite-volume Euler scheme are described. Unsteady transonic flows about an oscillating NACA 0012 airfoil are calculated. Computational results compare well with Euler solutions by the full boundary conditions on a body-fitted curvilinear grid and published experimental data. This study establishes the feasibility for computing unsteady fluid-structure interaction problems, where the use of a stationary rectilinear grid offers substantial advantages in saving computer time and program design since it does not require the generation and implementation of time-dependent body-fitted grids.

Numerical Simulation of Peen Forming

Peen forming is the most important process for the manufacturing of integral wing skin panels. In order to determine peening parameters, an equivalent deformation theory was proposed in which peening parameters are transformed into equivalent nodal forces. In this study, a linear elastic shell finite element method is adopted. Other related problems such as wing skin panel modeling and peening scheme analysis were also outlined. The method has been applied to manufacturing panels of a certain Chinese aircraft and was proved to be effective in reducing peening experimental tests and improving products quality.

PROCEDURE FOR COMPUTING THE POSSIBILITY AND FUZZY PROBABILITY OF FAILURE OF STRUCTURES

Traditionally, the calculation of reliability of fuzzy random structures is based on the well-known formulation of probability of fuzzy events. But sometimes the results of this formulation will not indicating the real state of safety of fuzzy-random structures. Based on the possibility theory, a computational procedure for the reliability analysis of fuzzy failure problems and random-fuzzy failure problems of mechanical structures that contain fuzzy variables were presented. A procedure for the analysis of structural reliability of problems of fuzzy failure criterion was also proposed. The failure possibility of fuzzy structures and possibility distribution of the probability of failure of fuzzy-random structures can be given by the proposed methods. It is shown that for the hybrid probabilistic and fuzzy reliability problems, the probability of failure should be suitably taken as a fuzzy variable in order to indicate the real safety of system objectively. Two examples illustrate the validity and rationality of the proposed methods.

COMPUTATION OF FIELD STRUCTURE AND AERODYNAMIC CHARACTERISTICS OF DELTA WINGS AT HIGH ANGLES OF ATTACK

A numerical investigation of the structure of the vortical flowfield over delta wings at high angles of attack in longitudinal and with small sideslip angle is presented. Three-dimensional Navier-Stokes numerical simulations were carried out to predict the complex leeward-side flowfield characteristics that are dominated by the effect of the breakdown of the leading-edge vortices. The methods that analyze the flowfield structure quantitatively were given by using flowfield data from the computational results. In the region before the vortex breakdown, the vortex axes are approximated as being straight line. As the angle of attack increases, the vortex axes are closer to the root chord, and farther away from the wing surface. Along the vortex axes, as the adverse pressure gradients occur, the axial velocity decreases, that is, A is negativee, so the vortex is unstable, and it is possible to breakdown. The occurrence of the breakdown results in the instability of lateral motion for a delta wing, and the lateral moment diverges after a small perturbation occurs at high angles of attack. However, after a critical angle of attack is reached the vortices breakdown completely at the wing apex, and the instability resulting from the vortex breakdown disappears.

ANALYSIS OF MICROBUCKLING FOR MONOMOLECULAR LAYERS ADHERING TO A SUBSTRATE

A two-dimensional linear spring model is established to study the microbuckling of a plane monomolecular layer adhering to a substrate. The model is for the layer subjected to a compressive load having an arbitrary angle with the chemical bond of the layer. The effects of the load angle, the strength of adhesion and the bending stiffness and shearing stiffness (the capability of resisting transverse bending and in-plane shearing) of the layer on the minimal buckling force and the critical buckling mode are discussed. It is found that the minimal buckling force increases with increasing load angle and, for a given bending stiffness, increases with increasing strength of adhesion and decreasing shearing stiffness. Furthermore, a critical condition under which the buckling of the layer can just occur is obtained, which is helpful to avoid buckling in an engineering application.

LARGE AMPLITUDE FREE VIBRATIONS OF LAMINATED COMPOSITE PLATES

This paper deals with large amplitude free flexural vibrations of laminated composite plates using a 9-node Heterosis degenerated isoparametric quadrilateral element, including the effects of transverse shear and rotary inertia. The nonlinear dynamic equations of the plates are formulated in von Karman's sense. Amplitude-frequency relationships are obtained through dynamic response history using; the Newmark numerical integration scheme. Detailed numerical results based on various parameters are presented for orthotropic laminated plates with different boundary conditions. The rectangular anti-symmetric cross-ply plates show the softening type of nonlinearity for initial small amplitudes. The displacement amplitudes decrease and nonlinear frequencies increase with the increment of time.

声-超声技术在碳-碳复合材料薄板损伤检测中的应用

根据声 -超声基本原理 ,以普通压电陶瓷传感器替代商用声发射超声换能器 ,采用传感器对称布置方式 ,对具有预冲击脱层损伤的碳 -碳复合材料薄板做了声 -超声测试。在对所得试验信号进行小波滤波后 ,应用了能量积分方法量化应力波因子 ,据此值可以区分材料有无损伤 。

基于约简遗传规划的线参数模型及在航空发动机起动建模中的应用(英文)

提出了一种新的约简遗传规划(PGP)算法和一种新的基于约简遗传规划的航空发动机起动动态线参数模型。这种模型采用遗传规划产生航空发动机起动模型的输入输出非线性模型集,并以二叉树结构表征函数项,运用正交最小二乘算法(OLS)估计二叉树分支(基本函数项)对于模型精度的贡献并去除复杂、冗余的函数项,从而加快遗传规划的收敛速度,最后通过 GP 进化可获得简单、可靠、准确的线参数非线性模型。发动机起动过程试车数据建模和与支持向量机的比较证明,这种方法可以产生适用性好、解析性强的线参数非线性模型,产生的模型可获得与支持向量机相当甚至更优的结果。

Highly maneuvering target tracking using multi-parameter fusion Singer model

An algorithm of highly maneuvering target tracking is proposed to solve the problem of large tracking error caused by strong maneuver. In this algorithm, a new estimator, named as multi-parameter fusion Singer (MF-Singer) model is derived based on the Singer model and the fuzzy reasoning method by using radial acceleration and velocity of the target, and applied to the problem of maneuvering target tracking in strong maneuvering environment and operating environment. The tracking performance of the MF-Singer model is evaluated and compared with other manuevering tracking models. It is shown that the MF-Singer model outperforms these algorithms in several examples.

INTERVAL ARITHMETIC AND STATIC INTERVAL FINITE ELEMENT METHOD

When the uncertainties of structures may be bounded in intervals, through some suitable discretization, interval finite element method can be constructed by combining the interval analysis with the traditional finite element method (FEM). The two parameters, median and deviation, were used to represent the uncertainties of interval variables. Based on the arithmetic rules of intervals, some properties and arithmetic rules of interval variables were demonstrated. Combining the procedure of interval analysis with FEM, a static linear interval finite element method was presented to solve the non-random uncertain structures. ne solving of the characteristic parameters of n-freedom uncertain displacement field of the static governing equation was transformed into 2 n-order linear equations. It is shown by a numerical example that the proposed method is practical and effective.