低温透平膨胀机用全动压气体轴承的设计与试验研究

介绍 15 0m3 /h制氧机用中压低温透平膨胀机 (PLK - 8 33× 2 / 2 0 - 5 )在转子转速12 16× 10 4r/min (比额定超速 14 % )工况下 ,获得了全动压气体轴承稳定运转的良好效果。还对螺旋槽止推轴承的性能进行了研究。试验结果表明 ,在国产高速低温透平膨胀机中改用全动压气体轴承是完全可行的。     

Numerical study on the spontaneous condensation flow in an air cryogenic turbo-expander using equilibrium and non-equilibrium models

The difficulty of data measurement in cryogenic environments and the complicated mechanism of nucleation process have restricted the design of wet type turbo-expander for cryogenic liquid plants. In this paper, equilibrium and non-equilibrium models are used to model the spontaneous condensation flow in a cryogenic turbo-expander along the main stream passage including nozzle, impeller and diffuser. The comparison shows a distinct difference of the predicted wetness fraction distribution along the streamline between the equilibrium model and the non-equilibrium model. In non-equilibrium model, the distributions of supercooling and nucleation rate along the length of turbo-expander are given for the analysis of flow characteristics. The comparison of outlet wetness fraction with the experimental data is also provided for verification and discussion. Both the effects of the rotation on nucleation and the effects of the nucleation on flow along suction side of the impeller are investigated.     

Simulation and experimental studies of the repaired composite laminate

Firstly the compress experiment is undertaken to investigate the efficiency of repaired panels in this paper, and then modeling of the mechanical behavior of the repaired composite panel under compressive static load is conducted by using of the finite element method. The effect of geometric non‐linearity on the stress‐strain response is considered in the numeric analysis. Fatherly, the user material subroutine (UMAT) is integrated with the ABAQUS package with the geometric non‐linearity effect for studying the damage initiation and its progression in the composite structure, and quadrilateral, linear, thick shell elements (S8R) are adopted. Finally, the predicted strain distribution, damage evolution and strength of the laminate are compared with the test results.     

Numerical and experimental analysis of residual stresses for fatigue crack growth

In this paper computational and experimental results are presented concerning residual stress effects on fatigue crack growth in a Compact Tension Shear (CTS) specimen under cyclic mode I loading. For a crack of constant length it is found that hardly any compressive residual stresses or crack closure effects are generated along the crack surfaces behind the crack tip through the considered cyclic mode I loading with a load ratio of R=0.1. Only if fatigue crack growth is modelled during the simulation of the cyclic loading process these well-known effects are found. On the other hand it is shown that they have hardly any influence on the residual stresses ahead of the crack tip and thus on further fatigue crack growth. For all cases considered the computational finite element results agree well with the experimental findings obtained through X-ray diffraction techniques.     

Numerical and experimental studies on aluminum sandwich plates of variable thickness

Abstract

The elastic flexural behavior of static deformation and free vibration of sandwich plates of variable thickness is investigated numerically and experimentally. In the analysis, the face plates are treated as Marguerre shells, and the core is assumed to be an antiplane core and to provide resistance to transverse shear and normal stresses only. Displacement continuity conditions are used at the interfaces between face plates and the core to derive the displacement field. Energy formulations are obtained and solved by the isoparametric finite element method. The numerical results are obtained to compare with the results in the existing literature and to show the effects of taper constant and face plate thickness on deflections and natural frequencies. Finally, experimental works based on the method of holographic interferometry are conducted to confirm the theoretical findings. Experimental and numerical data agree quite well in this work.     

Numerical predictions and experimental measurements of residual stresses in fatigue crack growth specimens

Controlling macro residual stress fields in a material while preserving a desired microstructure is often a challenging proposition. Processing techniques which induce or reduce residual stresses often also alter microstructural characteristics of the material through thermo-mechanical processes. A novel mechanical technique able to generate controlled residual stresses was developed. The method is based on a pin compression approach, and was used to produce well-controlled magnitudes and distributions of residual stresses in rectangular coupons and compact tension specimens typically used in fatigue crack growth testing. Residual stresses created through this method were first computationally modeled with finite element analysis, and then experimentally reproduced with various levels of pin compression. The magnitudes and distributions of residual stresses in experimental specimens were independently assessed with fracture mechanics methods and good correspondence was found between residual stresses produced using the pin compression and processing techniques. Fatigue crack growth data generated from specimens with low residual stresses, high residual stresses resulting from processing, and high residual stresses introduced through the new pin compression technique were compared and validated. The developed method is proposed to facilitate the acquisition and analysis of fatigue crack growth data generated in residual stresses, validate residual stress corrective models, and verify fatigue crack growth simulations and life predictions in the presence of residual stresses.     

Numerical studies on the effective shear modulus of particle reinforced composites with an inhomogeneous inter-phase

Numerical studies on the effective shear modulus of particle reinforced composites with an inhomogeneous inter-phase are performed. The influences of many parameters to the equivalent shear modulus of composites are carefully analyzed, including the inter-phase thickness, variation of interfacial properties, boundary conditions and volume fraction of particles, etc. Numerical results show that the Poisson ratio can be assumed as a constant across the whole inter-phase zone in the computation. The form of property variation across the inter-phase also greatly affects the effective shear modulus of composite. Numerical results predicted by the rigid boundary conditions are remarkably higher than those by the free boundary conditions and the exact solutions. The reasonability and exactness of the available models for predicting the effective shear modulus of composites are accessed by the numerical results in the present work.     

分舱段圆柱壳声散射数值和试验研究

基于有限元方法建立了填充不同介质的有限长分舱段圆柱壳声散射数值计算模型,仿真了填充空气-空气、空气-水、水-水三类两舱段圆柱壳声散射特性,并完成了三类两舱段圆柱壳体声散射试验,获取和分析了两舱段圆柱壳体声散射的时间角度谱和频率角度谱特性。利用物理声学方法分析了壳体表面、端面以及内部填充水介质对散射声场的影响,揭示了两舱段圆柱壳声散射频率角度谱中呈现的干涉条纹特征形成机理,为水下分舱段目标,如水下无人航行器的主动声呐探测和识别提供理论支撑。     

Integrated approach for prediction of stability limits for machining with large volumes of material removal

High-speed machining of thin-walled structures is widely used in the aeronautical industry. Higher spindle speed and machining feed rate, combined with a greater depth of cut, increases the removal rate and with it, productivity. The combination of higher spindle speed and depth of cut makes instabilities (chatter) a far more significant concern. Chatter causes reduced surface quality and accelerated tool wear. Since chatter is so prevalent, traditional cutting parameters and processes are frequently rendered ineffective and inaccurate. For the machine tool to reach its full utility, the chatter vibrations must be identified and avoided. In order to avoid chatter and implement optimum cutting parameters, the machine tool including all components and the work piece must be dynamically mapped to identify vibration characteristics. The aim of the presented work is to develop a model for the prediction of stability limits as a function of process parameters. The model consists of experimentally measured vibration properties of the spindle-tool, and finite element calculations of the work piece in (three) different stages of the process. Commercial software packages used for integration into the model prove to accomplish demands for functionality and performance. A reference geometry that is typical for an aircraft detail is used for evaluation of the prediction methodology. In order to validate the model, the stability limits predicted by the use of numerical simulation are compared with the results based on the experimental work.     

Numerical implementation of imperfect interfaces

One approach to characterizing interfacial stiffness is to introduce imperfect interfaces that allow displacement discontinuities whose magnitudes depend on interfacial traction and on properties of the interface or interphase region. This work implemented such imperfect interfaces into both finite element analysis and the material point method. The finite element approach defined imperfect interface elements that are compatible with static, linear finite element analysis. The material point method interfaces extended prior contact methods to include interfaces with arbitrary traction-displacement laws. The numerical methods were validated by comparison to new or existing stress transfer models for composites with imperfect interfaces. Some possible experiments for measuring the imperfect interface parameters needed for modeling are discussed.     

Numerical investigation on the failure phenomena of stiffened composite panels in post-buckling regime with discrete damages

The aim of this work is to investigate the final failure response of damaged composite stiffened panels in post buckling regime under compressive load, by using progressive failure analysis (PFA) methodology. The selected panel is characterized by T shaped stringers and it is representative of the upper skin panel, toward the wing tip, of the wing box of a typical regional aircraft. PFA methodology has been applied in order to predict in addition to the initiation of the local failure also its propagation up to the final collapse of the panel, in presence of local damage (barely visible impact damage, BVID) and in post-buckling regime. For this purpose, discrete damages have been considered in the skin of the panel. According to the indications contained in many guidelines finalized to the preliminary design of composite structures, a simplified design model of BVID has been considered in this work, in particular a hole 1/4 in. in diameter has been used to simulate this damage. The collapse load of the panel has been evaluated considering different locations of a single damage and also considering multi-damage maps (the latter are more representative of a real damage scenario). The results of PFA presented in this work illustrate the combined effect of the reduction of the panel stiffness and of the damage propagation, and the sensitivity of the buckling onset and of the residual strength of the panel with respect to different damage locations and damage density.     

Numerical analysis and experimental correlation of composite shell buckling and post-buckling

This paper deals with the buckling and post-buckling behaviour of carbon fibre reinforced plastic cylindrical shells under axial compression. The finite element analysis is used to investigate this problem and three different types of analysis are compared: eigenvalue analysis, non-linear Riks method and dynamic analysis. The effect of geometric imperfection shape and amplitude on critical loads is discussed. A numerical–experimental correlation is performed, using the results of experimental buckling tests. The geometric imperfections measured on the real specimens are accounted for in the finite element model. The results show the reliability of the method to follow the evolution of the cylinder shape from the buckling to the post-buckling field and good accuracy in reproducing the experimental post-buckling behaviour.     

Numerical and experimental study on morphing bi-stable composite laminates actuated by a heating method

This paper presents a new heating actuation method to induce the snap-through phenomenon of bi-stable laminates by manipulating the residual stress in laminates. The mechanism of the heating actuation method is analyzed, and the snap-through process is simulated by the finite element method. The heating actuation experiments on two types of laminates with different stacking sequences are performed. Good agreement is obtained between the experimental data and the finite element analysis (FEA). Subsequently, the heating actuation method is applied to several different bi-stable laminates and investigated by FEA. The FEA results show that this heating actuation method is an effective way to induce the snap-through phenomenon of bi-stable laminates of different thicknesses, sizes and shapes.     

Numerical simulation of crack problems using triangular finite elements with embedded interfaces

The aim of this paper is to propose numerical aspects for the modeling of discrete cracks in quasi-brittle materials using triangular finite elements with an embedded interface based on the formulation in [Computational Mechanics 27 (2001) 463]. The kinematics of the discontinuous displacement field and the variational formulation applied to a body with an internal discontinuity is given. The discontinuity is modeled by additional global degrees of freedom and the continuity of the displacement jumps across the element boundaries is enforced. To show the performance of the model, a single element test and two examples for mode-I dominated fracture, namely a tension test and a three-point bending beam, are discussed.     

Numerical and experimental investigation on lap shear fracture of Al/CFRP laminates

This paper presents a new approach to numerically investigate the lap shear fracture of a hybrid laminate made of Carbon Fibre Reinforced Plastic (CFRP) and metal foil plies (e.g. aluminium), validated by corresponding experiments. The numerical Finite Element (FE) model of the hybrid laminate, subjected to lap shear fracture, is composed of five laminas with alternating metal/CFRP layers with cohesive elements lying within Al/CFRP interface. In the FE model, individual CFRP laminas are assumed as an orthotropic homogenized continuum under plane stress, and aluminium facesheets are modelled as an elastic–plastic continuum. The Al/CFRP interface is represented via quadratic cohesive elements, the constitutive law of which is an exponentially decaying law representing the degrading behaviour of the interface (implemented as user element in ABAQUS). The numerical model captures the experimentally obtained results with minimal error, and predicts the failure modes successfully. The influence of specimen geometry (e.g. overlap length, total length, and total width) on lap shear fracture response is analyzed in detail in this study, too, in order to confirm the specimen design for the test, as there is still no corresponding test standard for hybrid laminates.     

Finite element simulation and experimental research on ZnO:Al by magnetron sputtering

Aluminum doped zinc oxide (ZnO:Al) thin films are suitable for the use as transparent conductive electrode in copper indium gallium selenide Cu(In,Ga)Se₂ thin film solar cells. The resistivity and film quality of ZnO:Al deposited on soda lime glass is nonuniform in magnetron sputtering process. According to the measurement results of magnetic field on the top of the target, obvious magnetic field distribution nonuniformity is observed along the vertical and horizontal directions respectively. With the longer distance between target and substrate, the magnetic field intensity becomes lower and flatter between the two magnet poles. Based on the simulation results by finite element analysis, it is verified the nonuniformity of magnetic field distribution influences the probability of Ar⁺ particles collision and the deposition of zinc oxide (ZnO) particles in different regions on substrate. The higher resistivity of ZnO:Al films is obtained where the magnetic field intensity is stronger.     

双环减速器辐射噪声数值仿真及试验研究

以中心输入式双环减速器为研究对象,综合考虑齿轮内部动态激励以及环板不平衡惯性激励,建立了减速器传动系统及结构系统的动力有限元分析模型,应用ANSYS软件对双环减速器进行固有模态及动态响应数值仿真。以振动位移作为边界条件,建立减速器箱体的声学边界元分析模型,在SYSNOISE软件中用直接边界元法计算了箱体表面声压及场点的辐射噪声。利用传动系统试验台对双环减速器进行振动加速度及辐射噪声测试,并与数值仿真结果对比分析,两者吻合良好。     

Numerical evaluation for debonding failure phenomenon of adhesively bonded joints at cryogenic temperatures

In the present study, material characteristics, such as inelastic constitutive behaviour and debonding failure, of an adhesively bonded joint (ABJ) at cryogenic temperature have been evaluated using a computational approach. The modified Bodner-Partom model (BP model) has been introduced to describe the material nonlinearities of ABJ. The Gurson-Tvergaard model (GT model) has also been implemented into the constitutive model in order to analyse the phenomenon of debonding failure. An ABAQUS user-defined subroutine UMAT is developed using a damage-coupled constitutive model based on an implicit formulation. The numerical results are compared with a series of lap shear tests of ABJ at cryogenic temperature in order to verify the proposed method.     

Numerical solution of elastic bodies in contact by FEM utilising equilibrium displacement fields

In this paper a new effective formulation of the computation of the statics of elastic bodies in contact is described and demonstrated. Line contact, plane strain, and negligible friction are assumed. The formulation is an extension of the standard finite element method (FEM). With the aim to utilise the Hertz theory directly, we use the exact solution of the elastic 2-dimensional halfspace loaded by Hertzian pressure distribution and enforce the contact condition by the method of Lagrange multipliers. In numerical examples we have focussed on the demonstration and evaluation of the accuracy of the new formulation for selected applications compared with the state of the art node-to-segment contact algorithm implemented in the software system ADINA. Proposed formulation is more accurate for problems where Hertz contact dominates the strain state, especially for small number of elements, whereas we obtained a fairly good agreement with ADINA for a more general bending problem.     

Numerical/experimental impact events on filament wound composite pressure vessel

Impacts on pressure vessels, produced by winding glass fibre with vinyl ester resin over a polyethylene liner, were numerically and experimentally investigated in the current work.Pressure vessels were experimentally tested under low velocity impact loads. Different locations and incident energies were tested in order to evaluate the induced damage and the capability of the developed numerical model.An advanced 3-D FE model was used for simulating the impact events. It is based on the combined use of interlaminar and intralaminar damage models. Puck and Hashin failure theories were used to evaluate the intralaminar damages (matrix cracking and fibre failure). Cohesive zone theory, by mean of cohesive elements, was used for modelling delamination onset and propagation.The experimental impact curves were accurately predicted by the numerical model for the different impact locations and energies. The overall damages, both intralaminar and interlaminar, were instead slightly over predicted for all the configurations.The model capabilities to simulate the low velocity impact events on the full scale composite structures were proved.