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.展开更多
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.展开更多
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.展开更多
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 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.展开更多
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.展开更多
A new multi-mesh contact algorithm for three-dimensional material point method is presented. The contact algorithm faithfully recovers the opposite acting forces between colliding bodies. Collision procedures between ...A new multi-mesh contact algorithm for three-dimensional material point method is presented. The contact algorithm faithfully recovers the opposite acting forces between colliding bodies. Collision procedures between regular bodies and/or rigid bodies are treated within the same framework. Multi-value of momentum and mass are defined on every node to describe the contact/sliding/separation procedure. Both normal and tangential velocities of each particle at the contact surface are calculated in respective individual mesh. A Coulomb friction is applied to describe the sliding or slipping between the contacting bodies. The efficiency of the contact algorithm is linearly related to the number of the contacting bodies because the overlapped nodes are labeled by sweeping the material particles of all bodies when the nodal momentum and mass are formed at every time step. Numerical simulation shows that our contact algorithm possesses high accuracy and low numerical energy dissipation, which is very important for solving collision problems.展开更多
The seismic stability of a cracked dam was examined in this study. Geometric nonlinearity and large deformations, as well as the contact condition at the crack site, were taken into consideration. The location of pene...The seismic stability of a cracked dam was examined in this study. Geometric nonlinearity and large deformations, as well as the contact condition at the crack site, were taken into consideration. The location of penetrated cracks was first identified using the concrete plastic-damage model based on the nonlinear finite element method （FEM）. Then, the hard contact algorithm was used to simulate the crack interaction in the normal direction, and the Coloumb friction model was used to simulate the crack interaction in the tangential direction. After verification of numerical models through a case study, the seismic stability of the Koyna Dam with two types of penetrated cracks is discussed in detail with different seismic peak accelerations, and the collapse processes of the cracked dam are also presented. The results show that the stability of the dam with two types of penetrated cracks can be ensured in an earthquake with a magnitude of the original Koyna earthquake, and the cracked dam has a large earthquake-resistant margin. The failure processes of the cracked dam in strong earthquakes can be divided into two stages： the sliding stage and the overturning stage. The sliding stage ends near the peak acceleration, and the top block slides a long distance along the crack before the collapse occurs. The maximum sliding displacement of the top block will decrease with an increasing friction coefficient at the crack site.展开更多
In this paper,the evaluations of metal ablation processes under high temperature,i.e.,the Al plate ablated by a laser and a heat carrier and the reactor pressure vessel ablated by a core melt,are studied by a novel pe...In this paper,the evaluations of metal ablation processes under high temperature,i.e.,the Al plate ablated by a laser and a heat carrier and the reactor pressure vessel ablated by a core melt,are studied by a novel peridynamic method.Above all,the peridynamic formulation for the heat conduction problem is obtained by Taylor’s expansion technique.Then,a simple and efficient moving boundary model in the peridynamic framework is proposed to handle the variable geometries,in which the ablated states of material points are described by an additional scalar field.Next,due to the automatic non-interpenetration properties of peridynamic method,a contact algorithm is established to determine the contact relationship between the ablated system and the additional heat carrier.In addition,the corresponding computational procedure is listed in detail.Finally,several numerical examples are carried out and the results verify the validity and accuracy of the present method.展开更多
A quasi-three-dimensional shell element model, which can beeffectively used to simulate the flanging and spring-backdeformation, is introduced into the independently developed CAEsoftware, KMAS, In this model, a doubl...A quasi-three-dimensional shell element model, which can beeffectively used to simulate the flanging and spring-backdeformation, is introduced into the independently developed CAEsoftware, KMAS, In this model, a double surface contact algorithm,which allows the gap between punch and die to change, and a spring-back treatment scheme based on finite element meshing are described.And then the flanging and spring-back deformations of the retractor'skickstand of a railcar made of stamped thick metal plate arenumerically simulated. The simulation results of flanging deformationare compared with those of international commercial software,PAM-STAMP, and experimental ones. Finally, a predicting scheme ofspring-back quantily for this problem is given.展开更多
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.展开更多
基金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.
基金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.
基金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.
基金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.
基金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.
基金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.
基金The project supported by the Science Foundation of Laboratory of Computational Physics,Science Foundation of China Academy of Engineering Physics,and National Natural Science Foundation of China under Grant Nos.10702010,10775018,10472052,and 10604010
文摘A new multi-mesh contact algorithm for three-dimensional material point method is presented. The contact algorithm faithfully recovers the opposite acting forces between colliding bodies. Collision procedures between regular bodies and/or rigid bodies are treated within the same framework. Multi-value of momentum and mass are defined on every node to describe the contact/sliding/separation procedure. Both normal and tangential velocities of each particle at the contact surface are calculated in respective individual mesh. A Coulomb friction is applied to describe the sliding or slipping between the contacting bodies. The efficiency of the contact algorithm is linearly related to the number of the contacting bodies because the overlapped nodes are labeled by sweeping the material particles of all bodies when the nodal momentum and mass are formed at every time step. Numerical simulation shows that our contact algorithm possesses high accuracy and low numerical energy dissipation, which is very important for solving collision problems.
基金supported by the National Basic Research Program of China (973 Program, Grant No. 2007CB714104)the National Natural Science Foundation of China (Grant No. 50779011)the Innovative Project for Graduate Students of Jiangsu Province (Grant No. CX10B_202Z)
文摘The seismic stability of a cracked dam was examined in this study. Geometric nonlinearity and large deformations, as well as the contact condition at the crack site, were taken into consideration. The location of penetrated cracks was first identified using the concrete plastic-damage model based on the nonlinear finite element method （FEM）. Then, the hard contact algorithm was used to simulate the crack interaction in the normal direction, and the Coloumb friction model was used to simulate the crack interaction in the tangential direction. After verification of numerical models through a case study, the seismic stability of the Koyna Dam with two types of penetrated cracks is discussed in detail with different seismic peak accelerations, and the collapse processes of the cracked dam are also presented. The results show that the stability of the dam with two types of penetrated cracks can be ensured in an earthquake with a magnitude of the original Koyna earthquake, and the cracked dam has a large earthquake-resistant margin. The failure processes of the cracked dam in strong earthquakes can be divided into two stages： the sliding stage and the overturning stage. The sliding stage ends near the peak acceleration, and the top block slides a long distance along the crack before the collapse occurs. The maximum sliding displacement of the top block will decrease with an increasing friction coefficient at the crack site.
基金supported by the National Natural Science Foundation of China(No.12102416).
文摘In this paper,the evaluations of metal ablation processes under high temperature,i.e.,the Al plate ablated by a laser and a heat carrier and the reactor pressure vessel ablated by a core melt,are studied by a novel peridynamic method.Above all,the peridynamic formulation for the heat conduction problem is obtained by Taylor’s expansion technique.Then,a simple and efficient moving boundary model in the peridynamic framework is proposed to handle the variable geometries,in which the ablated states of material points are described by an additional scalar field.Next,due to the automatic non-interpenetration properties of peridynamic method,a contact algorithm is established to determine the contact relationship between the ablated system and the additional heat carrier.In addition,the corresponding computational procedure is listed in detail.Finally,several numerical examples are carried out and the results verify the validity and accuracy of the present method.
基金the National Natural Science Foundation of China(No.19832020)the Ministry of Education of China.
文摘A quasi-three-dimensional shell element model, which can beeffectively used to simulate the flanging and spring-backdeformation, is introduced into the independently developed CAEsoftware, KMAS, In this model, a double surface contact algorithm,which allows the gap between punch and die to change, and a spring-back treatment scheme based on finite element meshing are described.And then the flanging and spring-back deformations of the retractor'skickstand of a railcar made of stamped thick metal plate arenumerically simulated. The simulation results of flanging deformationare compared with those of international commercial software,PAM-STAMP, and experimental ones. Finally, a predicting scheme ofspring-back quantily for this problem is given.
文摘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.
