Study of rectangular waveguides with grooves of irregular shape using the finite element method

Electromagnetic wave propagation through grooved waveguides is studied using the finite element method (FEM). The effect of grooves of irregular shape on TE₁₀, TE₂₀ mode frequencies and passband is studied. The variation in cutoff frequencies for TE₁₀, TE₂₀ mode and passband is observed.     

二维压电材料动强度因子的扩展有限元计算

采用扩展有限元求解二维弹性压电材料动断裂问题。扩展有限元的网格独立于裂纹,因此网格生成可大大地简化,且裂纹扩展时不需重构网格。采用相互作用积分技术计算动强度因子。比较了标准的力裂尖加强函数和力-电裂尖加强函数对动强度因子的影响,结果表明标准的力裂尖加强函数能有效地分析压电材料动断裂问题。分析了极化方向对动强度因子的影响。数值分析表明采用扩展有限元获得的动强度因子与其他数值方法解吻合得很好。     

A least‐squares coupling method between a finite element code and a discrete element code

This paper presents a coupling method between a discrete element code CeaMka3D and a finite element code Sem. The coupling is based on a least‐squares method, which adds terms of forces to finite element code and imposes the velocity at coupling particles. For each coupling face, a small linear system with a constant matrix is solved. This method remains conservative in energy and shows good results in applications. Copyright © 2014 John Wiley & Sons, Ltd.     

A high‐order approach for modelling transient wave propagation problems using the scaled boundary finite element method

A high‐order time‐domain approach for wave propagation in bounded and unbounded domains is proposed. It is based on the scaled boundary FEM, which excels in modelling unbounded domains and singularities. The dynamic stiffness matrices of bounded and unbounded domains are expressed as continued‐fraction expansions, which leads to accurate results with only about three terms per wavelength. An improved continued‐fraction approach for bounded domains is proposed, which yields numerically more robust time‐domain formulations. The coefficient matrices of the corresponding continued‐fraction expansion are determined recursively. The resulting solution is suitable for systems with many DOFs as it converges over the whole frequency range, even for high orders of expansion. A scheme for coupling the proposed improved high‐order time‐domain formulation for bounded domains with a high‐order transmitting boundary suggested previously is also proposed. In the time‐domain, the coupled model corresponds to equations of motion with symmetric, banded and frequency‐independent coefficient matrices, which can be solved efficiently using standard time‐integration schemes. Numerical examples for modal and time‐domain analysis are presented to demonstrate the increased robustness, efficiency and accuracy of the proposed method. Copyright © 2013 John Wiley & Sons, Ltd.     

An improved perturbation method for stochastic finite element model updating

In this paper, an improved perturbation method is developed for the statistical identification of structural parameters by using the measured modal parameters with randomness. On the basis of the first‐order perturbation method and sensitivity‐based finite element (FE) model updating, two recursive systems of equations are derived for estimating the first two moments of random structural parameters from the statistics of the measured modal parameters. Regularization technique is introduced to alleviate the ill‐conditioning in solving the equations. The numerical studies of stochastic FE model updating of a truss bridge are presented to verify the improved perturbation method under three different types of uncertainties, namely natural randomness, measurement noise, and the combination of the two. The results obtained using the perturbation method are in good agreement with, although less accurate than, those obtained using the Monte Carlo simulation (MCS) method. It is also revealed that neglecting the correlation of the measured modal parameters may result in an unreliable estimation of the covariance matrix of updating parameters. The statistically updated FE model enables structural design and analysis, damage detection, condition assessment, and evaluation in the framework of probability and statistics. Copyright © 2007 John Wiley & Sons, Ltd.     

A boundary‐perturbation finite element approach for shape optimization

A numerical method is proposed for the efficient solution of shape optimization problems, which combines the boundary perturbation technique and finite element analysis. The method is computationally efficient in that it requires a number of finite element analyses with a fixed geometry, as opposed to standard shape optimization which requires re‐analysis with varying geometry. The application of the method to general shape optimization is considered. In addition, a special optimization scheme is devised for a class of problems governed by linear partial differential equations. The performance of the method is illustrated via an example which involves acoustic wave scattering from an obstacle. Copyright © 2000 John Wiley & Sons, Ltd.     

Modelling cohesive crack growth using a two-step finite element-scaled boundary finite element coupled method

A two-step method, coupling the finite element method (FEM) and the scaled boundary finite element method (SBFEM), is developed in this paper for modelling cohesive crack growth in quasi-brittle normal-sized structures such as concrete beams. In the first step, the crack trajectory is fully automatically predicted by a recently-developed simple remeshing procedure using the SBFEM based on the linear elastic fracture mechanics theory. In the second step, interfacial finite elements with tension-softening constitutive laws are inserted into the crack path to model gradual energy dissipation in the fracture process zone, while the elastic bulk material is modelled by the SBFEM. The resultant nonlinear equation system is solved by a local arc-length controlled solver. Two concrete beams subjected to mode-I and mixed-mode fracture respectively are modelled to validate the proposed method. The numerical results demonstrate that this two-step SBFEM-FEM coupled method can predict both satisfactory crack trajectories and accurate load-displacement relations with a small number of degrees of freedom, even for crack growth problems with strong snap-back phenomenon. The effects of the tensile strength, the mode-I and mode-II fracture energies on the predicted load-displacement relations are also discussed.     

Fully-automatic modelling of cohesive crack growth using a finite element-scaled boundary finite element coupled method

This study develops a method coupling the finite element method (FEM) and the scaled boundary finite element method (SBFEM) for fully-automatic modelling of cohesive crack growth in quasi-brittle materials. The simple linear elastic fracture mechanics (LEFM)-based remeshing procedure developed previously is augmented by inserting nonlinear interface finite elements automatically. The constitutive law of these elements is modelled by the cohesive/fictitious crack model to simulate the fracture process zone, while the elastic bulk material is modelled by the SBFEM. The resultant nonlinear equation system is solved by a local arc-length controlled solver. The crack is assumed to grow when the mode-I stress intensity factor K_I vanishes in the direction determined by LEFM criteria. Other salient algorithms associated with the SBFEM, such as mapping state variables after remeshing and calculating K_I using a “shadow subdomain”, are also described. Two concrete beams subjected to mode-I and mixed-mode fracture respectively are modelled to validate the new method. The results show that this SBFEM-FEM coupled method is capable of fully-automatically predicting both satisfactory crack trajectories and accurate load-displacement relations with a small number of degrees of freedom, even for problems with strong snap-back. Parametric studies were carried out on the crack incremental length, the concrete tensile strength, and the mode-I and mode-II fracture energies. It is found that the K_I ? 0 criterion is objective with respect to the crack incremental length.     

Modelling crack propagation in reinforced concrete using a hybrid finite element–scaled boundary finite element method

A previously developed hybrid finite element–scaled boundary finite element method (FEM–SBFEM) is extended to model multiple cohesive crack propagation in reinforced concrete. This hybrid method can efficiently extract accurate stress intensity factors from the semi-analytical solutions of SBFEM and is also flexible in remeshing multiple cracks. Crack propagation in the concrete bulk is modelled by automatically inserted cohesive interface elements with nonlinear softening laws. The concrete–reinforcement interaction is also modelled by cohesive interface elements. The bond shear stress–slip relation of CEB-FIP Model Code 90 and an empirical confining stress–crack opening relation are used to characterise slip and split failure at the concrete–reinforcement interface, respectively. Three RC beams were simulated. The numerical results agreed well with both experimental and numerical results available in the literature. Parametric studies demonstrated the importance of modelling both slip and split failure mechanisms at the concrete–reinforcement interface.     

Developments in ultrasonic modeling with finite element analysis

An overview is given of finite element analysis and its application to the modeling of ultrasonic nondestructive evaluation phenomena. Following a discussion of the underlying weighted residual methodology, a mass-lumping technique is described which results in an efficient computer implementation for 2D geometries. Code predictions are compared with both analytical and experimental results, and data from studies of attenuation, anisotropy, defect interactions, and surface waves are given. Initial results from a full 3D formulation are also shown.     

A discontinuous Galerkin finite element method for dynamic and wave propagation problems in non‐linear solids and saturated porous media

A time‐discontinuous Galerkin finite element method (DGFEM) for dynamics and wave propagation in non‐linear solids and saturated porous media is presented. The main distinct characteristic of the proposed DGFEM is that the specific P3–P1 interpolation approximation, which uses piecewise cubic (Hermite's polynomial) and linear interpolations for both displacements and velocities, in the time domain is particularly proposed. Consequently, continuity of the displacement vector at each discrete time instant is exactly ensured, whereas discontinuity of the velocity vector at the discrete time levels still remains. The computational cost is then obviously saved, particularly in the materially non‐linear problems, as compared with that required for the existing DGFEM. Both the implicit and explicit algorithms are developed to solve the derived formulations for linear and materially non‐linear problems. Numerical results illustrate good performance of the present method in eliminating spurious numerical oscillations and in providing much more accurate solutions over the traditional Galerkin finite element method using the Newmark algorithm in the time domain. Copyright © 2003 John Wiley & Sons, Ltd.     

Control of wave propagation response using quasi crystals: A formulation based on spectral finite element

This paper presents wave propagation studies in quasi crystal structures and quasi crystal reinforced aluminium structures. The analysis is performed using frequency domain spectral finite element formulation. The analysis considers different 2-D decagonal and 3-D icosahedral quasi crystals. First, wave propagation analysis of quasi crystal structure alone is performed and the propagation of phonon and phason modes for different quasi crystals are studied. The study includes the propagation of axial and transverse wave responses in these quasi crystals. The study has found that the amplitude of the phason modes is very small compared to the phonon modes and the increase of the phason mode content (through increase in R) increases the phason mode amplitude, without affecting the phonon mode amplitudes. It is shown that the dominant axial phonon mode is non-dispersive and the dominant flexural phonon mode is dispersive. In the next study, the aluminium beam structure is reinforced with different quasi crystals in different configurations and the wave propagation of axial and transverse responses are studied. For all the combinations of quasi crystal aluminium beam combination, there is substantial suppression of responses both for the axial and the bending responses. Unsymmetrical configuration produces substantial non-dominant phonon modes which propagate dispersively. It is found that for a symmetric bi-morph configuration, the response is reduced significantly, about 68% and 75% for axial loading and 80% and 78% for flexural loading, respectively, for the 2-D decagonal quasi crystal and the 3-D icosahedral quasi crystal.     

Strain and stress computations in stochastic finite element methods

This paper focuses on the computation of statistical moments of strains and stresses in a random system model where uncertainty is modeled by a stochastic finite element method based on the polynomial chaos expansion. It identifies the cases where this objective can be achieved by analytical means using the orthogonality property of the chaos polynomials and those where it requires a numerical integration technique. To this effect, the applicability and efficiency of several numerical integration schemes are considered. These include the Gauss–Hermite quadrature with the direct tensor product—also known as the Kronecker product—Smolyak's approximation of such a tensor product, Monte Carlo sampling, and the Latin Hypercube sampling method. An algorithm for reducing the dimensionality of integration under a direct tensor product is also explored for optimizing the computational cost and complexity. The convergence rate and algorithmic complexity of all of these methods are discussed and illustrated with the non‐deterministic linear stress analysis of a plate. Copyright © 2007 John Wiley & Sons, Ltd.     

Mixed enrichment for the finite element method in heterogeneous media

Problems of multiple scales of interest or of locally nonsmooth solutions may often involve heterogeneous media. These problems are usually very demanding in terms of computations with the conventional finite element method. On the other hand, different enriched finite element methods such as the partition of unity, which proved to be very successful in treating similar problems, are developed and studied for homogeneous media. In this work, we present a new idea to extend the partition of unity finite element method to treat heterogeneous materials. The idea is studied in applications to wave scattering and heat transfer problems where significant advantages are noted over the standard finite element method. Although presented within the partition of unity context, the same enrichment idea can also be extended to other enriched methods to deal with heterogeneous materials. Copyright © 2014 John Wiley & Sons, Ltd.     

The perturbation method and the extended finite element method. An application to fracture mechanics problems

The extended finite element method has been successful in the numerical simulation of fracture mechanics problems. With this methodology, different to the conventional finite element method, discretization of the domain with a mesh adapted to the geometry of the discontinuity is not required. On the other hand, in traditional fracture mechanics all variables have been considered to be deterministic (uniquely defined by a given numerical value). However, the uncertainty associated with these variables (external loads, geometry and material properties, among others) it is well known. This paper presents a novel application of the perturbation method along with the extended finite element method to treat these uncertainties. The methodology has been implemented in a commercial software and results are compared with those obtained by means of a Monte Carlo simulation.     

Temporal‐spatial dispersion and stability analysis of finite element method in explicit elastodynamics

This work presents the temporal‐spatial (full) dispersion and stability analysis of plane square linear and biquadratic serendipity finite elements in explicit numerical solution of transient elastodynamic problems. Here, the central difference method, as an explicit time integrator, is exploited. The paper complements and extends the previous work on spatial/grid dispersion analysis of plane square biquadratic serendipity finite elements. We report on a computational strategy for temporal‐spatial dispersion relationships, where eigenfrequencies from grid/spatial dispersion analysis are adjusted to comply with the time integration method. Besides that, an ‘optimal’ lumped mass matrix for the studied finite element types is proposed and investigated. Based on the temporal‐spatial dispersion and stability analysis, relationships suggesting the ‘proper’ choice of mesh size and time step size from knowledge of the loading spectrum are presented. Copyright © 2015 John Wiley & Sons, Ltd.     

A scaled boundary finite element formulation for dynamic elastoplastic analysis

This study presents the development of the scaled boundary finite element method (SBFEM) to simulate elastoplastic stress wave propagation problems subjected to transient dynamic loadings. Material nonlinearity is considered by first reformulating the SBFEM to obtain an explicit form of shape functions for polygons with an arbitrary number of sides. The material constitutive matrix and the residual stress fields are then determined as analytical polynomial functions in the scaled boundary coordinates through a local least squares fit to evaluate the elastoplastic stiffness matrix and the residual load vector semianalytically. The treatment of the inertial force within the solution of the nonlinear system of equations is also presented within the SBFEM framework. The nonlinear equation system is solved using the unconditionally stable Newmark time integration algorithm. The proposed formulation is validated using several benchmark numerical examples.     

Generalized finite element method using proper orthogonal decomposition

A methodology is presented for generating enrichment functions in generalized finite element methods (GFEM) using experimental and/or simulated data. The approach is based on the proper orthogonal decomposition (POD) technique, which is used to generate low‐order representations of data that contain general information about the solution of partial differential equations. One of the main challenges in such enriched finite element methods is knowing how to choose, a priori, enrichment functions that capture the nature of the solution of the governing equations. POD produces low‐order subspaces, that are optimal in some norm, for approximating a given data set. For most problems, since the solution error in Galerkin methods is bounded by the error in the best approximation, it is expected that the optimal approximation properties of POD can be exploited to construct efficient enrichment functions. We demonstrate the potential of this approach through three numerical examples. Best‐approximation studies are conducted that reveal the advantages of using POD modes as enrichment functions in GFEM over a conventional POD basis. Copyright © 2009 John Wiley & Sons, Ltd.     

Modelling dynamic crack propagation using the scaled boundary finite element method

This study presents a novel application of the scaled boundary finite element method (SBFEM) to model dynamic crack propagation problems. Accurate dynamic stress intensity factors are extracted directly from the semi‐analytical solutions of SBFEM. They are then used in the dynamic fracture criteria to determine the crack‐tip position, velocity and propagation direction. A simple, yet flexible remeshing algorithm is used to accommodate crack propagation. Three dynamic crack propagation problems that include mode‐I and mix‐mode fracture are modelled. The results show good agreement with experimental and numerical results available in the literature. It is found that the developed method offers some advantages over conventional FEM in terms of accuracy, efficiency and ease of implementation. Copyright © 2011 John Wiley & Sons, Ltd.     

Analysis of deformed microstrip resonator using the finite element method

Microstrip resonator with shape deformation has been analysed using the finite element method (FEM). Keeping the front surface of the microstrip resonator fixed, the back surface is deformed. Five cases of the deformation of the back surface are considered: (i) the left vertical side of the back surface is shifted inward; (ii) both vertical sides of the back surface are shifted inward; (iii) the right vertical side of the back surface is shifted outward; (iv) both vertical sides of the back surface are shifted outward; (v) left vertical side of back surface is shifted inward and the right vertical side of it is shifted outward. Variation in cutoff frequency for the TM₀₁₁, TM₁₁₀, TM₁₁₁, TM₀₁₂, TM₁₁₂ and TM₂₁₀ is observed.