Based on existing algorithms, a newly developed contact search algorithm is proposed. The new algorithm consists of global search, local searching, local tracking and penetration calculation processes. It requires no ...Based on existing algorithms, a newly developed contact search algorithm is proposed. The new algorithm consists of global search, local searching, local tracking and penetration calculation processes. It requires no iteration steps. It can deal with not only general tool surfaces with vertical walls, but also tool surfaces meshed with elements having very poor aspect ratios. It is demonstrated that the FE code employing this new contact search algorithm becomes more reliable, efficient and accurate for sheet metal forming simulation than conventional ones.展开更多
A new node-pairs contact algorithm is proposed to deal with a composite material or bi-material interface crack face contact and friction problem （e.g., resistant coating and thermal barrier coatings） subjected to c...A new node-pairs contact algorithm is proposed to deal with a composite material or bi-material interface crack face contact and friction problem （e.g., resistant coating and thermal barrier coatings） subjected to complicated load conditions. To decrease the calculation scale and calculation errors, the local Lagrange multipliers are solved only on a pair of contact nodes using the Jacobi iteration method, and the constraint modification of the tangential multipliers are required. After the calculation of the present node-pairs Lagrange multiplier, it is turned to next contact node-pairs until all node-pairs have finished. Compared with an ordinary contact algorithm, the new local node-pairs contact algorithm is allowed a more precise element on the contact face without the stiffness matrix singularity. The stress intensity factors （SIFs） and the contact region of an infinite plate central crack are calculated and show good agreement with those in the literature. The contact zone near the crack tip as well as its influence on singularity of stress fields are studied. Furthermore, the frictional contacts are also considered and found to have a significant influence on the SIFs. The normalized mode-II stress intensity factors KII for the friction coefficient decrease by 16% when f changes from 1 to 0.展开更多
The normal stress distribution of the contact patch is the base for wheel/rail dynamic analysis as well as the wheel wear calculation.When the wheel polygonalization appears on the worn-type wheel tread of the high-sp...The normal stress distribution of the contact patch is the base for wheel/rail dynamic analysis as well as the wheel wear calculation.When the wheel polygonalization appears on the worn-type wheel tread of the high-speed trains,the lateral and longitudinal relative wheel/rail curvature around the contact patch will vary with the change of wheel circumferential position.In such a case,the commonly used Hertzian contact method cannot meet the calculation precision,while even the non-Hertzian method has relatively high precision,but its calculation efficiency cannot meet the demand of engineering application.A non-Hertzian method was introduced to analyze the influence of the wheel polygonalization on the contact patch.By studying the shape and the pressure distribution of the contact patches at different contact positions,it is found that when the wavelength of the wheel polygonalization is larger than 65 mm,the results of the improved semiHertzian method can agree well with the non-Hertzian method with reasonable computing time.Besides,to ensure that the calculation result can meet the accuracy requirement when the wheelset yaw motion occurs,the trace method is used to improve the semi-Hertzian method and verified by the contact method.展开更多
Peridynamics(PD)is an effective method for simulating the spontaneous initiation and propagation of tensile cracks in materials.However,it faces great challenges in simulating compression-shear cracking of geomaterial...Peridynamics(PD)is an effective method for simulating the spontaneous initiation and propagation of tensile cracks in materials.However,it faces great challenges in simulating compression-shear cracking of geomaterials due to the lack of efficient contact-friction models.This paper introduces an original contact-friction model that leverages twin mesh and potential function principles within PD to model rock cracking under tensile and compressive stresses.The contact detection algorithm,based on space segmentation axis-aligned bounding box(AABB)tree data structure,is used to address the significant challenge of highly efficient contact detection in compression and shear problems.In this method,the twin mesh and potential function are utilized to quantify contact detection and contact degree,as well as friction behavior.This is in contrast to the distance and circular contact area model,which lacks physical significance in the classical PD method.As demonstrated by the tests on specimens containing cracks,the proposed model can capture 8 types of secondary fractures,reduce the contact detection error by about 29%e56%,and increase the contact retrieval efficiency by over 1600 times compared to the classic PD models.This significantly enhances the capability of PD to simulate the initiation,expansion,and coalescence of intricate compression-shear cracks.展开更多
As the performance of the box-type multiple launch rocket system(BMLRS)improves,its mechanical structures,particularly the plane clearance design between the slider on the rocket and the guide inside the launch canist...As the performance of the box-type multiple launch rocket system(BMLRS)improves,its mechanical structures,particularly the plane clearance design between the slider on the rocket and the guide inside the launch canister,have grown increasingly complex.However,deficiencies still exist in the current launch modeling theory for BMLRS.In this study,a multi-rigid-flexible-body launch dynamics model coupling the launch platform and rocket was established using the multibody system transfer matrix method and the Newton-Euler formulation.Furthermore,considering the bending of the launch canister,a detection algorithm for slider-guide plane clearance contact was proposed.To quantify the contact force and friction effect between the slider and guide,the contact force model and modified Coulomb model were introduced.Both the modal and launch tests were conducted.Additionally,the modal convergence was verified.By comparing the modal experiments and simulation results,the maximum relative error of the eigenfrequency is 3.29%.thereby verifying the accuracy of the developed BMLRS dynamics model.Furthermore,the launch test validated the proposed plane clearance contact model.Moreover,the study investigated the influence of various model parameters on the dynamic characteristics of BMLRS,including launch canister bending stiffness,slider and guide material,slider-guide clearance,slider length and layout.This analysis of influencing factors provides a foundation for future optimization in BMLRS design.展开更多
Large deformation contact problems generally involve highly nonlinear behaviors,which are very time-consuming and may lead to convergence issues.The finite particle method(FPM)effectively separates pure deformation fr...Large deformation contact problems generally involve highly nonlinear behaviors,which are very time-consuming and may lead to convergence issues.The finite particle method(FPM)effectively separates pure deformation from total motion in large deformation problems.In addition,the decoupled procedures of the FPM make it suitable for parallel computing,which may provide an approach to solve time-consuming issues.In this study,a graphics processing unit(GPU)-based parallel algorithm is proposed for two-dimensional large deformation contact problems.The fundamentals of the FPM for planar solids are first briefly introduced,including the equations of motion of particles and the internal forces of quadrilateral elements.Subsequently,a linked-list data structure suitable for parallel processing is built,and parallel global and local search algorithms are presented for contact detection.The contact forces are then derived and directly exerted on particles.The proposed method is implemented with main solution procedures executed in parallel on a GPU.Two verification problems comprising large deformation frictional contacts are presented,and the accuracy of the proposed algorithm is validated.Furthermore,the algorithm’s performance is investigated via a large-scale contact problem,and the maximum speedups of total computational time and contact calculation reach 28.5 and 77.4,respectively,relative to commercial finite element software Abaqus/Explicit running on a single-core central processing unit(CPU).The contact calculation time percentage of the total calculation time is only 18%with the FPM,much smaller than that(50%)with Abaqus/Explicit,demonstrating the efficiency of the proposed method.展开更多
In this paper, we develop a general way to construct contact algorithms for contact dynamical systems. Such an algorithm requires the corresponding step transition map preserve the contact structure of the underlying ...In this paper, we develop a general way to construct contact algorithms for contact dynamical systems. Such an algorithm requires the corresponding step transition map preserve the contact structure of the underlying contact phase space. The constructions are based on the correspondence between the contact geometry of R2n+1 and the conic symplectic one of R2n+2 and therefore, the algorithms are derived naturally from the symplectic algorithms of Hamiltonian systems.展开更多
The grid-based multi-velocity field technique has become increasingly popular for simulating the Material Point Method(MPM)in contact problems.However,this traditional technique has some shortcomings,such as(1)early c...The grid-based multi-velocity field technique has become increasingly popular for simulating the Material Point Method(MPM)in contact problems.However,this traditional technique has some shortcomings,such as(1)early contact and contact penetration can occur when the contact conditions are unsuitable,and(2)the method is not available for contact problems involving rigid-nonrigid materials,which can cause numerical instability.This study presents a new hybrid contact approach for the MPM to address these limitations to simulate the soil and structure interactions.The approach combines the advantages of point-point and point-segment contacts to implement contact detection,satisfying the impenetrability condition and smoothing the corner contact problem.The proposed approach is first validated through a disk test on an inclined slope.Then,several typical cases,such as granular collapse,bearing capacity,and deformation of a flexible retaining wall,are simulated to demonstrate the robustness of the proposed approach compared with FEM or analytical solutions.Finally,the proposed method is used to simulate the impact of sand flow on a deformable structure.The results show that the proposed contact approach can well describe the phenomenon of soil-structure interaction problems.展开更多
Although the penalty algorithm is simple and direct in concept, it has a defect that the contact forces are badly dependent on the chosen penalty factor. An improved contact-impact algorithm for the explicit integrati...Although the penalty algorithm is simple and direct in concept, it has a defect that the contact forces are badly dependent on the chosen penalty factor. An improved contact-impact algorithm for the explicit integration FEM is proposed in the present paper. Based on the fact that bodies cannot penetrate into each other on the contact faces, a set of equations with the additional unknown contact forces on the slave nodes can be formed in a new system configuration. By solving these equations, the correct contact forces could be obtained without using the penalty factor.展开更多
Basetl on the finite element solution of the parametric varialional principle of elastic con/del problem, a corresponding parallel algorithm has been created bv utilizing the specialities of parallel computer and the ...Basetl on the finite element solution of the parametric varialional principle of elastic con/del problem, a corresponding parallel algorithm has been created bv utilizing the specialities of parallel computer and the architecture of concurrent processing in this paper. In this algorithm. the parallelisms have heen realized in the processes of creation and assembly of stiffness matrix, of the static condensation, of the solution of stresses and in many other aspects. The programme of this algorithm has been realized on ELXSI-6400 parallel computer of Xi'an Jiaotong University. The results of computation show that the computational time can be saved efficiently and it is an effective parallel algorithm for the analyses of contact problems.展开更多
A genetic algorithm （GA） was studied to simultaneously determine the thermal transport properties and the contact resistance of thin films deposited on a thick substrate. A pulsed photothermal reflectance （PPR） sy...A genetic algorithm （GA） was studied to simultaneously determine the thermal transport properties and the contact resistance of thin films deposited on a thick substrate. A pulsed photothermal reflectance （PPR） system was employed for the measurements. The GA was used to extract the thermal properties. Measurements were performed on SiO2 thin films of different thicknesses on silicon substrate. The results show that the GA accompanied with the PPR system is useful for the simultaneous determination of thermal properties of thin films on a substrate.展开更多
Three dimensional frictional contact problems are formulated as linear complementarity problems based on the parametric variational principle. Two aggregate-functionbased algorithms for solving complementarity problem...Three dimensional frictional contact problems are formulated as linear complementarity problems based on the parametric variational principle. Two aggregate-functionbased algorithms for solving complementarity problems are proposed. One is called the self-adjusting interior point algorithm, the other is called the aggregate function smoothing algorithm. Numerical experiment shows the efficiency of the proposed two algorithms.展开更多
基金the National Natural Science F oundation of China (5 9875 0 2 5 ) and Excellent Young Teacher Founda-tion of the Educational Departm ent of China
文摘Based on existing algorithms, a newly developed contact search algorithm is proposed. The new algorithm consists of global search, local searching, local tracking and penetration calculation processes. It requires no iteration steps. It can deal with not only general tool surfaces with vertical walls, but also tool surfaces meshed with elements having very poor aspect ratios. It is demonstrated that the FE code employing this new contact search algorithm becomes more reliable, efficient and accurate for sheet metal forming simulation than conventional ones.
基金supported by the National Basic Research Program of China(Grant No.2012CB026200)the National Natural Science Foundation of China(Grant No.50878048)
文摘A new node-pairs contact algorithm is proposed to deal with a composite material or bi-material interface crack face contact and friction problem （e.g., resistant coating and thermal barrier coatings） subjected to complicated load conditions. To decrease the calculation scale and calculation errors, the local Lagrange multipliers are solved only on a pair of contact nodes using the Jacobi iteration method, and the constraint modification of the tangential multipliers are required. After the calculation of the present node-pairs Lagrange multiplier, it is turned to next contact node-pairs until all node-pairs have finished. Compared with an ordinary contact algorithm, the new local node-pairs contact algorithm is allowed a more precise element on the contact face without the stiffness matrix singularity. The stress intensity factors （SIFs） and the contact region of an infinite plate central crack are calculated and show good agreement with those in the literature. The contact zone near the crack tip as well as its influence on singularity of stress fields are studied. Furthermore, the frictional contacts are also considered and found to have a significant influence on the SIFs. The normalized mode-II stress intensity factors KII for the friction coefficient decrease by 16% when f changes from 1 to 0.
基金the support of the National Natural Science Foundation of China through grant Nos.52102441 and U2034210Independent Research and Development Project of the State Key Laboratory of Traction Power through grant Nos.2019TPL-T18 and 2022TPL-T10Natural Science Foundation of Sichuan Province,through grant No.2025ZNSFSC0395。
文摘The normal stress distribution of the contact patch is the base for wheel/rail dynamic analysis as well as the wheel wear calculation.When the wheel polygonalization appears on the worn-type wheel tread of the high-speed trains,the lateral and longitudinal relative wheel/rail curvature around the contact patch will vary with the change of wheel circumferential position.In such a case,the commonly used Hertzian contact method cannot meet the calculation precision,while even the non-Hertzian method has relatively high precision,but its calculation efficiency cannot meet the demand of engineering application.A non-Hertzian method was introduced to analyze the influence of the wheel polygonalization on the contact patch.By studying the shape and the pressure distribution of the contact patches at different contact positions,it is found that when the wavelength of the wheel polygonalization is larger than 65 mm,the results of the improved semiHertzian method can agree well with the non-Hertzian method with reasonable computing time.Besides,to ensure that the calculation result can meet the accuracy requirement when the wheelset yaw motion occurs,the trace method is used to improve the semi-Hertzian method and verified by the contact method.
基金supported by the National Natural Science Foundation of China(Grant No.52278333)the China Scholarship Council(CSC)and the Science and Technology Department of Liaoning Province(Grant No.2024JH2/102500069).
文摘Peridynamics(PD)is an effective method for simulating the spontaneous initiation and propagation of tensile cracks in materials.However,it faces great challenges in simulating compression-shear cracking of geomaterials due to the lack of efficient contact-friction models.This paper introduces an original contact-friction model that leverages twin mesh and potential function principles within PD to model rock cracking under tensile and compressive stresses.The contact detection algorithm,based on space segmentation axis-aligned bounding box(AABB)tree data structure,is used to address the significant challenge of highly efficient contact detection in compression and shear problems.In this method,the twin mesh and potential function are utilized to quantify contact detection and contact degree,as well as friction behavior.This is in contrast to the distance and circular contact area model,which lacks physical significance in the classical PD method.As demonstrated by the tests on specimens containing cracks,the proposed model can capture 8 types of secondary fractures,reduce the contact detection error by about 29%e56%,and increase the contact retrieval efficiency by over 1600 times compared to the classic PD models.This significantly enhances the capability of PD to simulate the initiation,expansion,and coalescence of intricate compression-shear cracks.
基金supported by National Natural Science Foundation of China(Grant No.92266201).
文摘As the performance of the box-type multiple launch rocket system(BMLRS)improves,its mechanical structures,particularly the plane clearance design between the slider on the rocket and the guide inside the launch canister,have grown increasingly complex.However,deficiencies still exist in the current launch modeling theory for BMLRS.In this study,a multi-rigid-flexible-body launch dynamics model coupling the launch platform and rocket was established using the multibody system transfer matrix method and the Newton-Euler formulation.Furthermore,considering the bending of the launch canister,a detection algorithm for slider-guide plane clearance contact was proposed.To quantify the contact force and friction effect between the slider and guide,the contact force model and modified Coulomb model were introduced.Both the modal and launch tests were conducted.Additionally,the modal convergence was verified.By comparing the modal experiments and simulation results,the maximum relative error of the eigenfrequency is 3.29%.thereby verifying the accuracy of the developed BMLRS dynamics model.Furthermore,the launch test validated the proposed plane clearance contact model.Moreover,the study investigated the influence of various model parameters on the dynamic characteristics of BMLRS,including launch canister bending stiffness,slider and guide material,slider-guide clearance,slider length and layout.This analysis of influencing factors provides a foundation for future optimization in BMLRS design.
基金This work was supported by the National Key Research and Development Program of China[Grant No.2016YFC0800200]the National Natural Science Foundation of China[Grant Nos.51778568,51908492,and 52008366]+1 种基金Zhejiang Provincial Natural Science Foundation of China[Grant Nos.LQ21E080019 and LY21E080022]This work was also sup-ported by the Key Laboratory of Space Structures of Zhejiang Province(Zhejiang University)and the Center for Balance Architecture of Zhejiang University.
文摘Large deformation contact problems generally involve highly nonlinear behaviors,which are very time-consuming and may lead to convergence issues.The finite particle method(FPM)effectively separates pure deformation from total motion in large deformation problems.In addition,the decoupled procedures of the FPM make it suitable for parallel computing,which may provide an approach to solve time-consuming issues.In this study,a graphics processing unit(GPU)-based parallel algorithm is proposed for two-dimensional large deformation contact problems.The fundamentals of the FPM for planar solids are first briefly introduced,including the equations of motion of particles and the internal forces of quadrilateral elements.Subsequently,a linked-list data structure suitable for parallel processing is built,and parallel global and local search algorithms are presented for contact detection.The contact forces are then derived and directly exerted on particles.The proposed method is implemented with main solution procedures executed in parallel on a GPU.Two verification problems comprising large deformation frictional contacts are presented,and the accuracy of the proposed algorithm is validated.Furthermore,the algorithm’s performance is investigated via a large-scale contact problem,and the maximum speedups of total computational time and contact calculation reach 28.5 and 77.4,respectively,relative to commercial finite element software Abaqus/Explicit running on a single-core central processing unit(CPU).The contact calculation time percentage of the total calculation time is only 18%with the FPM,much smaller than that(50%)with Abaqus/Explicit,demonstrating the efficiency of the proposed method.
文摘In this paper, we develop a general way to construct contact algorithms for contact dynamical systems. Such an algorithm requires the corresponding step transition map preserve the contact structure of the underlying contact phase space. The constructions are based on the correspondence between the contact geometry of R2n+1 and the conic symplectic one of R2n+2 and therefore, the algorithms are derived naturally from the symplectic algorithms of Hamiltonian systems.
基金funding support from the National Nature Science Foundation of China(Grant No.52022060)the Key Laboratory of Impact and Safety Engineering(Ningbo University).
文摘The grid-based multi-velocity field technique has become increasingly popular for simulating the Material Point Method(MPM)in contact problems.However,this traditional technique has some shortcomings,such as(1)early contact and contact penetration can occur when the contact conditions are unsuitable,and(2)the method is not available for contact problems involving rigid-nonrigid materials,which can cause numerical instability.This study presents a new hybrid contact approach for the MPM to address these limitations to simulate the soil and structure interactions.The approach combines the advantages of point-point and point-segment contacts to implement contact detection,satisfying the impenetrability condition and smoothing the corner contact problem.The proposed approach is first validated through a disk test on an inclined slope.Then,several typical cases,such as granular collapse,bearing capacity,and deformation of a flexible retaining wall,are simulated to demonstrate the robustness of the proposed approach compared with FEM or analytical solutions.Finally,the proposed method is used to simulate the impact of sand flow on a deformable structure.The results show that the proposed contact approach can well describe the phenomenon of soil-structure interaction problems.
文摘Although the penalty algorithm is simple and direct in concept, it has a defect that the contact forces are badly dependent on the chosen penalty factor. An improved contact-impact algorithm for the explicit integration FEM is proposed in the present paper. Based on the fact that bodies cannot penetrate into each other on the contact faces, a set of equations with the additional unknown contact forces on the slave nodes can be formed in a new system configuration. By solving these equations, the correct contact forces could be obtained without using the penalty factor.
基金Supported by the National Funds of National Structutal Vibration & Strength Laboratory of Xi'an Jiaotong University
文摘Basetl on the finite element solution of the parametric varialional principle of elastic con/del problem, a corresponding parallel algorithm has been created bv utilizing the specialities of parallel computer and the architecture of concurrent processing in this paper. In this algorithm. the parallelisms have heen realized in the processes of creation and assembly of stiffness matrix, of the static condensation, of the solution of stresses and in many other aspects. The programme of this algorithm has been realized on ELXSI-6400 parallel computer of Xi'an Jiaotong University. The results of computation show that the computational time can be saved efficiently and it is an effective parallel algorithm for the analyses of contact problems.
基金the National Natural Science Foundation of China under Grant Nos. 59995550-5 , 90207003.
文摘A genetic algorithm （GA） was studied to simultaneously determine the thermal transport properties and the contact resistance of thin films deposited on a thick substrate. A pulsed photothermal reflectance （PPR） system was employed for the measurements. The GA was used to extract the thermal properties. Measurements were performed on SiO2 thin films of different thicknesses on silicon substrate. The results show that the GA accompanied with the PPR system is useful for the simultaneous determination of thermal properties of thin films on a substrate.
基金The project supported by the National Natural Science foundation of china(10225212,50178016.10302007)the National Kev Basic Research Special Foundation and the Ministry of Education of China
文摘Three dimensional frictional contact problems are formulated as linear complementarity problems based on the parametric variational principle. Two aggregate-functionbased algorithms for solving complementarity problems are proposed. One is called the self-adjusting interior point algorithm, the other is called the aggregate function smoothing algorithm. Numerical experiment shows the efficiency of the proposed two algorithms.
