Toward a multifrequency quasi‐static Ritz vector method for frequency‐dependent acoustic system application

Computational issues concerning the calculation of acoustic responses of a complex finite element (FE) model for various noise and vibration inputs have become prevalent. Such a model requires a significant amount of computation time because of repeated inversions of dynamic stiffness matrices. Thus, even state‐of‐the‐art computer hardware and software often face limitations where a model order reduction (MOR) scheme can help. The established MOR schemes such as Ritz vector or quasi‐static Ritz vector methods are efficient for general engineering systems, but these MOR methods become inaccurate for frequency response analyses in some acoustic systems with frequency‐dependent mass and stiffness matrices and force vectors (hereinafter frequency‐dependent acoustic systems). To cope with the inaccurate prediction by these methods for frequency‐dependent acoustic systems, this research presents and applies the multifrequency quasi‐static Ritz vector method. Unlike the Ritz vector or quasi‐static Ritz vector methods, the present multifrequency quasi‐static Ritz vector method employs direct Krylov subspace bases without an orthonormal procedure at multiple center frequencies. In comparison with the existing MOR scheme, a significant gain in computational efficiency is achieved, as well as enhanced accuracy. A comparison of these methods based on criteria such as efficiency, accuracy, and reliability was also conducted. Copyright © 2011 John Wiley & Sons, Ltd.     

高斯白噪声激励的阻尼耦合的两个杜芬-范德波振子的近似瞬态响应

研究受高斯白噪声作用下的非线性阻尼耦合的两个杜芬-范德波振子响应的近似瞬态概率密度。应用基于广义谐和函数的随机平均法,并将幅值的瞬态概率密度的近似解表示为拉盖尔正交基函数的级数和,其中系数是随时间变化的,用Galerkin法可得到幅值的近似瞬态概率密度,从而得到状态变量响应的近似瞬态概率密度,数值模拟结果表明该方法有很好的适用性及精度。     

Isogeometric shape optimization for quasi‐static processes

A framework to solve shape optimization problems for quasi‐static processes is developed and implemented numerically within the context of isogeometric analysis (IGA). Recent contributions in shape optimization within IGA have been limited to static or steady‐state loading conditions. In the present contribution, the formulation of shape optimization is extended to include time‐dependent loads and responses. A general objective functional is used to accommodate both structural shape optimization and passive control for mechanical problems. An adjoint sensitivity analysis is performed at the continuous level and subsequently discretized within the context of IGA. The methodology and its numerical implementation are tested using benchmark static problems of optimal shapes of orifices in plates under remote bi‐axial tension and pure shear. Under quasi‐static loading conditions, the method is validated using a passive control approach with an a priori known solution. Several applications of time‐dependent mechanical problems are shown to illustrate the capabilities of this approach. In particular, a problem is considered where an external load is allowed to move along the surface of a structure. The shape of the structure is modified in order to control the time‐dependent displacement of the point where the load is applied according to a pre‐specified target. Copyright © 2015 John Wiley & Sons, Ltd.     

A quasi‐linear reproducing kernel particle method

Reproducing kernel particle method (RKPM) has been applied to many large deformation problems. RKPM relies on polynomial reproducing conditions to yield desired accuracy and convergence properties but requires appropriate kernel support coverage of neighboring nodes to ensure kernel stability. This kernel stability condition is difficult to achieve for problems with large particle motion such as the fragment‐impact processes that commonly exist in extreme events. A new reproducing kernel formulation with ‘quasi‐linear’ reproducing conditions is introduced. In this approach, the first‐order polynomial reproducing conditions are approximately enforced to yield a nonsingular moment matrix. With proper error control of the first‐order completeness, nearly second‐order convergence rate in L2 norm can be achieved while maintaining kernel stability. The effectiveness of this quasi‐linear RKPM formulation is demonstrated by modeling several extremely large deformation and fragment‐impact problems. Copyright © 2016 John Wiley & Sons, Ltd.     

An upscale theory of particle simulation for two‐dimensional quasi‐static problems

An upscale theory of the particle simulation, which is based on the distinct element method, is presented for two‐dimensional quasi‐static problems. Since the present upscale theory is comprised of four similarity criteria between different length‐scale particle‐simulation models, it reveals the intrinsic relationship between the particle‐simulation solution obtained from a small length‐scale (e.g. a laboratory length‐scale) model and that obtained from a large length‐scale (e.g. a geological length‐scale) one. The present upscale theory is of significant theoretical value in the particle simulation of two‐dimensional systems, at least from the following two points of view. (1) If the mechanical response of a particle model of a small length‐scale is used to indirectly investigate that of a large length‐scale, then the present upscale theory provides the necessary conditions, under which the particle model of the small length‐scale needs to be satisfied so that a similarity between the mechanical responses of two different length‐scale particle models can be maintained. (2) If a particle model of a large length‐scale is used to directly investigate the mechanical response of the model, then the present upscale theory can be used to determine the necessary particle‐scale mechanical properties from the macroscopic mechanical properties that are obtained from either a laboratory test or an in situ measurement. The related simulation results from two typical examples of significantly different length‐scales (i.e. a metre‐scale and a kilometre‐scale) have demonstrated the usefulness and correctness of the proposed upscale theory for simulating different length‐scale problems in quasi‐static geological systems. Copyright © 2007 John Wiley & Sons, Ltd.     

A quasi‐interpolation method for solving stiff ordinary differential equations

Based on the idea of quasi‐interpolation and radial basis functions approximation, a numerical method is developed to quasi‐interpolate the forcing term of differential equations by using radial basis functions. A highly accurate approximation for the solution can then be obtained by solving the corresponding fundamental equation and a small size system of equations related to the initial or boundary conditions. This overcomes the ill‐conditioning problem resulting from using the radial basis functions as a global interpolant. Error estimation is given for a particular second‐order stiff differential equation with boundary layer. The result of computations indicates that the method can be applied to solve very stiff problems. With the use of multiquadric, a special class of radial basis functions, it has been shown that a reasonable choice for the optimal shape parameter is obtained by taking the same value of the shape parameter as the perturbed parameter contained in the stiff equation. Copyright © 2000 John Wiley & Sons, Ltd.     

Computational electro‐elasticity and magneto‐elasticity for quasi‐incompressible media immersed in free space

In this work, a mixed variational formulation to simulate quasi‐incompressible electro‐active or magneto‐active polymers immersed in the surrounding free space is presented. A novel domain decomposition is used to disconnect the primary coupled problem and the arbitrary free‐space mesh update problem. Exploiting this decomposition, we describe a block‐iterative approach to solving the linearised multiphysics problem, and a physically and geometrically based, three‐parameter method to update the free space mesh. Several application‐driven example problems are implemented to demonstrate the robustness of the mixed formulation for both electro‐elastic and magneto‐elastic problems involving both finite deformations and quasi‐incompressible media. Copyright © 2016 John Wiley & Sons, Ltd.     

A continuum damage mechanics method for fracture simulation of quasi‐brittle materials

This paper addresses a novel continuum damage‐based method for simulating failure process of quasi‐brittle materials starting from local damage initiation to final fracture. In the developed method, the preset characteristic length field is used to evaluate damage instead of element, which is used to reduce the spurious sensitivity. In addition, damage is only updated in the most dangerous location at a time for considering stress redistribution due to damage evolution, which is used to simulate competitive fracture process. As cases study, representative numerical simulations of two benchmark tests are given to verify the performance of the developed continuum damage‐based method together with a used damage model. The simulation results of the crack paths for two concrete specimens obtained from the developed method matched well with the corresponding experimental results. The results show that the developed continuum damage‐based method is effective and can be used to simulate damage and fracture process of brittle or quasi‐brittle materials. And the simulation results based on the developed method depend only the preset characteristic length field and not grid mesh.     

A 4‐node assumed strain quasi‐conforming shell element with 6 degrees of freedom

Quasi‐conforming formulations of 4‐node stress‐resultant shell elements are presented. The element formulations use interrelated displacement–rotation interpolations. The formulation also includes drilling degrees of freedom, which improves membrane behavior and allows the modeling of stiffened plates and shells. The proposed treatment for bending provides very good results in the 4‐node shell element. The stiffness matrices for the present elements are explicitly expressed and the stresses are taken accurately at the nodal points. Compared to elements using Gauss integration, where the stresses are most accurate at the integration points, the extrapolation procedure needed for post‐processing is eliminated in the present shell element. A lot of numerical tests were carried out for the validation of the present 4‐node shell element and the results are in good agreement with references. Copyright © 2003 John Wiley & Sons, Ltd.     

Multiscale domain decomposition analysis of quasi‐brittle heterogeneous materials

A hybrid multiscale framework is presented, which processes the material scales in a concurrent manner, borrowing features from hierarchical multiscale methods. The framework is used for the analysis of non‐linear heterogeneous materials and is capable of tackling strain localization and failure phenomena. Domain decomposition techniques, such as the ?nite element tearing and interconnecting method, are used to partition the material in a number of non‐overlapping domains and adaptive re?nement is performed at those domains that are affected by damage processes. This re?nement is performed in terms of material scale and ?nite element size. It is veri?ed that the results are independent of the chosen domain decomposition. Moreover, the multiscale analyses are validated with reference solutions obtained with a full ?ne‐scale solution procedure. Copyright © 2011 John Wiley & Sons, Ltd.     

The fracture process in quasi‐brittle materials simulated using a lattice dynamical model

The damage process in quasi‐brittle materials is characterized by the evolution of a micro‐crack field, followed by the joining of micro‐cracks, stress localization and crack instability. In network models, masses are lumped at nodal points which are interconnected by one‐dimensional elements with a bilinear constitutive relation, considering the energy consistency during the simulated process. In order to replicate the material imperfections, to render a realistic behaviour in damage localization, the model has not only random elastic and rupture properties, but also a geometric perturbation. In the present paper 2D plates with different levels of brittleness are simulated. The numerical results are presented in terms of global stress vs strain diagram, final network configuration, energy balance during the process and as geometric damage evolution. Therefore, the predictive potential of the lattice discrete element model to capture fracture processes in quasi‐brittle materials is demonstrated.     

A Fourier transformation‐based method for gradient‐enhanced modeling of fatigue

A key limitation of the most constitutive models that reproduce a degradation of quasi‐brittle materials is that they generally do not address issues related to fatigue. One reason is the huge computational costs to resolve each load cycle on the structural level. The goal of this paper is the development of a temporal integration scheme, which significantly increases the computational efficiency of the finite element method in comparison to conventional temporal integrations. The essential constituent of the fatigue model is an implicit gradient‐enhanced formulation of the damage rate. The evolution of the field variables is computed as a multiscale Fourier series in time. On a microchronological scale attributed to single cycles, the initial boundary value problem is approximated by linear BVPs with respect to the Fourier coefficients. Using the adaptive cycle jump concept, the obtained damage rates are transferred to a coarser macrochronological scale associated with the duration of material deterioration. The performance of the developed method is hence improved due to an efficient numerical treatment of the microchronological problem in combination with the cycle jump technique on the macrochronological scale. Validation examples demonstrate the convergence of the obtained solutions to the reference simulations while significantly reducing the computational costs.     

Incremental‐secant modulus iteration scheme and stress recovery for simulating cracking process in quasi‐brittle materials using XFEM

In this paper, an incremental‐secant modulus iteration scheme using the extended/generalized finite element method (XFEM) is proposed for the simulation of cracking process in quasi‐brittle materials described by cohesive crack models whose softening law is composed of linear segments. The leading term of the displacement asymptotic field at the tip of a cohesive crack (which ensures a displacement discontinuity normal to the cohesive crack face) is used as the enrichment function in the XFEM. The opening component of the same field is also used as the initial guess opening profile of a newly extended cohesive segment in the simulation of cohesive crack propagation. A statically admissible stress recovery (SAR) technique is extended to cohesive cracks with special treatment of non‐homogeneous boundary tractions. The application of locally normalized co‐ordinates to eliminate possible ill‐conditioning of SAR, and the influence of different weight functions on SAR are also studied. Several mode I cracking problems in quasi‐brittle materials with linear and bilinear softening laws are analysed to demonstrate the usefulness of the proposed scheme, as well as the characteristics of global responses and local fields obtained numerically by the XFEM. Copyright © 2006 John Wiley & Sons, Ltd.     

A quasi‐implicit characteristic–based penalty finite‐element method for incompressible laminar viscous flows

In this paper, a novel characteristic–based penalty (CBP) scheme for the finite‐element method (FEM) is proposed to solve 2‐dimensional incompressible laminar flow. This new CBP scheme employs the characteristic‐Galerkin method to stabilize the convective oscillation. To mitigate the incompressible constraint, the selective reduced integration (SRI) and the recently proposed selective node–based smoothed FEM (SNS‐FEM) are used for the 4‐node quadrilateral element (CBP‐Q4SRI) and the 3‐node triangular element (CBP‐T3SNS), respectively. Meanwhile, the reduced integration (RI) for Q4 element (CBP‐Q4RI) and NS‐FEM for T3 element (CBP‐T3NS) with CBP scheme are also investigated. The quasi‐implicit CBP scheme is applied to allow a large time step for sufficient large penalty parameters. Due to the absences of pressure degree of freedoms, the quasi‐implicit CBP‐FEM has higher efficiency than quasi‐implicit CBS‐FEM. In this paper, the CBP‐Q4SRI has been verified and validated with high accuracy, stability, and fast convergence. Unexpectedly, CBP‐Q4RI is of no instability, high accuracy, and even slightly faster convergence than CBP‐Q4SRI. For unstructured T3 elements, CBP‐T3SNS also shows high accuracy and good convergence but with pressure oscillation using a large penalty parameter; CBP‐T3NS produces oscillated wrong velocity and pressure results. In addition, the applicable ranges of penalty parameter for different proposed methods have been investigated.     

Thermo‐elasto‐plastic finite element analysis of quasi‐state processes in Eulerian reference frames

An Eulerian finite element formulation for quasi‐state one way coupled thermo‐elasto‐plastic systems is presented. The formulation is suitable for modeling material processes such as welding and laser surfacing. In an Eulerian frame, the solution field of a quasi‐state process becomes steady state for the heat transfer problem and static for the stress problem. A mixed small deformation displacement elasto‐plastic formulation is proposed. The formulation accounts for temperature dependent material properties and exhibits a robust convergence. Streamline upwind Petrov–Galerkin (SUPG) is used to remove spurious oscillations. Smoothing functions are introduced to relax the non‐differentiable evolution equations and allow for the use of gradient (stiffness) solution scheme via the Newton–Raphson method. A 3‐dimensional simulation of a laser surfacing process is presented to exemplify the formulation. Results from the Eulerian formulation are in good agreement with results from the conventional Lagrangian formulation. However, the Eulerian formulation is approximately 15 times faster than the Lagrangian. Copyright © 2001 John Wiley & Sons, Ltd.     

Improved Epitaxy of AlN Film for Deep‐Ultraviolet Light‐Emitting Diodes Enabled by Graphene

The growth of single‐crystal III‐nitride films with a low stress and dislocation density is crucial for the semiconductor industry. In particular, AlN‐derived deep‐ultraviolet light‐emitting diodes (DUV‐LEDs) have important applications in microelectronic technologies and environmental sciences but are still limited by large lattice and thermal mismatches between the epilayer and substrate. Here, the quasi‐van der Waals epitaxial (QvdWE) growth of high‐quality AlN films on graphene/sapphire substrates is reported and their application in high‐performance DUV‐LEDs is demonstrated. Guided by density functional theory calculations, it is found that pyrrolic nitrogen in graphene introduced by a plasma treatment greatly facilitates the AlN nucleation and enables fast growth of a mirror‐smooth single‐crystal film in a very short time of ≈0.5 h (≈50% decrease compared with the conventional process), thus leading to a largely reduced cost. Additionally, graphene effectively releases the biaxial stress (0.11 GPa) and reduces the dislocation density in the epilayer. The as‐fabricated DUV‐LED shows a low turn‐on voltage, good reliability, and high output power. This study may provide a revolutionary technology for the epitaxial growth of AlN films and provide opportunities for scalable applications of graphene films.     

Stress intensity factors for cracks emanating from a circular hole in an infinite quasi‐orthotropic plane

An infinite quasi‐orthotropic plane with a cracked circular hole under tensile loading at infinity is studied analytically. To this end, complex variable theory of Muskhelishvili is used. In addition, to obtain analytical functions, a new conformal mapping is proposed and expanded to series expressions. Stress intensity factors (SIFs) for two unequal cracks emanating from a circular hole are obtained. To validate the analytical SIFs in a quasi‐orthotropic plane, the results are compared with FEM and the results of isotropic plane. The SIFs for small cracks in a quasi‐orthotropic and an isotropic plane are different, because of difference between stress concentrations in points which cracks emanate from the hole. However, the results of quasi‐orthotropic plane converge to isotropic plane for the large cracks. Therefore, the SIFs of the large cracks in a quasi‐orthotropic plane can be replaced by the results of the center crack with equivalent length in an isotropic plane.     

Finite element analysis of quasi‐static indentation of woven fabric textile composites using different nose shape indenters

The finite element FE analysis of quasi‐static indentation event of various nose shape rigid indenters into woven fabric composite with carbon fiber as reinforcement has been performed and discussed in detail. It was found that indenter nose shape has large influence in terms of absorbed energy, indentation at failure and damage area. The FE software, ABAQUS^® was employed to simulate quasi‐static response of woven composite unit cell. Exhaustive parametric studies have been conducted with an aim to analyze the effect of change in indenter geometry on the indentation response of the woven composite unit cell. The developed FE model for the purpose of validation was compared with available experimental results and was found to be in reasonably good agreement. The failure morphologies, damage shape and damage size were evaluated, compared and deeply discussed for different nose shape indenters. Largest damaged areas were observed for flat and truncated indenters while the smallest for the conical one.     

Multiscale analysis method for thermo‐mechanical performance of periodic porous materials with interior surface radiation

This study develops a novel multiscale analysis method to predict thermo‐mechanical performance of periodic porous materials with interior surface radiation. In these materials, thermal radiation effect at microscale has an important impact on the macroscopic temperature and stress field, which is our particular interest in this paper. Firstly, the multiscale asymptotic expansions for computing the dynamic thermo‐mechanical coupling problem, which considers the mutual interaction between temperature and displacement field, are given successively. Then, the corresponding numerical algorithm based on the finite element‐difference method is brought forward in details. Finally, some numerical results are presented to verify the validity and relevancy of the proposed method by comparing it with a direct finite element analysis with detailed numerical models. The comparison shows that the new method is effective and valid for predicting the thermo‐mechanical performance and can capture the microstructure behavior of periodic porous materials exactly.s Copyright © 2015 John Wiley & Sons, Ltd.