Different material properties leads to different metal fracture behaviors. Even if the powder material is composed of plastic metal, the fracture still does not show macroscopic plastic deformation characteristics if ...Different material properties leads to different metal fracture behaviors. Even if the powder material is composed of plastic metal, the fracture still does not show macroscopic plastic deformation characteristics if the material contains a large number of voids. Eight node isoparametric elastic plastic finite element method was used to simulate the tensile process of sintered powder material. By setting a number of voids in the analyzed metal cuboid, the initial density was taken into consideration. The material properties of the three dimensional solid for the tensile simulation were defined with reference to the known pure iron material parameters. The load displacement curves during elongation were obtained with a universal testing machine, and then the simulated curves were compared with the experimental results. The factors that cause the stress concentration and strength decrease were analyzed according to the simulated equivalent von Mises stress distribution.展开更多
Base d on fluid velocity potential, an ALE finite element formulation for the analysi s of nonlinear sloshing problems has been developed. The ALE kinemat ical description is introduced to move the computational mesh...Base d on fluid velocity potential, an ALE finite element formulation for the analysi s of nonlinear sloshing problems has been developed. The ALE kinemat ical description is introduced to move the computational mesh independently of f luid motion, and the container fixed noninertial coordinate system is employed to establish the governing equations so that the mesh is needed to be updated in this coordinate system only. This leads to a very simple mesh moving algorithm which makes it easy to trace the motion of the moving boundaries and the free su rface without producing undesirable distortion of the computational mesh. The fi nite element method and finite difference method are used spacewise and timewise , respectively. A numerical example involving either forced horizontal oscillati on or forced pitching oscillation of the fluid filled container is presented to illustrate the effectiveness and the robustness of the method. In additi on, this work can be extended for the fluid structure interaction problems.展开更多
For studying the dynamic performance of subsea umbilical cable laying system and achieving the goal of cable tension and laying speed control,the rigid finite element method is used to discrete and transform the syste...For studying the dynamic performance of subsea umbilical cable laying system and achieving the goal of cable tension and laying speed control,the rigid finite element method is used to discrete and transform the system into a rigid-flexible coupling multi-body system which consists of rigid elements and spring-damping elements.The mathematical model of subsea umbilical cable laying system kinematic chain is presented with the second order Lagrange equation in the joint coordinate system,and dynamic modeling and simulation is performed with ADAMS.The dynamic analysis is conducted assuming the following three statuses:ideal laying,practical laying under wave disturbance,and practical laying with tension compensation.Results show that motion disturbances of the laying budge under sea waves,especially with heaving and pitching,will cause relatively serious fluctuations in cable tension and laying speed.Tension compensation,i.e.,active back tension torque control can restrict continuous tension increasing or decreasing effectively and rapidly,thus avoiding cable breach or buckling.展开更多
Numerical simulation of a two-dimensional nonlinear sloshing problem is preceded by the finite element method. Two theories are used. One is fully nonlinear theory; the other is time domain second order theory. A liqu...Numerical simulation of a two-dimensional nonlinear sloshing problem is preceded by the finite element method. Two theories are used. One is fully nonlinear theory; the other is time domain second order theory. A liquid sloshing in a rectangular container subjected to a horizontal excitation is simulated using these two theories. Numerical results are obtained and comparisons are made. It is found that a good agreement is obtained for the case of small amplitude oscillation. For the situation of large amplitude excitation, although the differences between using the two theories are obvious the second order solution can still exhibit typical nonlinear features of nonlinear wave.展开更多
The cross-flow(CF)vortex-induced vibration(VIV)of a deepwater steep wave riser(SWR)subjected to uniform or shear flow loads is investigated numerically.The model is based on a three-dimensional(3D)nonlinear elastic ro...The cross-flow(CF)vortex-induced vibration(VIV)of a deepwater steep wave riser(SWR)subjected to uniform or shear flow loads is investigated numerically.The model is based on a three-dimensional(3D)nonlinear elastic rod theory coupled with a wake oscillator model.In this numerical simulation,the nonlinear motion equations of the riser with large deformation features are established in a global coordinate system to avoid the transformation between global and local coordinate systems,and are discretized with the time-domain finite element method(FEM).A wakeoscillator model is employed to study the vortex shedding,and the lift force generated by the wake flow is described in a van der Pol equation.A Newmark-βiterative scheme is used to solve their coupling equation for the VIV response of the SWR.The developed model is validated against the existing experimental results for the VIV response of the top-tension riser(TTR).Then,the numerical simulations are executed to determine VIV characteristics of the SWR.The effects of both flow velocity and the spanwise length of the flow field on the drag coefficient in the inline(IL)direction and the lift coefficient in the CF direction are investigated systematically.The results indicate that compared with TTR,the low frequency and multi-modal vibration are the main components of the SWR due to the large deformation and flexible characteristics.For shear flow,the multi-frequency resonance dominates the VIV response of the SWR,especially at the hang-off segment.展开更多
The significant advantage of the complex resistivity method is to reflect the abnormal body through multi-parameters, but its inversion parameters are more than the resistivity tomography method. Therefore, how to eff...The significant advantage of the complex resistivity method is to reflect the abnormal body through multi-parameters, but its inversion parameters are more than the resistivity tomography method. Therefore, how to effectively invert these spectral parameters has become the focused area of the complex resistivity inversion. An optimized BP neural network (BPNN) approach based on Quantum Particle Swarm Optimization (QPSO) algorithm was presented, which was able to improve global search ability for complex resistivity multi-parameter nonlinear inversion. In the proposed method, the nonlinear weight adjustment strategy and mutation operator were used to enhance the optimization ability of QPSO algorithm. Implementation of proposed QPSO-BPNN was given, the network had 56 hidden neurons in two hidden layers (the first hidden layer has 46 neurons and the second hidden layer has 10 neurons) and it was trained on 48 datasets and tested on another 5 synthetic datasets. The training and test results show that BP neural network optimized by the QPSO algorithm performs better than the BP neural network without initial optimization on the inversion training and test models, and the mean square error distribution is better. At the same time, a double polarized anomalous bodies model was also used to verify the feasibility and effectiveness of the proposed method, the inversion results show that the QPSO-BP algorithm inversion clearly characterizes the anomalous boundaries and is closer to the values of the parameters.展开更多
Based on the actual measured well depth, azimuth and oblique angles, a novel interpolation method to obtain the well axis is developed. The initial stress of drill string at the reference state consistent with well ax...Based on the actual measured well depth, azimuth and oblique angles, a novel interpolation method to obtain the well axis is developed. The initial stress of drill string at the reference state consistent with well axis can be obtained from the curvature and the tortuosity of well axis. By using the principle of virtual work, the formula to compute the equivalent load vector of the initial stress was derived. In the derivation,the natural （curvilinear） coordinate system was adopted since both the curvature and the tortuosity were generally not zero. A set of displacement functions fully reflecting the rigid body modes was used. Some basic concepts in the finite element analysis of drill string were clarified. It is hoped that the proposed method would offer a theoretical basis for handling the geometric nonlinear problem of the drill string in a 3-D larg edisplacement wellbore.展开更多
文摘Different material properties leads to different metal fracture behaviors. Even if the powder material is composed of plastic metal, the fracture still does not show macroscopic plastic deformation characteristics if the material contains a large number of voids. Eight node isoparametric elastic plastic finite element method was used to simulate the tensile process of sintered powder material. By setting a number of voids in the analyzed metal cuboid, the initial density was taken into consideration. The material properties of the three dimensional solid for the tensile simulation were defined with reference to the known pure iron material parameters. The load displacement curves during elongation were obtained with a universal testing machine, and then the simulated curves were compared with the experimental results. The factors that cause the stress concentration and strength decrease were analyzed according to the simulated equivalent von Mises stress distribution.
文摘Base d on fluid velocity potential, an ALE finite element formulation for the analysi s of nonlinear sloshing problems has been developed. The ALE kinemat ical description is introduced to move the computational mesh independently of f luid motion, and the container fixed noninertial coordinate system is employed to establish the governing equations so that the mesh is needed to be updated in this coordinate system only. This leads to a very simple mesh moving algorithm which makes it easy to trace the motion of the moving boundaries and the free su rface without producing undesirable distortion of the computational mesh. The fi nite element method and finite difference method are used spacewise and timewise , respectively. A numerical example involving either forced horizontal oscillati on or forced pitching oscillation of the fluid filled container is presented to illustrate the effectiveness and the robustness of the method. In additi on, this work can be extended for the fluid structure interaction problems.
基金Supported by the National Science and Technology Major Project:Development of Large Oil-gas Fields&Coal-bed Methane(No.2011ZX05056)
文摘For studying the dynamic performance of subsea umbilical cable laying system and achieving the goal of cable tension and laying speed control,the rigid finite element method is used to discrete and transform the system into a rigid-flexible coupling multi-body system which consists of rigid elements and spring-damping elements.The mathematical model of subsea umbilical cable laying system kinematic chain is presented with the second order Lagrange equation in the joint coordinate system,and dynamic modeling and simulation is performed with ADAMS.The dynamic analysis is conducted assuming the following three statuses:ideal laying,practical laying under wave disturbance,and practical laying with tension compensation.Results show that motion disturbances of the laying budge under sea waves,especially with heaving and pitching,will cause relatively serious fluctuations in cable tension and laying speed.Tension compensation,i.e.,active back tension torque control can restrict continuous tension increasing or decreasing effectively and rapidly,thus avoiding cable breach or buckling.
文摘Numerical simulation of a two-dimensional nonlinear sloshing problem is preceded by the finite element method. Two theories are used. One is fully nonlinear theory; the other is time domain second order theory. A liquid sloshing in a rectangular container subjected to a horizontal excitation is simulated using these two theories. Numerical results are obtained and comparisons are made. It is found that a good agreement is obtained for the case of small amplitude oscillation. For the situation of large amplitude excitation, although the differences between using the two theories are obvious the second order solution can still exhibit typical nonlinear features of nonlinear wave.
基金financially supported by the National Natural Science Foundation of China(Grant Nos.52111530137 and 52025112)the Natural Science Found of Jiangsu Province(Grant No.BK20160556)the Jiangsu Provincial Higher Education Natural Science Research Major Project(Grant No.18KJA580003)。
文摘The cross-flow(CF)vortex-induced vibration(VIV)of a deepwater steep wave riser(SWR)subjected to uniform or shear flow loads is investigated numerically.The model is based on a three-dimensional(3D)nonlinear elastic rod theory coupled with a wake oscillator model.In this numerical simulation,the nonlinear motion equations of the riser with large deformation features are established in a global coordinate system to avoid the transformation between global and local coordinate systems,and are discretized with the time-domain finite element method(FEM).A wakeoscillator model is employed to study the vortex shedding,and the lift force generated by the wake flow is described in a van der Pol equation.A Newmark-βiterative scheme is used to solve their coupling equation for the VIV response of the SWR.The developed model is validated against the existing experimental results for the VIV response of the top-tension riser(TTR).Then,the numerical simulations are executed to determine VIV characteristics of the SWR.The effects of both flow velocity and the spanwise length of the flow field on the drag coefficient in the inline(IL)direction and the lift coefficient in the CF direction are investigated systematically.The results indicate that compared with TTR,the low frequency and multi-modal vibration are the main components of the SWR due to the large deformation and flexible characteristics.For shear flow,the multi-frequency resonance dominates the VIV response of the SWR,especially at the hang-off segment.
文摘The significant advantage of the complex resistivity method is to reflect the abnormal body through multi-parameters, but its inversion parameters are more than the resistivity tomography method. Therefore, how to effectively invert these spectral parameters has become the focused area of the complex resistivity inversion. An optimized BP neural network (BPNN) approach based on Quantum Particle Swarm Optimization (QPSO) algorithm was presented, which was able to improve global search ability for complex resistivity multi-parameter nonlinear inversion. In the proposed method, the nonlinear weight adjustment strategy and mutation operator were used to enhance the optimization ability of QPSO algorithm. Implementation of proposed QPSO-BPNN was given, the network had 56 hidden neurons in two hidden layers (the first hidden layer has 46 neurons and the second hidden layer has 10 neurons) and it was trained on 48 datasets and tested on another 5 synthetic datasets. The training and test results show that BP neural network optimized by the QPSO algorithm performs better than the BP neural network without initial optimization on the inversion training and test models, and the mean square error distribution is better. At the same time, a double polarized anomalous bodies model was also used to verify the feasibility and effectiveness of the proposed method, the inversion results show that the QPSO-BP algorithm inversion clearly characterizes the anomalous boundaries and is closer to the values of the parameters.
文摘Based on the actual measured well depth, azimuth and oblique angles, a novel interpolation method to obtain the well axis is developed. The initial stress of drill string at the reference state consistent with well axis can be obtained from the curvature and the tortuosity of well axis. By using the principle of virtual work, the formula to compute the equivalent load vector of the initial stress was derived. In the derivation,the natural （curvilinear） coordinate system was adopted since both the curvature and the tortuosity were generally not zero. A set of displacement functions fully reflecting the rigid body modes was used. Some basic concepts in the finite element analysis of drill string were clarified. It is hoped that the proposed method would offer a theoretical basis for handling the geometric nonlinear problem of the drill string in a 3-D larg edisplacement wellbore.
