首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 0 毫秒
1.
    
This paper presents a further improvement of the standard time‐domain BEM formulation for 2‐D elastodynamics. Linear‐time interpolation is assumed for both boundary displacements and tractions.As this assumption implies time continuity for the variables, a procedure to consider tractions time discontinuities must be worked out. The initial step of this procedure consists of adding to the basic BEM integral equation the integral equation for velocities: after discretization is accomplished, this is done only for nodes with prescribed values of displacements. Additional equations are then incorporated to the final system of algebraic equations, providing the means for the determination of the extra unknowns (represented by the tractions at the begining of each time interval). The integral equations are presented by employing the concept of finite part of integrals, and the kernels are integrated analitically on time: complete time‐integrated expressions are included in an appendix. Three numerical examples are presented in order to assess the accuracy of the proposed formulation. Copyright © 2000 John Wiley & Sons, Ltd.  相似文献   

2.
    
A simple, elegant approach is proposed to correct the error introduced by the truncation of the infinite boundary in the BEM modelling of two‐dimensional wave propagation problems in elastic half‐spaces. The proposed method exploits the knowledge of the far‐field asymptotic behaviour of the solution to adequately correct the BEM displacement system matrix for the truncated problem to account for the contribution of the omitted part of the boundary. The reciprocal theorem of elastodynamics is used for a convenient computation of this contribution involving the same boundary integrals that form the original BEM system. The method is formulated for a two‐dimensional homogeneous, isotropic, linearly elastic half‐space and its implementation in a frequency domain boundary element scheme is discussed in some detail. The formulation is then validated for a free Rayleigh pulse travelling on a half‐space and successfully tested for a benchmark problem with a known approximation to the analytical solution. Copyright © 2004 John Wiley & Sons, Ltd.  相似文献   

3.
A new direct BE formulation is proposed for the solution of elastodynamic problems. An analytical regularization procedure is devised via an integration‐by‐parts technique without introducing any hypothesis about either the discretization of the boundary geometry or the space–time interpolation of the elastic fields involved. The only requirement is for continuity in the displacement field. The regularization of the integral equations for static elasticity using the same approach, already available in the literature, is presented as a particular case of the general procedure introduced herein. The numerical implementation of this technique is discussed and two‐dimensional examples are presented. Copyright © 1999 John Wiley & Sons, Ltd.  相似文献   

4.
A new material damping model which is convenient for use in the time‐domain boundary element method (TDBEM) is presented and implemented in a proposed procedure. Since only fundamental solutions for linear elastic material are employed, the procedure has high efficiency and is easy to be integrated into current TDBEM codes. Analytical and numerical results for benchmark problems demonstrate that the accuracy of the proposed method is high. Copyright © 2001 John Wiley & Sons, Ltd.  相似文献   

5.
    
Evidence of numerical instabilities in two‐dimensional time domain direct boundary element methods is presented. The effects of numerical versus analytical integration of spatial integrals on stability are shown, and two new time‐stepping algorithms are introduced and compared to existing formulations. The so‐called new ‘direct half‐step’ scheme and the ‘epsilon’ scheme are shown to improve the numerical stability of direct boundary element methods. Copyright © 1999 John Wiley & Sons, Ltd.  相似文献   

6.
    
The boundary integral equations in 3‐d elastodynamics contain convolution integrals with respect to the time. They can be performed analytically or with the convolution quadrature method. The latter time‐stepping procedure's benefit is the usage of the Laplace‐transformed fundamental solution. Therefore, it is possible to apply this method also to problems where analytical time‐dependent fundamental solutions might not be known. To obtain a symmetric formulation, the second boundary integral equation has to be used which, unfortunately, requires special care in the numerical implementation since it involves hypersingular kernel functions. Therefore, a regularization for closed surfaces of the Laplace‐transformed elastodynamic kernel functions is presented which transforms the bilinear form of the hypersingular integral operator to a weakly singular one. Supplementarily, a weakly singular formulation of the Laplace‐transformed elastodynamic double layer potential is presented. This results in a time domain boundary element formulation involving at least only weakly singular integral kernels. Finally, numerical studies validate this approach with respect to different spatial and time discretizations. Further, a comparison with the wider used collocation method is presented. It is shown numerically that the presented formulation exhibits a good convergence rate and a more stable behavior compared with collocation methods. Copyright © 2008 John Wiley & Sons, Ltd.  相似文献   

7.
    
In this work a fast solver for large‐scale three‐dimensional elastodynamic crack problems is presented, implemented, and tested. The dual boundary element method in the Laplace transform domain is used for the accurate dynamic analysis of cracked bodies. The fast solution procedure is based on the use of hierarchical matrices for the representation of the collocation matrix for each computed value of the Laplace parameter. An ACA (adaptive cross approximation) algorithm is used for the population of the low rank blocks and its performance at varying Laplace parameters is investigated. A preconditioned GMRES is used for the solution of the resulting algebraic system of equations. The preconditioners are built exploiting the hierarchical arithmetic and taking full advantage of the hierarchical format. An original strategy, based on the computation of some local preconditioners only, is presented and tested to further speed up the overall analysis. The reported numerical results demonstrate the effectiveness of the technique for both uncracked and cracked solids and show significant reductions in terms of both memory storage and computational time. Copyright © 2010 John Wiley & Sons, Ltd.  相似文献   

8.
    
A domain‐decomposition algorithm has been developed to handle two‐phase flows with large deformation, breaking and fragmentation of the interface. The strategy couples a boundary element method with a Navier–Stokes solver combined with a level‐set technique for the tracking of the interface. The former is used in the fluid region where the interface can be modelled as a smooth surface. In the rest of the domain the field solver is applied. This results in an efficient and accurate method. In this paper, the features of the used strategy are described and the challenges connected with the coupling are deeply discussed. The numerical investigation highlighted the importance of a proper rational study when CFD methods are considered. In the present case, a crucial aspect is represented by the domain‐composition step, that is when the information from one solver to the other have to be properly reconstructed and made consistent with the receiver sub‐domain. Copyright © 2006 John Wiley & Sons, Ltd.  相似文献   

9.
    
In finite element formulations for poroelastic continua a representation of Biot's theory using the unknowns solid displacement and pore pressure is preferred. Such a formulation is possible either for quasi‐static problems or for dynamic problems if the inertia effects of the fluid are neglected. Contrary to these formulations a boundary element method (BEM) for the general case of Biot's theory in time domain has been published (Wave Propagation in Viscoelastic and Poroelastic Continua: A Boundary Element Approach. Lecture Notes in Applied Mechanics. Springer: Berlin, Heidelberg, New York, 2001.). If the advantages of both methods are required it is common practice to couple both methods. However, for such a coupled FE/BE procedure a BEM for the simplified dynamic Biot theory as used in FEM must be developed. Therefore, here, the fundamental solutions as well as a BE time stepping procedure is presented for the simplified dynamic theory where the inertia effects of the fluid are neglected. Further, a semi‐analytical one‐dimensional solution is presented to check the proposed BE formulation. Finally, wave propagation problems are studied using either the complete Biot theory as well as the simplified theory. These examples show that no significant differences occur for the selected material. Copyright © 2005 John Wiley & Sons, Ltd.  相似文献   

10.
研究了基于加速度时域信息构造残差量的结构损伤识别方法。针对该方法实施过程中可能产生的识别方程组病态问题,以及随之产生的求解结果对测量噪声敏感等问题,采用截断奇异值分解技术(TSVD)求解识别方程组。以IASC?ASCE的基准结构为算例,采用本文的方法,识别出两种情况的损伤位置,损伤程度的识别误差分别为3.8%、6.0%,验证了所研究的基于加速度时域信息的结构损伤识别方法是有效的。  相似文献   

11.
    
A D‐BEM approach, based on time‐weighting residuals, is developed for the solution of two‐dimensional scalar wave propagation problems. The modified basic equation of the D‐BEM formulation is generated by weighting, with respect to time, the basic D‐BEM equation, under the assumption of linear and cubic time variation for the potential and for the flux. A constant time‐weighting function is adopted. The time integration reduces the order of the time‐derivative that appears in the domain integral; as a consequence, the initial conditions are directly taken into account. An assessment of the potentialities of the proposed formulation is provided by the examples included at the end of the work. Copyright © 2009 John Wiley & Sons, Ltd.  相似文献   

12.
    
A precorrected fast Fourier transform (pFFT) accelerated boundary element method (BEM) for large‐scale transient elastodynamic analysis is developed and described in this paper. The frequency‐domain approach is used. To overcome the ‘wrap‐around’ problem associated with the discrete Fourier transform, the exponential window method (EWM) is employed and incorporated in the frequency‐domain BEM. An improved implementation scheme of the pFFT method based on polynomial interpolation technique is developed and applied to accelerate the elastodynamic BEM. This new scheme reduces the memory required to save the convolution matrix by a factor of 8. To further improve the efficiency of the code, a newly developed linear system solver based on the induced dimension reduction method is employed. Its performance is investigated and compared with that of the well‐known GMRES. The accuracy and computational efficiency of the method are evaluated and demonstrated by three examples: a classical benchmark, a plate subject to an impact loading and a porous cube with nearly half million DOFs subject to a step traction loading. Both analytical and experimental results are employed to validate the method. It has been found that the EWM can effectively resolve the wrap‐around problem and accurate time responses for an arbitrarily chosen time period can be obtained. Copyright © 2011 John Wiley & Sons, Ltd.  相似文献   

13.
    
An iterative coupling of finite element and boundary element methods for the time domain modelling of coupled fluid–solid systems is presented. While finite elements are used to model the solid, the adjacent fluid is represented by boundary elements. In order to perform the coupling of the two numerical methods, a successive renewal of the variables on the interface between the two subdomains is performed through an iterative procedure until the final convergence is achieved. In the case of local non‐linearities within the finite element subdomain, it is straightforward to perform the iterative coupling together with the iterations needed to solve the non‐linear system. In particular a more efficient and a more stable performance of the new coupling procedure is achieved by a special formulation that allows to use different time steps in each subdomain. Copyright © 2005 John Wiley & Sons, Ltd.  相似文献   

14.
The present paper describes a procedure that improves efficiency, stability and reduces artificial energy dissipation of the standard time-domain direct boundary element method (BEM) for acoustics and elastodynamics. Basically, the developed procedure modifies the boundary element convolution-related vector, being very easy to implement into existing codes. A stabilization parameter is introduced into the recent-in-time convolution operations and the operations related to the distant-in-time convolution contributions are approximated by matrix interpolations. As it is shown in the numerical examples presented at the end of the text, the proposed formulation substantially reduces the BEM computational cost, as well as its numerical instabilities.  相似文献   

15.
    
The time‐parallel framework for constructing parallel implicit time‐integration algorithms (PITA) is revisited in the specific context of linear structural dynamics and near‐real‐time computing. The concepts of decomposing the time‐domain in time‐slices whose boundaries define a coarse time‐grid, generating iteratively seed values of the solution on this coarse time‐grid, and using them to time‐advance the solution in each time‐slice with embarrassingly parallel time‐integrations are maintained. However, the Newton‐based corrections of the seed values, which so far have been computed in PITA and related approaches on the coarse time‐grid, are eliminated to avoid artificial resonance and numerical instability. Instead, the jumps of the solution on the coarse time‐grid are addressed by a projector which makes their propagation on the fine time‐grid computationally feasible while avoiding artificial resonance and numerical instability. The new PITA framework is demonstrated for a complex structural dynamics problem from the aircraft industry. Its potential for near‐real‐time computing is also highlighted with the solution of a relatively small‐scale problem on a Linux cluster system. Copyright © 2006 John Wiley & Sons, Ltd.  相似文献   

16.
In this paper a new time discontinuous Galerkin (TDG) formulation for non‐linear elastodynamics is presented. The new formulation embeds an energy correction which ensures truly energy decaying, thus allowing to achieve unconditional stability that, as shown in the paper, is not guaranteed by the classical TDG formulation. The resulting method is simple and easily implementable into existing finite element codes. Moreover, it inherits the desirable higher‐order accuracy and high‐frequency dissipation properties of the classical formulation. Numerical results illustrate the very good performance of the proposed formulation. Copyright © 2010 John Wiley & Sons, Ltd.  相似文献   

17.
    
A boundary element formulation is presented in this work for transient heat conduction analysis of three‐dimensional (3D) fiber‐reinforced materials. The cylindrical‐shaped fibers in a 3D matrix are idealized by a system of curvilinear line elements with a prescribed diameter. The variations in the temperature and flux fields in the circumferential direction are represented in terms of a trigonometric shape function together with a linear or quadratic variation in the longitudinal direction. This approach significantly reduces the modeling effort and the computing cost. The storage requirement for the convolution integrals is eliminated by adopting an accurate integration‐based convolution algorithm for the surface of the hole and the fibers as well as a fast convolution algorithm for the outer boundary recently developed by the present authors. Numerical examples are presented to demonstrate the accuracy and applicability of the proposed method of analysis of fiber‐reinforced materials. Copyright © 2007 John Wiley & Sons, Ltd.  相似文献   

18.
    
Two‐dimensional photonic crystal structures are analyzed by a recently developed hybrid technique combining the finite‐element time‐domain (FETD) method and the finite‐difference time‐domain (FDTD) method. This hybrid FETD/FDTD method uses the discontinuous Galerkin method as framework for domain decomposition. To the best of our knowledge, this is the first hybrid FETD/FDTD method that allows non‐conformal meshes between different FETD and FDTD subdomains. It is also highly parallelizable. These properties are very suitable for the computation of periodic structures with curved surfaces. Numerical examples for the computation of the scattering parameters of two‐dimensional photonic bandgap structures are presented as applications of the hybrid FETD/FDTD method. Numerical results demonstrate the efficiency and accuracy of the proposed hybrid method. Copyright © 2012 John Wiley & Sons, Ltd.  相似文献   

19.
A hypersingular time‐domain boundary element method (BEM) for transient elastodynamic crack analysis in two‐dimensional (2D), homogeneous, anisotropic, and linear elastic solids is presented in this paper. Stationary cracks in both infinite and finite anisotropic solids under impact loading are investigated. On the external boundary of the cracked solid the classical displacement boundary integral equations (BIEs) are used, while the hypersingular traction BIEs are applied to the crack‐faces. The temporal discretization is performed by a collocation method, while a Galerkin method is implemented for the spatial discretization. Both temporal and spatial integrations are carried out analytically. Special analytical techniques are developed to directly compute strongly singular and hypersingular integrals. Only the line integrals over an unit circle arising in the elastodynamic fundamental solutions need to be computed numerically by standard Gaussian quadrature. An explicit time‐stepping scheme is obtained to compute the unknown boundary data including the crack‐opening‐displacements (CODs). Special crack‐tip elements are adopted to ensure a direct and an accurate computation of the elastodynamic stress intensity factors from the CODs. Several numerical examples are given to show the accuracy and the efficiency of the present hypersingular time‐domain BEM. Copyright © 2008 John Wiley & Sons, Ltd.  相似文献   

20.
Two different boundary element methods (BEM) for crack analysis in two dimensional (2-D) antiplane, homogeneous, isotropic and linear elastic solids by considering frictional contact of the crack edges are presented. Hypersingular boundary integral equations (BIE) in time-domain (TD) and frequency domain (FD), with corresponding elastodynamic fundamental solutions are applied for this purpose. For evaluation of the hypersingular integrals involved in BIEs a special regularization process that converts the hypersingular integrals to regular integrals is applied. Simple regular formulas for their calculation are presented. For the problems solution while considering frictional contact of the crack edges a special iterative algorithm of Udzava's type is elaborated and used. Numerical results for crack opening, frictional contact forces and dynamic stress intensity factors (SIFs) are presented and discussed for a finite III-mode crack in an infinite domain subjected to a harmonic crack-face loading and considering crack edges frictional contact interaction using the TD and FD approaches.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号