应用桥联理论与有限元模型的航空复合材料结构失效分析

桥联理论是一个近年来发展迅速的细观力学理论,它仅需纤维和基体材料的力学性能即可预报复合材料的响应。考虑基体应力集中系数与界面脱粘的影响后,大幅提高了复合材料破坏与强度的预报精确度。本文通过商业有限元（FE）软件ABAQUS的二次开发功能,将桥联模型的最新发展编写入UGENS子程序,并对复合材料单层板、层合板以及航空航天工业常见的T型接头复杂结构进行了有限元强度分析,并与实验结果进行对比,结果吻合良好。     

基于能量法的跨声速风扇叶片气弹稳定性研究

基于能量法原理,采用弱耦合的方法对跨声速风扇叶片进行气弹稳定性分析;采用3维线性插值算法编程实现CSD/CFD数据交换。利用有限元法计算叶片的模态振型和固有振动频率;应用FLUENT动网格技术,对NASA R 67风扇叶片在不同模态振动下的非定常流场进行数值模拟,给出在前3阶模态振动下叶片表面的周期累积气动功和气动阻尼,并探讨了激波对跨声速风扇叶片气弹稳定性的影响。     

高阶有限元格式研究的一个新途径

 在作者提出的有限差分－有限元混合方法的基础上，将差分法中的三阶无波动、无自由参数的ENN格式较成功地推广应用到有限元法中，为发展高精度有限元格式提供了一个新途径。     

A novel surrogate modeling strategy of the mechanical properties of 3D braided composites

In this paper, a surrogate-based modeling methodology is developed and presented to predict the elastic properties of three dimensional (3D) four-directional braided composites. Using this approach, the prediction process becomes feasible with only a limited number of training points. The surrogate models constructed using Finite Element (FE) method and Diffuse Approximation, reduce the computational time and cost for preparing experimental samples. In the FE model, multiscale method is applied to couple the computations of elastic properties at microscale and mesoscale. Subsequently, a parametric study is performed to analyze the effects of the three design parameters on the elastic properties. Satisfactory results are obtained via the surrogate-based modeling predictions, which are compared with the experimental measurements. Moreover, the predictions obtained from surrogate models concur well with the FE predictions. This study orients a new direction for predicting the mechanical properties based on surrogate models which can effectively reduce the sample preparation cost and computational efforts.     

制造误差对气体静压圆柱轴承静态特性的FEM分析

在考虑轴承的制造误差情况下,采用Galerkin有限元方法（FEM）对小孔节流型气体静压圆柱轴承的静态特性,即承载能力和刚度,进行了理论计算和分析。所考虑的制造误差为典型的几何形状误差和位置误差。计算结果表明,几何形状误差和位置误差对该种轴承的静态特性有着不同程度的影响。     

Micro-mechanical model for predicting the elasto-plastic behavior of composites based on secant formulation method

To predict the Elasto-Plastic Behaviors (EPBs) of aligned inclusions reinforced composites, this paper develops an interpolative Mori-Tanaka/Double-Inclusion (MT-DI) homogenization model with the secant formulation, and gives the numerical implementation algorithms of the developed MT-DI model with the secant formulation. The Finite Element (FE) homogenization method is implemented to provide the “exact” EPBs of the composites and thus validate the MT-DI model with the secant formulation. The MT-DI model with the 2nd-order secant formulation is validated to provide the more accurate predictions, while the MT-DI model with the 1st-order secant formulation always gives the stiffer predictions. The results show that using the macro-stress and macro-strain as the inputs, the MT-DI model with the secant formulation gives the identical predictions. The predictions of the MT-DI model with the secant formulation vary between those of the MT and DI models with the secant formulation. Meanwhile, the MT-DI model with the secant formulation does not predict the accurate EPBs for the phases of the composites.     

Multi-layered plate finite element models with node-dependent kinematics for smart structures with piezoelectric components

This article presents a type of plate Finite Element (FE) models with adaptive mathematical refinement capabilities for modeling laminated smart structures with piezoelectric layers or distributed patches. The p-version shape functions are used in combination with the higher-order Layer-Wise (LW) kinematics adopting hierarchical Legendre polynomials. Node-Dependent Kinematics (NDK) is employed to implement local LW models in the regions with piezoelectric components and simulate the global substrate structure with the Equivalent Single-Layer (ESL) approach. Through the proposed NDK FE models, the electro-mechanical behavior of smart structures can be predicted with high fidelity and numerical efficiency, and various patch configurations can be conveniently modeled through one set of mesh grids. Moreover, the effectiveness and efficiency of the NDK FE approach are assessed through numerical examples and its application is demonstrated.     

Integration of nanoindentation and finite element method for interpretable tensile properties:A cross-scale calculation method of uneven joints

Nanoindentation testing and its Reverse Analysis Method(RAM) show great potential in understanding the tensile properties of metallic alloys with various microstructures. Nevertheless,the tensile properties of heterogeneous materials such as nickel-based superalloy welded joints have not been well interpreted by combining the microstructures and nanoindentation results, due to their diverse and complex microscopic zones, which throws shade on the properties of separated zones in the material. He...     

有大预紧效应的多叶径向箔片轴承的分析

针对有大预紧效应多叶径向箔片轴承的结构特点,通过引入辅助分析部件建立箔片轴承的非线性接触有限元仿真模型。在分析过程中通过改变作为刚性体的辅助部件和作为柔性体轴承部件之间的几何位置来模拟大预紧多叶箔片径向轴承的装配过程并最终求得箔片轴承中接触应力应变。将处理轴承间隙压力分布的控制方程-雷诺方程和非线性的接触求解过程耦合起来迭代求解,第一次对有大预紧效应的多叶箔片径向轴承进行了完全气弹耦合数值分析。      

燕尾形榫联接修形的数值研究

将边界元法用于榫联接修形的数值分析研究。对分析具有摩擦接触问题的边界元法的理论及实现方法作了讨论;将边界元分析应用于航空发动机轮盘榫槽与叶片样头联接的计算,并与激光散斑法的实测结果进行比较;对修形的榫联接进行了计算分析,初步揭示了修形榫联接接触边界上力学参数的变化规律,与试验结果相符合。     

金属基复合材料的蠕变力学研究进展

通过对国内外近10年来在金属基复合材料蠕变力学领域的研究状况的分析，着重从理论、试验和有限元数值模拟三方面评述已经取得的成果和现实需要解决的问题。最后对我国在该领域未来的研究方向和需要开展的工作提出了建议。     

Dynamic parametric modeling-based model updating strategy of aeroengine casings

For accurate Finite Element (FE) modeling for the structural dynamics of aeroengine casings, Parametric Modeling-based Model Updating Strategy (PM-MUS) is proposed based on efficient FE parametric modeling and model updating techniques regarding uncorrelated/correlated mode shapes. Casings structure is parametrically modeled by simplifying initial structural FE model and equivalently simulating mechanical characteristics. Uncorrelated modes between FE model and experiment are reasonably handled by adopting an objective function to recognize correct correlated modes pairs. The parametrized FE model is updated to effectively describe structural dynamic characteristics in respect of testing data. The model updating technology is firstly validated by the detailed FE model updating of one fixed–fixed beam structure in light of correlated/uncorrelated mode shapes and measured mode data. The PM-MUS is applied to the FE parametrized model updating of an aeroengine stator system (casings) which is constructed by the proposed parametric modeling approach. As revealed in this study, (A) the updated models by the proposed updating strategy and dynamic test data is accurate, and (B) the uncorrelated modes like close modes can be effectively handled and precisely identify the FE model mode associated the corresponding experimental mode, and (C) parametric modeling can enhance the dynamic modeling updating of complex structure in the accuracy of mode matching. The efforts of this study provide an efficient dynamic model updating strategy (PM-MUS) for aeroengine casings by parametric modeling and experimental test data regarding uncorrelated modes.     

A dynamic modeling approach for nonlinear vibration analysis of the L-type pipeline system with clamps

There exists a lot of research on the nonlinear vibration of the pipeline system with different boundary conditions. To the best of our knowledge, little research on the actual constraint of the clamp has been performed. In this paper, according to hysteresis loops of the clamp obtained from experimental test, the simplified bilinear stiffness and damping model is proposed. Then the Finite Element (FE) model of L-type pipeline system with clamps is established using Timoshenko beam theory in combination with aforementioned stiffness-damping model. Both hammering and shaker tests verify the FE model via the comparisons of natural frequencies and vibration responses. The results show that the maximum errors of natural frequencies and vibration responses are about 8.31% and 17.6%, respectively. The proposed model can simulate the dynamic characteristics of the L-type pipeline system with clamps well, which is helpful to provide some guidance for the early design stage of pipeline in aero-engine.     

Hierarchical model updating strategy of complex assembled structures with uncorrelated dynamic modes

In structural simulation and design, an accurate computational model directly determines the effectiveness of performance evaluation. To establish a high-fidelity dynamic model of a complex assembled structure, a Hierarchical Model Updating Strategy(HMUS) is developed for Finite Element(FE) model updating with regard to uncorrelated modes. The principle of HMUS is first elaborated by integrating hierarchical modeling concept, model updating technology with proper uncorrelated mode treatment, and...     

Characteristics and mechanisms of turboshaft engine axial compressor casing containment

To investigate the containment characteristics and mechanisms of axial compressor blade and casing in turboshaft engine, experimental and simulation research is conducted on Titanium alloy axial compressor blades and stainless steel simulator casings in this paper. Experiments for four thicknesses (from 0.8 mm to 1.4 mm) of casings are presented on high-speed spin tester. Perforation, ricochet with and without failure of the casings are obtained in test results. Three obvious bulges or dishing region are observed, petaling failure occurs in the first bulge or the third deformation region. Parabolic and elongated dimples are observed at the fracture surface. Finite Element (FE) models with calibrated Johnson-Cook material behavior law are built and analyzed by using explicit dynamic software for a better understanding on the containment behavior. Good agreement is obtained between the experimental observations and numerical predictions. The evolution of the impact force, energy absorption, temperature increase and the cracks’ propagation are analyzed. Three force peaks occur in the impact process. Energy analysis reveals that penetration condition of ricochet with failure leads to most internal energy of the casing.     

基于能量有限元法和虚拟模态综合法的高频冲击响应分析方法

为了分析结构受到高频冲击载荷激励后的瞬态响应,提出了一种基于能量有限元法（EFEM）和虚拟模态综合（VMSS）法的高频冲击响应分析方法。通过能量有限元法进行高频稳态分析,获得结构频响函数（FRF）的频段平均值,然后结合虚拟模态综合法得到虚拟模态振型系数,最后通过Duhamel积分获得结构在高频冲击载荷作用下的瞬态响应。对一简支梁模型进行算例分析,将本文方法的结果与传统有限元法（FEM）和统计能量法（SEA）的分析结果进行对比,验证了所提方法的有效性,也表明该方法具有模型简单、分析速度快等优点。     

Nonlinear vibration response characteristics of a dual-rotor-bearing system with squeeze film damper

Rolling bearing and Squeeze Film Damper (SFD) are used in rotor support structures, and most researches on the nonlinear rotor-bearing system are focused on the simple rotor-bearing systems. This work emphasizes the comparative analysis of the influence of SFD on the nonlinear dynamic behavior of the dual-rotor system supported by rolling bearings. Firstly, a reduced dynamic model is established by combining the Finite Element (FE) method and the free-interface method of component mode synthesis. The proposed model is verified by comparing the natural characteristics obtained from an FE model with those from the experiment. Then, the steady-state vibration responses of the system with or without SFD are solved by the numerical integration method. The influences of the ball bearing clearance, unbalance, centralizing spring stiffness and oil film clearance of SFD on the nonlinear steady-state vibration responses of the dual-rotor system are analyzed. Results show that SFD can effectively suppress the amplitude jump of the dual rotor system sustaining two rotors unbalance excitations. As the ball bearing clearance or unbalance increases, the amplitude jump phenomenon becomes more obvious, the resonance hysteresis phenomenon strengthens or weakens, the resonant peaks shift to the left or the right, respectively. SFD with unreasonable parameters will aggravate the system vibration, the smaller the oil film clearance, the better the damping performance of the SFD, the larger the centralizing spring stiffness is, the larger resonance amplitudes are.     

Nonlinear characteristic investigation of magnetorheological damper-rotor system with local nonlinearity

Magnetorheological (MR) dampers show superior performance in reducing rotor vibration, but their high nonlinearity will cause nonsynchronous response, resulting in fatigue and instability of rotors. Herein, we are devoted to the investigation of the nonlinear characteristics of MR damper mounted on a flexible rotor. First, Reynolds equations with bilinear constitutive equations of MR fluid are employed to derive nonlinear oil film forces. Then, the Finite Element (FE) model of rotor system is developed, where the local nonlinear support forces produced by MR damper and its coupling effects with the rotor are considered. A hybrid numerical method is proposed to solve the nonlinear FE motion equations of the MR damper-rotor system. To validate the proposed model, a rotor test bench with two dual-coil MR dampers is constructed, upon which experimental studies on the dynamic characteristics of MR damper-rotor system are carried out. The effects of different system parameters, including rotational speed, excitation current and amount of unbalance, on nonlinear dynamic behaviors of MR damper-rotor system are evaluated. The results show that the system may appear chaos, jumping, and other complex nonlinear phenomena, and the level of the nonlinearity can be effectively alleviated by applying suitable excitation current and oil supply pressure.     

Three-dimensional temperature prediction in cylindrical turning with large-chamfer insert based on a modified slip-line field approach

Chamfered inserts have found broader applications in metal cutting process especially in high-performance machining of hard-to-cut materials for their excellent edge resistance and cutting toughness. However, excessive heat generation and resulting high cutting temperature eventually cause severe tool wear and poor surface integrity, which simultaneously limits the optimal selection of machining parameters. In the present study, an analytical thermal–mechanical model is proposed for the prediction of the three-dimensional (3-D) temperature field in cylindrical turning with chamfered round insert based on a modified slip-line field approach. First, an innovative discretization method is introduced in a general 3-D coordinate system to provide a comprehensive demonstration of the irregular cutting geometry and heat generation zones. Then, a plasticity-theory-based slip-line field model is developed and employed to determine the intensities and geometries of every elementary heat sources in Primary Deformation Zones (PDZ), Secondary Deformation Zones (SDZ) and Dead Metal Zones (DMZ). At last, a 3-D analytical model is suggested to calculate the temperature increases caused by the entire heat sources and associated images. The maximum cutting temperature region predicted is found existing upon the chip-tool contact area rather than the tool edge. Moreover, the rationalities of cutting parameters employed are analyzed along with theoretical material removal rates and ensuing maximum cutting temperatures. The results indicate that the cutting conditions with large depth of cut and high cutting speed are more desirable than those with high feed rates. The proposed models are respectively verified through a series of 3-D Finite Element (FE) simulations and dry cutting experiments of Inconel 718 with chamfered round insert. Satisfactory agreement has been reached between the predictions and simulations as well as the measurements, which confirms the correctness and effectiveness of the presented analytical model.     

An integrated preforming-performance model for high-fidelity performance analysis of cured woven composite part with non-orthogonal yarn angles

Preforming process would change yarn angle and yarn orientation, its influence on the material properties and material orientations needs to be considered in the performance analysis.However, most current performance models fail to account for the preforming effect. An integrated performance model accounting for the impact of preforming has been developed. In this integrated model, part geometry, yarn angle and orientation after preforming of multiple prepreg layers are predicted by Finite Eleme...