In this paper,a substructure method of three-dimensional semi-analytic boundary element is established.The seismic scattering by three-dimensional topography of a hill can be analyzed by the method in frequency domain...In this paper,a substructure method of three-dimensional semi-analytic boundary element is established.The seismic scattering by three-dimensional topography of a hill can be analyzed by the method in frequency domain.Using this method,the computational effort and storage space are reduced considerably.Finally,analytic results are given.展开更多
In this paper,we present the applications of Boundary Element Method(BEM) to simulate the electro-mechanical coupling responses of Micro-Electro-Mechanical systems(MEMS). The algorithm is programmed in our research gr...In this paper,we present the applications of Boundary Element Method(BEM) to simulate the electro-mechanical coupling responses of Micro-Electro-Mechanical systems(MEMS). The algorithm is programmed in our research group based on BEM modeling for electrostatics and elastostatics.Good agreement is shown while the simulation results of the pull-in voltages are compared with the theoretical/experimental ones for some examples.展开更多
This paper is an attempt to solve the soil-pile interaction problems using the boundary element method（BEM）.A computer package called PGroupN,which deals mainly with the analysis of the pile group problem,is employe...This paper is an attempt to solve the soil-pile interaction problems using the boundary element method（BEM）.A computer package called PGroupN,which deals mainly with the analysis of the pile group problem,is employed in this study.Parametric studies are carried out to assess the impacts of the pile diameter,pile length,ratio of spacing to diameter and the thickness of soil stratum.The external load is applied incrementally and,at each increment,a check is made that the stress state at the pile-soil interfaces does not violate the yield criteria.This is achieved by specifying the limited stresses of the soil for the axial pile shaft capacity and end-bearing resistance.The elements of the pile-soil interface yielded can take no additional load,and any increase in load is therefore redistributed between the remaining elements until all elements have failed.Thus,by successive application of loading increments,the entire load-displacement relationship for the pile group is determined.It is found that as the applied load reaches the ultimate bearing capacity of the pile group,all the piles will share the same amount of load.An exception to this case is for the center pile in a group of 9 piles embedded in clay,which is not consistent with the behaviors of the other piles in the group even if the load reaches the ultimate state.For the 4 piles group embedded in clay,the maximum load carried by the base does not exceed 8% of the load carried by each pile with different diameters.This low percentage ascertains that the piles embedded in cohesive soils carry most of the load throughout their shafts.展开更多
This paper presents an elasto-viscoplastic consistent tangent operator (CTO) based boundary element formulation, and application for calculation of path-domain independentJ integrals (extension of the classicalJ integ...This paper presents an elasto-viscoplastic consistent tangent operator (CTO) based boundary element formulation, and application for calculation of path-domain independentJ integrals (extension of the classicalJ integrals) in nonlinear crack analysis. When viscoplastic deformation happens, the effective stresses around the crack tip in the nonlinear region is allowed to exceed the loading surface, and the pure plastic theory is not suitable for this situation. The concept of consistency employed in the solution of increment viscoplastic problem, plays a crucial role in preserving the quadratic rate asymptotic convergence of iteractive schemes based on Newton's method. Therefore, this paper investigates the viscoplastic crack problem, and presents an implicit viscoplastic algorithm using the CTO concept in a boundary element framework for path-domain independentJ integrals. Applications are presented with two numerical examples for viscoplastic crack problems andJ integrals.展开更多
In high seismicity areas, it is important to consider kinematic effects to properly design pile foundations.Kinematic effects are due to the interaction between pile and soil deformations induced by seismic waves. One...In high seismicity areas, it is important to consider kinematic effects to properly design pile foundations.Kinematic effects are due to the interaction between pile and soil deformations induced by seismic waves. One of the effect is the arise of significant strains in weak soils that induce bending moments on piles. These moments can be significant in presence of a high stiffness contrast in a soil deposit. The single pile kinematic interaction problem is generally solved with beam on dynamic Winkler foundation approaches(BDWF) or using continuous models. In this work, a new boundary element method(BEM)based computer code(KIN SP) is presented where the kinematic analysis is preceded by a free-field response analysis. The analysis results of this method, in terms of bending moments at the pile-head and at the interface of a two-layered soil, are influenced by many factors including the soil-pile interface discretization. A parametric study is presented with the aim to suggest the minimum number of boundary elements to guarantee the accuracy of a BEM solution, for typical pile-soil relative stiffness values as a function of the pile diameter, the location of the interface of a two-layered soil and of the stiffness contrast. KIN SP results have been compared with simplified solutions in literature and with those obtained using a quasi-three-dimensional(3D) finite element code.展开更多
A general algorithm is applied to the regularization of nearly singular integrals in the boundary element method of planar potential problems. For linear elements, the strongly singular and hypersingular integrals of ...A general algorithm is applied to the regularization of nearly singular integrals in the boundary element method of planar potential problems. For linear elements, the strongly singular and hypersingular integrals of the interior points very close to boundary were categorized into two forms. The factor leading to the singularity was transformed out of the integral representations with integration by parts, so non-singular regularized formulas were presented for the two forms of integrals. Furthermore, quadratic elements are used in addition to linear ones. The quadratic element very close to the internal point can be divided into two linear ones, so that the algorithm is still valid. Numerical examples demonstrate the effectiveness and accuracy of this algorithm. Especially for problems with curved boundaries, the combination of quadratic elements and linear elements can give more accurate results.展开更多
A radial integral boundary element method(BEM)is used to simulate the phase change problem with a mushy zone in this paper.Three phases,including the solid phase,the liquid phase,and the mushy zone,are considered in t...A radial integral boundary element method(BEM)is used to simulate the phase change problem with a mushy zone in this paper.Three phases,including the solid phase,the liquid phase,and the mushy zone,are considered in the phase change problem.First,according to the continuity conditions of temperature and its gradient on the liquid-mushy interface,the mushy zone and the liquid phase in the simulation can be considered as a whole part,namely,the non-solid phase,and the change of latent heat is approximated by heat source which is dependent on temperature.Then,the precise integration BEM is used to obtain the differential equations in the solid phase zone and the non-solid phase zone,respectively.Moreover,an iterative predictor-corrector precise integration method(PIM)is needed to solve the differential equations and obtain the temperature field and the heat flux on the boundary.According to an energy balance equation and the velocity of the interface between the solid phase and the mushy zone,the front-tracking method is used to track the move of the interface.The interface between the liquid phase and the mushy zone is obtained by interpolation of the temperature field.Finally,four numerical examples are provided to assess the performance of the proposed numerical method.展开更多
Both the orthotropy and the stress concentration are common issues in modem structural engineering. This paper introduces the boundary element method （BEM） into the elastic and elastoplastic analyses for 2D orthotro...Both the orthotropy and the stress concentration are common issues in modem structural engineering. This paper introduces the boundary element method （BEM） into the elastic and elastoplastic analyses for 2D orthotropic media with stress concentration. The discretized boundary element formulations are established, and the stress formulae as well as the fundamental solutions are derived in matrix notations. The numerical procedures are proposed to analyze both elastic and elastoplastic problems of 2D orthotropic me- dia with stress concentration. To obtain more precise stress values with fewer elements, the quadratic isoparametric element formulation is adopted in the boundary discretization and numerical procedures. Numerical examples show that there are significant stress concentrations and different elastoplastic behaviors in some orthotropic media, and some of the computational results are compared with other solutions. Good agreements are also observed, which demonstrates the efficiency and reliability of the present BEM in the stress concentration analysis for orthotropic media.展开更多
In this paper, a method of transforming volume integrals to boundary integrals is given for complicated loadings such as a i(y)x i and b i(x)y i . In the present method the volume in...In this paper, a method of transforming volume integrals to boundary integrals is given for complicated loadings such as a i(y)x i and b i(x)y i . In the present method the volume integrals are approximately transformed to boundary integrals.展开更多
For the discretization of higher order elements, this paper presents a modified integral domain method to remove the irregular frequencies inherited in the integral equation of wave diffraction and radiation from a su...For the discretization of higher order elements, this paper presents a modified integral domain method to remove the irregular frequencies inherited in the integral equation of wave diffraction and radiation from a surface-piercing body. The set of over-determined linear equations obtained from the method is modified into a normal set of linear equations by superposing a set of linear equations with zero solutions. Numerical experiments have also been carried out to find the optimum choice of the size of the auxiliary domain and the discretization on it.展开更多
Cavity resonance noise of passenger car tires is generated by interacting excitation between a tire structure and the fill gas (air), and generally lies in a frequency range of 200?250 Hz. As such, this noise is stron...Cavity resonance noise of passenger car tires is generated by interacting excitation between a tire structure and the fill gas (air), and generally lies in a frequency range of 200?250 Hz. As such, this noise is strongly perceived and may be a serious source of driver annoyance. Thus, many studies regarding the cavity noise mechanism and its reduction have already been conducted. In this work, a vibro-acoustic coupled analysis was conducted between a tire structure and air cavity. Using this analysis, we can more accurately simulate the tire noise performance in the region of the cavity resonance frequency. An analysis of the effects of variation of tire contour design factors was conducted, using design-of-experiments methods. Finally, a multi-objective optimization was performed using in-house codes to reduce the cavity noise level while minimizing the loss of other performances, such as diminished ride comfort and handling caused by the variations of contour. As a result of this optimization, an optimized contour shape was derived, which satisfied the multi-objective performances.展开更多
基金This project was sponsored by the Earthquake Science Foundation, China
文摘In this paper,a substructure method of three-dimensional semi-analytic boundary element is established.The seismic scattering by three-dimensional topography of a hill can be analyzed by the method in frequency domain.Using this method,the computational effort and storage space are reduced considerably.Finally,analytic results are given.
基金The project supported by the 973 Program (G1999033108)the National Natural Science Foundation of China (10125211)
文摘In this paper,we present the applications of Boundary Element Method(BEM) to simulate the electro-mechanical coupling responses of Micro-Electro-Mechanical systems(MEMS). The algorithm is programmed in our research group based on BEM modeling for electrostatics and elastostatics.Good agreement is shown while the simulation results of the pull-in voltages are compared with the theoretical/experimental ones for some examples.
文摘This paper is an attempt to solve the soil-pile interaction problems using the boundary element method（BEM）.A computer package called PGroupN,which deals mainly with the analysis of the pile group problem,is employed in this study.Parametric studies are carried out to assess the impacts of the pile diameter,pile length,ratio of spacing to diameter and the thickness of soil stratum.The external load is applied incrementally and,at each increment,a check is made that the stress state at the pile-soil interfaces does not violate the yield criteria.This is achieved by specifying the limited stresses of the soil for the axial pile shaft capacity and end-bearing resistance.The elements of the pile-soil interface yielded can take no additional load,and any increase in load is therefore redistributed between the remaining elements until all elements have failed.Thus,by successive application of loading increments,the entire load-displacement relationship for the pile group is determined.It is found that as the applied load reaches the ultimate bearing capacity of the pile group,all the piles will share the same amount of load.An exception to this case is for the center pile in a group of 9 piles embedded in clay,which is not consistent with the behaviors of the other piles in the group even if the load reaches the ultimate state.For the 4 piles group embedded in clay,the maximum load carried by the base does not exceed 8% of the load carried by each pile with different diameters.This low percentage ascertains that the piles embedded in cohesive soils carry most of the load throughout their shafts.
基金The project supported by National Natural Science Foundation of China(9713008)Zhejiang Natural Science Foundation Special Funds No. RC.9601
文摘This paper presents an elasto-viscoplastic consistent tangent operator (CTO) based boundary element formulation, and application for calculation of path-domain independentJ integrals (extension of the classicalJ integrals) in nonlinear crack analysis. When viscoplastic deformation happens, the effective stresses around the crack tip in the nonlinear region is allowed to exceed the loading surface, and the pure plastic theory is not suitable for this situation. The concept of consistency employed in the solution of increment viscoplastic problem, plays a crucial role in preserving the quadratic rate asymptotic convergence of iteractive schemes based on Newton's method. Therefore, this paper investigates the viscoplastic crack problem, and presents an implicit viscoplastic algorithm using the CTO concept in a boundary element framework for path-domain independentJ integrals. Applications are presented with two numerical examples for viscoplastic crack problems andJ integrals.
文摘In high seismicity areas, it is important to consider kinematic effects to properly design pile foundations.Kinematic effects are due to the interaction between pile and soil deformations induced by seismic waves. One of the effect is the arise of significant strains in weak soils that induce bending moments on piles. These moments can be significant in presence of a high stiffness contrast in a soil deposit. The single pile kinematic interaction problem is generally solved with beam on dynamic Winkler foundation approaches(BDWF) or using continuous models. In this work, a new boundary element method(BEM)based computer code(KIN SP) is presented where the kinematic analysis is preceded by a free-field response analysis. The analysis results of this method, in terms of bending moments at the pile-head and at the interface of a two-layered soil, are influenced by many factors including the soil-pile interface discretization. A parametric study is presented with the aim to suggest the minimum number of boundary elements to guarantee the accuracy of a BEM solution, for typical pile-soil relative stiffness values as a function of the pile diameter, the location of the interface of a two-layered soil and of the stiffness contrast. KIN SP results have been compared with simplified solutions in literature and with those obtained using a quasi-three-dimensional(3D) finite element code.
文摘A general algorithm is applied to the regularization of nearly singular integrals in the boundary element method of planar potential problems. For linear elements, the strongly singular and hypersingular integrals of the interior points very close to boundary were categorized into two forms. The factor leading to the singularity was transformed out of the integral representations with integration by parts, so non-singular regularized formulas were presented for the two forms of integrals. Furthermore, quadratic elements are used in addition to linear ones. The quadratic element very close to the internal point can be divided into two linear ones, so that the algorithm is still valid. Numerical examples demonstrate the effectiveness and accuracy of this algorithm. Especially for problems with curved boundaries, the combination of quadratic elements and linear elements can give more accurate results.
基金the National Natural Science Foundation of China(No.11672064)。
文摘A radial integral boundary element method(BEM)is used to simulate the phase change problem with a mushy zone in this paper.Three phases,including the solid phase,the liquid phase,and the mushy zone,are considered in the phase change problem.First,according to the continuity conditions of temperature and its gradient on the liquid-mushy interface,the mushy zone and the liquid phase in the simulation can be considered as a whole part,namely,the non-solid phase,and the change of latent heat is approximated by heat source which is dependent on temperature.Then,the precise integration BEM is used to obtain the differential equations in the solid phase zone and the non-solid phase zone,respectively.Moreover,an iterative predictor-corrector precise integration method(PIM)is needed to solve the differential equations and obtain the temperature field and the heat flux on the boundary.According to an energy balance equation and the velocity of the interface between the solid phase and the mushy zone,the front-tracking method is used to track the move of the interface.The interface between the liquid phase and the mushy zone is obtained by interpolation of the temperature field.Finally,four numerical examples are provided to assess the performance of the proposed numerical method.
基金The project supported by the Basic Research Foundation of Tsinghua University,the National Foundation for Excellent Doctoral Thesis(200025)the National Natural Science Foundation of China(19902007).
文摘Both the orthotropy and the stress concentration are common issues in modem structural engineering. This paper introduces the boundary element method （BEM） into the elastic and elastoplastic analyses for 2D orthotropic media with stress concentration. The discretized boundary element formulations are established, and the stress formulae as well as the fundamental solutions are derived in matrix notations. The numerical procedures are proposed to analyze both elastic and elastoplastic problems of 2D orthotropic me- dia with stress concentration. To obtain more precise stress values with fewer elements, the quadratic isoparametric element formulation is adopted in the boundary discretization and numerical procedures. Numerical examples show that there are significant stress concentrations and different elastoplastic behaviors in some orthotropic media, and some of the computational results are compared with other solutions. Good agreements are also observed, which demonstrates the efficiency and reliability of the present BEM in the stress concentration analysis for orthotropic media.
文摘In this paper, a method of transforming volume integrals to boundary integrals is given for complicated loadings such as a i(y)x i and b i(x)y i . In the present method the volume integrals are approximately transformed to boundary integrals.
基金This work is a part of the research project financially supported by the National Natural Science Foundation of China
文摘For the discretization of higher order elements, this paper presents a modified integral domain method to remove the irregular frequencies inherited in the integral equation of wave diffraction and radiation from a surface-piercing body. The set of over-determined linear equations obtained from the method is modified into a normal set of linear equations by superposing a set of linear equations with zero solutions. Numerical experiments have also been carried out to find the optimum choice of the size of the auxiliary domain and the discretization on it.
文摘Cavity resonance noise of passenger car tires is generated by interacting excitation between a tire structure and the fill gas (air), and generally lies in a frequency range of 200?250 Hz. As such, this noise is strongly perceived and may be a serious source of driver annoyance. Thus, many studies regarding the cavity noise mechanism and its reduction have already been conducted. In this work, a vibro-acoustic coupled analysis was conducted between a tire structure and air cavity. Using this analysis, we can more accurately simulate the tire noise performance in the region of the cavity resonance frequency. An analysis of the effects of variation of tire contour design factors was conducted, using design-of-experiments methods. Finally, a multi-objective optimization was performed using in-house codes to reduce the cavity noise level while minimizing the loss of other performances, such as diminished ride comfort and handling caused by the variations of contour. As a result of this optimization, an optimized contour shape was derived, which satisfied the multi-objective performances.
