共查询到20条相似文献,搜索用时 742 毫秒
1.
A variant of the boundary element method, called the boundary contour method (BCM), offers a further reduction in dimensionality. Consequently, boundary contour analysis of two-dimensional problems does not require any numerical integration at all. In another development, a boundary contour implementation of a regularized hypersingular boundary integral equation (HBIE) using quadratic elements and end-node collocation was proposed and the technique is termed the hypersingular boundary contour method (HBCM). As reported in that work, the approach requires highly refined meshes in order to numerically enforce the stress continuity across boundary contour elements. This continuity requirement is very crucial since the regularized HBIE is only valid at collocation points where the stress tensor is continuous, while the computed stress at the endpoints of a boundary contour element, which is a non-conforming element, is generally not. This paper presents a new implementation of the HBCM for which the regularized HBIE is collocated at the mid-node of a boundary contour element. As the computed stress tensor is continuous at these mid-nodes, there is no need for unusually refined meshes. Some numerical tests herein show that, for the same mesh density, the HBCM using mid-node collocation has a comparable accuracy as the BCM. 相似文献
2.
Summary The computation of interior displacements and stresses with the boundary element method (BEM) often requires the evaluation of nearly singular integrals. These integrals arise from the singular behaviour of the kernel functions in the Somigliana identity and the Somigliana stress identity. Treating them numerically in a standard way leads to inaccuracy near the boundary. This effect is always present in the calculation of field variables near the boundary and is called boundary layer effect. In this paper regularization procedures are proposed which consist of an indirect evaluation of singular integrals and a special coordinate transformation. The proposed procedures eliminate the boundary layer effect for both, the calculation of displacements and stresses. In a numerical example of elastostatics the developed strategies are shown to work. Due to the generality of the proposed procedures they can be extended to any standard boundary element formulation for problems with bounded domains. 相似文献
3.
4.
An experimental-numerical hybrid method for the stress separation in photoelasticity is proposed in this study. In the proposed
method, boundary conditions for a local finite element model, that is, tractions along boundaries are inversely determined
from photoelastic fringes. Two algorithms are proposed for determining the boundary condition. One is a linear algorithm in
which the tractions are obtained by the method of linear least-squares from both principal stress difference and principal
direction. Another is the nonlinear algorithm in which the tractions are determined only from the principal stress difference.
After determining the boundary conditions for the local finite element model, the stresses can be obtained by finite element
direct analysis. The effectiveness is demonstrated by applying the proposed method to a perforated plate under tension and
contact problems. Results show that the boundary conditions of the local finite element model can be determined from the photoelastic
fringes and then the individual stresses can be obtained by the proposed method. Furthermore, the stresses can be evaluated
even if the boundary condition is complicated such as at the contact surface. It is expected that the proposed method can
be powerful tool for stress analysis. 相似文献
5.
对于体积力问题,可将体积力的域内积分转换为边界上的积分运算[1],这将减少域内划分单元,使边界元法处理体积力时更简秉方便,本文在文[1]和[2]的基础上,提出用单元的解析积分取代数值积分,计算离心体力问题,推导了相应的直线单元即常单元和线性单元的解析式,算例表明,本文方法数值精度高,尤其在造近边界的内点时,效果迹较好。 相似文献
6.
本文通过引入拉氏因子,以带线性约束最小二乘逼近的Golub方法,很好地处理了结点应力拟合的边值问题.利用三角形快速插值原理求出了等值线,并以CAD为环境绘制了单元图和等值线图,实现了有限元分析中高精度等值线的自动处理. 相似文献
7.
A direct formulation of the boundary element method using a complex variable numerical approach is presented for the time-dependent inelastic stresses in edge notched and cracked creeping metallic structural components subject to high temperature gradients. Particular attention is focused on the numerical evaluation of energy rate contour integrals in single edge cracked specimens in tension. The constitutive models used in the numerical calculation are internal state variable creep–plasticity or elastic power law creep model. Numerical results are compared with solution obtained from other methods for different loading rates. 相似文献
8.
Summary A boundary element formulation is presented for the solution of the equations of fully coupled thermoelasticity for materials of arbitrary degree of anisotropy. By employing the fundamental solutions of anisotropic elastostatics and stationary heat conduction, a system of equations with time-independent matrices is obtained. Since the fundamental solutions are uncoupled and time-independent, a domain integral remains in the representation formula which contains the time-dependence as well as the thermoelastic coupling. This domain integral is transformed to the boundary by means of the dual reciprocity method. By taking this approach, the use of dynamic fundamental solutions is avoided, which enables an efficient calculation of system matrices. In addition, the solution of transient processes as well as, free and forced vibration analysis becomes straightforward and can be carried out with standard time-stepping schemes and eigensystem solvers. Another important advantage of the present formulation is its versatility, since it includes a number of simplified thermoelastic theories, viz. the theory of thermal stresses, coupled and uncoupled quasi-static thermoelasticity, and stationary thermoelasticity. The accuracy of the new thermoelastic boundary element method is demonstrated by a number of example problems.
Support by the Deutsche Forschungsgemeinschaft (DFG) of the Graduate Collegium Modelling and discretization methods for continua and fluids (GKKS) at the University of Stuttgart is gratefully acknowledged. 相似文献
9.
Stochastic boundary element method in elasticity 总被引:1,自引:0,他引:1
The stochastic boundary element method is developed to analyze elasticity problems with random material and/or geometrical
parameters and randomly perturbed boundaries. Based on the first-order Taylor series expansion, the boundary integration equations
concerning the mean and deviation of the displacements are derived, respectively. It is found that the randomness of material
parameters is equivalent to a random body force, so the mean and covariance matrices of unknown boundary displacements and
tractions can be obtained. Furthermore, the mean and covariance of displacements and stresses at inner points can also be
obtained. Numerical examples show that the proposed stochastic boundary element method gives satisfactory solutions, as compared
with those obtained by theoretical analysis or other numerical methods.
The project supported by the National Natural Science Foundation of China and the State Education Commission Foundation of
China 相似文献
10.
11.
12.
The flow about submerged, fully cavitating axisymmetric bodies at both zero and non-zero angle of attack is considered in this paper. A cavity closure model that relates the point of detachment, the angle that the separating streamline makes with the body and the cavity length is described. The direct boundary element method is used to solve the potential flow problem and to determine the cavity shape. A momentum integral boundary layer solver is included in the formulation so that shear stresses can be incorporated into the drag calculations. The numerical predictions based on the proposed closure model are compared with water tunnel measurements and photographs. 相似文献
13.
This paper presents a new strategy of using the radial integration boundary element method (RIBEM) to solve non-homogeneous heat conduction and thermoelasticity problems. In the method, the evaluation of the radial in-tegral which is used to transform domain integrals to equivalent boundary integrals is carried out on the basis of elemental nodes. As a result, the computational time spent in evaluating domain integrals can be saved considerably in comparison with the conventional RIBEM. Three numerical examples are given to demonstrate the correctness and computational efficiency of the proposed approach. 相似文献
14.
A parallel fast multipole BEM and its applications to large-scale analysis of 3-D fiber-reinforced composites 总被引:2,自引:0,他引:2
In this paper, an adaptive boundary element method (BEM) is presented for solving 3-D elasticity problems. The numerical scheme is accelerated by the new version of fast multipole method (FMM) and parallelized on distributed memory architectures. The resulting solver is applied to the study of representative volume element (RVE) for short fiber-reinforced composites with complex inclusion geometry. Numerical examples performed on a 32-processor cluster show that the proposed method is both accurate and efficient, and can solve problems of large size that are challenging to existing state-of-the-art domain methods.The project supported by the National Natural Science Foundation of China (10472051).The English text was polished by Yunming Chen. 相似文献
15.
A mixed boundary element and finite element numerical algorithm for the simultaneous prediction of the electric fields, viscous flow fields, thermal fields and surface deformation of electrically conducting droplets in an electrostatic field is described in this paper. The boundary element method is used for the computation of the electric potential distribution. This allows the boundary conditions at infinity to be directly incorporated into the boundary integral formulation, thereby obviating the need for discretization at infinity. The surface deformation is determined by solving the normal stress balance equation using the weighted residuals method. The fluid flow and thermal fields are calculated using the mixed finite element method. The computational algorithm for the simultaneous prediction of surface deformation and fluid flow involves two iterative loops, one for the electric field and surface deformation and the other for the surface tension driven viscous flows. The two loops are coupled through the droplet surface shapes for viscous fluid flow calculations and viscous stresses for updating the droplet shapes. Computing the surface deformation in a separate loop permits the freedom of applying different types of elements without complicating procedures for the internal flow and thermal calculations. Tests indicate that the quadratic, cubic spline and spectral boundary elements all give approximately the same accuracy for free surface calculations; however, the quadratic elements are preferred as they are easier to implement and also require less computing time. Linear elements, however, are less accurate. Numerical simulations are carried out for the simultaneous solution of free surface shapes and internal fluid flow and temperature distributions in droplets in electric fields under both microgravity and earthbound conditions. Results show that laser heating may induce a non-uniform temperature distribution in the droplets. This non-uniform thermal field results in a variation of surface tension along the surface of the droplet, which in turn produces a recirculating fluid flow in the droplet. The viscous stresses cause additional surface deformation by squeezing the surface areas above and below the equator plane. 相似文献
16.
平均源边界节点法ASBNM是一种最近提出的边界型无网格法。该方法仅使用边界节点不涉及任何单元和积分的概念,具有方法简单和程序设计容易等特点。但是,对于依赖于边界积分方程的边界型无网格法,关键问题是如何准确高效地估计影响矩阵的对角元。本文提出直接计算影响矩阵对角元的方法,是已有ASBNM法的改进,将对角元的计算转化为一个纯几何问题,因此适用于任何二维边值问题。数值算例证明了本文方法的有效性和准确性。 相似文献
17.
A stress function-based approach is proposed to analyze the free-edge interlaminar stresses of piezo-bonded symmetric laminates. The proposed method satisfies the traction free boundary conditions, as well as surface free conditions. The symmetric laminated structure was excited under electric fields that can generate induced strain, resulting in pure extension in the laminated plate. The governing equations were obtained by taking the principle of complementary virtual work. To verify the proposed method, cross-ply, angle-ply and quasi-isotropic laminates were analyzed. The stress concentrations predicted by the present method were compared with those analyzed by the finite element method. The results show that the stress function-based analysis of piezo-bonded laminated composite structures is an efficient and accurate method for the initial design stage of piezo-bonded composite structures. 相似文献
18.
The present paper presents a general treatment of the transient thermoelastic stresses in a rotating nonhomogeneous anisotropic
solid under compressive initial stress. The system of fundamental equations is solved by means of a boundary element method
(BEM) and the numerical calculations are carried out for the temperature, displacement components and stress components. The
results indicate that the effects of inhomogeneity and initial stress are very pronounced. 相似文献
19.
Xiangqiao Yan 《Mechanics Research Communications》2004,31(2):263
In this paper, a numerical analysis of perpendicular cracks under general in-plane loading is performed by using a hybrid displacement discontinuity method which consists of the non-singular displacement discontinuity element presented by Crouch and Starfied and the crack tip displacement discontinuity elements by the author. In the boundary element implementation the left or the right crack tip displacement discontinuity element is placed locally at corresponding left or right crack tip on top of the ordinary non-singular displacement discontinuity elements that cover the entire crack surface and the other boundary. The present numerical results show that the numerical approach is simple, yet very accurate for calculating numerically stress intensity factors for perpendicular cracks under general in-plane loading. 相似文献
20.
A boundary element method is proposed for studying periodic shallow water problems. The numerical model is based on the shallow water equation. The key feature of this method is that the boundary integral equations are derived using the weighted residual method and the fundamental solutions for shallow water wave problems are obtained by solving the simultaneous singular equations. The accuracy of this method is studied for the wave reflection problem in a rectangular tank. As a result of this test, it has been shown that the number of element divisions and the distribution of nodes are significant to the accuracy. For numerical examples of external problems, the wave diffraction problems due to single cylindrical, double cylindrical and plate obstructions are analysed and compared with the exact and other numerical solutions. Relatively accurate solutions are obtained. 相似文献