Ribonucleic Acid(RNA)contact prediction holds great significance for modeling RNA 3D structures and further understanding RNA biological functions.The rapid growth of RNA sequencing data has driven the development of ...Ribonucleic Acid(RNA)contact prediction holds great significance for modeling RNA 3D structures and further understanding RNA biological functions.The rapid growth of RNA sequencing data has driven the development of diverse computational methods for RNA contact prediction,and a benchmark evaluation of these methods remains essential.In this work,we first classified RNA contact prediction methods into statistical inference-based and neural networkbased ones.We then evaluated eight state-of-the-art methods on three test sets:a sequencediverse set,a structurally non-redundant set and a CASP RNA targets set.Our evaluation shows that for identifying non-local and long-range contacts,neural network-based methods outperform statistical inference-based ones,with SPOT-RNA-2D achieving the best performance,followed by CoCoNet and RNAcontact.However,for identifying the long-range tertiary contacts,which are vital for stabilizing RNA tertiary structure,statistical inference-based methods exhibit superior performance with GREMLIN emerging as the top performer.This work provides a comprehensive benchmarking of RNA contact prediction methods,highlighting their strengths and limitations to guide further methodological improvements and applications in RNA structure modeling.展开更多
This paper investigates the frictional adhesive contact of a rigid,electrically/magnetically conductive spherical indenter sliding past a multiferroic coating deposed onto a rigid substrate,based on the hybrid element...This paper investigates the frictional adhesive contact of a rigid,electrically/magnetically conductive spherical indenter sliding past a multiferroic coating deposed onto a rigid substrate,based on the hybrid element method.The adhesion behavior is described based on the Maugis-Dugdale model.The adhesion-driven conjugate gradient method is employed to calculate the distribution of unknown pressures,while the discrete convolution-fast Fourier transform is utilized to compute the deformations,surface electric and magnetic potentials as well as the subsurface stresses,electric displacements,and magnetic inductions.The goal of this study is to investigate the influences of adhesion parameter,friction coefficient,coating thickness,and surface electric and magnetic charge densities on contact behaviors,such as contact area and pressures,electric and magnetic potentials,and subsurface stresses.展开更多
This study aims to address the feasibility of planned islanding operation and to investigate the effect of unplanned islanding using the master-slave islanding method for controlling the distributed generation units d...This study aims to address the feasibility of planned islanding operation and to investigate the effect of unplanned islanding using the master-slave islanding method for controlling the distributed generation units during grid-connected and islanding operation. Neplan desktop power simulation tool was used for the modelling and simulation of a realistic MV network with four different distributed generation technologies (diesel, gas, hydro and wind) along with their excitation and governor control systems, while an exponential model was used to represent the loads in the network. The dynamic and steady state behavior of the four distributed generation technologies were investigated during grid-connected operation and two transition modes to the islanding situation, planned and unplanned. The obtained results that validated through various case studies have shown that a suitable planned islanding transition could provide support to critical loads at the event of electricity utility outages.展开更多
Contact bounce of relay, which is the main cause of electric abrasion and material erosion, is inevitable. By using the mode expansion form, the dynamic behavior of two different reed systems for aerospace relays is a...Contact bounce of relay, which is the main cause of electric abrasion and material erosion, is inevitable. By using the mode expansion form, the dynamic behavior of two different reed systems for aerospace relays is analyzed. The dynamic model uses Euler-Bernoulli beam theory for cantilever beam, in which the driving force (or driving moment) of the electromagnetic system is taken into account, and the contact force between moving contact and stationary contact is simulated by the Kelvin-Voigt vis-coelastic...展开更多
The impact dynamics of a flexible multibody system is investigated. By using a partition method, the system is divided into two parts, the local impact region and the region away from the impact. The two parts are con...The impact dynamics of a flexible multibody system is investigated. By using a partition method, the system is divided into two parts, the local impact region and the region away from the impact. The two parts are connected by specific boundary conditions, and the system after partition is equivalent to the original system. According to the rigid-flexible coupling dynamic theory of multibody system, system＇s rigid-flexible coupling dynamic equations without impact are derived. A local impulse method for establishing the initial impact conditions is proposed. It satisfies the compatibility con- ditions for contact constraints and the actual physical situation of the impact process of flexible bodies. Based on the contact constraint method, system＇s impact dynamic equa- tions are derived in a differential-algebraic form. The contact/separation criterion and the algorithm are given. An impact dynamic simulation is given. The results show that system＇s dynamic behaviors including the energy, the deformations, the displacements, and the impact force during the impact process change dramatically. The impact makes great effects on the global dynamics of the system during and after impact.展开更多
A kind of improved contact frictional model on basis of traditional Coulomb Friction model is adopted.Corresponding contact element is also given.The contact algorithm on basis of augmented Lagrange method is introduc...A kind of improved contact frictional model on basis of traditional Coulomb Friction model is adopted.Corresponding contact element is also given.The contact algorithm on basis of augmented Lagrange method is introduced and successfully applied to complex contact friction problem.Test example and actual engineering case all show that the algorithm of the model is efficient and computation results agree well with general rules.展开更多
This paper investigates the effects of charge parameters of the underwater contact explosion based on the axisymmetric smoothed particle hydrodynamics （SPH） method. The dynamic boundary particle is proposed to impro...This paper investigates the effects of charge parameters of the underwater contact explosion based on the axisymmetric smoothed particle hydrodynamics （SPH） method. The dynamic boundary particle is proposed to improve the pressure fluctuation and numerical accuracy near the symmetric axis. An in-depth study is carried out over the influence of charge shapes and detonation modes on the near-field loads in terms of the peak pressure and impulse of shock waves. For different charge shapes, the cylindrical charge with different length-diameter ratios may cause strong directivity of peak pressure and impulse in the near field. Compared with spherical charge, the peak pressure of cylindrical charge may be either weakened or enhanced in different directions. Within a certain range, the greater the length-diameter ratio is, the more obvious the effect will be. The weakened ratio near the detonation end may reach 25% approximately, while the enhanced ratio may reach around 20% in the opposite direction. However, the impulse in different directions seems to be uniform. For different detonation modes, compared with point-source explosion, the peak pressure of plane-source explosion is enhanced by about 5%. Besides, the impulse of plane-source explosion is enhanced by around 5% near the detonation end, but close to those of the point-source explosion in other directions. Based on the material constitutive relation in the axisymmetric coordinates, a simple case of underwater contact explosion is simulated to verify the above conclusions, showing that the charge parameters of underwater contact explosion should not be ignored.展开更多
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.展开更多
An experimental method of measuring the real contact area of transparent blocks based on the principle of total internal reflection is presented, intending to support the investigation of friction characteristics, hea...An experimental method of measuring the real contact area of transparent blocks based on the principle of total internal reflection is presented, intending to support the investigation of friction characteristics, heat conduction, and energy dissipation at the contact interface. A laser sheet illuminates the contact interface, and the transmitted laser sheet is projected onto a screen. Then the contact information is acquired from the screen by a camera. An improved Otsu method is proposed to process the data of experimental images. It can compute the threshold of the overall image and filter out all the pixels one by one. Through analyzing the experimental results, we describe the relationship between the real contact area and the positive pressure during a continuous loading process, at different loading rates, with the polymethyl methacrylate（PMMA）material. A hysteresis phenomenon in the relationship between the real contact area and the positive pressure is found and explained.展开更多
The simulation of real contact area between materials is foundationally important for the contact mechanics of mechanical structures. The Greenwood and Williamson(GW) model and the Majumdar(MB) model are the basic mod...The simulation of real contact area between materials is foundationally important for the contact mechanics of mechanical structures. The Greenwood and Williamson(GW) model and the Majumdar(MB) model are the basic models in this field, which are widely accepted and proven to be valid in many experiments and engineering. Although the contact models have evolved considerably in recent years, the verifications of the models are most based on the indirect methods such as electrical conductivity and contact stiffness, because of the lack of effective methods to directly measure the variation of contact surface. In this paper, the total reflection(TR) method is introduced into the verification of contact models.An experiment system based on TR method is constructed to measure the real contact area of two PMMA specimens. The comparison analysis between the results of experiment and models suggests that the experiment result has the same trend with simulation, the MB model has better agreement with the experimental result because this method can take into account the variation of radius and the merging of asperities, while the GW model has a huge deviation because of the dependence on resolution and the lack of considering the variation of radius and asperity's merging process. Taking the interaction of asperities into account could give a better result that is closer to the experiment. Our results and analysis prove that the experimental methods in this paper could be used as a more direct and valid method to quantitatively measure the real contact area and to verify the contact models.展开更多
Two modeling methods of the root insert for wind turbine blade are presented,i.e.,the local mesh optimization method(LMOM)and the global modeling method(GMM).Based on the optimized mesh of the local model for the meta...Two modeling methods of the root insert for wind turbine blade are presented,i.e.,the local mesh optimization method(LMOM)and the global modeling method(GMM).Based on the optimized mesh of the local model for the metal contact interface,LMOM is proposed to analyze the load path and stress distribution characteristics,while GMM is used to calculate and analyze the stress distribution characteristics of the resin layer established between the bushing and composite layers of root insert.To validate the GMM,a tension test is carried out.The result successfully shows that the shear strain expresses a similar strain distribution tendency with the GMM′s results.展开更多
There are two separate traditional approaches to model contact problems: continuum and atomistic theory. Continuum theory is successfully used in many domains, but when the scale of the model comes to nanometer, conti...There are two separate traditional approaches to model contact problems: continuum and atomistic theory. Continuum theory is successfully used in many domains, but when the scale of the model comes to nanometer, continuum approximation meets challenges. Atomistic theory can catch the detailed behaviors of an individual atom by using molecular dynamics (MD) or quantum mechanics, although accurately, it is usually time-consuming. A multiscale method coupled MD and finite element (FE) is presented. To mesh the FE region automatically, an adaptive method based on the strain energy gradient is introduced to the multiscale method to constitute an adaptive multiscale method. Utilizing the proposed method, adhesive contacts between a rigid cylinder and an elastic substrate are studied, and the results are compared with full MD simulations. The process of FE meshes refinement shows that adaptive multiscale method can make FE mesh generation more flexible. Comparison of the displacements of boundary atoms in the overlap region with the results from full MD simulations indicates that adaptive multiscale method can transfer displacements effectively. Displacements of atoms and FE nodes on the center line of the multiscale model agree well with that of atoms in full MD simulations, which shows the continuity in the overlap region. Furthermore, the Von Mises stress contours and contact force distributions in the contact region are almost same as full MD simulations. The method presented combines multiscale method and adaptive technique, and can provide a more effective way to multiscale method and to the investigation on nanoscale contact problems.展开更多
The discontinuous dynamical problem of multi-point contact and collision in multi-body system has always been a hot and difficult issue in this field.Based on the Gauss’principle of least constraint,a unified optimiz...The discontinuous dynamical problem of multi-point contact and collision in multi-body system has always been a hot and difficult issue in this field.Based on the Gauss’principle of least constraint,a unified optimization model for multibody system dynamics with multi-point contact and collision is established.The paper presents the study of the numerical solution scheme,in which particle swarm optimization method is used to deal with the corresponding optimization model.The article also presents the comparison of the Gauss optimization method(GOM)and the hybrid linear complementarity method(i.e.combining differential algebraic equations(DAEs)and linear complementarity problems(LCP)),commonly used to solve the dynamic contact problem of multibody systems with bilateral constraints.The results illustrate that,the GOM has the same advantage of dynamical modelling with LCP and when the redundant constraint exists,the GOM always has a unique solution and so no additional processing is needed,whereas the corresponding DAE-LCP method may have singular cases with multiple solutions or no solutions.Using numerical examples,the GOM is verified to effectively solve the dynamics of multibody systems with redundant unilateral and bilateral constraints without additional redundancy processing.The GOM can also be applied to the optimal control of systems in the future and combined with the parameter optimization of systems to handle dynamic problems.The work given provides the dynamics and control of the complex system with a new train of thought and method.展开更多
Several effective numerical methods for solving the elasto-plastic contact problems with friction are pres- ented.First,a direct substitution method is employed to impose the contact constraint conditions on condensed...Several effective numerical methods for solving the elasto-plastic contact problems with friction are pres- ented.First,a direct substitution method is employed to impose the contact constraint conditions on condensed finite ele- ment equations,thus resulting in a reduction by half in the dimension of final governing equations.Second,an algorithm composed of contact condition probes and elasto-plastic iterations is utilized to solve the governing equation,which distinguishes two kinds of nonlinearities,and makes the solution unique.In addition,Positive-Negative Sequence Modifica- tion Method is used to condense the finite element equations of each substructure and an analytical integration is intro- duced to determine the elasto-plastic status after each time step or each iteration,hence the computational efficiency is en- hanced to a great extent.Finally,several test and practical examples are pressented showing the validity and versatility of these methods and algorithms.展开更多
In this paper,an efficien formulation based on the Lagrangian method is presented to investigate the contact–impact problems of f exible multi-body systems.Generally,the penalty method and the Hertz contact law are t...In this paper,an efficien formulation based on the Lagrangian method is presented to investigate the contact–impact problems of f exible multi-body systems.Generally,the penalty method and the Hertz contact law are the most commonly used methods in engineering applications.However,these methods are highly dependent on various non-physical parameters,which have great effects on the simulation results.Moreover,a tremendous number of degrees of freedom in the contact–impact problems will influenc thenumericalefficien ysignificantl.Withtheconsideration of these two problems,a formulation combining the component mode synthesis method and the Lagrangian method is presented to investigate the contact–impact problems in fl xible multi-body system numerically.Meanwhile,the finit element meshing laws of the contact bodies will be studied preliminarily.A numerical example with experimental verificatio will certify the reliability of the presented formulationincontact–impactanalysis.Furthermore,aseries of numerical investigations explain how great the influenc of the finit element meshing has on the simulation results.Finally the limitations of the element size in different regions are summarized to satisfy both the accuracy and efficien y.展开更多
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..展开更多
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.展开更多
The spatial and temporal evolution of real contact area of contact interface with loads is a challenge.It is generally believed that there is a positive linear correlation between real contact area and normal load.How...The spatial and temporal evolution of real contact area of contact interface with loads is a challenge.It is generally believed that there is a positive linear correlation between real contact area and normal load.However,with the development of measuring instruments and methods,some scholars have found that the growth rate of real contact area will slow down with the increase of normal load under certain conditions,such as large-scale interface contact with small roughness surface,which is called the nonlinear phenomenon of real contact area.At present,there is no unified conclusion on the explanation of this phenomenon.We set up an experimental apparatus based on the total reflection principle to verify this phenomenon and analyze its mechanism.An image processing method is proposed,which can be used to quantitative analysis micro contact behaviors on macro contact phenomenon.The weighted superposition method is used to identify micro contact spots,to calculate the real contact area,and the color superimposed image is used to identify micro contact behaviors.Based on this method,the spatiotemporal evolution mechanism of real contact area nonlinear phenomena is quantitatively analyzed.Furthermore,the influence of nonlinear phenomenon of real contact area on the whole loading and unloading process is analyzed experimentally.It is found that the effects of fluid between contact interface,normal load amplitude and initial contact state on contact behavior cannot be ignored in large-scale interface contact with small roughness surface.展开更多
A comparison of direct integration methods is madeand their efficiency is investigated for impact problems.New-mark,Wilson-θ,Central Difference and Houbolt Methodsare used as direct integration methods.Impact analysi...A comparison of direct integration methods is madeand their efficiency is investigated for impact problems.New-mark,Wilson-θ,Central Difference and Houbolt Methodsare used as direct integration methods.Impact analysisincludes that of elastic and large deformation based uponupdated Lagrangian including buckling check.The resultsshow that the direct integration methods give differentresults in different contact-impact cases.展开更多
A hybrid method is established by combining photoelastic experiment and finite element analysis.The method is used to evaluate contact stress distribution on dry friction interfaces,such as the contact interfaces betw...A hybrid method is established by combining photoelastic experiment and finite element analysis.The method is used to evaluate contact stress distribution on dry friction interfaces,such as the contact interfaces between shrouds of fan blades and turbine blades.The photoelastic stress frozen experiment method is used to decide the displacement boundary conditions of numerical calculation.Higher accuracy and efficiency of solving problems are improved by the method.Technical difficulty and high cost of experiment are also avoided by the method.Good agreement of the stress distribution by using the hybrid method and experiment is obtained.展开更多
基金supported by grants from the National Natural Science Foundation of China(Grant No.12205223 to YLT,Grant No.12375038 to ZJT and Grant No.11605125 to YZS)the Department of Education of Hubei Province(Grant No.Q20221705 to YLT)。
文摘Ribonucleic Acid(RNA)contact prediction holds great significance for modeling RNA 3D structures and further understanding RNA biological functions.The rapid growth of RNA sequencing data has driven the development of diverse computational methods for RNA contact prediction,and a benchmark evaluation of these methods remains essential.In this work,we first classified RNA contact prediction methods into statistical inference-based and neural networkbased ones.We then evaluated eight state-of-the-art methods on three test sets:a sequencediverse set,a structurally non-redundant set and a CASP RNA targets set.Our evaluation shows that for identifying non-local and long-range contacts,neural network-based methods outperform statistical inference-based ones,with SPOT-RNA-2D achieving the best performance,followed by CoCoNet and RNAcontact.However,for identifying the long-range tertiary contacts,which are vital for stabilizing RNA tertiary structure,statistical inference-based methods exhibit superior performance with GREMLIN emerging as the top performer.This work provides a comprehensive benchmarking of RNA contact prediction methods,highlighting their strengths and limitations to guide further methodological improvements and applications in RNA structure modeling.
基金support from the National Natural Science Foundation of China(12102085)the Postdoctoral Science Foundation of China(2023M730504)+2 种基金the Sichuan Province Regional Innovation and Cooperation Project(2024YFHZ0210)supported by the European Union-NextGenerationEU through the Italian Ministry of University and Research under the following programs:(NM)PRIN2022(Projects of Relevant National Interest)grant no.2022SJ8HTC-Electroactive Gripper for Micro-Object Manipulation(ELFIN)(NM)PRIN2022 PNRR(Projects of Relevant National Interest)grant no.P2022MAZHX-Tribological Modeling for Sustainable Design of Industrial Frictional Interfaces(TRIBOSCORE).
文摘This paper investigates the frictional adhesive contact of a rigid,electrically/magnetically conductive spherical indenter sliding past a multiferroic coating deposed onto a rigid substrate,based on the hybrid element method.The adhesion behavior is described based on the Maugis-Dugdale model.The adhesion-driven conjugate gradient method is employed to calculate the distribution of unknown pressures,while the discrete convolution-fast Fourier transform is utilized to compute the deformations,surface electric and magnetic potentials as well as the subsurface stresses,electric displacements,and magnetic inductions.The goal of this study is to investigate the influences of adhesion parameter,friction coefficient,coating thickness,and surface electric and magnetic charge densities on contact behaviors,such as contact area and pressures,electric and magnetic potentials,and subsurface stresses.
文摘This study aims to address the feasibility of planned islanding operation and to investigate the effect of unplanned islanding using the master-slave islanding method for controlling the distributed generation units during grid-connected and islanding operation. Neplan desktop power simulation tool was used for the modelling and simulation of a realistic MV network with four different distributed generation technologies (diesel, gas, hydro and wind) along with their excitation and governor control systems, while an exponential model was used to represent the loads in the network. The dynamic and steady state behavior of the four distributed generation technologies were investigated during grid-connected operation and two transition modes to the islanding situation, planned and unplanned. The obtained results that validated through various case studies have shown that a suitable planned islanding transition could provide support to critical loads at the event of electricity utility outages.
文摘Contact bounce of relay, which is the main cause of electric abrasion and material erosion, is inevitable. By using the mode expansion form, the dynamic behavior of two different reed systems for aerospace relays is analyzed. The dynamic model uses Euler-Bernoulli beam theory for cantilever beam, in which the driving force (or driving moment) of the electromagnetic system is taken into account, and the contact force between moving contact and stationary contact is simulated by the Kelvin-Voigt vis-coelastic...
基金supported by the National Natural Science Foundation of China(Nos.11132007,11272155,and 10772085)the Fundamental Research Funds for the Central Universities(No.30920130112009)the 333 Project of Jiangsu Province of China(No.BRA2011172)
文摘The impact dynamics of a flexible multibody system is investigated. By using a partition method, the system is divided into two parts, the local impact region and the region away from the impact. The two parts are connected by specific boundary conditions, and the system after partition is equivalent to the original system. According to the rigid-flexible coupling dynamic theory of multibody system, system＇s rigid-flexible coupling dynamic equations without impact are derived. A local impulse method for establishing the initial impact conditions is proposed. It satisfies the compatibility con- ditions for contact constraints and the actual physical situation of the impact process of flexible bodies. Based on the contact constraint method, system＇s impact dynamic equa- tions are derived in a differential-algebraic form. The contact/separation criterion and the algorithm are given. An impact dynamic simulation is given. The results show that system＇s dynamic behaviors including the energy, the deformations, the displacements, and the impact force during the impact process change dramatically. The impact makes great effects on the global dynamics of the system during and after impact.
文摘A kind of improved contact frictional model on basis of traditional Coulomb Friction model is adopted.Corresponding contact element is also given.The contact algorithm on basis of augmented Lagrange method is introduced and successfully applied to complex contact friction problem.Test example and actual engineering case all show that the algorithm of the model is efficient and computation results agree well with general rules.
基金supported by the National Natural Science Foundation of China(No.51379039)the Excellent Young Scientists Fund(No.51222904)
文摘This paper investigates the effects of charge parameters of the underwater contact explosion based on the axisymmetric smoothed particle hydrodynamics （SPH） method. The dynamic boundary particle is proposed to improve the pressure fluctuation and numerical accuracy near the symmetric axis. An in-depth study is carried out over the influence of charge shapes and detonation modes on the near-field loads in terms of the peak pressure and impulse of shock waves. For different charge shapes, the cylindrical charge with different length-diameter ratios may cause strong directivity of peak pressure and impulse in the near field. Compared with spherical charge, the peak pressure of cylindrical charge may be either weakened or enhanced in different directions. Within a certain range, the greater the length-diameter ratio is, the more obvious the effect will be. The weakened ratio near the detonation end may reach 25% approximately, while the enhanced ratio may reach around 20% in the opposite direction. However, the impulse in different directions seems to be uniform. For different detonation modes, compared with point-source explosion, the peak pressure of plane-source explosion is enhanced by about 5%. Besides, the impulse of plane-source explosion is enhanced by around 5% near the detonation end, but close to those of the point-source explosion in other directions. Based on the material constitutive relation in the axisymmetric coordinates, a simple case of underwater contact explosion is simulated to verify the above conclusions, showing that the charge parameters of underwater contact explosion should not be ignored.
基金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.
基金Project supported by the National Natural Science Foundation of China(Grant No.11272171)the Beijing Natural Science Foundation,China(Grant No.3132030)the Education Ministry Doctoral Fund of China(Grant No.20120002110070)
文摘An experimental method of measuring the real contact area of transparent blocks based on the principle of total internal reflection is presented, intending to support the investigation of friction characteristics, heat conduction, and energy dissipation at the contact interface. A laser sheet illuminates the contact interface, and the transmitted laser sheet is projected onto a screen. Then the contact information is acquired from the screen by a camera. An improved Otsu method is proposed to process the data of experimental images. It can compute the threshold of the overall image and filter out all the pixels one by one. Through analyzing the experimental results, we describe the relationship between the real contact area and the positive pressure during a continuous loading process, at different loading rates, with the polymethyl methacrylate（PMMA）material. A hysteresis phenomenon in the relationship between the real contact area and the positive pressure is found and explained.
基金Project supported by the National Natural Science Foundation of China(Grant No.11872033)the Beijing Natural Science Foundation,China(Grant No.3172017)
文摘The simulation of real contact area between materials is foundationally important for the contact mechanics of mechanical structures. The Greenwood and Williamson(GW) model and the Majumdar(MB) model are the basic models in this field, which are widely accepted and proven to be valid in many experiments and engineering. Although the contact models have evolved considerably in recent years, the verifications of the models are most based on the indirect methods such as electrical conductivity and contact stiffness, because of the lack of effective methods to directly measure the variation of contact surface. In this paper, the total reflection(TR) method is introduced into the verification of contact models.An experiment system based on TR method is constructed to measure the real contact area of two PMMA specimens. The comparison analysis between the results of experiment and models suggests that the experiment result has the same trend with simulation, the MB model has better agreement with the experimental result because this method can take into account the variation of radius and the merging of asperities, while the GW model has a huge deviation because of the dependence on resolution and the lack of considering the variation of radius and asperity's merging process. Taking the interaction of asperities into account could give a better result that is closer to the experiment. Our results and analysis prove that the experimental methods in this paper could be used as a more direct and valid method to quantitatively measure the real contact area and to verify the contact models.
基金supported jointly by the National Basic Research Program of China("973"Program)(No2014CB046200)the National Science Foundation of Jiangsu Province(No.BK2014059)+1 种基金the Priority Academic Program Development of Jiangsu Higher Education Institutionsthe National Natural Science Foundation of China(No.11172135)
文摘Two modeling methods of the root insert for wind turbine blade are presented,i.e.,the local mesh optimization method(LMOM)and the global modeling method(GMM).Based on the optimized mesh of the local model for the metal contact interface,LMOM is proposed to analyze the load path and stress distribution characteristics,while GMM is used to calculate and analyze the stress distribution characteristics of the resin layer established between the bushing and composite layers of root insert.To validate the GMM,a tension test is carried out.The result successfully shows that the shear strain expresses a similar strain distribution tendency with the GMM′s results.
基金supported by National Natural Science Foundation of China (Grant Nos. 51205313, 50975232)Northwestern Polytechnical University Foundation for Fundamental Research of China (Grant No.JC20110249)
文摘There are two separate traditional approaches to model contact problems: continuum and atomistic theory. Continuum theory is successfully used in many domains, but when the scale of the model comes to nanometer, continuum approximation meets challenges. Atomistic theory can catch the detailed behaviors of an individual atom by using molecular dynamics (MD) or quantum mechanics, although accurately, it is usually time-consuming. A multiscale method coupled MD and finite element (FE) is presented. To mesh the FE region automatically, an adaptive method based on the strain energy gradient is introduced to the multiscale method to constitute an adaptive multiscale method. Utilizing the proposed method, adhesive contacts between a rigid cylinder and an elastic substrate are studied, and the results are compared with full MD simulations. The process of FE meshes refinement shows that adaptive multiscale method can make FE mesh generation more flexible. Comparison of the displacements of boundary atoms in the overlap region with the results from full MD simulations indicates that adaptive multiscale method can transfer displacements effectively. Displacements of atoms and FE nodes on the center line of the multiscale model agree well with that of atoms in full MD simulations, which shows the continuity in the overlap region. Furthermore, the Von Mises stress contours and contact force distributions in the contact region are almost same as full MD simulations. The method presented combines multiscale method and adaptive technique, and can provide a more effective way to multiscale method and to the investigation on nanoscale contact problems.
基金This study was funded by the National Natural Science Foundation of China(Grant 11272167).
文摘The discontinuous dynamical problem of multi-point contact and collision in multi-body system has always been a hot and difficult issue in this field.Based on the Gauss’principle of least constraint,a unified optimization model for multibody system dynamics with multi-point contact and collision is established.The paper presents the study of the numerical solution scheme,in which particle swarm optimization method is used to deal with the corresponding optimization model.The article also presents the comparison of the Gauss optimization method(GOM)and the hybrid linear complementarity method(i.e.combining differential algebraic equations(DAEs)and linear complementarity problems(LCP)),commonly used to solve the dynamic contact problem of multibody systems with bilateral constraints.The results illustrate that,the GOM has the same advantage of dynamical modelling with LCP and when the redundant constraint exists,the GOM always has a unique solution and so no additional processing is needed,whereas the corresponding DAE-LCP method may have singular cases with multiple solutions or no solutions.Using numerical examples,the GOM is verified to effectively solve the dynamics of multibody systems with redundant unilateral and bilateral constraints without additional redundancy processing.The GOM can also be applied to the optimal control of systems in the future and combined with the parameter optimization of systems to handle dynamic problems.The work given provides the dynamics and control of the complex system with a new train of thought and method.
基金The Project Supported by National Natural Science Foundation of China
文摘Several effective numerical methods for solving the elasto-plastic contact problems with friction are pres- ented.First,a direct substitution method is employed to impose the contact constraint conditions on condensed finite ele- ment equations,thus resulting in a reduction by half in the dimension of final governing equations.Second,an algorithm composed of contact condition probes and elasto-plastic iterations is utilized to solve the governing equation,which distinguishes two kinds of nonlinearities,and makes the solution unique.In addition,Positive-Negative Sequence Modifica- tion Method is used to condense the finite element equations of each substructure and an analytical integration is intro- duced to determine the elasto-plastic status after each time step or each iteration,hence the computational efficiency is en- hanced to a great extent.Finally,several test and practical examples are pressented showing the validity and versatility of these methods and algorithms.
基金supported by the National Science Foundation of China (Grants 11132007,11272203)
文摘In this paper,an efficien formulation based on the Lagrangian method is presented to investigate the contact–impact problems of f exible multi-body systems.Generally,the penalty method and the Hertz contact law are the most commonly used methods in engineering applications.However,these methods are highly dependent on various non-physical parameters,which have great effects on the simulation results.Moreover,a tremendous number of degrees of freedom in the contact–impact problems will influenc thenumericalefficien ysignificantl.Withtheconsideration of these two problems,a formulation combining the component mode synthesis method and the Lagrangian method is presented to investigate the contact–impact problems in fl xible multi-body system numerically.Meanwhile,the finit element meshing laws of the contact bodies will be studied preliminarily.A numerical example with experimental verificatio will certify the reliability of the presented formulationincontact–impactanalysis.Furthermore,aseries of numerical investigations explain how great the influenc of the finit element meshing has on the simulation results.Finally the limitations of the element size in different regions are summarized to satisfy both the accuracy and efficien y.
文摘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..
基金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 Foundation of China(Grant No.11872033)the Beijing Natural Science Foundation,China(Grant No.3172017).
文摘The spatial and temporal evolution of real contact area of contact interface with loads is a challenge.It is generally believed that there is a positive linear correlation between real contact area and normal load.However,with the development of measuring instruments and methods,some scholars have found that the growth rate of real contact area will slow down with the increase of normal load under certain conditions,such as large-scale interface contact with small roughness surface,which is called the nonlinear phenomenon of real contact area.At present,there is no unified conclusion on the explanation of this phenomenon.We set up an experimental apparatus based on the total reflection principle to verify this phenomenon and analyze its mechanism.An image processing method is proposed,which can be used to quantitative analysis micro contact behaviors on macro contact phenomenon.The weighted superposition method is used to identify micro contact spots,to calculate the real contact area,and the color superimposed image is used to identify micro contact behaviors.Based on this method,the spatiotemporal evolution mechanism of real contact area nonlinear phenomena is quantitatively analyzed.Furthermore,the influence of nonlinear phenomenon of real contact area on the whole loading and unloading process is analyzed experimentally.It is found that the effects of fluid between contact interface,normal load amplitude and initial contact state on contact behavior cannot be ignored in large-scale interface contact with small roughness surface.
文摘A comparison of direct integration methods is madeand their efficiency is investigated for impact problems.New-mark,Wilson-θ,Central Difference and Houbolt Methodsare used as direct integration methods.Impact analysisincludes that of elastic and large deformation based uponupdated Lagrangian including buckling check.The resultsshow that the direct integration methods give differentresults in different contact-impact cases.
基金Youth Foundation of Beijing Polytechnic University
文摘A hybrid method is established by combining photoelastic experiment and finite element analysis.The method is used to evaluate contact stress distribution on dry friction interfaces,such as the contact interfaces between shrouds of fan blades and turbine blades.The photoelastic stress frozen experiment method is used to decide the displacement boundary conditions of numerical calculation.Higher accuracy and efficiency of solving problems are improved by the method.Technical difficulty and high cost of experiment are also avoided by the method.Good agreement of the stress distribution by using the hybrid method and experiment is obtained.
