The various factors influencing the settlement of composite foundation have been more completely studied through numerical simulation. The influence on the settlement of composite foundation of the geometry and mechan...The various factors influencing the settlement of composite foundation have been more completely studied through numerical simulation. The influence on the settlement of composite foundation of the geometry and mechanical properties of the pile, soil, cushion, and the interface between pile and soil have been investigated through computer simulation, in which the contact elements with zero thickness are used. Some valuable conclusions for the settlement of composite foundation have been obtained: (1) The method using the contact element of zero thickness is successful when used in the simulation of the settlement of composite foundation; (2) Among the factors influencing the settlement of composite foundation, the compression modulus of the soil is the largest, and the cohesion of the soil is the second largest; (3) The effects on settlement of the internal friction angle of the soil, the elastic modulus, the radius, and the length of the pile, and the elestic modulus of the cushion are also more obvious.展开更多
A uniform transient temperature field model of electrical contacts operation was found by analyzing the process of closing arc constriction resistance Joule heat ~ breaking arc. Essential parameters of Ag/La2NiO4 elec...A uniform transient temperature field model of electrical contacts operation was found by analyzing the process of closing arc constriction resistance Joule heat ~ breaking arc. Essential parameters of Ag/La2NiO4 electrical contact material for transient temperature field calculation were obtained through tests of electrical contact experimental instrument under 18 V DC in different cur- rents, other correlation experiments, and calculation anal- ysis. The finite element method was applied to solve the transient temperature field, and the features and distribution of the transient temperature field were obtained. The condition of material erosion and mass transfer can be forecasted by those calculation results. It is beneficial to research about the lifetime of Ag/La2NiO4 electrical material.展开更多
An improved version of the regular boundary element method, the artificial boundary node approach, is derived. A simple contact algorithm is designed and implemented into the direct boundary element, regular boundary ...An improved version of the regular boundary element method, the artificial boundary node approach, is derived. A simple contact algorithm is designed and implemented into the direct boundary element, regular boundary element and artificial boundary node approaches. The exisiting and derived approaches are tested using some case studies. The results of the artificial boundary node approach are compared with those of the existing boundary element program, the regular element approach, ANSYS and analytical solution whenever possible. The results show the effectiveness of the artificial boundary node approach for a wider range of boundary offsets.展开更多
The evolution of leakage is studied using detailed contact finite element analysis. The distribution of stress at the gasket is analyzed using a contact condition based on slide-line elements using ABAQUS, a commercia...The evolution of leakage is studied using detailed contact finite element analysis. The distribution of stress at the gasket is analyzed using a contact condition based on slide-line elements using ABAQUS, a commercial finite element code, Slide-line elements also take into account pressure penetration as contact that is lost between flange and gasket. Results are presented for a particular flange, a raised face flange sealed by a mild steel gasket. A comparison of the results from the gasket contact analysis and the contact conditions specified by the ASME Boiler and Pressure Vessel Code, Sections VIII, Division 1 shows that the conditions specified in the ASME Code predict leakage relatively accurately.展开更多
A scheme of boundary element method for moving contact of two-dimensional elastic bodies using conforming discretization is presented. Both the displacement and the traction boundary conditions are satisfied on the co...A scheme of boundary element method for moving contact of two-dimensional elastic bodies using conforming discretization is presented. Both the displacement and the traction boundary conditions are satisfied on the contacting region in the sense of discretization. An algorithm to deal with the moving of the contact boundary on a larger possible contact region is presented. The algorithm is generalized to rolling contact problem as well. Some numerical examples of moving and rolling contact of 2D elastic bodies with or without friction, including the bodies with a hole-type defect, are given to show the effectiveness and the accuracy of the presented schemes.展开更多
Each surface of roughness has different shape of asperity which is modeled with various shapes of analytical models. In this paper, the differences among various models of shape of asperity investigate using the Finit...Each surface of roughness has different shape of asperity which is modeled with various shapes of analytical models. In this paper, the differences among various models of shape of asperity investigate using the Finite Element Method (FEM) and various analytical models. The contact stresses in rough surfaces are calculated analytically using various asperity shape models. Finite element analysis is also carried out assuming three types of material properties namely, the linear, the elastic-perfect plastic and the elastic-nonlinear hardening. The analytical results are compared with the results obtained by the finite element method. The results illustrate for using a deterministic approach which the numerical models are suitable. In hertz model, the result of force is very big in interface of causing deformation plastic, while Model Zhao has almost same result with FEM nonlinear property model. It is observed that the results obtained from Zhao’s model are generally in a better agreement with the results obtained from various finite element models especially in elastic-plastic and plastic zones, hence it may be concluded that Zhao’s model can be used for analyzing the rough surfaces in contact mechanics.展开更多
Ceramics are good alternative to metal as bearing couple materials because of their better wear resistance. A Finite Element(FE) study was performed to investigate the contact mechanics and stress distribution of Cera...Ceramics are good alternative to metal as bearing couple materials because of their better wear resistance. A Finite Element(FE) study was performed to investigate the contact mechanics and stress distribution of Ceramic-on-Ceramic (COC) hip resurfacingprostheses. It was focused in particular on a parametric study to examine the effects of radial clearance, loading,alumina coating on the implants, bone quality, and fixation of cup-bone interface. It was found that a reduction in the radialclearance had the most significant effect on the predicted contact pressure distribution among all of the parameters considered inthis study. It was determined that there was a significant influence of non-metallic materials, such as the bone underneath thebearing components, on the predicted contact mechanics. Stress shielding within the bone tissue was found to be a major concernwhen regarding the use of ceramic as an alternative to metallic resurfacing prostheses. Therefore, using alumina implantswith a metal backing was found to be the best design for ceramic resurfacing prostheses in this study. The loading, bone quality,and acetabular cup fixation conditions were found to have only minor effects on the predicted contact pressure distribution alongthe bearing surfaces.展开更多
Vertical rigidity of the space self adaptive 530 high rigidity mill is calculated by applying the boundary element method (BEM) of three dimension elastic contact problem,which can update the existed deforming s...Vertical rigidity of the space self adaptive 530 high rigidity mill is calculated by applying the boundary element method (BEM) of three dimension elastic contact problem,which can update the existed deforming separation calculating theory and corresponding methods of material mechanics,elastic mechanics and finite element method.The method has less hypotheses and stronger synthesis in contact type calculating model.The advantages of the method are high calculating rate,high calculating accuracy,etc..展开更多
The discrete element method(DEM) has been extensively adopted to investigate many complex geotechnical related problems due to its capability to incorporate the discontinuous nature of granular materials. In particula...The discrete element method(DEM) has been extensively adopted to investigate many complex geotechnical related problems due to its capability to incorporate the discontinuous nature of granular materials. In particular, when simulating large deformations or distortion of soil(e.g. cavity expansion),DEM can be very effective as other numerical solutions may experience convergence problems. Cavity expansion theory has widespread applications in geotechnical engineering, particularly to the problems concerning in situ testing, pile installation and so forth. In addition, the behaviour of geomaterials in a macro-level is utterly determined by microscopic properties, highlighting the importance of contact models. Despite the fact that there are numerous contact models proposed to mimic the realistic behaviour of granular materials, there are lack of studies on the effects of these contact models on the soil response.Hence, in this study, a series of three-dimensional numerical simulations with different contact constitutive models was conducted to simulate the response of sandy soils during cylindrical cavity expansion. In this numerical investigation, three contact models, i.e. linear contact model, rolling resistance contact model,and Hertz contact model, are considered. It should be noted that the former two models are linear based models, providing linearly elastic and frictional plasticity behaviours, whereas the latter one consists of nonlinear formulation based on an approximation of the theory of Mindlin and Deresiewicz. To examine the effects of these contact models, several cylindrical cavities were created and expanded gradually from an initial radius of 0.055 m to a final radius of 0.1 m. The numerical predictions confirm that the calibrated contact models produced similar results regarding the variations of cavity pressure, radial stress, deviatoric stress, volumetric strain, as well as the soil radial displacement. However, the linear contact model may result in inaccurate predictions when highly angular soil particles are involved. In addition, considering the excessive soil displacement induced by the pile installation(i.e. cavity expansion), a minimum distance of11 a(a is the cavity radius) is recommend for practicing engineers to avoid the potential damages to the existing piles and adjacent structures.展开更多
The thrust contact between the Togo and the Dahomeyan formations in Ghana is a lithotectonic boundary that exists between two major Precambrian formations which are of importance to geologists owing to the fact that P...The thrust contact between the Togo and the Dahomeyan formations in Ghana is a lithotectonic boundary that exists between two major Precambrian formations which are of importance to geologists owing to the fact that Precambrian rocks in Ghana host almost all economic minerals and metals. The lithostratigraphy of the Togo-Dahomyan thrust contact rocks from a borehole in Kwabenya near Accra (the Capital of Ghana) has been studied and major crustal chemical elements assayed using Instrumental Neutron Activation Analysis (INAA) techniques. The results have revealed elemental compositions and the mineralogical make up of the lithostratigraphic units of these contact rocks and the general geology of the Togo and the Dahomeyan formations. The profile shows a thrust contact that exists between the Togo formations which are underlying the Dahomeyan formations. The Togo formations here are made mainly of quartzite in fresh schist and the Dahomeyan made of gneisses. In between these two major geologic formations are the rocks of the contact which are intercalation of quartzite and schist. The rocks are felsic with an average felsic index (F) of 85.74 and feldspar rich with K- (orthoclase) feldspars dominating the Dahomeyan rocks. Iron and titanium oxides are depleted with depth from 52 m depth below surface downward and potassium oxide was enriched with depth from 42 m below surface downward. Major mineral forming elements such as aluminum and calcium had varied levels in the Togo and in the Dahomeyan rocks.展开更多
A simple interface element for analyzing contact friction problems is developed. Taking nodal displacements and contact stresses as unknowns, this element can simulate frictional slippage, decoupling and re-bonding of...A simple interface element for analyzing contact friction problems is developed. Taking nodal displacements and contact stresses as unknowns, this element can simulate frictional slippage, decoupling and re-bonding of two bodies initially mating or having gaps at a common interface. The method is based on the Finite Element Method and load incremental theory. The geometric and static constraint conditions on contact surfaces are treated as additional conditions and are included in stiffness equations. This simple element has the advantages of easy implementation into standard finite element programs and fast speed for convergence as well as high accuracy for stress distribution in interface. Undesirable stress oscillations are also investigated whenever large stress gradients exist over the contact surfaces. Exact integration or the conventional Gauss integration scheme used to evaluate the interpolation function matrix of the interface element is found to be the source of the oscillations. Eigenmode analysis demonstrates that the stress behavior of an interface element can be improved by using the Newton-Cotes integration scheme. Finally, the test example of a strip footing problem is presented.展开更多
An algorithm is presented for controlling two-dimensional motion contact bodies with conforming discretization. Since a kind of special boundary element is utilized in the algorithm, the displacement compatibility and...An algorithm is presented for controlling two-dimensional motion contact bodies with conforming discretization. Since a kind of special boundary element is utilized in the algorithm, the displacement compatibility and traction equilibrium conditions at nodes can be satisfied simultaneously in arbitrary locations of the contact interface. In addition, a method is also proposed in which the contact boundary location can be moved flexibly on the possible contact boundary. This method is effective to deal with moving and rolling contact problems on a possible larger moving or rolling contact region. Numerical examples show effectiveness of the presented scheme.展开更多
The adaptive element techniques of contact problem are studied by means of penalty method, and the error estimators are discussed. Based on error estimators, algorithm of the adaptive element techniques is developed, ...The adaptive element techniques of contact problem are studied by means of penalty method, and the error estimators are discussed. Based on error estimators, algorithm of the adaptive element techniques is developed, then the Gauss - Newton iterations are used which allow the nonlinear problem to be transformed into a sequence of linear sub- problems then easily solved. In addition, the algorithm can be applied into the simulation of de -bonding of fiber - reinforced composites.展开更多
The boundary element method in framework is given to evaluate three dimensional frictional contact problems. Elasto plastic material behavior is taken into account by mean of an initial stress formulation and Von Mis...The boundary element method in framework is given to evaluate three dimensional frictional contact problems. Elasto plastic material behavior is taken into account by mean of an initial stress formulation and Von Mises yield criterion. The amount of tangential traction at contact surface is limited by Coulomb's friction law and constant shear rule. From some numerical results of a plate rolling problem, it is demonstrated here that the BEM can be used to efficiently and accurately analyze this class of forming problems.展开更多
Through defining slide yield function and floating potential function of thermo-contact surface, the complementary equation of thermo-contact boundary has been reached, the fundamental equations to solve 3D thermo-con...Through defining slide yield function and floating potential function of thermo-contact surface, the complementary equation of thermo-contact boundary has been reached, the fundamental equations to solve 3D thermo-contact coupled problem have been listed. On this foundation, the finite element equation and definite solution condition of contact heat transfer have been given out. Based on virtual work principle and contact element technology, the finite element equation of 3D elastic contact system has been deduced under the effect of thermal stress. The pseudo load brought by contact gap have been introduced into this equation in order to reflect the contact state change. During iteration, once contact rigidity matrix is formed, it won’t change, which will make calculation reduce greatly.展开更多
The formulation of boundary element method for handling contact problems with friction and the technique for high speed contact analysis are presented. This formulation is based on the idea of modifying the length of...The formulation of boundary element method for handling contact problems with friction and the technique for high speed contact analysis are presented. This formulation is based on the idea of modifying the length of contact elements without altering the total number of elements. The high precision of solution and high speed analysis are verified according to the results of conventional method and analysis method.展开更多
A hydromechanical interface element is proposed for the consideration of the hydraulic-mechanical coupling effect along the interface.The fully coupled governing equations and the relevant finite element formulations ...A hydromechanical interface element is proposed for the consideration of the hydraulic-mechanical coupling effect along the interface.The fully coupled governing equations and the relevant finite element formulations are derived in detail for the interface element.All the involved matrices are of the same form as those of a solid element,which makes the incorporation of the model into a finite element program straightforward.Three examples are then numerically simulated using the interface element.Reasonable results confirm the correctness of the proposed model and motivate its application in hydromechanical contact problems in the future.展开更多
A new contact searching algorithm for contact-impact systems is proposed in this paper.In terms of the cell structure and the linked-list,this algo- rithm solves the problem of sorting and searching contacts in three ...A new contact searching algorithm for contact-impact systems is proposed in this paper.In terms of the cell structure and the linked-list,this algo- rithm solves the problem of sorting and searching contacts in three dimensions by transforming it to a retrieving process from two one-dimensional arrays,and binary searching is no longer required.Using this algorithm, the cost of contact searching is reduced to the order of O(N)instead of O(Nlog_2N)for traditional ones,where N is the node number in the system.Moreover,this algorithm can handle contact systems with arbitrary mesh layouts.Due to the simplicity of this algorithm it can be easily implemented in a dynamic explicit finite element program.Our numerical experi- mental result shows that this algorithm is reliable arid efficient for contact searching of three dimensional systems.展开更多
A new nonlinear force model based on experimental data is proposed to replace the classical Hertzian contact model to solve the fractional index nonlinearity in a ball bearing system. Firstly, the radial force and the...A new nonlinear force model based on experimental data is proposed to replace the classical Hertzian contact model to solve the fractional index nonlinearity in a ball bearing system. Firstly, the radial force and the radial deformation are measured by statics experiments, and the data are fitted respectively by using the Hertzian contact model and the cubic polynomial model. Then~ the two models are compared with the approximation formula appearing in Aeroengine Design Manual. In consequence, the two models are equivalent in an allowable deformation range. After that, the relationship of contact force and contact deformation for single rolling element between the races is cal- culated based on statics equilibrium to obtain the two kinds of nonlinear dynamic models in a rigid-rotor ball bearing system. Finally~ the displacement response and frequency spectrum for the two system models are compared quantitatively at different rotational speeds, and then the structures of frequency-amplitude curves over a wide speed range are compared qualitatively under different levels of radial clearance, amplitude of excitation, and mass of supporting rotor. The results demonstrate that the cubic polynomial model can take place of the Hertzian contact model in a range of deformation.展开更多
文摘The various factors influencing the settlement of composite foundation have been more completely studied through numerical simulation. The influence on the settlement of composite foundation of the geometry and mechanical properties of the pile, soil, cushion, and the interface between pile and soil have been investigated through computer simulation, in which the contact elements with zero thickness are used. Some valuable conclusions for the settlement of composite foundation have been obtained: (1) The method using the contact element of zero thickness is successful when used in the simulation of the settlement of composite foundation; (2) Among the factors influencing the settlement of composite foundation, the compression modulus of the soil is the largest, and the cohesion of the soil is the second largest; (3) The effects on settlement of the internal friction angle of the soil, the elastic modulus, the radius, and the length of the pile, and the elestic modulus of the cushion are also more obvious.
基金financially supported by the National Science Foundation of China-Yunnan United Foundation(No.U0837601)the National Natural Science Foundation of China(No.51267007)the Natural Science Foundation of Yunnan Province(No.2010CD126,No.2012FB195)
文摘A uniform transient temperature field model of electrical contacts operation was found by analyzing the process of closing arc constriction resistance Joule heat ~ breaking arc. Essential parameters of Ag/La2NiO4 electrical contact material for transient temperature field calculation were obtained through tests of electrical contact experimental instrument under 18 V DC in different cur- rents, other correlation experiments, and calculation anal- ysis. The finite element method was applied to solve the transient temperature field, and the features and distribution of the transient temperature field were obtained. The condition of material erosion and mass transfer can be forecasted by those calculation results. It is beneficial to research about the lifetime of Ag/La2NiO4 electrical material.
文摘An improved version of the regular boundary element method, the artificial boundary node approach, is derived. A simple contact algorithm is designed and implemented into the direct boundary element, regular boundary element and artificial boundary node approaches. The exisiting and derived approaches are tested using some case studies. The results of the artificial boundary node approach are compared with those of the existing boundary element program, the regular element approach, ANSYS and analytical solution whenever possible. The results show the effectiveness of the artificial boundary node approach for a wider range of boundary offsets.
文摘The evolution of leakage is studied using detailed contact finite element analysis. The distribution of stress at the gasket is analyzed using a contact condition based on slide-line elements using ABAQUS, a commercial finite element code, Slide-line elements also take into account pressure penetration as contact that is lost between flange and gasket. Results are presented for a particular flange, a raised face flange sealed by a mild steel gasket. A comparison of the results from the gasket contact analysis and the contact conditions specified by the ASME Boiler and Pressure Vessel Code, Sections VIII, Division 1 shows that the conditions specified in the ASME Code predict leakage relatively accurately.
基金The project supported by the National Natural Science Foundation of China (19772025)
文摘A scheme of boundary element method for moving contact of two-dimensional elastic bodies using conforming discretization is presented. Both the displacement and the traction boundary conditions are satisfied on the contacting region in the sense of discretization. An algorithm to deal with the moving of the contact boundary on a larger possible contact region is presented. The algorithm is generalized to rolling contact problem as well. Some numerical examples of moving and rolling contact of 2D elastic bodies with or without friction, including the bodies with a hole-type defect, are given to show the effectiveness and the accuracy of the presented schemes.
文摘Each surface of roughness has different shape of asperity which is modeled with various shapes of analytical models. In this paper, the differences among various models of shape of asperity investigate using the Finite Element Method (FEM) and various analytical models. The contact stresses in rough surfaces are calculated analytically using various asperity shape models. Finite element analysis is also carried out assuming three types of material properties namely, the linear, the elastic-perfect plastic and the elastic-nonlinear hardening. The analytical results are compared with the results obtained by the finite element method. The results illustrate for using a deterministic approach which the numerical models are suitable. In hertz model, the result of force is very big in interface of causing deformation plastic, while Model Zhao has almost same result with FEM nonlinear property model. It is observed that the results obtained from Zhao’s model are generally in a better agreement with the results obtained from various finite element models especially in elastic-plastic and plastic zones, hence it may be concluded that Zhao’s model can be used for analyzing the rough surfaces in contact mechanics.
文摘Ceramics are good alternative to metal as bearing couple materials because of their better wear resistance. A Finite Element(FE) study was performed to investigate the contact mechanics and stress distribution of Ceramic-on-Ceramic (COC) hip resurfacingprostheses. It was focused in particular on a parametric study to examine the effects of radial clearance, loading,alumina coating on the implants, bone quality, and fixation of cup-bone interface. It was found that a reduction in the radialclearance had the most significant effect on the predicted contact pressure distribution among all of the parameters considered inthis study. It was determined that there was a significant influence of non-metallic materials, such as the bone underneath thebearing components, on the predicted contact mechanics. Stress shielding within the bone tissue was found to be a major concernwhen regarding the use of ceramic as an alternative to metallic resurfacing prostheses. Therefore, using alumina implantswith a metal backing was found to be the best design for ceramic resurfacing prostheses in this study. The loading, bone quality,and acetabular cup fixation conditions were found to have only minor effects on the predicted contact pressure distribution alongthe bearing surfaces.
文摘Vertical rigidity of the space self adaptive 530 high rigidity mill is calculated by applying the boundary element method (BEM) of three dimension elastic contact problem,which can update the existed deforming separation calculating theory and corresponding methods of material mechanics,elastic mechanics and finite element method.The method has less hypotheses and stronger synthesis in contact type calculating model.The advantages of the method are high calculating rate,high calculating accuracy,etc..
文摘The discrete element method(DEM) has been extensively adopted to investigate many complex geotechnical related problems due to its capability to incorporate the discontinuous nature of granular materials. In particular, when simulating large deformations or distortion of soil(e.g. cavity expansion),DEM can be very effective as other numerical solutions may experience convergence problems. Cavity expansion theory has widespread applications in geotechnical engineering, particularly to the problems concerning in situ testing, pile installation and so forth. In addition, the behaviour of geomaterials in a macro-level is utterly determined by microscopic properties, highlighting the importance of contact models. Despite the fact that there are numerous contact models proposed to mimic the realistic behaviour of granular materials, there are lack of studies on the effects of these contact models on the soil response.Hence, in this study, a series of three-dimensional numerical simulations with different contact constitutive models was conducted to simulate the response of sandy soils during cylindrical cavity expansion. In this numerical investigation, three contact models, i.e. linear contact model, rolling resistance contact model,and Hertz contact model, are considered. It should be noted that the former two models are linear based models, providing linearly elastic and frictional plasticity behaviours, whereas the latter one consists of nonlinear formulation based on an approximation of the theory of Mindlin and Deresiewicz. To examine the effects of these contact models, several cylindrical cavities were created and expanded gradually from an initial radius of 0.055 m to a final radius of 0.1 m. The numerical predictions confirm that the calibrated contact models produced similar results regarding the variations of cavity pressure, radial stress, deviatoric stress, volumetric strain, as well as the soil radial displacement. However, the linear contact model may result in inaccurate predictions when highly angular soil particles are involved. In addition, considering the excessive soil displacement induced by the pile installation(i.e. cavity expansion), a minimum distance of11 a(a is the cavity radius) is recommend for practicing engineers to avoid the potential damages to the existing piles and adjacent structures.
文摘The thrust contact between the Togo and the Dahomeyan formations in Ghana is a lithotectonic boundary that exists between two major Precambrian formations which are of importance to geologists owing to the fact that Precambrian rocks in Ghana host almost all economic minerals and metals. The lithostratigraphy of the Togo-Dahomyan thrust contact rocks from a borehole in Kwabenya near Accra (the Capital of Ghana) has been studied and major crustal chemical elements assayed using Instrumental Neutron Activation Analysis (INAA) techniques. The results have revealed elemental compositions and the mineralogical make up of the lithostratigraphic units of these contact rocks and the general geology of the Togo and the Dahomeyan formations. The profile shows a thrust contact that exists between the Togo formations which are underlying the Dahomeyan formations. The Togo formations here are made mainly of quartzite in fresh schist and the Dahomeyan made of gneisses. In between these two major geologic formations are the rocks of the contact which are intercalation of quartzite and schist. The rocks are felsic with an average felsic index (F) of 85.74 and feldspar rich with K- (orthoclase) feldspars dominating the Dahomeyan rocks. Iron and titanium oxides are depleted with depth from 52 m depth below surface downward and potassium oxide was enriched with depth from 42 m below surface downward. Major mineral forming elements such as aluminum and calcium had varied levels in the Togo and in the Dahomeyan rocks.
文摘A simple interface element for analyzing contact friction problems is developed. Taking nodal displacements and contact stresses as unknowns, this element can simulate frictional slippage, decoupling and re-bonding of two bodies initially mating or having gaps at a common interface. The method is based on the Finite Element Method and load incremental theory. The geometric and static constraint conditions on contact surfaces are treated as additional conditions and are included in stiffness equations. This simple element has the advantages of easy implementation into standard finite element programs and fast speed for convergence as well as high accuracy for stress distribution in interface. Undesirable stress oscillations are also investigated whenever large stress gradients exist over the contact surfaces. Exact integration or the conventional Gauss integration scheme used to evaluate the interpolation function matrix of the interface element is found to be the source of the oscillations. Eigenmode analysis demonstrates that the stress behavior of an interface element can be improved by using the Newton-Cotes integration scheme. Finally, the test example of a strip footing problem is presented.
基金This work was financially supported by the National Nature Science Foundation of China(No.19902001), the National Excellent Yout
文摘An algorithm is presented for controlling two-dimensional motion contact bodies with conforming discretization. Since a kind of special boundary element is utilized in the algorithm, the displacement compatibility and traction equilibrium conditions at nodes can be satisfied simultaneously in arbitrary locations of the contact interface. In addition, a method is also proposed in which the contact boundary location can be moved flexibly on the possible contact boundary. This method is effective to deal with moving and rolling contact problems on a possible larger moving or rolling contact region. Numerical examples show effectiveness of the presented scheme.
文摘The adaptive element techniques of contact problem are studied by means of penalty method, and the error estimators are discussed. Based on error estimators, algorithm of the adaptive element techniques is developed, then the Gauss - Newton iterations are used which allow the nonlinear problem to be transformed into a sequence of linear sub- problems then easily solved. In addition, the algorithm can be applied into the simulation of de -bonding of fiber - reinforced composites.
文摘The boundary element method in framework is given to evaluate three dimensional frictional contact problems. Elasto plastic material behavior is taken into account by mean of an initial stress formulation and Von Mises yield criterion. The amount of tangential traction at contact surface is limited by Coulomb's friction law and constant shear rule. From some numerical results of a plate rolling problem, it is demonstrated here that the BEM can be used to efficiently and accurately analyze this class of forming problems.
文摘Through defining slide yield function and floating potential function of thermo-contact surface, the complementary equation of thermo-contact boundary has been reached, the fundamental equations to solve 3D thermo-contact coupled problem have been listed. On this foundation, the finite element equation and definite solution condition of contact heat transfer have been given out. Based on virtual work principle and contact element technology, the finite element equation of 3D elastic contact system has been deduced under the effect of thermal stress. The pseudo load brought by contact gap have been introduced into this equation in order to reflect the contact state change. During iteration, once contact rigidity matrix is formed, it won’t change, which will make calculation reduce greatly.
文摘The formulation of boundary element method for handling contact problems with friction and the technique for high speed contact analysis are presented. This formulation is based on the idea of modifying the length of contact elements without altering the total number of elements. The high precision of solution and high speed analysis are verified according to the results of conventional method and analysis method.
基金supported by the Innovation Plan for Postgraduate Students sponsored by the Education Department of Jiangsu Province,China (CX08B 107Z)
文摘A hydromechanical interface element is proposed for the consideration of the hydraulic-mechanical coupling effect along the interface.The fully coupled governing equations and the relevant finite element formulations are derived in detail for the interface element.All the involved matrices are of the same form as those of a solid element,which makes the incorporation of the model into a finite element program straightforward.Three examples are then numerically simulated using the interface element.Reasonable results confirm the correctness of the proposed model and motivate its application in hydromechanical contact problems in the future.
基金The project supported by the National Natural Science Foundation of China(59875045)and the State Key Laboratory of Automobile Safety and Energy Saving(K9705)
文摘A new contact searching algorithm for contact-impact systems is proposed in this paper.In terms of the cell structure and the linked-list,this algo- rithm solves the problem of sorting and searching contacts in three dimensions by transforming it to a retrieving process from two one-dimensional arrays,and binary searching is no longer required.Using this algorithm, the cost of contact searching is reduced to the order of O(N)instead of O(Nlog_2N)for traditional ones,where N is the node number in the system.Moreover,this algorithm can handle contact systems with arbitrary mesh layouts.Due to the simplicity of this algorithm it can be easily implemented in a dynamic explicit finite element program.Our numerical experi- mental result shows that this algorithm is reliable arid efficient for contact searching of three dimensional systems.
基金Project supported by the National Basic Research Program(973 Program)of China(No.2015CB057400)the National Natural Science Foundation of China(No.11602070)+1 种基金the China Postdoctoral Science Foundation(No.2016M590277)the Heilongjiang Postdoctoral Financial Assistance(No.LBH-Z16067)
文摘A new nonlinear force model based on experimental data is proposed to replace the classical Hertzian contact model to solve the fractional index nonlinearity in a ball bearing system. Firstly, the radial force and the radial deformation are measured by statics experiments, and the data are fitted respectively by using the Hertzian contact model and the cubic polynomial model. Then~ the two models are compared with the approximation formula appearing in Aeroengine Design Manual. In consequence, the two models are equivalent in an allowable deformation range. After that, the relationship of contact force and contact deformation for single rolling element between the races is cal- culated based on statics equilibrium to obtain the two kinds of nonlinear dynamic models in a rigid-rotor ball bearing system. Finally~ the displacement response and frequency spectrum for the two system models are compared quantitatively at different rotational speeds, and then the structures of frequency-amplitude curves over a wide speed range are compared qualitatively under different levels of radial clearance, amplitude of excitation, and mass of supporting rotor. The results demonstrate that the cubic polynomial model can take place of the Hertzian contact model in a range of deformation.