Time-Domain Analysis of Body Freedom Flutter Based on 6DOF Equation

The reduced weight and improved efficiency of modern aeronautical structures result in a decreasing separation of frequency ranges of rigid and elastic modes.Particularly,a high-aspect-ratio flexible flying wing is prone to body freedomflutter(BFF),which is a result of coupling of the rigid body short-periodmodewith 1st wing bendingmode.Accurate prediction of the BFF characteristics is helpful to reflect the attitude changes of the vehicle intuitively and design the active flutter suppression control law.Instead of using the rigid body mode,this work simulates the rigid bodymotion of the model by using the six-degree-of-freedom(6DOF)equation.A dynamicmesh generation strategy particularly suitable for BFF simulation of free flying aircraft is developed.An accurate Computational Fluid Dynamics/Computational Structural Dynamics/six-degree-of-freedom equation(CFD/CSD/6DOF)-based BFF prediction method is proposed.Firstly,the time-domain CFD/CSD method is used to calculate the static equilibrium state of the model.Based on this state,the CFD/CSD/6DOF equation is solved in time domain to evaluate the structural response of themodel.Then combinedwith the variable stiffnessmethod,the critical flutter point of the model is obtained.This method is applied to the BFF calculation of a flyingwing model.The calculation results of the BFF characteristics of the model agree well with those fromthe modalmethod andNastran software.Finally,the method is used to analyze the influence factors of BFF.The analysis results show that the flutter speed can be improved by either releasing plunge constraint or moving the center ofmass forward or increasing the pitch inertia.

Dynamic Coupling Analysis of Semisubmersible Platform Float-over Method for Docking Case

In this paper,the multi-body coupled dynamic characteristics of a semisubmersible platform and an HYSY 229 barge were investigated.First,coupled hydrodynamic analysis of the HYSY 229 barge and the semisubmersible platform was performed.Relevant hydrodynamic parameters were obtained using the retardation function method of three-dimensional frequency-domain potential flow theory.The results of the hydrodynamic analysis were highly consistent with the test findings,verifying the accuracy of the multifloating hydrodynamic coupling analysis,and key hydrodynamic parameters were solved for different water depths and the coupling effect.According to the obtained results,the hydrodynamic influence was the largest in shallow waters when the coupling effect was considered.Furthermore,the coupled motion equation combined with viscous damping,fender system,and mooring system was established,and the hydrodynamics,floating body motion,and dynamic response of the fender system were analyzed.Motion analysis revealed good agreement among the surge,sway,and yaw motions of the two floating bodies.However,when the wave period reached 10 s,the motion of the two floating bodies showed severe shock,and a relative motion was also observed.Therefore,excessive constraints should be added between the two floating bodies during construction to ensure construction safety.The numerical analysis and model test results of the semisubmersible platform and HYSY 229 barge at a water depth of 42 m and sea conditions of 0°,45°,and 90° were in good agreement,and the error was less than 5%.The maximum movement of the HYSY 229 barge reached 2.61 m in the sway direction,whereas that of the semisubmersible platform was 2.11 m.During construction,excessive constraints should be added between the two floating bodies to limit their relative movement and ensure construction safety.

Interval Finite Element Analysis of Wing Flutter

The influences of uncertainties in structural parameters on the flutter speed of wing are studied. On the basis of the deterministic flutter analysis model of wing, the uncertainties in structural parameters are considered and described by interval numbers. By virtue of first-order Taylor series expansion, the lower and upper bound curves of the transient decay rate coefficient versus wind velocity are given. So the interval estimation of the flutter critical wind speed of wing can be obtained, which is more reasonable than the point esti- mation obtained by the deterministic flutter analysis and provides the basis for the further non-probabilistic interval reliability analysis of wing flutter. The flow chart for interval fmite element model of flutter analysis of wing is given. The proposed interval finite element model and the stochastic finite element model for wing flutter analysis are compared by the examples of a three degrees of freedom airfoil and fuselage and a 15° sweptback wing, and the results have shown the effectiveness and feasibility of the presented model. The prominent advantage of the proposed interval finite element model is that only the bounds of uncertain parameters are required, and the probabilistic distribution densities or other statistical characteristics are not needed.

Structural design and modal behaviors analysis of a new swept baffled inflatable wing

Inflatable wing has significant application value in the design of loitering munitions because of its advantages such as lightweight and foldability.However,due to the flexible characteristics,aeroelastic behaviors of inflatable wings such as flutter are nonnegligible in flight.By designing a certain angle between the inflatable beam and the wing span,the structural dynamic and even the aeroelastic performance of the inflatable wing can be effectively improved.Based on the analysis of the mechanical and geometric characteristics of the inflatable structure,a new inflatable wing with sweep arranged inflatable beams is proposed,and the main design variables and methods are analyzed.For purpose of investigating the aeroelastic performance of the swept baffled inflatable wing,the modal behaviors by considering the wet mode are studied.In consideration of the deficiencies of the traditional wet modal analysis method,by introducing the influence on the additional stiffness of flow field,an added massstiffness method is proposed in this paper,and the advantages are verified by ground vibration experiments.On this basis,the effects of baffles sweep angle,pressure,and boundary conditions on the modal parameters and aeroelastic performance of inflatable wing are analyzed.The results show that the aeroelastic performance of the inflatable wing can be designed by changing the baffles sweep angle,which is enlightened for the aeroelastic tailoring design on inflatable wings.

Parameter sensitivities analysis for classical flutter speed of a horizontal axis wind turbine blade

The parameter sensitivities affecting the flutter speed of the NREL （National Renewable Energy Laboratory） 5-MW baseline HAWT （horizontal axis wind turbine） blades are analyzed. An aeroelastic model, which comprises an aerodynamic part to calculate the aerodynamic loads and a structural part to determine the structural dynamic responses, is established to describe the classical flutter of the blades. For the aerodynamic part, Theodorsen unsteady aerodynamics model is used. For the structural part, Lagrange’s equation is employed. The flutter speed is determined by introducing “V–g” method to the aeroelastic model, which converts the issue of classical flutter speed determination into an eigenvalue problem. Furthermore, the time domain aeroelastic response of the wind turbine blade section is obtained with employing Runge-Kutta method. The results show that four cases （i.e., reducing the blade torsional stiffness, moving the center of gravity or the elastic axis towards the trailing edge of the section, and placing the turbine in high air density area） will decrease the flutter speed. Therefore, the judicious selection of the four parameters （the torsional stiffness, the chordwise position of the center of gravity, the elastic axis position and air density） can increase the relative inflow speed at the blade section associated with the onset of flutter.

Match Point Solution for Robust Flutter Analysis in Constant-Mach Prediction

This paper presents a method for robust flutter computation which uses flight altitude as the perturbation variable in order to obtain a match point solution. The air density and sound speed of standard atmosphere model are approximated as the polynomial function of altitude, such that the flight altitude becomes the single perturbation variable that describes the aeroelastic system. The uncertainties of generalized stiffness and damping are considered and the uncertain aeroelastic system can be formulated as linear fractional transformation （LFT） representation which is suitable for/.t analysis framework. Finally, the match point solutions of robust flutter margins can be computed with structured singular value （SSV） theory. The robust flutter analysis method provided in this paper is suitable for constant-Mach flight flutter test and provides valuable reference for flight envelope expansion.

Isotopic analysis based on terahertz spectrum

As a new promising detection technology in the terahertz research field,the terahertz time-domain spec-troscopy(THz-TDS)has very broad application potential in many fields because its advantage on the characteristic spectrum,wide spectrum and non-destructive analysis of interested substances.In this paper,the terahertz absorption spectra of gases mixed with 12 CO and 13 CO in the spec-trum range of 0.5–2.5 THz are measured by terahertz time-domain spectroscopy for the first time.Several isotopo-logues can be clearly distinguished based on the difference in their rotational energies and the consequent terahertz spectrum.The experimental results show that 12 CO and 13 CO have obvious characteristic absorption peaks in the spectrum range of 0.5–2.5 THz due to the difference in rotational energy,and the rotational constant B can be calculated according to the experimental values to distin-guish the two gaseous isotopologues.The frequency posi-tions of the characteristic absorption peak measured by this experiment and the rotation constant B calculated accord-ing to the experimental values are compared with those previous theoretical calculations and experimental results,and they are in good agreement.This result lays a foun-dation for developing more sophisticated terahertz instru-ments to the detection of different isotopologues.

Time-Domain Analysis of Underground Station-Layered Soil Interaction Based on High-Order Doubly Asymptotic Transmitting Boundary

Based on the modified scale boundary finite element method and continued fraction solution,a high-order doubly asymptotic transmitting boundary(DATB)is derived and extended to the simulation of vector wave propagation in complex layered soils.The high-order DATB converges rapidly to the exact solution throughout the entire frequency range and its formulation is local in the time domain,possessing high accuracy and good efficiency.Combining with finite element method,a coupled model is constructed for time-domain analysis of underground station-layered soil interaction.The coupled model is divided into the near and far field by the truncated boundary,of which the near field is modelled by FEM while the far field is modelled by the high-order DATB.The coupled model is implemented in an open source finite element software,OpenSees,in which the DATB is employed as a super element.Numerical examples demonstrate that results of the coupled model are stable,accurate and efficient compared with those of the extended mesh model and the viscous-spring boundary model.Besides,it has also shown the fitness for long-time seismic response analysis of underground station-layered soil interaction.Therefore,it is believed that the coupled model could provide a new approach for seismic analysis of underground station-layered soil interaction and could be further developed for engineering.

Three-dimensional nonlinear flutter analysis of long-span suspension bridges during erection

In this work, the aerodynamic stability of the Yichang Suspension Bridge over Yangtze River during erection was determined by three dimensional nonlinear flutter analysis, in which the nonlinearities of structural dynamic characteristics and aeroelastic forces caused by large deformation are fully considered. An interesting result obtained was that the bridge was more stable when the stiffening girders were erected in a non symmetrical manner as opposed to the traditional symmetrical erection schedule. It was also found that the severe decrease in the aerodynamic stability was due to the nonlinear effects. Therefore, the nonlinear factors should be considered accurately in aerodynamic stability analysis of long span suspension bridges during erection.

Brillouin optical time-domain analysis for geotechnical monitoring

In this paper, we show some recent experimental applications of Brillouin optical time-domain analysis （BOTDA） based sensors for geotechnical monitoring. In particular, how these sensors can be applied to detecting early movements of soil slopes by the direct embedding of suitable fiber cables in the ground is presented. Furthermore, the same technology can be used to realize innovative inclinometers, as well as smart foundation anchors.

A Recursive Method of Time-Frequency Analysis for the Signal Processing of Flutter Test with Progression Variable Speed

Focused on the non-statlonarity and real-time analysis of signal in flutter test with progression variable speed （FTPVS）, a new method of recursive time-frequency analysis is presented. The time-varying system is tracked on-line by building a time-varying parameter model, and then the relevant parameter spectrum can be obtained. The feasibility and advantages of the method are examined by digital simulation. The results of FTPVS at low-speed wind-tunnel promise the engineering application perspective of the method.

Flutter Analysis of Aircraft Wing Using Equivalent-Plate Models with Orthogonal Polynomials

A method of equivalent simplification,using equivalent-plate models(EPMs),is developed.It is to achieve goals of rapid modeling and effective analysis in structural dynamics and flutter analysis of complex wing structures.It is on the assumption that the wing structures discussed are composed of skin,beams and ribs,and the different plate units(such as skin,beam web,rib web)are not distinguished in modeling,which is to avoid the complex pre-processing and make it more generalized.Taking the effect of transverse shear deformation into consideration,the equivalence is based on the first-order shear deformation theory,and it can import the model files of MSC/NASTRAN and process the information to accomplish the equivalent modeling.The Ritz method is applied with the Legendre polynomials,which is used to define the geometry,structure and displacements of the wing.Particularly,the selection of Legendre polynomials as trial functions brings good accuracy to the modeling and can avoid the ill-conditions.This is in contrast to the EPM method based on the classical plate theory.Through vibration and flutter analysis,the results obtained by using EPM agree well with those obtained by the finite element method,which indicates the accuracy and effectiveness in vibration and flutter analysis of the EPM method.

Adsorption and Transfer Properties for Toluene and pDichlorobenzene in Dense CO_2/Silica Gel Packed Bed System by Moment Method and Time-Domain Analysis

Based on chromatographic theory, the moment method and the time-domain fitting analysis were applied to measure and evaluate the adsorption equilibrium constant and mass transfer properties (axial dispersion coefficient and effective intra-particle diffusivity) for toluene and p-dichlorobenzene on silica gel adsorbent in the subcritical and supercritical CO2. An apparatus based on supercritical fluid chromatography was established and the experiments were performed at temperatures of 298.15-318.15 K and pressures of 7.5-17.8 MPa. The two methods have been compared. The results show that for the systems studied here the moment method can give reasonable values for both adsorption equilibrium constant and mass transfer properties, but the time-domain analysis only can obtain the adsorption equilibrium constant. The dependence of adsorption equilibrium constant and mass transfer properties on temperature and pressure was investigated.

Numerical analysis of multicomponent responses of surface-hole transient electromagnetic method

We calculate the multicomponent responses of surface-hole transient electromagnetic method. The methods and models are unsuitable as geoelectric models of conductive surrounding rocks because they are based on regular local targets. We also propose a calculation and analysis scheme based on numerical simulations of the subsurface transient electromagnetic fields. In the modeling of the electromagnetic fields, the forward modeling simulations are performed by using the finite-difference time-domain method and the discrete image method, which combines the Gaver–Stehfest inverse Laplace transform with the Prony method to solve the initial electromagnetic fields. The precision in the iterative computations is ensured by using the transmission boundary conditions. For the response analysis, we customize geoelectric models consisting of near-borehole targets and conductive wall rocks and implement forward modeling simulations. The observed electric fields are converted into induced electromotive force responses using multicomponent observation devices. By comparing the transient electric fields and multicomponent responses under different conditions, we suggest that the multicomponent-induced electromotive force responses are related to the horizontal and vertical gradient variations of the transient electric field at different times. The characteristics of the response are determined by the varying the subsurface transient electromagnetic fields, i.e., diffusion, attenuation and distortion, under different conditions as well as the electromagnetic fields at the observation positions. The calculation and analysis scheme of the response consider the surrounding rocks and the anomalous field of the local targets. It therefore can account for the geological data better than conventional transient field response analysis of local targets.

Distributed Fiber Optic Monitoring and Stability Analysis of a Model Slope under Surcharge Loading

In the discipline of geotechnical engineering, fiber optic sensor based distributed monitoring has played an increasingly important role over the past few decades. Compared with conventional sensors, fiber optic sensors have a variety of exclusive advantages, such as smaller size, higher precision, and better corrosion resistance. These innovative monitoring technologies have been successfully applied for performance monitoring of geo-structures and early warning of potential geo- hazards around the world. In order to investigate their ability to monitor slope stability problems, a medium-sized model of soil nailed slope has been constructed in laboratory. The fully distributed Brillouin optical time-domain analysis （BOTDA） sensing technology was employed to measure the horizontal strain distributions inside the model slope. During model construction, a specially designed strain sensing fiber was buried in the soil mass. Afterward, the surcharge loading was applied on the slope crest in stages using hydraulic jacks and a reaction frame. During testing, an NBX-6o5o BOTDA sensing interrogator was used to collect the fiber optic sensing data. The test results have been analyzed in detail, which shows that the fiber optic sensors can capture the progressive deformation and failure pattern of the model slope. The limit equilibrium analyses were also conducted to obtain the factors ofsafety of the slope under different surface loadings. It is found that the characteristic maximum strains can reflect the stability of the model slope and an empirical relationship was obtained, This study verified the effectiveness of the distributed BOTDA sensing technology in performance monitoring of slope.

Assessment of right ventricular afterload by pressure waveform analysis in acute pulmonary hypertension

AIM: To characterize hydraulic right ventricle (RV) af- terload by pulmonary arterial pressure waveform analy- sis in an acute pulmonary hypertension (PH) model. METHODS: Pulmonary artery (PA) flow and pressure were recorded in six anesthetized sheep. Acute iso- baric PH was induced by phenylephrine (active) and PA mechanical constriction (passive). We estimated the amplitude of the forward and reflected pressure waves according to the inflection point. In most cases the in- flection pressure was smooth, thus the inflection point was defined as the time at which the first derivative ofpulmonary arterial pressure reached its first minimum. We calculated the input and characteristic (Z C , time- domain Li method) impedances, the capacitance index (stroke volume/pulse pressure), the augmentation index (AI) (reflected pressure/pulse pressure), the frac- tional pulse pressure (pulse pressure/mean pressure) and the wasted energy generated by the RV due to wave reflection during ejection (E W ). RESULTS: Pulse pressure, fractional pulse pressure, AI and Z C increased and capacitance index decreased during passive PH with respect to control (P < 0.05). In contrast, Z C and the capacitance index did not change and E W and the AI decreased during active PH. Pulse pressure correlated with E W and Z C and the AI was cor- related with E W (r > 0.6, P < 0.05). CONCLUSION: PA pressure waveform analysis al- lows the quantification of the dynamic RV afterload. Prospective clinical studies will be necessary to validate this time-domain approach to evaluate the dynamic RV afterload in chronic PH.

Thermally Induced Vibration Analysis of Flexible Beams Based on Isogeometric Analysis

Spacecraft flexible appendages may experience thermally induced vibrations(TIV)under sudden heating loads,which in consequence will be unable to complete their intended missions.Isogeometric analysis(IGA)utilizes,in an isoparametric concept,the same high order and high continuity non-uniform rational B-splines(NURBS)to represent both the geometry and the physical field of the structure.Compared to the traditional Lagrange polynomial based finite element method where only C0-continuity across elements can be achieved,IGA is geometrically exact and naturally fulfills the C1-continuity requirement of Euler–Bernoulli(EB)beam elements,therefore,does not need extra rotational degrees-of-freedom.In this paper,we present a thermally induced vibration analysis framework based on the isogeometric method where thermal and structural behaviors are coupled.We fully exploited the higher order,higher continuous and geometric exactness of the NURBS basis with both benchmarks and sophisticated problems.In particular,we studied the thermally induced vibrations of the Hubble Space Telescope(HST)solar panel where main factors influencing thermal flutters are studied,and where possible improvements of the analytical reference methods are discussed.Additionally,thermally induced vibrations of the thin-walled lenticular tubes are studied and two new configurations of the tube are proposed to effectively suppress the thermally induced vibrations.Numerical examples of both benchmarks and sophisticated problems confirm the accuracy and efficiency of the isogeometric analysis framework for thermally induced vibration analysis of space structures.

Analysis of photonic crystal on LED extraction efficiency

Because of its very low light extraction efficiency(LEE),LED is limited to be widely used under the condition of the internal quantum efficiency which up to 90%.In order to fullfill the design of a more efficient GaN-based blue light LED,the model including deeply etched surface photonic crystals(PhCs)LED is discussed using mode analysis method from light waveguide theory.The distributions of all order modes in GaN layer are obtained by the effective index approximation.The light extraction efficiencies are also calculated by finite-difference time-domain method(FDTD).The emulated results fully coincide with the former analysis.Because the manufacture of the surface photonic crystal is feasible,the work can be very meaningful to design and manufacture the high efficiency GaN-based blue light LED in factory for a large amount.

开口断面钢-混结合梁悬索桥颤振特性及颤振形态研究

为准确评估开口断面双边箱钢-混结合梁悬索桥颤振性能和多模态参与效应,以典型的钢-混结合梁悬索桥——怀化洞庭溪沅水特大桥为背景,开展风洞模型试验及数值分析。制作加劲梁节段缩尺模型并进行风洞试验,测试原桥梁结构的颤振性能,分析调整阻尼比、采用气动措施及结构措施对结构颤振特性的影响,采用三维和二维两种颤振分析方法分析设置中央扣结构的颤振模态及多模态参与效应。结果表明:该钢-混结合梁悬索桥存在较明显的颤振起振点,且颤振是以扭转为主的弯扭耦合振动,同时其扭转和竖向振动存在较大的相位差(相差近180°);节段模型风洞试验与三维和二维颤振分析一致得到设置中央扣后,频率较高的正对称扭转模态先于反对称模态发生颤振,且三维与二维颤振分析结果相差不大。

PC-based artif icial neural network inversion for airborne time-domain electromagnetic data

Traditionally, airborne time-domain electromagnetic （ATEM） data are inverted to derive the earth model by iteration. However, the data are often highly correlated among channels and consequently cause ill-posed and over-determined problems in the inversion. The correlation complicates the mapping relation between the ATEM data and the earth parameters and thus increases the inversion complexity. To obviate this, we adopt principal component analysis to transform ATEM data into orthogonal principal components （PCs） to reduce the correlations and the data dimensionality and simultaneously suppress the unrelated noise. In this paper, we use an artificial neural network （ANN） to approach the PCs mapping relation with the earth model parameters, avoiding the calculation of Jacobian derivatives. The PC-based ANN algorithm is applied to synthetic data for layered models compared with data-based ANN for airborne time-domain electromagnetic inversion. The results demonstrate the PC-based ANN advantages of simpler network structure, less training steps, and better inversion results over data-based ANN, especially for contaminated data. Furthermore, the PC-based ANN algorithm effectiveness is examined by the inversion of the pseudo 2D model and comparison with data-based ANN and Zhody＇s methods. The results indicate that PC-based ANN inversion can achieve a better agreement with the true model and also proved that PC-based ANN is feasible to invert large ATEM datasets.