SIMULATION OF COMPOSITE NON-LINEAR MECHANICAL BEHAVIOR OF CMCS BY FEM-BASED MULTI- SCALE APPROACH

The non-linear behavior of continuous fiber reinforced C/SiC ceramic matrix composites （CMCs） under tensile loading is modeled by three-dimensional representative volume element （RVE） models of the composite. The theoretical background of the multi-scale approach solved by the finite element method （FEM） is recalled first- ly. Then the geometric characters of three kinds of damage mechanisms, i.e. micro matrix cracks, fiber/matrix interface debonding and fiber fracture, are studied. Three kinds of RVE are proposed to model the microstructure of C/SiC with above damage mechanisms respectively. The matrix cracking is modeled by critical matrix strain en- ergy （CMSE） principle while a maximum shear stress criterion is used for modeling fiber/matrix interface debond- ing. The behavior of fiber fracture is modeled by the famous Weibull statistic theory. A numerical example of con- tinuous fiber reinforced C/SiC composite under tensile loading is performed. The results show that the stress/ strain curve predicted by the developed model agrees with experimental data.     

SIMULATION OF LARGE CIVIL AIRCRAFT LANDING RESPONSE UNDER CONSIDERATIONS OF FLEXIBLE AIRFRAME

To evaluate the landing response of the large civil aircraft in the conceptual design phase , a method for simulating aircraft landing is given.The model for the shock absorber is investigated.The flexible airframe model is established using finite element model ( FEM ) to analyze its modes.Then , the whole aircraft model with flexible airframe is made for the multibody simulation.Tail-down , two-point , three-point and sideslip landing scenarios are studied.The influence on the landing performance considering mode superposition of the flexible airframe is analyzed.Both longitudinal and spanwise positions of the main landing gear are changed to research the influence on the landing performance.Results show that the method is feasible.The shock absorber axial force of the main landing gear with the flexible airframe is smaller than that of rigid airframe.The number of mode superposition and the position of main landing gear can influence the landing response.     

Lattice Boltzmann Flux Solver： An Efficient Approach for Numerical Simulation of Fluid Flows

A lattice Boltzmann flux solver （LBFS） is presented for simulation of fluid flows. Like the conventional computational fluid dynamics （CFD） solvers, the new solver also applies the finite volume method to discretize the governing differential equations, but the numerical flux at the cell interface is not evaluated by the smooth function approximation or Riemann solvers. Instead, it is evaluated from local solution of lattice Boltzmann equation （LBE） at cell interface. Two versions of LBFS are presented in this paper. One is to locally apply one-dimensional compressible lattice Boltzmann （LB） model along the normal direction to the cell interface for simulation of compressible inviscid flows with shock waves. The other is to locally apply multi-dimensional LB model at cell interface for simulation of incompressible viscous and inviscid flows. The present solver removes the drawbacks of conventional lattice Boltzmann method （LBM） such as limitation to uniform mesh, tie-up of mesh spacing and time interval, limitation to viscous flows. Numerical examples show that the present solver can be well applied to simulate fluid flows with non-uniform mesh and curved boundary.     

Identification of Time-Varying Modal Parameters for Thermo- Elastic Structure Subject to Unsteady Heating

A time-varying modal parameter identification method combined with Bayesian information criterion （BIC） and grey correlation analysis （GCA） is presented for a kind of thermo-elastic structures with sparse natural frequencies and subject to an unsteady temperature field. To demonstrate the method, the thermo-elastic structure to be identified is taken as a simply-supported beam with an axially movable boundary and subject to both random excitation and an unsteady temperature field, and the dynamic outputs of the beam are first simulated as the meas- ured data for the identification. Then, an improved time-varying autoregressive （TVAR） model is generated from the simulated input and output of the system. The time-varying coefficients of the TVAR model are expahded as a finite set of time basis functions that facilitate the time-varying coefficients to be time invariant. According to the BIC for preliminarily determining the scope of the order number, the grey system theory is introduced to determine the order of TVAR and the dimension of the basis functions simultaneously via the absolute grey correlation degree （AGCD）. Finally, the time-varying instantaneous frequencies of the system are estimated by using the recursive least squares method. The identified results are capable of tracking the slow time-varying natural frequencies with high accuracy no matter for noise-free or noisy estimation.     

Dynamical Model Updating Based on Modal Tests with Changed Structure

A new approach to modifying the stiffness and mass matrices of finite element models is presented to improve the calculation precision. By measuring the mode frequencies and shapes of both of the original and the new structures with changed stiffness and mass, the stiffness and mass matrices of the finite element model can be updated through matrices calculation and solving algebra equations. Taking a multi-freedom model as an example, the relation between the number of the modes and the correction precision of stiffness and mass matrix elements is researched. The facility and precision of the method are totally confirmed especially when the modeling error is known limited to a definite local range. The feasibility of the approach is proven by an effective engineering application to the model updating of a wing piece used in flutter test.     

Modeling and Simulation for Refueling Boom and Receiver in Coupled Mode

In coupled mode,the major problem of boom refueling system is undesirable nozzle loads.An automated load alleviation system(ALAS)is needed to alleviate nozzle loads.In order to simulate dynamic of the system and to validate ALAS,dynamic model is developed.Two models are established,which are the static model and the moving model,named after the two relative states between the fixed boom and the extension boom.Kane method is employed as main method considering system′s multi-body characteristics.D′Alembert′s principle is used to calculate nozzle loads.Simulation is conducted to research the effects of position disturbance and velocity disturbance on nozzle loads.Results indicate that position disturbance plays a more significant role in inducing nozzle loads.A fuzzy control law based ALAS is validated using the formulated model.It is concluded that this model can simulate system dynamic and validate ALAS.     

SPH-BASED NUMERICAL ANALYSIS FOR GRANULAR MATERIAL MODEL OF SAND AND REGOLITH

The development of an innovative drilling system makes it necessary to model the material of extra-ter- restrial soil, i.e. regolith or lunar soil. And then the drilling process is numerically simulated. Since real regolith is very scarce and costly, sand properties are investigated by means of experiments and sand model is constructed with the capability of partly reflecting these properties. With the use of smooth particle hydrodynamics （SPH） method in LS-DYNA, data analysis and modeling process are presented to reach the following achievements：（1） Develop a general model approach of sand using SPH method; （2） Compare SPH results with experimental data to validate it; （3） Adapt the sand model to lunar soil while fitting the simulation of dual-reciprocating drilling （DRD） process.     

EXPERIMENTAL MODELING OF MIDDLE-RISE BUILDING VIA AMBIENT MODAL IDENTIFICATION

Experimental modeling of a middle-rise office building via ambient modal identification is presented. A 200-DOF (Dimension of freedom) test model is designed to correlate with finite element mode. A newly developed frequency-spatial domain decomposition ( FSDD ) technique is used to identify modal characteristics of the full-size building by using ambient response measurements. In the interested frequency ranges of 0～4.5 Hz and 0～ 6.5 Hz altogether 9 bending and torsion modes are identified. As one of the major focuses of the project, the accurate damping estimation is conducted based on FSDD. The identified modal frequencies and mode shapes are utilized for finite element model tuning. Excellent agreement has been achieved with respect to the final tuned finite element (FE) model up to 9 modes.     

WAKE GEOMETRY CALCULATIONS FOR TILT-ROTOR USING VISCOUS VORTEX METHOD

A tilt-rotor unsteady flow analytical method has been developed based upon viscous vortex-particle meth- od. In this method, the vorticity field is divided into small assembled vortex particles. Vortex motion and diffusion are obtained by solving the velocity-vorticity-formed incompressible Navier-Stokes equations using a grid-free La- grangian simulation method. Generation of the newly vortex particles is calculated by using the Weissinger-L lifting surface model. Furthermore, in order to significantly improve computational efficiency, a fast multiple method （FMM） is introduced into the calculation of induced velocity and its gradient. Finally, the joint vertical experimen- tal （JVX） tilt-rotor is taken as numerical examples to analyze. The wake geometry and downwash are investigated for both hover and airplane modes. The proposed method for tilt-rotor flow analysis is verified by comparing its re- sults with those available measured data. Comparison indicates that the current method can accurately capture the complicated tilt-rotor wake variation and be suitable for aerodynamic interaction simulation in complex environ- ments. Additionally, the aerodynamic interactional characteristics of dual-rotor wake are discussed in different ro- tor distance. Results show that there are significant differences on interactional characteristics between hover mode and airplane mode.     

STRAIGHTFORWARD MULTI-SCALE BOUNDARY ELEMENT METHOD FOR GLOBAL/LOCAL MECHANICAL ANALYSIS OF ELASTIC HETEROGENEOUS MATERIAL

A straightforward multi-scale boundary element method is proposed for global and local mechanical analysis of heterogeneous material.The method is more accurate and convenient than finite element based multi-scale method.The formulations of this method are derived by combining the homogenization approach and the fundamental equations of boundary element method.The solution gives the convenient formulations to compute global elastic constants and the local stress field.Finally,two numerical examples of porous material are presented to prove the accuracy and the efficiency of the proposed method.The results show that the method does not require the iteration to obtain the solution of the displacement in micro level.     

HIGH CYCLE FATIGUE RELIABILITY ANALYSIS ON ROTOR HUB BASED ON APPROXIMATION TECHNIQUE

A high cycle fatigue reliability analysis approach to helicopter rotor hub is proposed under working load spectrum. Automatic calculation for the approach is implemented through writing the calculating programs. In the system, the modification of geometric model of rotor hub is controlled by several parameters, and finite element method and S-N curve method are then employed to solve the fatigue life by automatically assigned parameters. A database between assigned parameters and fatigue life is obtained via Latin Hypercube Sampling （LHS） on toler- ance zone of rotor hub. Different data-fitting technologies are used and compared to determine a highest-precision approximation for this database. The parameters are assumed to he independent of each other and follow normal distributions. Fatigue reliability is then computed by the Monte Carlo （MC） method and the mean-value first order second moment （MFOSM） method. Results show that the approach has high efficiency and precision, and is suit- able for engineering application.     

Aircraft Electric Anti-skid Braking System Based on Fuzzy-PID Controller with Parameter Self-adjustment Feature

The principle of electric braking system is analyzed and an anti-skid braking system based on the slip rate control is proposed. The fuzzy-PID controller with parameter self-adjustment feature is designed for the anti-skid braking system. The dynamic model of aircraft ground braking is established in the simulation environment of MATLAB/SIMULINK, and simulation results of dry runway and wet runway are presented. The results show that the fuzzy-PID controller with parameter self-adjustment feature for the electric anti-skid braking system keeps working in the state of stability and the brake efficiencies are increased to 93% on dry runway and 82% on wet runway respectively.     

Construction of Crack Perturbation Model and Forward Semi-analytical Model of Attached Eddy Current Sensor

Fatigue damage monitoring is critical metallic structure health monitoring of aircraft.The sensor should be high sensitive,easy to be integrated into structure and well adaptable for poor working conditions.Therefore,an attached eddy current sensor with flexible plane is put forward and its characteristics are analyzed.By extracting material′s conductivity as the crack features,forward semi-analytical model is established and parameter optimizations are carried out.Crack perturbation model of attached eddy current sensor is constructed,and perturbation voltages of sensing channels under three-dimension structural crack are obtained.To verify the sensor′s performance,monitoring experiment on crack extension is conducted under condition of 3 MHz frequency.The validation experimental results show that perturbation model of 2A12-T4 aluminum alloy agrees well with experiment results,and perturbation model errors of four sensing channels are within 25%.The attached eddy current sensor is capable of testing the crack nondestructively and measuring the crack extension quantitatively with the accuracy of 1mm.     

Three-Dimentional Surface Light Irradiance Reconstruction in Bioluminescence Tomography

Reconstruction of 3D surface irradiance distribution using multiple views captured by charged coupled device(CCD)camera is the basis of solving the light source in bioluminescence tomography(BLT).A simple and convenient mapping technique based on the pin-hole imaging model and Lambert′s cosine law was presented to establish the relationship between gray levels and irradiance intensities.Compared with previous integrating sphere camera calibration used in BLT,the proposed method can effectively avoid heavy burden of simulation experiment to obtain the corresponding relationship of gray levels and irradiance intensities.The accuracy and feasibility of the proposed method are validated with no more than 1mm location error by different types of phantom experiments.The mapping approach is also applicable to other noncontact optical imaging system.     

Full Wavefield Analysis for Damage Assessment in Composite Materials

In this paper damage assessment based on guided elastic wave propagation phenomenon is presented.Guided waves are generated by piezoelectric transducer and registered by scanning laser doppler vibrometer(SLDV).Signal processing is based on the analysis of full wavefield measurements gathered from dense mesh of measurement points spanned over area of investigated samples.Full wavefield measurement approach allows creation of animations presenting the guided wave propagation in the structure.Moreover such approach is suitable for analysis of interaction of guided waves with discontinuities located in structure.In the research attention is paid especially on analysis of phenomenon of S0/A0′guided wave mode conversion due to interaction with investigated discontinuities-teflon inserts and impact damage.The presented work is related to glass fibre reinforced polymer(GFRP)samples.In the research,auxiliary non-destructive testing(NDT)method is also utilized.The aim of this method is to indicate the depth of discontinuity,and to prove that delamination was created in the case of impact damage.Auxiliary method is based on terahertz spectroscopy(THz)where the analysis of propagation of electromagnetic waves in the terahertz band is conducted.THz spectroscopy method can be utilized for damage assessment in the dielectric materials like GFRP.     

Fast Numerical Solution to Forward Kinematics of General Stewart Mechanism Using Quaternion

The forward kinematics of the general Stewart mechanism is studied and a fast numerical method is presented.Quaternion is utilized to model the forward kinematics and the equations are merely a system of quadratic ones.The numerical method is a nice simplification of the Newton-Raphson method when applied to this system.A simulation of the movement control of the Stewart mechanism is accomplished,confirming the effectiveness of the proposed algorithm in real-time conditions.     

Analysis for Transmission of Composite Structure with Graphene Using Equivalent Circuit Model

A simple analytical method is presented to analyze the transmission of electromagnetic plane waves through multilayer stacked composite two-dimensional （2D） structures at microwave frequencies. Unlike the traditional structure, high impedance surface with graphene sheet is proposed. The structure includes graphene and thin metal patches and meshes. Simple analytical formulas are introduced for the surface impedance of graphene and for the grid impedance of electrically dense arrays of metal square patches or strips. The result of transmission proper- ties is based on the dynamic tunable model of the high impedance surface, which considers the surface conductivity of graphene layer. The transmission coefficient obtained by using the equivalent circuit method is validated against full-wave numerical simulations. The considered equivalent circuit method can be useful in the design of graphene tunable planar devices.     

An Empirical Method for Prediction of Hypersonic Rarefied Flow-Field Structure

Numerical simulations are presented about the effects of gas rarefaction on hypersonic flow field.Due to the extremely difficult experiment,limited wind-tunnel conditions and high cost,most problems in rarefied flow regime are investigated through numerical methods,in which the direct simulation Monte-Carlo(DSMC)method is widely adopted.And the unstructured DSMC method is employed here.Flows around a vertical plate at a given velocity 7 500 m/s are simulated.For gas rarefaction is judged by the free-stream Knudsen number(Kn),two vital factors are considered:molecular number density and the plate′s length.Cases in which Kn varies from 0.035 to13.36 are simulated.Flow characters in the whole rarefied regime are described,and flow-field structure affected by Knis analyzed.Then,the dimensionless position D*of a certain velocity in the stagnation line is chosen as the marker of flow field to measure its variation.Through flow-field tracing and least-square numerical method analyzing,it is proved that hypersonic rarefied flow field expands outward linearly with the increase of 1/2Kn.An empirical method is proposed,which can be used for the prediction of the hypersonic flow-field structure at a given inflow velocity,especially the shock wave position.     

Rotor Airload and Acoustics Prediction Based on CFD/CSD Coupling Method

An advanced airload and noise prediction method based on computational fluid dynamics/computational structural dynamics(CFD/CSD)coupling for helicopter rotor has been developed in this paper.In the present method,Navier-Stokes equation is applied as the governing equation,and a moving overset grid system is generated in order to account for the blade motions in rotation,flapping and pitching.The blade structural analysis is based on 14-DOF Euler beam model,and the finite element discretization is conducted on Hamilton′s variational principle and moderate deflection theory.Aerodynamic noise is calculated by Farassat 1 Aformula derived from FW-H equation.Using the developed method,numerical example of UH-60 Ais performed for aeroelastic loads calculation in a low-speed forward flight,and the calculated results are compared with both those from isolated CFD method and available experimental data.Then,rotor noise is emphatically calculated by CFD/CSD coupling method and compared with the isolated CFD method.The results show that the aerodynamic loads calculated from CFD/CSD method are more satisfactory than those from isolated CFD method,and the exclusion of blade structural deformation in rotor noise calculation may cause inaccurate results in low-speed forward flight state.     

Attitude Analysis in Process Conflict for C919 Aircraft Manufacturing

Based on the option prioritization in graph model for conflict resolution of two decision makers(DMs),new logical and matrix representations of four stability concepts for DMs′attitude are proposed.The logical representation of attitude is defined,and converted to the matrix form in order to develop a decision support system(DSS)efficiently.Compared with existing definitions of DMs′attitude based on states,the proposed definitions of attitude based on options are convenient and more effective to generate preferences since that of states can be significantly larger than that of options in a large conflict.In addition,it is easier to obtain the information of the prioritization of option statements than to obtain preference of states for users.The proposed representations are applied to the process conflict during aircraft manufacturing to demonstrate the efficiency of the new approach.