Search results

In this work, an SFEM is proposed for solving acoustic problems by redistributing the entries in the mass matrix to “tune” the balance between “stiffness” and “mass” of discrete equation systems, aiming to minimize the dispersion error. The paper aims to discuss this issue.

This is done by simply shifting the four integration points’ locations when computing the entries of the mass matrix in the scheme of SFEM, while ensuring the mass conservation. The proposed method is devised for bilinear quadratic elements.

The balance between “stiffness” and “mass” of discrete equation systems is critically important in simulating wave propagation problems such as acoustics. A formula is also derived for possibly the best mass redistribution in terms of minimizing dispersion error reduction. Both theoretical and numerical examples demonstrate that the present method possesses distinct advantages compared with the conventional SFEM using the same quadrilateral mesh.

After introducing the mass-redistribution technique, the magnitude of the leading relative dispersion error (the quadratic term) of MR-SFEM is bounded by (5/8), which is much smaller than that of original SFEM models with traditional mass matrix (13/4) and consistence mass matrix (2). Owing to properly turning the balancing between stiffness and mass, the MR-SFEM achieves higher accuracy and much better natural eigenfrequencies prediction than the original SFEM does.

This paper aims to provide a comprehensive review of uncertainty quantification methods supported by evidence-based comparison studies. Uncertainties are widely encountered in engineering practice, arising from such diverse sources as heterogeneity of materials, variability in measurement, lack of data and ambiguity in knowledge. Academia and industries have long been researching for uncertainty quantification (UQ) methods to quantitatively account for the effects of various input uncertainties on the system response. Despite the rich literature of relevant research, UQ is not an easy subject for novice researchers/practitioners, where many different methods and techniques coexist with inconsistent input/output requirements and analysis schemes.

This confusing status significantly hampers the research progress and practical application of UQ methods in engineering. In the context of engineering analysis, the research efforts of UQ are most focused in two largely separate research fields: structural reliability analysis (SRA) and stochastic finite element method (SFEM). This paper provides a state-of-the-art review of SRA and SFEM, covering both technology and application aspects. Moreover, unlike standard survey papers that focus primarily on description and explanation, a thorough and rigorous comparative study is performed to test all UQ methods reviewed in the paper on a common set of reprehensive examples.

Over 20 uncertainty quantification methods in the fields of structural reliability analysis and stochastic finite element methods are reviewed and rigorously tested on carefully designed numerical examples. They include FORM/SORM, importance sampling, subset simulation, response surface method, surrogate methods, polynomial chaos expansion, perturbation method, stochastic collocation method, etc. The review and comparison tests comment and conclude not only on accuracy and efficiency of each method but also their applicability in different types of uncertainty propagation problems.

The research fields of structural reliability analysis and stochastic finite element methods have largely been developed separately, although both tackle uncertainty quantification in engineering problems. For the first time, all major uncertainty quantification methods in both fields are reviewed and rigorously tested on a common set of examples. Critical opinions and concluding remarks are drawn from the rigorous comparative study, providing objective evidence-based information for further research and practical applications.

Since the most of structures and structural components suffers from cyclic loadings, the study on the fatigue failure due to the crack growth has a great importance. The purpose of this paper is to present a three-dimensional fatigue crack growth simulation of embedded cracks using s-version finite element method (SFEM). Using the numerical results, the validity of the fitness-for-service (FFS) code evaluation method is verified.

In this paper, three-dimensional fatigue crack propagation analysis of embedded cracks is performed using the SFEM. SFEM is a numerical analysis method in which the shape of the structure is represented by a global mesh, and cracks are modeled by local meshes independently. The independent global and local meshes are superimposed to obtain the displacement solution of the problem simultaneously.

The fatigue crack growth of arbitrary shape of cracks is slow compared to that of the simplified circular crack and the crack approximated based on the FFS code of the Japan Society of Mechanical Engineers (JSME). The results tell us that the FFS code of JSME can provide a conservative evaluation of the fatigue crack growth and the residual life time.

This paper presents a three-dimensional fatigue crack growth simulation of embedded cracks using SFEM. Using this method, it is possible to apply mixed mode loads to complex shaped cracks that are closer to realistic conditions.

The purpose of the study is to present a frequency domain spectral finite element model (SFEM) based on fast Fourier transform (FFT) for wave propagation analysis of smart laminated composite beams with embedded delamination. For generating and sensing high-frequency elastic waves in composite beams, piezoelectric materials such as lead zirconate titanate (PZT) are used because they can act as both actuators and sensors. The present model is used to investigate the effects of parametric variation of delamination configuration on the propagation of fundamental anti-symmetric wave mode in piezoelectric composite beams.

The spectral element is derived from the exact solution of the governing equation of motion in frequency domain, obtained through fast Fourier transformation of the time domain equation. The beam is divided into two sublaminates (delamination region) and two base laminates (integral regions). The delamination region is modeled by assuming constant and continuous cross-sectional rotation at the interfaces between the base laminate and sublaminates. The governing differential equation of motion for delaminated composite beam with piezoelectric lamina is obtained using Hamilton’s principle by introducing an electrical potential function.

A detailed study of the wave response at the sensor shows that the A0 mode can be used for delamination detection in a wide region and is more suitable for detecting small delamination. It is observed that the amplitude and time of arrival of the reflected A0 wave from a delamination are strongly dependent on the size, position of the delamination and the stacking sequence. The degraded material properties because of the loss of stiffness and density in damaged area differently alter the S0 and A0 wave response and the group speed. The present method provides a potential technique for researchers to accurately model delaminations in piezoelectric composite beam structures. The delamination position can be identified if the time of flight of a reflected wave from delamination and the wave propagation speed of A0 (or S0) mode is known.

Spectral finite element modeling of delaminated composite beams with piezoelectric layers has not been reported in the literature yet. The spectral element developed is validated by comparing the present results with those available in the literature. The spectral element developed is then used to investigate the wave propagation characteristics and interaction with delamination in the piezoelectric composite beam.

To generalize the traditional 2nd order stochastic perturbation technique for input random variables and fields and to demonstrate for flow problems.

The methodology is based on an n‐th order expansion (perturbation) for input random parameters and state functions around their expected value to recover probabilistic moments of the response. A finite element formulation permits stochastic simulations on irregular meshes for practical applications.

The methodology permits approximation of expected values and covariances of quantities such as the fluid pressure and flow velocity using both symbolic and discrete FEM computations. It is applied to inviscid irrotational flow, Poiseulle flow and viscous Couette flow with randomly perturbed boundary conditions, channel height and fluid viscosity to illustrate the scheme.

The focus of the present work is on the basic concepts as a foundation for extension to engineering applications. The formulation for the viscous incompressible problem can be implemented by extending a 3D viscous primitive variable finite element code as outlined in the paper. For the case where the physical parameters are temperature dependent this will necessitate solution of highly non‐linear stochastic differential equations.

Techniques presented here provide an efficient approach for numerical analyses of heat transfer and fluid flow problems, where input design parameters and/or physical quantities may have small random fluctuations. Such an analysis provides a basis for stochastic computational reliability analysis.

The mathematical formulation and computational implementation of the generalized perturbation‐based stochastic finite element method (SFEM) is the main contribution of the paper.

To develop a new method for estimation of damage configuration in composite laminate structure using acoustic wave propagation signal and a reduction‐prediction neural network to deal with high dimensional spectral data.

A reduction‐prediction network, which is a combination of an independent component analysis (ICA) and a multi‐layer perceptron (MLP) neural network, is proposed to quantify the damage state related to transverse matrix cracking in composite laminates using acoustic wave propagation model. Given the Fourier spectral response of the damaged structure under frequency band‐selective excitation, the problem is posed as a parameter estimation problem. The parameters are the stiffness degradation factors, location and approximate size of the stiffness‐degraded zone. A micro‐mechanics model based on damage evolution criteria is incorporated in a spectral finite element model (SFEM) for beam type structure to study the effect of transverse matrix crack density on the acoustic wave response. Spectral data generated by using this model is used in training and testing the network. The ICA network called as the reduction network, reduces the dimensionality of the broad‐band spectral data for training and testing and sends its output as input to the MLP network. The MLP network, in turn, predicts the damage parameters.

Numerical demonstration shows that the developed network can efficiently handle high dimensional spectral data and estimate the damage state, damage location and size accurately.

Only numerical validation based on a damage model is reported in absence of experimental data. Uncertainties during actual online health monitoring may produce errors in the network output. Fault‐tolerance issues are not attempted. The method needs to be tested using measured spectral data using multiple sensors and wide variety of damages.

The developed network and estimation methodology can be employed in practical structural monitoring system, such as for monitoring critical composite structure components in aircrafts, spacecrafts and marine vehicles.

A new method is reported in the paper, which employs the previous works of the authors on SFEM and neural network. The paper addresses the important problem of high data dimensionality, which is of significant importance from practical engineering application viewpoint.

To provide an explicit representation for wide‐sense stationary stochastic fields which can be used in stochastic finite element modelling to describe random material properties.

This method represents wide‐sense stationary stochastic fields in terms of multiple Fourier series and a vector of mutually uncorrelated random variables, which are obtained by minimizing the mean‐squared error of a characteristic equation and solving a standard algebraic eigenvalue problem. The result can be treated as a semi‐analytic solution of the Karhunen‐Loève expansion.

According to the Karhunen‐Loève theorem, a second‐order stochastic field can be decomposed into a random part and a deterministic part. Owing to the harmonic essence of wide‐sense stationary stochastic fields, the decomposition can be effectively obtained with the assistance of multiple Fourier series.

The proposed explicit representation of wide‐sense stationary stochastic fields is accurate, efficient and independent of the real shape of the random structure in consideration. Therefore, it can be readily applied in a variety of stochastic finite element formulations to describe random material properties.

This paper discloses the connection between the spectral representation theory of wide‐sense stationary stochastic fields and the Karhunen‐Loève theorem of general second‐order stochastic fields, and obtains a Fourier‐Karhunen‐Loève representation for the former stochastic fields.

The simulation of eddy currents in laminated iron cores by the finite element method (FEM) is of great interest in the design of electrical devices. Modeling each laminate by finite elements leads to extremely large nonlinear systems of equations impossible to solve with present computer resources reasonably. The purpose of this study is to show that the multiscale finite element method (MSFEM) overcomes this difficulty.

A new MSFEM approach for eddy currents of laminated nonlinear iron cores in three dimensions based on the magnetic vector potential is presented. How to construct the MSFEM approach in principal is shown. The MSFEM with the Biot–Savart field in the frequency domain, a higher-order approach, the time stepping method and with the harmonic balance method are introduced and studied.

Various simulations demonstrate the feasibility, efficiency and versatility of the new MSFEM.

The novel MSFEM solves true three-dimensional eddy current problems in laminated iron cores taking into account of the edge effect.

The main purpose of the paper is to propose a new approach to stochastic computational modeling of interface defects in fiber‐reinforced composites. Interface defects with random radius and total number at the fiber‐matrix interface are modeled as an interphase between original composite components with the thickness obeying all the discontinuities and material parameters of this new, fictitious material are obtained by modified spatial averaging method. Such a model is used in the stochastic finite element analysis of composites in their original configuration. Next, the probabilistic moments of global effective properties of the entire composite are estimated, thanks to the traditional Monte Carlo simulation method implementation. Numerical experiments show that introduction of the interface defects results in significant increase of randomness level of the composite displacements and the homogenized elastic characteristics. Computer programs implemented can find their applications in digital image‐based analysis and the reliability analyses for fiber‐reinforced composites.

The purpose of this paper is to discuss the characteristics of several stochastic simulation methods applied in computation issue of structure health monitoring (SHM).

On the basis of the previous studies, this research focusses on four promising methods: transitional Markov chain Monte Carlo (TMCMC), slice sampling, slice-Metropolis-Hasting (M-H), and TMCMC-slice algorithm. The slice-M-H is the improved slice sampling algorithm, and the TMCMC-slice is the improved TMCMC algorithm. The performances of the parameters samples generated by these four algorithms are evaluated using two examples: one is the numerical example of a cantilever plate; another is the plate experiment simulating one part of the mechanical structure.

Both the numerical example and experiment show that, identification accuracy of slice-M-H is higher than that of slice sampling; and the identification accuracy of TMCMC-slice is higher than that of TMCMC. In general, the identification accuracy of the methods based on slice (slice sampling and slice-M-H) is higher than that of the methods based on TMCMC (TMCMC and TMCMC-slice).

The stochastic simulation methods evaluated in this paper are mainly two categories of representative methods: one introduces the intermediate probability density functions, and another one is the auxiliary variable approach. This paper provides important references about the stochastic simulation methods to solve the ill-conditioned computation issue, which is commonly encountered in SHM.