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.     

飞机整体壁板广义槽分层识别方法及其实现算法(英文)

To automatically obtain a machining area in numerical control （NC） programming, a data model of generalized pocket is established by analyzing aircraft integral panel characteristics, and a feature recognition approach is proposed. First, by reference to the practical slice-machining process of an aircraft integral panel, both the part and the blank are sliced in the Z-axis direction; hence a feature profile is created according to the slicing planes and the contours are formed by the intersection of the slicing planes with the part and its blank. Second, the auxiliary features of the generalized pocket are also determined based on the face type and the position, to correct the profile of the pocket. Finally, the generalized pocket feature relationship tree is constructed by matching the vertical relationships among the features. Machining feature information produced by using this method can be directly used to calculate the cutter path. The validity and practicability of the method is verified by NC programming for aircraft panels.     

Residual Strength of Stiffened LY12CZ Aluminum Alloy Panels with Widespread Fatigue Damage

Experimental and analytical investigations on the residual strength of the stiffened LY12CZ aluminum alloy panels with widespread fatigue damage （WFD） are conducted. Nine stiffened LY12CZ aluminum alloy panels with three different types of damage are tested for residual strength. Each specimen is pre-cracked at rivet holes by saw cuts and subjected to a monotonically increasing tensile load until failure is occurred and the failure load is recorded. The stress intensity factors at the tips of the lead crack and the adjacent WFD cracks of the stiffened aluminum alloy panels are calculated by compounding approach and finite element method （FEM） respectively. The residual strength of the stiffened panels with WFD is evaluated by the engineering method with plastic zone linkup criterion and the FEM with apparent fracture toughness criterion respectively. The predicted residual strength agrees well with the experiment results. It indicates that in engineering practice these methods can be used for residual strength evaluation with the acceptable accuracy. It can be seen from this research that WFD can significantly reduce the residual strength and the critical crack length of the stiffened panels with WFD. The effect of WFD crack length on residual strength is also studied.     

Effect of Coriolis and centrifugal forces on flow and heat transfer at high rotation number and high density ratio in non orthogonally internal cooling channel

Numerical predictions of three-dimensional flow and heat transfer are performed for a two-pass square channel with 45° staggered ribs in non-orthogonally mode-rotation using the second moment closure model. At Reynolds number of 25,000, the rotation numbers studied were 0,0.24, 0.35 and 1.00. The density ratios were 0.13, 0.23 and 0.50. The results show that at high buoyancy parameter and high rotation number with a low density ratio, the flow in the first passage is governed by the secondary flow induced by the rotation whereas the secondary flow induced by the skewed ribs was almost distorted. As a result the heat transfer rate is enhanced on both co-trailing and co-leading sides compared to low and medium rotation number. In contrast, for the second passage, the rotation slightly reduces the heat transfer rate on co-leading side at high rotation number with a low density ratio and degrades it significantly on both co-trailing and co-leading sides at high buoyancy parameter compared to the stationary, low and medium rotation numbers. The numerical results are in fair agreement with available experimental data in the bend region and the second passage, while in the first passage were overestimated at low and medium rotation numbers.     

Material parameter modeling and solution technique using birth–death element for notched metallic panel repaired with bonded composite patch

This paper seeks to outline a novel three-layer model and a new birth–death element solution technique to evaluate static strength of notched metallic panel repaired with bonded composite patch and to optimize material parameters. The higher order 3D, 8-node isotropic solid element and 8-node anisotropic layered solid element with three degrees of freedom per node are respectively implemented to model substrate panel, adhesive layer and composite patch to establish three-layer model of repaired panel. The new solving technique based on birth–death element is developed to allow solution of the stress pattern of repaired panel for identifying failure mode. The new model and its solution are used to model failure mode and residual strength of repaired panel, and the obtained results have a good agreement with the experimental findings. Finally, the influences of material parameter of adhesive layer and composite patch on the residual strength of repaired panel are investigated for optimizing material properties to meet operational and environmental constraints.     

孔间距和前倾角对进出口扩张形气膜孔流量系数的影响(英文)

Discharge coefficients of 3-in-1 hole of three inclination angles and three spacing between holes are presented which described the discharge behavior of a row of holes. The inlet and outlet of the 3-in-1 hole both have a 15° lateral expansion. The flow conditions considered are mainstream turbulence intensities and density ratios of secondary flow to mainstream. The momentum flux ratios varied in the range from 1 to 4. The comparison is made of the discharge coefficients of three shaped holes to find an optimal hole with low flow loss. The results show that the discharge coefficients of 3-in- 1 hole are highest in three shaped holes and therefore this article is focused on the measurements of discharge coefficients of 3-in-1 hole for various geometries and aerodynamic parameters. The measured results of 3-in-1 hole indicate that turbulence intensities, density ratios and momentum flux ratios have weak influence on discharge coefficients for inclination angle of 20°. The high turbulence intensity yields the small discharge coefficients for inclination angle of 45° and 90°. The increased both momentum flux ratios and density ratios lead to the increased discharge coefficients for inclination angle of 45° and 90°. The increased inclination angle causes the rapidly increased discharge coefficients. There is a weak dependence of discharge coefficients on hole pitches.     

Wing aeroelasticity analysis based on an integral boundary-layer method coupled with Euler solver

An interactive boundary-layer method, which solves the unsteady flow, is developed for aeroelastic computation in the time domain. The coupled method combines the Euler solver with the integral boundary-layer solver(Euler/BL) in a ‘‘semi-inverse" manner to compute flows with the inviscid and viscous interaction. Unsteady boundary conditions on moving surfaces are taken into account by utilizing the approximate small-perturbation method without moving the computational grids. The steady and unsteady flow calculations for the LANN wing are presented. The wing tip displacement of high Reynolds number aero-structural dynamics(HIRENASD) Project is simulated under different angles of attack. The flutter-boundary predictions for the AGARD445.6 wing are provided. The results of the interactive boundary-layer method are compared with those of the Euler method and experimental data. The study shows that viscous effects are significant for these cases and the further data analysis confirms the validity and practicability of the coupled method.     

A vision-based navigation approach with multiple radial shape marks for indoor aircraft locating

Since GPS signals are unavailable for indoor navigation, current research mainly focuses on vision-based locating with a single mark. An obvious disadvantage with this approach is that locating will fail when the mark cannot be seen. The use of multiple marks can solve this problem. However, the extra process to design and identify different marks will significantly increase system complexity. In this paper, a novel vision-based locating method is proposed by using marks with feature points arranged in a radial shape. The feature points of the marks consist of inner points and outer points. The positions of the inner points are the same in all marks, while the positions of the outer points are different in different marks. Unlike traditional camera locating methods （the PnP methods）, the proposed method can calculate the camera location and the positions of the outer points simultaneously. Then the calculation results of the positions of the outer points are used to identify the mark. This method can make navigation with multiple marks more efficient. Simulations and real world experiments are carried out, and their results show that the proposed method is fast, accurate and robust to noise.     

Flow and Heat Transfer Characteristics in Rotating Two-pass Channels Cooled by Superheated Steam

In a modern gas turbine,using superheated steam to cool the vane and blade for internal convection cooling is a promising alternative to traditional compressor air.However,further investigations of steam cooling need to be performed.In this paper,the three-dimensional flow and heat transfer characteristics of steam are numerically investigated in two-pass square channels with 45° ribbed walls under stationary and rotating conditions.The investigated rotation numbers are 0 and 0.24.The simulation is carried out by solving the Reynolds averaged Navier-Stokes equations employing the Reynolds stress turbulence model,especially considering two additional terms for Coriolis and rotational buoyancy forces caused by the rotating effect.For comparison,calculations for the air-cooled channels are done first at a Reynolds number of 25 000 and inlet coolant-to-wall density ratio of 0.13.The results are compared with the experiment data.Then the flow and heat transfer in steam-cooled channels are analyzed under the same operating conditions.The results indicate that the superheated steam has better heat transfer performance than air.Due to the combined effect of rotation,skewed ribs and 180° sharp turn,the secondary flow pattern in steam-cooled rotating two-pass channels is quite complex.This complex secondary flow pattern leads to strong anisotropic turbulence and high level of anisotropy of Reynolds stresses,which have a significant impact on the local heat transfer coefficient distributions.     

Cathode design investigation based on iterative correction of predicted profile errors in electrochemical machining of compressor blades

Electrochemical machining (ECM) is an effective and economical manufacturing method for machining hard-to-cut metal materials that are often used in the aerospace field. Cathode design is very complicated in ECM and is a core problem influencing machining accuracy, especially for complex profiles such as compressor blades in aero engines. A new cathode design method based on iterative correction of predicted profile errors in blade ECM is proposed in this paper. A math-ematical model is first built according to the ECM shaping law, and a simulation is then carried out using ANSYS software. A dynamic forming process is obtained and machining gap distributions at different stages are analyzed. Additionally, the simulation deviation between the prediction profile and model is improved by the new method through correcting the initial cathode profile. Further-more, validation experiments are conducted using cathodes designed before and after the simulation correction. Machining accuracy for the optimal cathode is improved markedly compared with that for the initial cathode. The experimental results illustrate the suitability of the new method and that it can also be applied to other complex engine components such as diffusers.     

带筋整体壁板预应力喷丸成形数值模拟及变形预测

带筋整体壁板特别是高筋整体壁板具有优异的结构效率、减重效益和密封效果，在航空航天等领域备受青睐。预应力喷丸成形是大中型、长寿命、高性能复杂带筋整体壁板的一种不可多得的有效成形手段。数值模拟是促进带筋整体壁板预应力喷丸成形技术研究、进步与应用的一种极具潜力的途径。针对带筋整体壁板预应力喷丸成形数值模拟，建立了基于响应面函数的多弹丸撞击有限元模型、基于应变中性层内移的反弯曲应力场法模拟模型、预应力喷丸成形RBF神经网络预测模型，实现了喷丸成形应力场法数值模拟、预应力喷丸成形较高精度数值模拟与变形预测，为带筋整体壁板预应力喷丸成形技术研究和实际应用提供了一种更为便捷、高效、经济的途径。     

整体壁板激光辅助预应力成形

整体壁板成形是大飞机制造的关键难题之一。本文探索了激光辅助预应力成形方法应用于整体壁板这一类筋条加强结构工件的可行性。该方法的原理是：首先对整体壁板施加弹性载荷,然后用激光束扫描并加热其弹性变形能集中区域,通过降低该区域材料的屈服极限促使弹性变形能转化为塑性功,实现整体壁板的塑性成形。对网格筋条加强结构试件进行了激光辅助预应力成形试验和有限元分析以及时效成形试验,验证了激光辅助预应力成形方法的可行性,并讨论了激光扫描道次与工件成形效率的关系。     

A progressive approach to predict shot peening process parameters for forming integral panel of Al7050-T7451

In this paper, a progressive approach to predict the multiple shot peening process parameters for complex integral panel is proposed. Firstly, the invariable parameters in the forming process including shot size, mass flow, peening distance and peening angle are determined according to the empirical and machine type. Then, the optimal value of air pressure for the whole shot peening is selected by the experimental data. Finally, the feeding speed for every shot peening path is predicted by regression equation. The integral panel part with thickness from 2 mm to 5 mm and curvature radius from 3200 mm to 16000 mm is taken as a research object, and four experiments are conducted. In order to design specimens for acquiring the forming data, one experiment is conducted to compare the curvature radius of the plate and stringer-structural specimens, which were peened along the middle of the two stringers. The most striking finding of this experiment is that the outer shape error range is below 3.9%, so the plate specimens can be used in predicting feeding speed of the integral panel. The second experiment is performed and results show that when the coverage reaches the limit of 80%, the minimum feeding speed is 50 mm/s. By this feeding speed, the forming curvature radius of the specimens with different thickness from the third experiment is measured and compared with the research object, and the optimal air pressure is 0.15 MPa. Then, the plate specimens with thickness from 2 mm to 5 mm are peened in the fourth experiment, and the measured curvature radius data are used to calculate the feeding speed of different shot peening path by regressive analysis method. The algorithm is validated by forming a test part and the average deviation is 0.496 mm. It is shown that the approach can realize the forming of the integral panel precisely.     

Design and Optimization of Press Bend Forming Path for Producing Aircraft Integral Panels with Compound Curvatures

In order to find out the optimal press bend forming path in fabricating aircraft integral panels, this article proposes a new method on the basis of the authors' previous work. It is composed of the finite element method (FEM) equivalent model, the surface curvature analysis, the artificial neural network response surface and the genetic algorithm. The method begins with analyzing the objective's shape curvature to determine the bending position. Then it optimizes the punch travel at each bending position by the following steps: (1) Establish a multi-step press bend forming FEM equivalent model, with which the FEM experiments designed with the Taguchi method are performed. (2) Construct a back-propagation (BP) neural network response surface with the data from the FEM experiments. (3) Use the genetic algorithm to optimize the neural network response surface as the objective function. Finally, this method is verified by press bending a complicated double-curvature grid-type stiffened panel and bears out its effectiveness and intrinsic worth in designing the press bend forming path.     

A toolpath strategy for improving geometric accuracy in double-sided incremental sheet forming

The double-sided incremental forming (DSIF) improved the process flexibility compared to other incremental sheet forming (ISF) processes. Despite the flexible nature, it faces the challenge of low geometric precision like ISF variants. In this work, two strategies are used to overcome this. First, a novel method is employed to determine the optimal support tool location for improving geometric precision. In this method, the toolpath oriented the tools to each other systematically in the circumferential direction. Besides, it squeezed the sheet by the same amount at the point of interest. The impacts of various support tool positions in the circumferential direction are evaluated for geometric precision. The results demonstrate that the support tool should support the master tool within 10° to its local normal in the circumferential direction to improve the geometric accuracy. Second, a two-stage process reduced the geometric error of the part by incrementally accommodating the springback error by artificially increasing the step size for the second stage. With the optimal support tool position and two-stage DSIF, the geometric precision of the part has improved significantly. The proposed method is compared to the best DSIF toolpath strategies for geometric accuracy, surface roughness, forming time, and sheet thickness fluctuations using grey relational analysis (GRA). It outperforms the other toolpath strategies including single-stage DSIF, accumulative double-sided incremental forming (ADSIF), and two-stage mixed double sided incremental forming (MDSIF). Our approach can improve geometric precision in complex parts by successfully employing the support tool and managing the springback incrementally.     

整体加筋壁板的破损安全特性与断裂控制分析

用有限元和断裂力学方法分析大型整体机翼下壁板的破损安全特性。对一个九桁条铝合金整体加筋板,裂纹从断裂筋条下向两侧均匀扩展的开裂形式,进行了应力强度因子和剩余强度计算,并和另一个具有相同构形的铆接加筋板的结果进行对比。当裂纹长度在两倍桁条间距以内时两种加筋板的剩余强度水平相当。但是当整体壁板中的裂纹穿过外侧相邻桁条时,整体桁条的止裂能力逐渐减弱。研究对整体壁板蒙皮胶接止裂条的断裂控制措施及有限元建模分析方法,描述止裂条胶层局部脱粘的处理方法和迭代过程。止裂条材料分别选用了钛合金和单向层压复合材料。计算结果说明胶接止裂措施能显著提高整体加筋板的破损安全性能。     

大型机翼整体壁板系统化喷丸成形技术

大型机翼整体壁板是现代大型飞机重要的大型承力整体结构件并且通常直接构成飞机的气动外形。喷丸成形是现代大型轻质高强铝合金整体壁板件成形制造的首选技术方法，但如何实现大型机翼整体壁板的精确喷丸成形一直是现代航空制造技术领域的一个难点问题。针对这一工程问题，本文采用系统化的方法，将影响大型机翼整体壁板喷丸成形精度的因素分解为壁板平面板坯误差、成形参数设计准确度、成形参数控制精度、环境因素。针对这些因素，采用基于变形位能最小的板坯优化设计来减小由板坯导致的成形误差；采用数据拟合、人工神经网络以及解析模型计算相结合的喷丸成形参数综合设计方法来提高喷丸参数设计的精度和效率；建立了板坯修正模型以修正环境温度、喷丸设备参数波动等因素对成形件形状和尺寸的影响；对于从喷丸设备上下线后仍存在的外形贴模误差，则采用手提喷丸机进行局部的渐进式校形喷丸至外形贴模。壁板喷丸成形的工程实践表明，本文所提出的系统化方法能够有效提高大型机翼整体壁板喷丸成形的精度和效率，并可满足工业生产的需求。     

含离散源损伤复合材料加筋板的压缩特性

 通过对含有离散源损伤的复合材料加筋板的压缩实验和有限元模拟,研究了离散源损伤的损伤扩展与破坏特性以及对复合材料加筋板剩余强度的影响。 结果表明，复合材料加筋板的离散源损伤用穿透蒙皮切断桁条的切口来模拟是合适的，蒙皮上的穿透切口前端有很高的应变集中，桁条被切断导致加筋板传力路线改变；离散源损伤使复合材料加筋板的压缩强度下降显著，达到60％左右； 基于Hashin失效准则的渐进损伤有限元数值模拟方法，可以有效地模拟含切口加筋板的宏观损伤扩展和破坏过程，计算结果与实验值吻合较好。     

Improving the process forming limit considering forming defects in the transitional region in local loading forming of Ti-alloy rib-web components

The isothermal local loading forming technology provides a feasible way to form Ti-alloy large-scale rib-web components in aerospace and aviation fields.However,the local load ing process forming limit is restricted by forming defects in the transitional region.In this work,the feasibility of controlling forming defects and improving the process forming limit by adjusting die parameters is explored through finite element (FE) simulation.It is found that the common cavum and folding defects in the transitional region are significantly influenced by the fillet radii of left rib and middle rib,respectively.The cavum and folding defects can be effectively controlled by increas ing the fillet radii of left rib and middle rib,respectively.The process forming limits considering forming defects in the transitional region are determined by the stepwise searching method under various die parameters.Moreover,the relationship between the process forming limit and die parameters is developed through the response surface methodology (RSM).The developed RSM models suggest that increasing the fillet radii of left and middle ribs is effective to improve the pro cess forming limit during local loading forming of rib-web components.The results will provide technical basis for the design of die parameters and the reduction amount,which is of great impor tance to control forming defects and improve the process forming limit in local loading forming of Ti-alloy large-scale rib-web components.     

Analysis of the effect of stiffener profile on buckling strength in composite isogrid stiffened shell under axial loading

In this paper, the buckling behavior of thin-walled GFRP cylindrical shells with triangular lattice patterned reinforcements formed by helical and circumferential ribs under axial force is analyzed. In this analysis, various models of composite isogrid stiffened cylindrical shells with outer diameter of 150 millimeters, shell thickness of 0.5 millimeters and height of 280 millimeters, stiffened with 6 helical and 2 circumferential ribs and all with the same material properties of shell and ribs are used. Ribs have constant section areas but different shapes and cross section profiles. The effects of these differences on buckling strengths of structures under axial load are studied. For analysis and modeling of structures, Finite Element Analysis method and ANSYS software were used. The results (elastic buckling load) for each model were derived and based on these results, ratio of buckling strengths to weight parameters were calculated for each model and were compared to results obtained from other models. The effect of profile of the ribs on the buckling of shells under axial loading can be concluded from the results. Results showed that stiffening the shells increased the buckling load from 10% up to 36% while decreased the buckling load to weight ratio to 42% up to 52% of an unstiffened shell.