Core-drilling kinematic modeling and analysis of Jiaolong submersible manipulator

The complicated topographies of the deep sea pose significant challenges for the core drilling with the Jiaolong submersible manipulator.To address this problem,we proposed a core-drilling kinematic model and evaluated the core-drilling behavior of the submersible manipulator by comprehensively considering the uncertain posture of the Jiaolong submersible.First,we established a forward kinematic model for the core-drilling task in deep sea,which satisfied the requirement of gravitational-direction core drilling.Based on the forward kinematic equations,we then built a double-redundancy inverse kinematic model,which was able to determine the required motion trajectories of six active joints according to the desired core-drilling trajectory.The core-drilling workspaces and the motions of the Jiaolong submersible manipulator were assessed with several calculation examples.The established forward and inverse kinematic models are constructed with clear analytic equations,and thus are directly applicable to the Jiaolong submersible manipulator-based core-drilling task.

Improved Stiffness Modeling for An Exechon‑Like Parallel Kinematic Machine(PKM)and Its Application

Hole drilling or contour milling for the large and complex workpieces such as automobile panels and aircraft fuselages makes a high combined demand on machining accuracy,stiffness and workspace of machining equipment.Therefore,a 5-DOF(degrees of freedom)parallel kinematic machine(PKM)with redundant constraints is proposed.Based on the kinematics analysis of the parallel mechanism using intermediate variables,the kinematics problems of the PKM are solved through equivalent kinematics model.The structural stiffness matrix method is adopted to model the stiffness of the parallel mechanism of the PKM,where the stiffness of each joint and branch component is obtained by stiffness formula and finite element analysis.And the stiffness model of the parallel mechanism is improved by correction coefficient matrix,each element of which is constructed as a polynomial function of three independent end variables of the parallel mechanism.The terminal stiffness matrices obtained by simulation result are used to determine the coefficients of polynomial function by least square fitting to describe the correction coefficient over the workspace of the parallel mechanism quantitatively.The experiment results prove that the modification method can greatly improve the stiffness model of the parallel mechanism.To enhance the machining accuracy of the PKM,the proposed kinematics model and the improved stiffness model are utilized to optimize the working stiffness of parallel machine by searching the best relative position of parallel machine and workpiece.A plate workpiece taken as example is examined in the case study section,which demonstrates the effectiveness of optimization method.

Evolution of convex structure during counter-rotating electrochemical machining based on kinematic modeling

Counter-rotating electrochemical machining(CRECM)is a novel electrochemical machining(ECM)method,which can be used to machine convex structures with complex shapes on the outer surface of casings.In this study,the evolution of the convex structure during CRECM is studied.The complex motion form of CRECM is replaced by an equivalent kinematic model,in which the movement of the cathode tool is realized by matrix equations.The trajectory of the cathode tool center satisfies the Archimedes spiral equation,and the feed depth in adjacent cycles is a constant.The simulation results show that the variations of five quality indexes for the convex structure:as machining time increases,the height increases linearly,and the width reduces linearly,the fillets at the top and root fit the rational function,and the inclination angle of the convex satisfies the exponential function.The current density distributions with different rotation angles is investigated.Owing to the differential distribution of current density on workpiece surface,the convex is manufactured with the cathode window transferring into and out of the processing area.Experimental results agree very well with the simulation,which indicates that the proposed model is effective for prediction the evolution of the convex structure in CRECM.

Dynamic Modeling and Investigation of Maneuver Characteristics of A Deep-Sea Manned Submarine Vehicle

A deep-sea Manned Submarine Vehicle （MSV） is usually required to move at a low forward speed and a low rotational speed when it executes investigation tasks. In this condition, the motion is in large drift angles, and the maneuverability hydrodynamic forces cannot be expressed properly in the conventional mathematical model of submersible motion. In this paper, firstly, a general equation of MSV with six-freedom motion is presented, and the numerical simulation of descent/ascent motion and helix motion is conducted to reveal the general maneuver characteristics of MSV. Secondly, according to the data of captive model tests of large drift angles of MSV, the regression analysis of position hydrodynamic forces and rotation hydrodynamic forces is carried out, and the results of regression analysis of maneuverability hydrody- namic characteristics are analyzed to reveal the special maneuver characteristics. Thirdly, a special new mathematical model of MSV with the whole range of drift angles motion is presented, which can be used to predict hydrodynamic performance of motion in the 0° - 180° range of drift angles. The results are applied to the design of maneuverability hydrodynamic forces, development of control system and simulator of a practical MSV.

浮动基多机协调吊运系统的工作空间分析

At present,the cranes used at sea have several shortcomings in terms of flexibility,efficiency,and safety.Therefore,a floating multi-robot coordinated lifting system is proposed to fulfill the offshore lifting requirements.First,the structure of the lifting system is established according to the lifting task,the kinematic model of the system is developed by using the D–H coordinate transformation,and the dynamic model is developed based on rigid-body dynamics and hydrodynamics.Then,the static and dynamic workspace of the lifting system are analyzed,and the solving steps of the workspace are given by using the Monte–Carlo method.The effect of the load mass and the maximum allowable tension of the cable on the workspace is examined by simulation.Results show that the lifting system has limited carrying capacity and a data reference for selecting the structural parameters by analyzing the factors affecting the workspace.Findings provide a basis for further research on the optimal design of structural parameters and the determination of safe configurations of the lifting system.

Magnetostrictive and Kinematic Model Considering the Dynamic Hysteresis and Energy Loss for GMA

Due to the influence of magnetic hysteresis and energy loss inherent in giant magnetostrictive materials （GMM）, output displacement accuracy of giant magnetostrictive actuator （GMA） can not meet the precision and ultra precision machining. Using a GMM rod as the core driving element, a GMA which may be used in the field of precision and ultra precision drive engineering is designed through modular design method. Based on the Armstrong theory and elastic Gibbs free energy theory, a nonlinear magnetostriction model which considers magnetic hysteresis and energy loss characteristics is established. Moreover, the mechanical system differential equation model for GMA is established by utilizing D＇Alembert＇s principle. Experimental results show that the model can preferably predict magnetization property, magnetic potential orientation, energy loss for GMM. It is also able to describe magnetostrictive elongation and output displacement of GMA. Research results will provide a theoretical basis for solving the dynamic magnetic hysteresis, energy loss and working precision for GMA fundamentally.

A Kinematic Thermal Model for Descending Slabs with Velocity Boundary Layers:A Case Study for the Tonga Subducting Slab

For the purpose of investigating the influence of metastable olivine（MO） phase transformations on both deep seismicity and stagnation of slabs,we constructed a 2-dimensional finite element thermal model for a 120 Ma-old 50°dipping oceanic lithosphere descending at 10 cm/yr with velocity boundary layers,which would mitigate the interference of constant velocity field for the slab. The resulting temperatures show that most of intermediate and deep earthquakes occurring within the Tonga slab are occurring inside the 800℃and 1200℃isotherm,respectively.The elevation of olivine transformation near～410 km and respective persistence of metastable olivine and spinel within the transition zone and beneath 660 km would thus result in bimodal positive,zonal,negative density anomalies,respectively.These results together with the resulting pressure anomalies may reflect the stress pattern of the Tonga slab：（i） slab pull force exerts above a depth of～230 km;（ii） MO existence changes the buoyancy force within the transition zone and facilitates slab stagnation at a depth of 660 km;（iii） as the subducting materials accumulated over 660 km,deepest earthquakes occur due to MO transformation;（iv） a flattened‘slab’ may penetrate into the lower mantle due to the density increment of Sp transformation.

Parameter Identification and Application of Slippage Kinematics for Tracked Mobile Robots

A new parameter identification method is proposed to solve the slippage problem when tracked mobile robots execute turning motions.Such motion is divided into two states in this paper:pivot turning and coupled turning between angular velocity and linear velocity.In the processing of pivot turning,the slippage parameters could be obtained by measuring the end point in a square path.In the process of coupled turning,the slippage parameters could be calculated by measuring the perimeter of a circular path and the linear distance between the start and end points.The identification results showed that slippage parameters were affected by velocity.Therefore,a fuzzy rule base was established with the basis on the identification data,and a fuzzy controller was applied to motion control and dead reckoning.This method effectively compensated for errors resulting in unequal tension between the left and right tracks,structural dimensions and slippage.The results demonstrated that the accuracy of robot positioning and control could be substantially improved on a rigid floor.

Three-Dimensional Kinematics Simulation of Robot Fighting Platform in Virtual Environment

A method of 3 D kinematics simulation of robot fighting platform （RFP） in virtual environment is proposed with the aim of enhancing vision telepresence. Based on the theory of space coordinate transformation, kinematics equat!ons of RFP are formulated; followed by applying a method of modeling using 3DMAX software to build an RFP＇s 3D geometric model before a 3D kinematics simulation system of RFP is completed based on virtual reality technology and Open Inventor VC ＋ ＋. Test results have indicated that this system can perform RFP＇s kinematics simulation in virtual environment. It can also imitate RFP＇s motion states and environmental features well. Moreover, not only can better real-time performances and interactions be achieved but also operator＇s vision telepresence be enhanced, therefore this approach may help lay the foundation for the realization of RFP＇s teleoperation with vision telepresence.

Preliminary results on kinematic model of tectonic blocks derived from high precision GPS observations in Southwest China

In the paper, the kinematic model of tectonic blocks in southwest China is studied based on the precision GPS observations carried out under the major subject of 'Studies on Current Crustal Movement and Geodynamics' which belongs to the State Climbing Project. It is believed that at present, the data of high precision GPS observation may provide convincing information related to the horizontal movement of tectonic blocks in the Chinese mainland. The preliminary results obtained from the kinematic model have given some direct evidences for the research of dynamic mechanism of crustal deformation in the Chinese mainland and on the basis of which, the kinematic characteristics and their relations to the seismicity and seismic risk in the reobserved region are analysed. The preliminary observation results are encouraging.

A Generic Kinematic Model for Three Main Types of Five-axis Machine Tools

Material removal is one of the most used processes in manufacturing. Five-axis CNC machines are believed to be the best tools in sculptured surface machining. In this study, a generic and unified kinematic model was developed as a viable alternative to the particular solutions that are only applicable to individual machine configurations. This versatile model is then used to verify the feasibility of the two rotational joints within the kinematic chain of three main types of a five-axis machine-tool. This versatile model is very useful applied to the design of five-axis machine tools.

Kinematic source model for simulation of near-fault ground motion field using explicit finite element method

This paper briefly reviews the characteristics and major processes of the explicit finite element method in modeling the near-fault ground motion field. The emphasis is on the finite element-related problems in the finite fault source modeling. A modified kinematic source model is presented, in which vibration with some high frequency components is introduced into the traditional slip time function to ensure that the source and ground motion include sufficient high frequency components. The model presented is verified through a simple modeling example. It is shown that the predicted near-fault ground motion field exhibits similar characteristics to those observed in strong motion records, such as the hanging wall effect, vertical effect, fling step effect and velocity pulse effect, etc.

Anisotropic Constitutive Modeling of Compressible Biological Tissue

The anisotropic continuum stored energy density (ACSED) functional is applied for accurate constitutive modeling of biological tissues and finite element implementation without the isochoric—volumetric split, the anisotropic—isotropic split, or the anisotropic invariant split. Related stress and elasticity tensors in the reference and current configurations are worked out. A new kinematic model is derived based on the tangent Poisson’s ratio as a cubic polynomial function of stretch. The ACSED model, along with the kinematic model, accurately fits uniaxial extension test data for compressible human skin, bovine articular cartilage, and human aorta samples.

A Method for Automatically Establishing a Mathematical Model of Kinematic Analysis to a Planar Mechanism with a Multiple Joint and Prismatic Pair

A method for automatically establishing a mathematical model of kinematic analysis to a planar mechanism with multiple joint and prismatic pair is presented. The breadth （ or depth ） first search spanning tree can be obtained based on an adjacency matrix of the mechanism. Then the kinematic chain （or mechanism）＇s basic loops can be obtained. On the basis of these basic loops, a mathematical model of kinematic analysis can be established and solved automatically. In the sense of a calculative mechanism, structural analysis of the kinematic chain relates to the kinematic analysis of a mechanism. Thus, an effective way is supplied to the given mechanism＇s kinematic analysis for automatic modeling and solving, and a method is supplied to the structural type to optimize kinematic synthesis.

Accurate parameter estimation of systematic odometry errors for two-wheel differential mobile robots

Odometry using incremental wheel encoder sensors provides the relative robot pose estimation.However,the odometry suffers from the accumulation of kinematic modeling errors of wheels as the robot's travel distance increases.Therefore,the systematic errors need to be calibrated.The University of Michigan Benchmark(UMBmark)method is a widely used calibration scheme of the systematic errors in two-wheel differential mobile robots.In this paper,the accurate parameter estimation of systematic errors is proposed by extending the conventional method.The contributions of this paper can be summarized as two issues.The first contribution is to present new calibration equations that reduce the systematic odometry errors.The new equations were derived to overcome the limitation of conventional schemes.The second contribution is to propose the design guideline of the test track for calibration experiments.The calibration performance can be improved by appropriate design of the test track.The simulations and experimental results show that the accurate parameter estimation can be implemented by the proposed method.

Study on Dynamic Characteristics of Hydraulic Pumping Unit on Offshore Platform

A new technology of offshore oil rod pumping production is developed for offshore heavy oil recovery. A new type of miniature hydraulic pumping unit with long-stroke, low pumping speed and compact structure is designed based on the spatial characteristics of offshore platforms. By combining the strengths of sinusoidal velocity curve and trapezoidal velocity curve, a kinematical model of the acceleration, the velocity and displacement of the pumping unit＇s hanging point is established. The results show that the pumping unit has good kinematic characteristics of smooth motion and small dynamic load. The multi-degree-of-freedom dynamic model of the single-well pumping unit is established. The first and second order natural frequencies of the sucker rod string subsystem and the pumping unit subsystem are studied. The results show that the first and the second order natural frequencies among the pumping rod string, pumping unit-platform subsystem and the dynamic excitation have differences over 5 times from each other, indicating that resonance phenomenon will not appear during the operation and the dynamic requirements for field use are met in the system.

Algebraic Method‑Based Point‑to‑Point Trajectory Planning of an Under‑Constrained Cable‑Suspended Parallel Robot with Variable Angle and Height Cable Mast

To avoid impacts and vibrations during the processes of acceleration and deceleration while possessing flexible working ways for cable-suspended parallel robots(CSPRs),point-to-point trajectory planning demands an under-constrained cable-suspended parallel robot(UCPR)with variable angle and height cable mast as described in this paper.The end-effector of the UCPR with three cables can achieve three translational degrees of freedom(DOFs).The inverse kinematic and dynamic modeling of the UCPR considering the angle and height of cable mast are completed.The motion trajectory of the end-effector comprising six segments is given.The connection points of the trajectory segments(except for point P3 in the X direction)are devised to have zero instantaneous velocities,which ensure that the acceleration has continuity and the planned acceleration curve achieves smooth transition.The trajectory is respectively planned using three algebraic methods,including fifth degree polynomial,cycloid trajectory,and double-S velocity curve.The results indicate that the trajectory planned by fifth degree polynomial method is much closer to the given trajectory of the end-effector.Numerical simulation and experiments are accomplished for the given trajectory based on fifth degree polynomial planning.At the points where the velocity suddenly changes,the length and tension variation curves of the planned and unplanned three cables are compared and analyzed.The OptiTrack motion capture system is adopted to track the end-effector of the UCPR during the experiment.The effectiveness and feasibility of fifth degree polynomial planning are validated.

Practical approach to determining traffic congestion propagation boundary due to traffic incidents

A practical approach for predicting the congestion boundary due to traffic incidents was proposed. Based on the kinematic wave theory and Van Aerde single-regime flow model, a model for estimating the congestion propagation speed for the basic road segment was developed. Historical traffic flow data were used to analyze the time variant characteristics of the urban traffic flow for each road type. Then, the saturation flow rate was used for analyzing the impact of the traffic incident on the traversing traffic flow at the congestion area. The base congestion propagation speed for each road type was calculated based on field data, which were provided by the remote traffic microwave sensors(RTMS), floating car data(FCD) system and screen line survey. According to a comparative analysis of the congestion propagation speed, it is found that the expressway, major arterial, minor arterial and collector are decreasingly influenced by the traffic incident. Subsequently, the impact of turning movements at intersections on the congestion propagation was considered. The turning ratio was adopted to represent the impact of turning movements, and afterward the corresponding propagation pattern at intersections was analyzed. Finally, an implementation system was designed on a geographic information system(GIS) platform to display the characteristics of the congestion propagation over the network. The validation results show that the proposed approach is able to capture the congestion propagation properties in the actual road network.

Modeling and predicting low-speed vehicle emissions as a function of driving kinematics

An instantaneous emission model was developed to model and predict the real driving emissions of the low-speed vehicles. The emission database used in the model was measured by using portable emission measurement system （PEMS） under actual traffic conditions in the rural area, and the characteristics of the emission data were determined in relation to the driving kinematics （speed and acceleration） of the low-speed vehicle. The input of the emission model is driving cycle, and the model requires instantaneous vehicle speed and acceleration levels as input variables and uses them to interpolate the pollutant emission rate maps to calculate the transient pollutant emission rates, which will be accumulated to calculate the total emissions released during the whole driving cycle. And the vehicle fuel consumption was determined through the carbon balance method. The model predicted the emissions and fuel consumption of an in-use low-speed vehicle type model, which agreed well with the measured data.

Ground-motion simulation for the MW6.1 Ludian earthquake on 3 August 2014 using the stochastic finite-fault method

The stochastic finite-fault simulation method was applied to synthesize the horizontal ground acceleration seismograms produced by the MW6.1 Ludian earthquake on August 3,2014.For this purpose,we produced first a total of 200 kinematic source models for the Ludian event,which are characterized by the heterogeneous slip on the conjugated ruptured fault and the slip-dependent spreading of the rupture front.The results indicated that the heterogeneous slip and the spatial extent of the ruptured fault play dominant roles in the spatial distribution of ground motions in the near-fault area.The peak ground accelerations(PGAs)and 5%-damped pseudospectral accelerations(PSAs)at periods shorter than 0.5 s estimated on the resulting synthetics generally match well with the observations at stations with Joyner-Boore distances(RJB)greater than 20 km.The synthetic PGVs and PSAs at periods of 0.5 s and 0.75 s are in good agreement with predicted medians by the Yu14 model(Yu et al.,2014).However,the synthetic results are generally much lower than the predicted medians by BSSA14 model(Boore et al.,2014).Moreover,the ground motion variability caused by the randomness in the source rupture process was evaluated by these synthetics.The standard deviations of PSAs on the base-10 logarithmic scale,Sigma[log10(PSA)],are closely dependent on either the spectral period or the RJB.The Sigma[log10(PSA)]remains a constant approximately 0.55 at periods shorter than 0.1 s,and then increase continuously up to^0.13 as the period increases from 0.1 to 2.0 s.The Sigma[log10(PSA)]values at periods of 0.1‒2.0 s show the downward tendency as the RJB values increase.However,the Sigma[log10(PSA)]​values at periods shorter than 0.1 s decrease as the RJB values increase up to^50 km,and then increase with the increasing RJB.Furthermore,we found that the ground-motion variability shows the significant dependence on the azimuth.