Galerkin based generalized ANOVA for the solution of stochastic steady state diffusion problems

This work presents a novel approach, referred here as Galerkin based generalized analysis of variance decomposition (GG-ANOVA), for the solution of stochastic steady state diffusion problems. The proposed approach utilizes generalized ANOVA (G-ANOVA) expansion to represent the unknown stochastic response and Galerkin projection to decompose the stochastic differential equation into a set of coupled differential equations. The coupled set of partial differential equations obtained are solved using finite difference method and homotopy algorithm. Implementation of the proposed approach for solving stochastic steady state diffusion problems has been illustrated with three numerical examples. For all the examples, results obtained are in excellent agreement with the benchmark solutions. Additionally, for the second and third problems, results obtained have also been compared with those obtained using polynomial chaos expansion (PCE) and conventional G-ANOVA. It is observed that the proposed approach yields highly accurate result outperforming both PCE and G-ANOVA. Moreover, computational time required using GG-ANOVA is in close proximity of G-ANOVA and less as compared to PCE.     

Vibration modeling of large repetitive sandwich structures with viscoelastic core

A double scale asymptotic method (DSAM) is proposed for vibration modeling of large repetitive sandwich structures with a viscoelastic core. The method decomposes the initial nonlinear vibration problem into two small linear ones. The first one is defined on few basic cells while the second is a differential global amplitude equation with complex coefficients. Their numerical computations permit determination of the damping properties as well as pass and stop bands avoiding the direct computation on the whole structure. Viscoelastic frequency dependent core with fractional and anelastic displacement field models are considered. The resulting nonclassical problems are solved by asymptotic numerical method coupled with automatic differentiation. Based on the presented method, a large reduction of the needed computational time and memory is obtained. The accuracy and efficiency of the proposed method are validated with comparisons to the direct simulations by discretization of the whole structure using asymptotic numerical method coupled with automatic differentiation.     

Topology optimization of acoustic–structure interaction problems using a mixed finite element formulation

The paper presents a gradient‐based topology optimization formulation that allows to solve acoustic–structure (vibro‐acoustic) interaction problems without explicit boundary interface representation. In acoustic–structure interaction problems, the pressure and displacement fields are governed by Helmholtz equation and the elasticity equation, respectively. Normally, the two separate fields are coupled by surface‐coupling integrals, however, such a formulation does not allow for free material re‐distribution in connection with topology optimization schemes since the boundaries are not explicitly given during the optimization process. In this paper we circumvent the explicit boundary representation by using a mixed finite element formulation with displacements and pressure as primary variables (a u /p‐formulation). The Helmholtz equation is obtained as a special case of the mixed formulation for the elastic shear modulus equating to zero. Hence, by spatial variation of the mass density, shear and bulk moduli we are able to solve the coupled problem by the mixed formulation. Using this modelling approach, the topology optimization procedure is simply implemented as a standard density approach. Several two‐dimensional acoustic–structure problems are optimized in order to verify the proposed method. Copyright © 2006 John Wiley & Sons, Ltd.     

基于能量有限元法的功能梯度梁振动分析

功能梯度材料(Functionally Graded Material,FGM)由于其优良的结构性能和重要的应用价值,近些年来得到了广泛的研究和关注。采用能量有限元法对功能梯度梁和耦合梁的弯曲振动特性进行研究,推导了功能梯度材料梁的能量密度控制方程、能量有限元矩阵方程以及耦合梁的能量有限元方程,从而得到梁中的能量密度和能量流。以一简支功能梯度梁为例,分别采用该方法和传统有限元法计算了梁弯曲振动时的能量密度,通过对比验证了能量有限元法求解的准确性。在此基础上进一步对耦合功能梯度梁结构的能量密度和能量流进行了求解,得到其能量分布特征。该研究为基于能量有限元法分析复杂功能梯度材料结构的振动特性提供了理论基础。     

Finite element model for the coupled radiative transfer equation and diffusion approximation

In this paper a coupled radiative transfer equation and diffusion approximation model for light propagation in tissues is proposed. The light propagation is modelled with the radiative transfer equation in sub‐domains in which the assumptions of the diffusion approximation are not valid. The diffusion approximation is used elsewhere in the domain. The two equations are coupled through their boundary conditions and they are solved simultaneously using the finite element method. The method is tested with simulations. The results of the proposed approach are compared with finite element solutions of the radiative transfer equation, the diffusion approximation and a previously proposed hybrid model. The results show that the method improves the accuracy of the forward solution for diffuse optical tomography compared to the conventional diffusion model. Copyright © 2005 John Wiley & Sons, Ltd.     

Consistent arbitrary Lagrangian Eulerian formulation for large deformation thermo-mechanical analysis

Arbitrary Lagrangian Eulerian (ALE) method is widely used for simulation of large deformation problems, such as metal forming. However, in many such applications, modeling of the heat generation and transfer in conjunction with the stress analysis is necessary. In this work, a fully coupled dynamic ALE formulation is developed. The ALE form of energy balance equation is derived, and is coupled with the dynamic, rate dependent ALE stress analysis. The proposed formulation is used for simulation of a few thermo-mechanical problems. The effectiveness and efficiency of the ALE method is verified by comparing the results of this simulation with available experimental and numerical results.     

基于等效线性化方法的非线性振动能量采集器功率分析

输出功率分析是振动能量采集器结构设计、参数选择的重要依据.有鉴于传统功率分析方法比较复杂,以一类电磁式振动能量采集器为对象,提出了一种新的非线性振动能量采集器功率分析方法,其核心是将非线性振动方程进行等效线性化处理,并借助传递函数开展功率优化分析.首先,考虑非线性磁力以及基尔霍夫电流定律,建立了通用化的采集器1.5自由...     

Flexural-Torsional Buckling and Vibration Analysis of Composite Beams

In this paper the general flexural-torsional buckling and vibration problems of composite Euler-Bernoulli beams of arbitrarily shaped cross section are solved using a boundary element method. The general character of the proposed method is verified from the formulation of all basic equations with respect to an arbitrary coordinate system, which is not restricted to the principal one. The composite beam consists of materials in contact each of which can surround a finite number of inclusions. It is subjected to a compressive centrally applied load together with arbitrarily transverse and/or torsional distributed or concentrated loading, while its edges are restrained by the most general linear boundary conditions. The resulting problems are (i) the flexural-torsional buckling problem, which is described by three coupled ordinary differential equations and (ii) the flexural-torsional vibration problem, which is described by three coupled partial differential equations. Both problems are solved employing a boundary integral equation approach. Besides the effectiveness and accuracy of the developed method, a significant advantage is that the method can treat composite beams of both thin and thick walled cross sections taking into account the warping along the thickness of the walls. The proposed method overcomes the shortcoming of possible thin tube theory (TTT) solution, which its utilization has been proven to be prohibitive even in thin walled homogeneous sections. Example problems of composite beams are analysed, subjected to compressive or vibratory loading, to illustrate the method and demonstrate its efficiency and wherever possible its accuracy. Moreover, useful conclusions are drawn from the buckling and dynamic response of the beam.     

Bi-directional evolutionary structural optimization for design-dependent fluid pressure loading problems

This article presents an evolutionary topology optimization method for compliance minimization of structures under design-dependent pressure loads. In traditional density based topology optimization methods, intermediate values of densities for the solid elements arise along the iterations. Extra boundary parametrization schemes are demanded when these methods are applied to pressure loading problems. An alternative methodology is suggested in this article for handling this type of load. With an extended bi-directional evolutionary structural optimization method associated with a partially coupled fluid–structure formulation, pressure loads are modelled with hydrostatic fluid finite elements. Due to the discrete nature of the method, the problem is solved without any need of pressure load surfaces parametrization. Furthermore, the introduction of a separate fluid domain allows the algorithm to model non-constant pressure fields with Laplace's equation. Three benchmark examples are explored in order to show the achievements of the proposed method.     

A meshless method for generalized linear or nonlinear Poisson-type problems

This paper presents a meshless method, based on coupling virtual boundary collocation method (VBCM) with the radial basis functions (RBF) and the analog equation method (AEM), to analyze generalized linear or nonlinear Poisson-type problems. In this method, the AEM is used to construct equivalent equations to the original differential equation so that a simpler fundamental solution of the Laplacian operator, instead of other complicated ones which are needed in conventional BEM, can be employed. While global RBF is used to approximate fictitious body force which appears when the analog equation method is introduced, and VBCM are utilized to solve homogeneous solution based on the superposition principle. As a result, a new meshless method is developed for solving nonlinear Poisson-type problems. Finally, some numerical experiments are implemented to verify the efficiency of the proposed method and numerical results are in good agreement with the analytical ones. It appears that the proposed meshless method is very effective for nonlinear Poisson-type problems.     

The time-marching method of fundamental solutions for wave equations

In this paper, a meshless numerical algorithm is developed for the solution of multi-dimensional wave equations with complicated domains. The proposed numerical method, which is truly meshless and quadrature-free, is based on the Houbolt finite difference (FD) scheme, the method of the particular solutions (MPS) and the method of fundamental solutions (MFS). The wave equation is transformed into a Poisson-type equation with a time-dependent loading after the time domain is discretized by the Houbolt FD scheme. The Houbolt method is used to avoid the difficult problem of dealing with time evolution and the initial conditions to form the linear algebraic system. The MPS and MFS are then coupled to analyze the governing Poisson equation at each time step. In this paper we consider six numerical examples, namely, the problem of two-dimensional membrane vibrations, the wave propagation in a two-dimensional irregular domain, the wave propagation in an L-shaped geometry and wave vibration problems in the three-dimensional irregular domain, etc. Numerical validations of the robustness and the accuracy of the proposed method have proven that the meshless numerical model is a highly accurate and efficient tool for solving multi-dimensional wave equations with irregular geometries and even with non-smooth boundaries.     

A coupled thermo-mechanical bond-based peridynamics for simulating thermal cracking in rocks

A coupled thermo-mechanical bond-based peridynamical (TM-BB-PD) method is developed to simulate thermal cracking processes in rocks. The coupled thermo-mechanical model consists of two parts. In the first part, temperature distribution of the system is modeled based on the heat conduction equation. In the second part, the mechanical deformation caused by temperature change is calculated to investigate thermal fracture problems. The multi-rate explicit time integration scheme is proposed to overcome the multi-scale time problem in coupled thermo-mechanical systems. Two benchmark examples, i.e., steady-state heat conduction and transient heat conduction with deformation problem, are performed to illustrate the correctness and accuracy of the proposed coupled numerical method in dealing with thermo-mechanical problems. Moreover, two kinds of numerical convergence for peridynamics, i.e., m- and \(\delta \)-convergences, are tested. The thermal cracking behaviors in rocks are also investigated using the proposed coupled numerical method. The present numerical results are in good agreement with the previous numerical and experimental data. Effects of PD material point distributions and nonlocal ratios on thermal cracking patterns are also studied. It can be found from the numerical results that thermal crack growth paths do not increases with changes of PD material point spacing when the nonlocal ratio is larger than 4. The present numerical results also indicate that thermal crack growth paths are slightly affected by the arrangements of PD material points. Moreover, influences of thermal expansion coefficients and inhomogeneous properties on thermal cracking patterns are investigated, and the corresponding thermal fracture mechanism is analyzed in simulations. Finally, a LdB granite specimen with a borehole in the heated experiment is taken as an application example to examine applicability and usefulness of the proposed numerical method. Numerical results are in good agreement with the previous experimental and numerical results. Meanwhile, it can be found from the numerical results that the coupled TM-BB-PD has the capacity to capture phenomena of temperature jumps across cracks, which cannot be captured in the previous numerical simulations.     

A perturbational approach to magneto-thermal problems of a deformed sphere levitated in a magnetic field

This paper presents a perturbation method for the solution of the electromagnetic and thermal problems of a deformed sphere levitated in an alternating magnetic field. The analytical solutions of the electromagnetic field distribution, the Joule heat generation, the magnetic lifting force and the temperature field are obtained based on a linear perturbation theory. The Maxwell equations are first simplified in terms of the vector potential and then solved by the method of separation of variables. The time-averaged Joule-heat source is calculated and coupled to the Fourier heat-conduction equation. The coupled equation is solved for temperature distributions within the deformed sphere by a combined approach of series expansion and variation of parameters. Both asymptotic and numerical analyses are provided. The total power absorption and temperature field for both single and multiple coils are also discussed.     

A robust formulation for solving fluid-structure interaction problems

 In this article, a coupled finite element and boundary element approach for the acoustic radiation and scattering from submerged elastic bodies of arbitrary shape is presented. An alternative to the direct boundary element method for acoustics is proposed. By taking an auxiliary source surface inside the radiating boundary and following the usual discretization and integration procedures employed in the boundary element method, both the singularities of the integrands and the nonuniqueness problems do not arise. In addition, the difficulty of slope discontinuity also can be overcome. This procedure is formulated in a similar fashion of wave superposition method, except that the direct boundary integral equations are adopted. The proposed formulation employ the surface Helmholtz integral equation and its normal gradient like that adopted in the Burton–Miller approach, but do not employ any coupling constant. Typical examples are presented that demonstrate the efficiency of the proposed technique. Received 9 April 2000     

Solution of viscoelastic scattering problems in linear acoustics using hp Boundary/Finite Element Method

The interaction of acoustic waves with submerged structures remains one of the most difficult and challenging problems in underwater acoustics. Many techniques such as coupled Boundary Element (BE)/Finite Element (FE) or coupled Infinite Element (IE)/Finite Element approximations have evolved. In the present work, we focus on the steady‐state formulation only, and study a general coupled hp‐adaptive BE/FE method. A particular emphasis is placed on an a posteriori error estimation for the viscoelastic scattering problems. The highlights of the proposed methodology are as follows: (1) The exterior Helmholtz equation and the Sommerfeld radiation condition are replaced with an equivalent Burton–Miller (BM) boundary integral equation on the surface of the scatterer. (2) The BM equation is coupled to the steady‐state form of viscoelasticity equations modelling the behaviour of the structure. (3) The viscoelasticity equations are approximated using hp‐adaptive FE isoparametric discretizations with order of approximation p⩾5 in order to avoid the ‘locking’ phenomenon. (4) A compatible hp superparametric discretization is used to approximate the BM integral equation. (5) Both the FE and BE approximations are based on a weak form of the equations, and the Galerkin method, allowing for a full convergence analysis. (6) An a posteriori error estimate for the coupled problem of a residual type is derived, allowing for estimating the error in pressure on the wet surface of the scatterer. (7) An adaptive scheme, an extension of the Texas Three Step Adaptive Strategy is used to manipulate the mesh size h and the order of approximation p so as to approximately minimize the number of degrees of freedom required to produce a solution with a specified accuracy. The use of this hp‐scheme may exhibit exponential convergence rates. Several numerical experiments illustrate the methodology. These include detailed convergence studies for the problem of scattering of a plane acoustic wave on a viscoelastic sphere, and adaptive solutions of viscoelastic scattering problems for a series of MOCK0 models. Copyright © 1999 John Wiley & Sons, Ltd.     

Three-dimensional fracture analysis of CT specimens with a ductile damage model based on endochronic plasticity theory

For many fracture problems of practical engineering importance, the three-dimensional effects are significant and a three-dimensional analysis for the problems is thus required. In this paper, an endochronic theory coupled with anisotropic damage is first established, which is actually an elasto-plastic damage theory coupled with isotropic-nonlinear kinematic hardening. The ductile damage evolution equation is derived from the orthogonality rule with a new intrinsic time scale introduced especially for damage evolution. Then, a three-dimensional finite element program incorporating the endochronic damage model is formulated and emploved to analyze the widely used CT fracture specimen. Two failure criteria are proposed for the prediction of crack initiation direction and crack initiation load. From the analysis, significant three-dimensional effects are observed and the crack is estimated to initiate first at the middle of the crack front line. Experiments have been conducted to verify the proposed theory and the results are found to compare well with the theoretical values.     

A simple solution procedure to 3D-piezoelectric problems: Isotropic BEM coupled with a point collocation method

This paper presents a simple strategy allowing to adapt well established isotropic BEM approach for the solution of coupled problems with anisotropic material parameters. The method which is illustrated on the case of a piezoelectric material is based on the partition of the primary fields into complementary and particular parts. The complementary fields solve the isotropic form of the partial differential equation while particular fields are obtained by a point collocation of the strong form equation. Using the local radial point interpolation method, the effectiveness and accuracy of approach is demonstrated on some examples allowing a comparison with literature results.     

A pure boundary element method approach for solving hypersingular boundary integral equations

This paper is concerned with discretization and numerical solution of a regularized version of the hypersingular boundary integral equation (HBIE) for the two-dimensional Laplace equation. This HBIE contains the primary unknown, as well as its gradient, on the boundary of a body. Traditionally, this equation has been solved by combining the boundary element method (BEM) together with tangential differentiation of the interpolated primary variable on the boundary. The present paper avoids this tangential differentiation. Instead, a “pure” BEM method is proposed for solving this class of problems. Dirichlet, Neumann and mixed problems are addressed in this paper, and some numerical examples are included in it.     

Temperature-constrained topology optimization of thermo-mechanical coupled problems

Temperature-constrained topology optimization for thermo-mechanical coupled problems under a design-dependent temperature field considering the thermal expansion effect remains an open problem. A temperature-constrained topology optimization method is proposed for thermo-mechanical coupled problems. In this article, the temperature values at the heat sources are constrained. The numerical results reveal that the temperature constraints play an important role in topology optimization of thermo-mechanical coupled problems. The optimized structure obtained by the presented method not only has certain strength but also decreases the temperature significantly compared with structures obtained by other methods without considering temperature constraints. The proposed method is applied to the design of the cooling system of a battery package. Numerical examples verify the efficiency of the presented method.