Under the frame of multibody dynamics, the contact dynamics of elasto-plastic spatial thin beams is numerically studied by using the spatial thin beam elements of absolute nodal coordinate formulation（ANCF）. The int...Under the frame of multibody dynamics, the contact dynamics of elasto-plastic spatial thin beams is numerically studied by using the spatial thin beam elements of absolute nodal coordinate formulation（ANCF）. The internal force of the elasto-plastic spatial thin beam element is derived under the assumption that the plastic strain of the beam element depends only on its longitudinal deformation.A new body-fixed local coordinate system is introduced into the spatial thin beam element of ANCF for efficient contact detection in the contact dynamics simulation. The linear isotropic hardening constitutive law is used to describe the elasto-plastic deformation of beam material, and the classical return mapping algorithm is adopted to evaluate the plastic strains. A multi-zone contact approach of thin beams previously proposed by the authors is also introduced to detect the multiple contact zones of beams accurately, and the penalty method is used to compute the normal contact force of thin beams in contact. Four numerical examples are given to demonstrate the applicability and effectiveness of the proposed elasto-plastic spatial thin beam element of ANCF for flexible multibody system dynamics.展开更多
In several previous studies,it was reported that a supported pipe with small geometric imperfections would lose stability when the internal flow velocity became sufficiently high.Recently,however,it has become clear t...In several previous studies,it was reported that a supported pipe with small geometric imperfections would lose stability when the internal flow velocity became sufficiently high.Recently,however,it has become clear that this conclusion may be at best incomplete.A reevaluation of the problem is undertaken here by essentially considering the flow-induced static deformation of a pipe.With the aid of the absolute nodal coordinate formulation(ANCF)and the extended Lagrange equations for dynamical systems containing non-material volumes,the nonlinear governing equations of a pipe with three different geometric imperfections are introduced and formulated.Based on extensive numerical calculations,the static equilibrium configuration,the stability,and the nonlinear dynamics of the considered pipe system are determined and analyzed.The results show that for a supported pipe with the geometric imperfection of a half sinusoidal wave,the dynamical system could not lose stability even if the flow velocity reaches an extremely high value of 40.However,for a supported pipe with the geometric imperfection of one or one and a half sinusoidal waves,the first-mode buckling instability would take place at high flow velocity.Moreover,based on a further parametric analysis,the effects of the amplitude of the geometric imperfection and the aspect ratio of the pipe on the static deformation,the critical flow velocity for buckling instability,and the nonlinear responses of the supported pipes with geometric imperfections are analyzed.展开更多
As a subsequent work of previous studies of authors, a new parallel computation approach is proposed to simulate the coupled dynamics of a rigid-flexible multibody system and compressible fluid. In this approach, the ...As a subsequent work of previous studies of authors, a new parallel computation approach is proposed to simulate the coupled dynamics of a rigid-flexible multibody system and compressible fluid. In this approach, the smoothed particle hydrodynamics(SPH) method is used to model the compressible fluid, the natural coordinate formulation(NCF) and absolute nodal coordinate formulation(ANCF) are used to model the rigid and flexible bodies, respectively. In order to model the compressible fluid properly and efficiently via SPH method, three measures are taken as follows. The first is to use the Riemann solver to cope with the fluid compressibility, the second is to define virtual particles of SPH to model the dynamic interaction between the fluid and the multibody system, and the third is to impose the boundary conditions of periodical inflow and outflow to reduce the number of SPH particles involved in the computation process. Afterwards, a parallel computation strategy is proposed based on the graphics processing unit(GPU) to detect the neighboring SPH particles and to solve the dynamic equations of SPH particles in order to improve the computation efficiency. Meanwhile, the generalized-alpha algorithm is used to solve the dynamic equations of the multibody system. Finally, four case studies are given to validate the proposed parallel computation approach.展开更多
In the present study,the dynamics of the tendon system of a tension-leg platform(TLP)is investigated through the absolute nodal coordinate formulation(ANCF).Based on the energy conversion principle,the stiffness,gener...In the present study,the dynamics of the tendon system of a tension-leg platform(TLP)is investigated through the absolute nodal coordinate formulation(ANCF).Based on the energy conversion principle,the stiffness,generalized elastic force,external load and mass matrices of the element are deduced to perform the element assembling by using the finite element method.Then the motion equation of the tendon/riser is established.In this study,the TLP in the International Ship Structures Committee(ISSC)model under the first and second wave forces is considered as the case study.The simulation is performed in the MATLAB environment.Moreover,the accuracy and reliability of the programs are verified for cases of beam model with theoretical solutions.It is found that the motion response of tendons is affected by the TLP movement and environmental load,simultaneously.Then,the motion response is calculated using the SESAM software and exported as the boundary of ANCF tendons.Finally,the static and dynamic characteristics of the four tendons of ISSC TLP are analyzed systematically by the ANCF method.Performed analysis proves the effectiveness and feasibility of the ANCF method.It is concluded that the proposed method is a powerful scheme for calculating the dynamics of tendon/riser in the field of ocean engineering.展开更多
The recently developed hard-magnetic soft(HMS)materials can play a significant role in the actuation and control of medical devices,soft robots,flexible electronics,etc.To regulate the mechanical behaviors of the cant...The recently developed hard-magnetic soft(HMS)materials can play a significant role in the actuation and control of medical devices,soft robots,flexible electronics,etc.To regulate the mechanical behaviors of the cantilevered pipe conveying fluid,the present work introduces a segment made of the HMS material located somewhere along the pipe length.Based on the absolute node coordinate formulation(ANCF),the governing equations of the pipe conveying fluid with an HMS segment are derived by the generalized Lagrange equation.By solving the derived equations with numerical methods,the static deformation,linear vibration characteristic,and nonlinear dynamic response of the pipe are analyzed.The result of the static deformation of the pipe shows that when the HMS segment is located in the middle of the pipe,the downstream portion of the pipe centerline will keep a straight shape,providing that the pipe is stable with a relatively low flow velocity.Therefore,it is possible to precisely regulate the ejection direction of the fluid flow by changing the magnetic and fluid parameters.It is also found that the intensity and direction of the external magnetic field greatly affect the stability and dynamic response of the pipe with an HMS segment.In most cases,the magnetic actuation increases the critical flow velocity for the flutter instability of the pipe system and suppresses the vibration amplitude of the pipe.展开更多
Creatures with longer bodies in nature like snakes and eels moving in water commonly generate a large swaying of their bodies or tails,with the purpose of producing significant frictions and collisions between body an...Creatures with longer bodies in nature like snakes and eels moving in water commonly generate a large swaying of their bodies or tails,with the purpose of producing significant frictions and collisions between body and fluid to provide the power of consecutive forward force.This swaying can be idealized by considering oscillations of a soft beam immersed in water when waves of vibration travel down at a constant speed.The present study employs a kind of large deformations induced by nonlinear vibrations of a soft pipe conveying fluid to design an underwater bio-inspired snake robot that consists of a rigid head and a soft tail.When the head is fixed,experiments show that a second mode vibration of the tail in water occurs as the internal flow velocity is beyond a critical value.Then the corresponding theoretical model based on the absolute nodal coordinate formulation(ANCF)is established to describe nonlinear vibrations of the tail.As the head is free,the theoretical modeling is combined with the computational fluid dynamics(CFD)analysis to construct a fluid-structure interaction(FSI)simulation model.The swimming speed and swaying shape of the snake robot are obtained through the FSI simulation model.They are in good agreement with experimental results.Most importantly,it is demonstrated that the propulsion speed can be improved by 21%for the robot with vibrations of the tail compared with that without oscillations in the pure jet mode.This research provides a new thought to design driving devices by using nonlinear flow-induced vibrations.展开更多
为了探究子午线轮胎面内振动特性,建立了旋转超弹性厚壁REF模型。并在具有随动坐标系的绝对节点坐标法(absolute nodal coordinate formulation,ANCF)框架下,提出平面内旋转环扇单元,对旋转超弹性厚壁弹性基环模型(ring on the elastic ...为了探究子午线轮胎面内振动特性,建立了旋转超弹性厚壁REF模型。并在具有随动坐标系的绝对节点坐标法(absolute nodal coordinate formulation,ANCF)框架下,提出平面内旋转环扇单元,对旋转超弹性厚壁弹性基环模型(ring on the elastic foundation,REF)进行离散。考虑到非线性本构关系和旋转运动中的广义惯性力,建立了该模型的非线性运动微分方程以及线性化运动微分方程。通过与现有文献中的试验结果进行对比,验证了上述模型的有效性,并分析了弹性基参数、环体厚度和旋转角速度对子午线轮胎固有特性的影响规律。展开更多
基金supported in part by the National Natural Science Foundation of China (Grants 11290151 and 11221202)supported in part by the Beijing Higher Education Young Elite Teacher Project (Grant YETP1201)
文摘Under the frame of multibody dynamics, the contact dynamics of elasto-plastic spatial thin beams is numerically studied by using the spatial thin beam elements of absolute nodal coordinate formulation（ANCF）. The internal force of the elasto-plastic spatial thin beam element is derived under the assumption that the plastic strain of the beam element depends only on its longitudinal deformation.A new body-fixed local coordinate system is introduced into the spatial thin beam element of ANCF for efficient contact detection in the contact dynamics simulation. The linear isotropic hardening constitutive law is used to describe the elasto-plastic deformation of beam material, and the classical return mapping algorithm is adopted to evaluate the plastic strains. A multi-zone contact approach of thin beams previously proposed by the authors is also introduced to detect the multiple contact zones of beams accurately, and the penalty method is used to compute the normal contact force of thin beams in contact. Four numerical examples are given to demonstrate the applicability and effectiveness of the proposed elasto-plastic spatial thin beam element of ANCF for flexible multibody system dynamics.
基金supported by the National Natural Science Foundation of China(Nos.11972167,12072119)the Alexander von Humboldt Foundation。
文摘In several previous studies,it was reported that a supported pipe with small geometric imperfections would lose stability when the internal flow velocity became sufficiently high.Recently,however,it has become clear that this conclusion may be at best incomplete.A reevaluation of the problem is undertaken here by essentially considering the flow-induced static deformation of a pipe.With the aid of the absolute nodal coordinate formulation(ANCF)and the extended Lagrange equations for dynamical systems containing non-material volumes,the nonlinear governing equations of a pipe with three different geometric imperfections are introduced and formulated.Based on extensive numerical calculations,the static equilibrium configuration,the stability,and the nonlinear dynamics of the considered pipe system are determined and analyzed.The results show that for a supported pipe with the geometric imperfection of a half sinusoidal wave,the dynamical system could not lose stability even if the flow velocity reaches an extremely high value of 40.However,for a supported pipe with the geometric imperfection of one or one and a half sinusoidal waves,the first-mode buckling instability would take place at high flow velocity.Moreover,based on a further parametric analysis,the effects of the amplitude of the geometric imperfection and the aspect ratio of the pipe on the static deformation,the critical flow velocity for buckling instability,and the nonlinear responses of the supported pipes with geometric imperfections are analyzed.
基金supported by the 111 China Project(Grant No.B16003)the National Natural Science Foundation of China(Grant Nos.11290151,11702022,and 11221202)
文摘As a subsequent work of previous studies of authors, a new parallel computation approach is proposed to simulate the coupled dynamics of a rigid-flexible multibody system and compressible fluid. In this approach, the smoothed particle hydrodynamics(SPH) method is used to model the compressible fluid, the natural coordinate formulation(NCF) and absolute nodal coordinate formulation(ANCF) are used to model the rigid and flexible bodies, respectively. In order to model the compressible fluid properly and efficiently via SPH method, three measures are taken as follows. The first is to use the Riemann solver to cope with the fluid compressibility, the second is to define virtual particles of SPH to model the dynamic interaction between the fluid and the multibody system, and the third is to impose the boundary conditions of periodical inflow and outflow to reduce the number of SPH particles involved in the computation process. Afterwards, a parallel computation strategy is proposed based on the graphics processing unit(GPU) to detect the neighboring SPH particles and to solve the dynamic equations of SPH particles in order to improve the computation efficiency. Meanwhile, the generalized-alpha algorithm is used to solve the dynamic equations of the multibody system. Finally, four case studies are given to validate the proposed parallel computation approach.
基金financially supported by the National Natural Science Foundation of China(Grant Nos.51879047 and 51890915)the Engineering Development Program of Deepwater Semi-submersible Production Storage and Unloading Platform of China(Grant No.SSBQ-2020-HN-03-03)the Natural Science Foundation of Heilongjiang Province of China(Grant No.E2017029)。
文摘In the present study,the dynamics of the tendon system of a tension-leg platform(TLP)is investigated through the absolute nodal coordinate formulation(ANCF).Based on the energy conversion principle,the stiffness,generalized elastic force,external load and mass matrices of the element are deduced to perform the element assembling by using the finite element method.Then the motion equation of the tendon/riser is established.In this study,the TLP in the International Ship Structures Committee(ISSC)model under the first and second wave forces is considered as the case study.The simulation is performed in the MATLAB environment.Moreover,the accuracy and reliability of the programs are verified for cases of beam model with theoretical solutions.It is found that the motion response of tendons is affected by the TLP movement and environmental load,simultaneously.Then,the motion response is calculated using the SESAM software and exported as the boundary of ANCF tendons.Finally,the static and dynamic characteristics of the four tendons of ISSC TLP are analyzed systematically by the ANCF method.Performed analysis proves the effectiveness and feasibility of the ANCF method.It is concluded that the proposed method is a powerful scheme for calculating the dynamics of tendon/riser in the field of ocean engineering.
基金supported by the National Natural Science Foundation of China(Nos.11972167 and 12072119)the China National Postdoctoral Program for Innovative Talents(No.BX20220118)+1 种基金the China Postdoctoral Science Foundation(No.2021M701306)the Third Batch Postdoctoral Program for the Innovative Talents in Hubei Province of China。
文摘The recently developed hard-magnetic soft(HMS)materials can play a significant role in the actuation and control of medical devices,soft robots,flexible electronics,etc.To regulate the mechanical behaviors of the cantilevered pipe conveying fluid,the present work introduces a segment made of the HMS material located somewhere along the pipe length.Based on the absolute node coordinate formulation(ANCF),the governing equations of the pipe conveying fluid with an HMS segment are derived by the generalized Lagrange equation.By solving the derived equations with numerical methods,the static deformation,linear vibration characteristic,and nonlinear dynamic response of the pipe are analyzed.The result of the static deformation of the pipe shows that when the HMS segment is located in the middle of the pipe,the downstream portion of the pipe centerline will keep a straight shape,providing that the pipe is stable with a relatively low flow velocity.Therefore,it is possible to precisely regulate the ejection direction of the fluid flow by changing the magnetic and fluid parameters.It is also found that the intensity and direction of the external magnetic field greatly affect the stability and dynamic response of the pipe with an HMS segment.In most cases,the magnetic actuation increases the critical flow velocity for the flutter instability of the pipe system and suppresses the vibration amplitude of the pipe.
基金the National Natural Science Foundation of China(No.12072119)。
文摘Creatures with longer bodies in nature like snakes and eels moving in water commonly generate a large swaying of their bodies or tails,with the purpose of producing significant frictions and collisions between body and fluid to provide the power of consecutive forward force.This swaying can be idealized by considering oscillations of a soft beam immersed in water when waves of vibration travel down at a constant speed.The present study employs a kind of large deformations induced by nonlinear vibrations of a soft pipe conveying fluid to design an underwater bio-inspired snake robot that consists of a rigid head and a soft tail.When the head is fixed,experiments show that a second mode vibration of the tail in water occurs as the internal flow velocity is beyond a critical value.Then the corresponding theoretical model based on the absolute nodal coordinate formulation(ANCF)is established to describe nonlinear vibrations of the tail.As the head is free,the theoretical modeling is combined with the computational fluid dynamics(CFD)analysis to construct a fluid-structure interaction(FSI)simulation model.The swimming speed and swaying shape of the snake robot are obtained through the FSI simulation model.They are in good agreement with experimental results.Most importantly,it is demonstrated that the propulsion speed can be improved by 21%for the robot with vibrations of the tail compared with that without oscillations in the pure jet mode.This research provides a new thought to design driving devices by using nonlinear flow-induced vibrations.
文摘为了探究子午线轮胎面内振动特性,建立了旋转超弹性厚壁REF模型。并在具有随动坐标系的绝对节点坐标法(absolute nodal coordinate formulation,ANCF)框架下,提出平面内旋转环扇单元,对旋转超弹性厚壁弹性基环模型(ring on the elastic foundation,REF)进行离散。考虑到非线性本构关系和旋转运动中的广义惯性力,建立了该模型的非线性运动微分方程以及线性化运动微分方程。通过与现有文献中的试验结果进行对比,验证了上述模型的有效性,并分析了弹性基参数、环体厚度和旋转角速度对子午线轮胎固有特性的影响规律。
