共查询到20条相似文献,搜索用时 15 毫秒
1.
S. Hamzehei Javaran S. Shojaee 《International journal for numerical methods in engineering》2017,112(13):2067-2086
In this paper, new spherical Hankel shape functions are used to reformulate boundary element method for 2‐dimensional elastostatic and elastodynamic problems. To this end, the dual reciprocity boundary element method is reconsidered by using new spherical Hankel shape functions to approximate the state variables (displacements and tractions) of Navier's differential equation. Using enrichment of a class of radial basis functions (RBFs), called spherical Hankel RBFs hereafter, the interpolation functions of a Hankel boundary element framework has been derived. For this purpose, polynomial terms are added to the functional expansion that only uses spherical Hankel RBF in the approximation. In addition to polynomial function fields, the participation of spherical Bessel function fields has also increased robustness and efficiency in the interpolation. It is very interesting that there is no Runge phenomenon in equispaced Hankel macroelements, unlike equispaced classic Lagrange ones. Several numerical examples are provided to demonstrate the effectiveness, robustness and accuracy of the proposed Hankel shape functions and in comparison with the classic Lagrange ones, they show much more accurate and stable results. 相似文献
2.
Yijun Liu 《International journal for numerical methods in engineering》2006,65(6):863-881
A new fast multipole boundary element method (BEM) is presented in this paper for large‐scale analysis of two‐dimensional (2‐D) elastostatic problems based on the direct boundary integral equation (BIE) formulation. In this new formulation, the fundamental solution for 2‐D elasticity is written in a complex form using the two complex potential functions in 2‐D elasticity. In this way, the multipole and local expansions for 2‐D elasticity BIE are directly linked to those for 2‐D potential problems. Furthermore, their translations (moment to moment, moment to local, and local to local) turn out to be exactly the same as those in the 2‐D potential case. This formulation is thus very compact and more efficient than other fast multipole approaches for 2‐D elastostatic problems using Taylor series expansions of the fundamental solution in its original form. Several numerical examples are presented to study the accuracy and efficiency of the developed fast multipole BEM formulation and code. BEM models with more than one million equations have been solved successfully on a laptop computer. These results clearly demonstrate the potential of the developed fast multipole BEM for solving large‐scale 2‐D elastostatic problems. Copyright © 2005 John Wiley & Sons, Ltd. 相似文献
3.
This paper presents a new boundary element application for free vibration analysis of 2D elastic structures. The dual reciprocity
method is applied using four compact supported radial basis functions for approximating the domain inertia terms. The eigen-problem
of displacement is then solved considering the traction contribution by means of static condensation. The formulation is also
extended to consider additional internal nodes to improve accuracy. Three numerical problems are studied to demonstrate the
validity and accuracy of the developed formulation. The results are compared to those obtained from analytical and other numerical
solutions. A parametric study is set up to demonstrate the effect of the compact support radius on the final results and on
the sparsity of system matrices. 相似文献
4.
M.E. Biancolini A. Chiappa F. Giorgetti C. Groth U. Cella P. Salvini 《International journal for numerical methods in engineering》2018,115(12):1411-1429
A common issue in multiphysics analysis regards a reliable way for loose couplings, because the same object is modeled using different mesh refinements, each one suited for a proper field of physics. Output data originating from a simulation environment are transferred as input data to a different model to run a new analysis. It is strongly desirable that such information transfers in a conservative way in terms of general balance. This paper faces the problem of pressure mapping between widely dissimilar meshes. The proposed procedure yields two steps: pressure interpolation by means of radial basis functions and fuzzy subset correction. The first step is pointwise interpolation that exploits a series of basis functions. The second step applies to the outcome of the first one to reestablish load balance between the two models through the introduction of a smooth correction field. Practical tests from the aeronautical field allow validating the method. 相似文献
5.
6.
The Boundary Element Method is a very effective method for solving linear differential equations. To use it also in the consideration of non-linear problems some different procedures were developed, among them the dual reciprocity method and the particular integral method. Both procedures use interpolation conditions for the approximation with radial basis functions. In this paper a method is presented which avoids problems connected with interpolation by means of quasi-interpolation. It is possible to solve differential equations of the kind Δmu=p(u) approximately; the application to two non-linear problems of plate theory yield good results. Hints to a theoretical examination of the method including sufficient conditions for feasibility and convergence are given. © 1997 by John Wiley & Sons, Ltd. 相似文献
7.
Hideaki Kaneko Peter A. Padilla 《International journal for numerical methods in engineering》1999,45(4):491-495
In this note, we make a few comments concerning the paper of Hughes and Akin (Int. J. Numer. Meth. Engng., 15 , 733–751 (1980)). Our primary goal is to demonstrate that the rate of convergence of numerical solutions of the finite element method with singular basis functions depends upon the location of additional collocation points associated with the singular elements. Copyright © 1999 John Wiley & Sons, Ltd. 相似文献
8.
9.
Adam Ghoneim 《International journal for numerical methods in engineering》2016,107(10):813-852
A new technique for sharp‐interface modeling of dendritic solidification is proposed using a meshfree interface finite element method such that the liquid–solid interface is represented implicitly and allowed to arbitrarily intersect the finite elements. At the interface‐embedded elements, meshfree interface points without connectivity are imposed directly at the zero level set while meshfree interpolants are constructed using radial basis functions. This ensures both the partition of unity and the Kronecker delta properties are satisfied allowing for precise and easy imposition of Dirichlet boundary conditions at the interface. The constructed meshfree interpolants are also used for solving a variational level set equation based on the Ginzburg–Landau energy functional minimization such that reinitialization is completely eliminated and fast marching algorithms for interfacial velocity extension are not necessary resulting in an efficient algorithm with excellent volume conservation. The meshfree interface finite element method is used for modeling dendritic solidification in a pure melt where it is found suitable in handling the complex interfacial dynamics often encountered in dendritic growth. Mathematical formulation and implementation followed by numerical results and analysis will be presented and discussed. Copyright © 2016 John Wiley & Sons, Ltd. 相似文献
10.
A local RBF collocation method for band structure computations of 2D solid/fluid and fluid/solid phononic crystals
下载免费PDF全文

Hui Zheng Chuanzeng Zhang Yuesheng Wang Wen Chen Jan Sladek Vladimir Sladek 《International journal for numerical methods in engineering》2017,110(5):467-500
In this paper, an efficient local radial basis function collocation method (LRBFCM) is presented for computing the band structures of the two‐dimensional (2D) solid/fluid and fluid/solid phononic crystals. Both systems of solid scatterers embedded in a fluid matrix (solid/fluid phononic crystals) and fluid scatterers embedded in a solid matrix (fluid/solid phononic crystals) are investigated. The solid–fluid interactions are taken into account by properly formulating and treating the continuity/equilibrium conditions on the solid–fluid interfaces, which require an accurate computation of the normal derivatives of the displacements and the pressure on the fluid–solid interfaces and the unit‐cell boundaries. The developed LRBFCM for the mixed wave propagation problems in 2D solid/fluid and fluid/solid phononic crystals is validated by the corresponding results obtained by the finite element method (FEM). To the best knowledge of the authors, the present LRBFCM has yet not been applied to the band structure computations of 2D solid/fluid and fluid/solid phononic crystals. For different lattice forms, scatterers' shapes, acoustic impedance ratios, and material combinations (solid scatterers in fluid matrix or fluid scatterers in solid matrix), numerical results are presented and discussed to reveal the efficiency and the accuracy of the developed LRBFCM for calculating the band structures of 2D solid/fluid and fluid/solid phononic crystals. A comparison of the present numerical results with that of the FEM shows that the present LRBFCM is much more efficient than the FEM for the band structure computations of the considered 2D solid/fluid and fluid/solid phononic crystals. Copyright © 2016 John Wiley & Sons, Ltd. 相似文献
11.
M. L. Bittencourt 《International journal for numerical methods in engineering》2005,63(11):1530-1558
This paper presents nodal and modal shape functions for triangle and tetrahedron finite elements. The functions are constructed based on the fully tensorial expansions of one‐dimensional polynomials expressed in barycentric co‐ordinates. The nodal functions obtained from the application of the tensorial procedure are the standard h‐Lagrange shape functions presented in the literature. The modal shape functions use Jacobi polynomials and have a natural global C0 inter‐element continuity. An efficient Gauss–Jacobi numerical integration procedure is also presented to decrease the number of points for the consistent integration of the element matrices. An example illustrates the approximation properties of the modal functions. Copyright © 2005 John Wiley & Sons, Ltd. 相似文献
12.
Hui Wang Qing‐Hua Qin Detdexa Arounsavat 《International journal for numerical methods in engineering》2007,69(6):1262-1277
This investigation provides a hybrid Trefftz finite element approach for analysing minimal surface problems. The approach is based on combining Trefftz finite element formulation with radial basis functions (RBF) and the analogue equation method (AEM). In this method, use of the analogue equation approach avoids the difficulty of treating the non‐linear terms appearing in the soap bubble equation, making it possible to solve non‐linear problems with the Trefftz method. Global RBF is used to approximate the inhomogeneous term induced from non‐linear functions and other loading terms. Finally, some numerical experiments are implemented to verify the efficiency of this method. Copyright © 2006 John Wiley & Sons, Ltd. 相似文献
13.
Andrea A. Mammoli 《International journal for numerical methods in engineering》2002,55(9):1115-1128
The boundary integral equation that results from the application of the reciprocity theorem to non‐linear or non‐homogeneous differential equations generally contains a domain integral. While methods exist for the meshless evaluation of these integrals, mesh‐based domain integration is generally more accurate and can be performed more quickly with the application of fast multipole methods. However, polygonalization of complex multiply‐connected geometries can become a costly task, especially in three‐dimensional analyses. In this paper, a method that allows a mesh‐based integration in complex domains, while retaining a simple mesh structure, is described. Although the technique is intended for the numerical solution of more complex differential equations, such as the Navier–Stokes equations, for simplicity the method is applied to the solution of a Poisson equation, in domains of varying complexity. It is shown that the error introduced by the auxiliary domain subtraction method is comparable to the discretization error, while the complexity of the mesh is significantly reduced. The behaviour of the error in the boundary solution observed with the application of the new method is analogous to the behaviour observed with conventional cell‐based domain integration. Copyright © 2002 John Wiley & Sons, Ltd. 相似文献
14.
This paper is concerned with the development of a meshless local approach based on the finite collocation method for solving Cauchy problems of 2-D elliptic PDEs in annulus domains. In the proposed approach, besides the collocation of unknown solution, the governing equation is also enforced in the local domains. Moreover, to improve the accuracy, the method considers auxiliary points in local subdomains and imposes the governing PDE operator at these points, without changing the global system size. Localization property of the method reduces the ill-conditioning of the problem and makes it efficient for Cauchy problem. To show the efficiency of the method, four test problems containing Laplace, Poisson, Helmholtz and modified Helmholtz equations are given. A numerical comparison with traditional local RBF method is given in the first test problem. 相似文献
15.
D.F. Gilhooley J.R. Xiao R.C. Batra M.A. McCarthy J.W. Gillespie Jr. 《Computational Materials Science》2008,41(4):467-481
The meshless local Petrov–Galerkin (MLPG) method is used for analysing two-dimensional (2D) static and dynamic deformations of functionally graded materials (FGMs) with material response modelled as either linear elastic or as linear viscoelastic. The multiquadric radial basis function (RBF) is employed to approximate the trial solution. Results are computed with two different choices of test functions, namely a fourth-order spline weight function, and a Heaviside step function, each having a compact support. No background mesh is used to numerically evaluate integrals appearing in the weak formulation of the problem, thus the method is truly meshless. A benefit of using RBFs is that they possess the Kronecker delta property; thus it is easy to satisfy essential boundary conditions. For five problems, the computed results are found to match well with those either from their analytical solutions or numerical solutions of other researchers who employed different algorithms. For a dynamic problem, the Laplace-transform technique is utilised. The numerical examples illustrate that displacements and stress distributions in a structure made of an FGM differ considerably from those at the corresponding points in the same structure made of a homogeneous material. Thus, the inhomogeneity in material properties can be exploited to optimise stress distribution, minimise deflection and reduce the maximum stress. 相似文献
16.
A parameterization level set method is presented for structural shape and topology optimization of compliant mechanisms involving large displacements. A level set model is established mathematically as the Hamilton–Jacobi equation to capture the motion of the free boundary of a continuum structure. The structural design boundary is thus described implicitly as the zero level set of a level set scalar function of higher dimension. The radial basis function with compact support is then applied to interpolate the level set function, leading to a relaxation and separation of the temporal and spatial discretizations related to the original partial differential equation. In doing so, the more difficult shape and topology optimization problem is now fully parameterized into a relatively easier size optimization of generalized expansion coefficients. As a result, the optimization is changed into a numerical process of implementing a series of motions of the implicit level set function via an existing efficient convex programming method. With the concept of the shape derivative, the geometrical non‐linearity is included in the rigorous design sensitivity analysis to appropriately capture the large displacements of compliant mechanisms. Several numerical benchmark examples illustrate the effectiveness of the present level set method, in particular, its capability of generating new holes inside the material domain. The proposed method not only retains the favorable features of the implicit free boundary representation but also overcomes several unfavorable numerical considerations relevant to the explicit scheme, the reinitialization procedure, and the velocity extension algorithm in the conventional level set method. Copyright © 2008 John Wiley & Sons, Ltd. 相似文献
17.
Maryam Mohammadi Reza Mokhtari Hamid Panahipour 《Engineering Analysis with Boundary Elements》2013,37(12):1642-1652
The paper introduces a Galerkin method in the reproducing kernel Hilbert space. It is implemented as a meshless method based on spatial trial spaces spanned by the Newton basis functions in the “native” Hilbert space of the reproducing kernel. For the time-dependent PDEs it leads to a system of ordinary differential equations. The method is used for solving the 2D nonlinear coupled Burgers' equations having Dirichlet and mixed boundary conditions. The numerical solutions for different values of Reynolds number (Re) are compared with analytical solutions as well as other numerical methods. It is shown that the proposed method is efficient, accurate and stable for flow with reasonably high Re in the case of Dirichlet boundary conditions. 相似文献
18.
Przemyslaw Litewka Jerzy Rakowski 《International journal for numerical methods in engineering》2001,52(3):273-286
The purpose of this paper is to analyse free vibrations of arches with influence of shear and axial forces taken into account. Arches with various depth of cross‐section and various types of supports are considered. In the calculations, the curved finite element elaborated by the authors is adopted. It is the plane two‐node, six‐degree‐of‐freedom arch element with constant curvature. Its application to the static analysis yields the exact results, coinciding with the analytical ones. This feature results from the use of the exact shape functions in derivation of the element stiffness matrix. In the free vibration analysis the consistent mass matrix is used. It is obtained on the base of the same functions. Their coefficients contain the influences of shear flexibility and compressibility of the arch. The numerical results are compared with the results obtained for the simple diagonal mass matrix representing the lumped mass model. The natural frequencies are also compared with the ones for the continuous arches for which the analytically determined frequencies are known. The advantage of the paper is a thorough analysis of selected examples, where the influences of shear forces, axial forces as well as the rotary and tangential inertia on the natural frequencies are examined. Copyright © 2001 John Wiley & Sons, Ltd. 相似文献
19.
The simulation of macrosegregation as a consequence of solidification of a binary Al-4.5%Cu alloy in a 2-dimensional rectangular enclosure is tackled in the present paper. Coupled volume-averaged governing equations for mass, energy, momentum and species transfer are considered. The phase properties are resolved from the Lever solidification rule, the mushy zone is modeled by the Darcy law and the liquid phase is assumed to behave like an incompressible Newtonian fluid. Double diffusive effects in the melt are modeled by the thermal and solutal Boussinesq hypothesis. The physical model is solved by the novel Local Radial Basis Function Collocation Method (LRBFCM). The involved physical relevant fields are represented on overlapping 5-noded sub-domains through collocation by using multiquadrics Radial Basis Functions (RBF). The involved first and second derivatives of the fields are calculated from the respective derivatives of the RBFs. The fields are solved through explicit time stepping. The pressure-velocity coupling is calculated through a local pressure correction scheme. The evolution of the solidification process is presented through temperature, velocity, liquid fraction and species concentration histories in four sampling points. The fully solidified state is analyzed through final macrosegregation map in three vertical and three horizontal cross-sections. The results are compared with the classical Finite Volume Method (FVM). A surprisingly good agreement of the numerical solution of both methods is shown and therefore the results can be used as a reference for future verification studies. The advantages of the represented meshless approach are its simplicity, accuracy, similar coding in 2D and 3D, and straightforward applicability in non-uniform node arrangements. The paper probably for the first time shows an application of a meshless method in such a highly non-linear and multi-physics problem. 相似文献
20.
Yanan Liu Yinghua Liu Keqin Ding 《International journal for numerical methods in engineering》2017,112(10):1295-1322
In this paper, a coupling technique is developed for the combination of the wavelet‐Galerkin method (WGM) with the finite element method (FEM). In this coupled method, the WGM and FEM are respectively used in different sub‐domains. The WGM sub‐domain and the FEM sub‐domain are connected by a transition region that is described by B‐spline basis functions. The basis functions of WGM and FEM are modified in the transition region to ensure the basic polynomial reconstruction condition and the compatibility of displacements and their gradients. In addition, the elements of FEM and WGM are not necessary to conform to the transition region. This newly developed coupled method is applied to the stress analysis of 2D and 3D elasticity problems in order to investigate its performance and study parameters. Numerical results show that the present coupled method is accurate and stable. The new method has a promising potential for practical applications and can be easily extended for coupling of other numerical methods. Copyright © 2017 John Wiley & Sons, Ltd. 相似文献