Based on the field velocity method,a novel approach for simulating unsteady pitching and plunging motion of an airfoil is presented in this paper.Responses to pitching and plunging motions of the airfoil are simulated...Based on the field velocity method,a novel approach for simulating unsteady pitching and plunging motion of an airfoil is presented in this paper.Responses to pitching and plunging motions of the airfoil are simulated under different conditions.The obtained results are compared with those of moving grid method and good agreement is achieved.In the conventional field velocity method,the Euler solver is usually used to simulate the movement of the airfoil.However,when viscous effect is considered,unsteady Navier-Stokes equations have to be solved and the viscous flux correction must be taken into account.In this work,the viscous flux correction is introduced into the conventional field velocity method when non-uniform grid speed distribution is occurred.Numerical experiments for the flow around NACA0012 airfoil showed that the proposed approach can well simulate viscous moving boundary flow problems.展开更多
There is a common difficulty in elastic-plastic impact codes such as EPIC[2,3] NONSAP[4], etc.. Most of these codes use the simple linear functions usually taken from static problem to represent the displacement compo...There is a common difficulty in elastic-plastic impact codes such as EPIC[2,3] NONSAP[4], etc.. Most of these codes use the simple linear functions usually taken from static problem to represent the displacement components. In such finite element formulation, the stress components are constant in each element and they are discontinuous in any two neighboring elements. Therefore, the bases of using the virtual work principle in such elements are unreliable. In this paper, we introduce a new method, namely, the compatible stress iterative method, to eliminate the above-said difficulty. The calculated examples show that the calculation using the new method in dynamic finite element analysis of high velocity impact is valid and stable, and the element stiffness can be somewhat reduced.展开更多
We apply the method of guidance by a required velocity for solving the optimal control problem over spacecraft’s reorientation from known initial attitude into a required final attitude.We suppose that attitude contr...We apply the method of guidance by a required velocity for solving the optimal control problem over spacecraft’s reorientation from known initial attitude into a required final attitude.We suppose that attitude control is carried out by impulse jet engines.For optimization of fuel consumption,the controlling moments are calculated and formed according to the method of free trajectories together with principle of iterative control using the quaternions for generating commands to actuators.Optimal solution corresponds to the principle“acceleration-free rotation-separate corrections-free rotation-braking”.Rotation along a hitting trajectory is supported by insignificant correction of the uncontrolled motion at discrete instants between segments of acceleration and braking.Various strategies of forming the correction impulses during stage of free motion are suggested.Improving accuracy of achievement of spacecraft's final position is reached by terminal control using information about current attitude and angular velocity measurements for determining an instant of beginning of braking(condition for start of braking based on actual motion parameters is formulated in analytical form).The described method is universal and invariant relative to moments of inertia.Developed laws of attitude control concern the algorithms with prognostic model,the synthesized control modes are invariant with respect to both external perturbations and parametric errors.Results of mathematical modeling are presented that demonstrate practical feasibility and high efficiency of designed algorithms.展开更多
The data-driven phenomenological models based on deformation measurements have been widely utilized to predict the slope failure time(SFT).The observational and model uncertainties could lead the predicted SFT calcula...The data-driven phenomenological models based on deformation measurements have been widely utilized to predict the slope failure time(SFT).The observational and model uncertainties could lead the predicted SFT calculated from the phenomenological models to deviate from the actual SFT.Currently,very limited study has been conducted on how to evaluate the effect of such uncertainties on SFT prediction.In this paper,a comprehensive slope failure database was compiled.A Bayesian machine learning(BML)-based method was developed to learn the model and observational uncertainties involved in SFT prediction,through which the probabilistic distribution of the SFT can be obtained.This method was illustrated in detail with an example.Verification studies show that the BML-based method is superior to the traditional inverse velocity method(INVM)and the maximum likelihood method for predicting SFT.The proposed method in this study provides an effective tool for SFT prediction.展开更多
Cylindrical Cam Mechanism which is one of the best eq uipments to accomplish an accurate motion transmission is widely used in the fie lds of industries, such as machine tool exchangers, textile machinery and automa t...Cylindrical Cam Mechanism which is one of the best eq uipments to accomplish an accurate motion transmission is widely used in the fie lds of industries, such as machine tool exchangers, textile machinery and automa tic transfer equipments. This paper proposes a new approach for the shape design and manufacturing of the cylindrical cam. The design approach uses the relative velocity concept and the manufacturing approach uses the inverse kinematics concept. For the shape desig n, the contact points between the cam and the follower roller are calculated bas ed on relative velocity of which the direction is on the common tangential line, and then the whole shape of cam is determined from transformation of the coordi nate system. For the manufacturing procedures, the location and the orientation of cutter path can be allocated corresponding to the designed shape data. The in tegral NC code for multi-axis CNC machining center is created using the inverse kinematics concept from the data of the location and the orientation of cutter path. As the advantages of the proposed approach, the machine tool is designed t o having an alternative size in fabricating the general cam, while the tool must be fitted to diameter size of the follower in the conventional approach. Finally, CAD/CAM program, "Cylindrical DAM", is developed on C++ lan guage. This program can perform shape design, manufacturing and kinematics simul ation, which can make integral NC code for multi-axis CNC machining center. The proposed method can be applied easily on fields of industries.展开更多
This paper employs a velocity plus displacement（V＋D）-based equivalent force control（EFC） method to solve the velocity/displacement difference equation in a real-time substructure test. This method uses type 2 fee...This paper employs a velocity plus displacement（V＋D）-based equivalent force control（EFC） method to solve the velocity/displacement difference equation in a real-time substructure test. This method uses type 2 feedback control loops to replace mathematical iteration to solve the nonlinear dynamic equation. A spectral radius analysis of the amplification matrix shows that the type 2 EFC-explicit, Newmark-β method has beneficial numerical characteristics for this method. Its stability limit of Ω = 2 remains unchanged regardless of the system damping because the velocity is achieved with very high accuracy during simulation. In contrast, the stability limits of the central difference method using direct velocity prediction and the EFC-average acceleration method with linear interpolation are shown to decrease with an increase in system damping. In fact, the EFC-average acceleration method is shown to change from unconditionally stable to conditionally stable. We also show that if an over-damped system with a damping ratio of 1.05 is considered, the stability limit is reduced to Ω =1.45. Finally, the results from an experiment with a single-degree-of-freedom structure installed with a magneto-rheological（MR） damper are presented. The results demonstrate that the proposed method is able to follow both displacement and velocity commands with moderate accuracy, resulting in improved test performance and accuracy for structures that are sensitive to both velocity and displacement inputs. Although the findings of the study are promising, additional test data and several further improvements will be required to draw general conclusions.展开更多
A time splitting technique is common to many free surface ocean models. The different truncation errors in the equations of the internal and external modes require a numerical adjustment to make sure that algorithms c...A time splitting technique is common to many free surface ocean models. The different truncation errors in the equations of the internal and external modes require a numerical adjustment to make sure that algorithms correctly satisfy continuity equations and conserve tracers quantities. The princeton ocean model （POM） has applied a simple method of adjusting the vertical mean of internal velocities to external velocities at each internal time step. However, due to the Asselin time filter method adopted to prevent the numerical instability, the method of velocity adjustment used in POM can no longer guarantee the satisfaction of the continuity equation in the internal mode, though a special treatment is used to relate the surface elevation of the internal mode with that of the external mode. The error is proved to be a second-order term of the coefficient in the Asselin filter. One influence of this error in the numerical model is the failure of the kinetic boundary condition at the sea floor. By a regional experiment and a quasi-global experiment, the magni- tudes of this error are evaluated, and several sensitivity tests of this error are performed. The characteristic of this error is analyzed and two alternative algorithms are suggested to reduce the error.展开更多
An improved algorithm for velocity field of general configurations ispresentd for low-order panel method based on the internal Dirichlet boundary condi-tion. A direct calculating method for the velocity distribution b...An improved algorithm for velocity field of general configurations ispresentd for low-order panel method based on the internal Dirichlet boundary condi-tion. A direct calculating method for the velocity distribution by means of a limit pro-cess combining with analytic evaluation of higher-order singular integrals instead of theconventional method of doublet strength gradient is devised in order to avoid the diffi-culty of edge extrapolation of doublet strength. The problem of substantialunderpredictions of the induced drag coefficient obtained from the VSAERO analysisdisappears for the present improved algorithm. Illustrative calculations for several testcases such as swept back wing, swept forward wing and wing-body combination showthat the accuracy of results may be improved and is competitive with high-order panelmethod. In addition, the present direct integral method can be used to evaluate the ve-locity distribution for external flow field correctly, where the method of gradient cannot be used at all.展开更多
Based on the seismic observation report data provided by the Xinjiang Digital Seismic Network from 2009 to 2014,we calculate the wave velocity ratio and its background value for medium and small earthquakes by using t...Based on the seismic observation report data provided by the Xinjiang Digital Seismic Network from 2009 to 2014,we calculate the wave velocity ratio and its background value for medium and small earthquakes by using the multi-station method in Tianshan,Xinjiang.This paper analyzes the variation of the wave velocity ratio disturbance value to highlight the abnormal,and also back-traces 7 moderate earthquakes at the research area.The results show that:(1)the background value of the wave velocity ratio is almost 1.70,the wave velocity ratio obviously decreases in the middle-eastern part of Tianshan and the region near the Puchang fault;(2)the wave velocity ratio disturbance value is mostly low in the epicenter before four earthquakes of M≥5.0 from 2011 to 2013 in the study area;(3)before 7 moderate strong earthquakes,the earthquake events with low value of the wave velocity ratio account for over 60% of corresponding total events near the epicenters,and the low value of the wave velocity ratio is relatively obvious before moderate earthquakes.展开更多
Great interest has been aroused on deeply-situated Osaka clay since Kobe Earthquake in 1994. In this paper is presented the analysis on the results of a series of lab tests on Osaka clay situated from 100 m to 1500 m ...Great interest has been aroused on deeply-situated Osaka clay since Kobe Earthquake in 1994. In this paper is presented the analysis on the results of a series of lab tests on Osaka clay situated from 100 m to 1500 m under the ground. The wave velocity method, bender element method, LDT and the formula derived by the authors are used, focus is put on the pre-failure mechanical behavior of the clay. The analysis shows that, (i) pore-pressure coefficient B is less than 1.0, (ii) the relationship between shear modulus and Poisson's ratio is not in agreement with that reported before, (iii) the modulus measured with LDT is still less than that measured with bender element method, and (iv) there are two threshold strains, within which the clay can be considered as elastic, and both of them are larger than that reported before.展开更多
Predicting the failure time of a landslide is considered as challenging work in the field of landslide research,and inverse velocity is proved to be an effective and convenient method.The onset of acceleration(OOA)has...Predicting the failure time of a landslide is considered as challenging work in the field of landslide research,and inverse velocity is proved to be an effective and convenient method.The onset of acceleration(OOA)has a crucial effect on the prediction failure time from the inverse velocity method.However,a simple method to identify OOA points is lacked,and most of the identifications rely on expert experience.Therefore,this study presents an application of a simple framework developed to identify the OOA by analyzing monitoring velocity data in three steps,including selection of the absolute value of velocity,reliable area identification and OOA identification.A new parameter based on exponential moving average(EMA)is developed to identify the landslide OOA.The framework is applied to three historical case studies to test its practicability and effectiveness.The forecasting results show a good correspondence between the accuracy rate and the coefficient of determination(R2).The predicted failure time according to the linear extrapolation starting from the identified OOA points is acceptable with a high R2 and high accuracy.展开更多
This paper focus on the accuracy enhancement of parallel kinematics machine through kinematics calibration. In the calibration processing, well-structured identification Jacobian matrix construction and end-effector p...This paper focus on the accuracy enhancement of parallel kinematics machine through kinematics calibration. In the calibration processing, well-structured identification Jacobian matrix construction and end-effector position and orientation measurement are two main difficulties. In this paper, the identification Jacobian matrix is constructed easily by numerical calculation utilizing the unit virtual velocity method. The generalized distance errors model is presented for avoiding measuring the position and orientation directly which is difficult to be measured. At last, a measurement tool is given for acquiring the data points in the calibration processing. Experimental studies confirmed the effectiveness of method. It is also shown in the paper that the proposed approach can be applied to other typed parallel manipulators.展开更多
A photonic lantem is studied which is formed by seven multimode fibers inserted into a pure silica capillary tube. The core of the tapered end has a uniform refractive index because the polymer claddings are removed b...A photonic lantem is studied which is formed by seven multimode fibers inserted into a pure silica capillary tube. The core of the tapered end has a uniform refractive index because the polymer claddings are removed before the fibers are inserted. Consequently, the light distribution is also uniform.Two theories describing a slowly varying waveguide and multimode coupling are used to analyze the photonic lantern. The transmission loss decreases as the length of the tapered part increases. For a device with a taper length of 3.4 cm, the loss is about 1.06 dB on average for light propagating through the taper from an inserted fiber to the tapered end and 0.99 dB in the reverse direction. For a device with a taper length of 0.7 cm, the two loss values are 2.63 dB and 2.53 dB, respectively. The results show that it is possible to achieve a uniform light distribution with the tapered end and a low-loss transmission in the device if parameters related to the lantern are reasonably defined.展开更多
The alignment accuracy of the strap-down inertial navigation system （SINS） of airborne weapon is greatly degraded by the dynamic wing flexure of the aircraft. An adaptive Kalman filter uses innovation sequences base...The alignment accuracy of the strap-down inertial navigation system （SINS） of airborne weapon is greatly degraded by the dynamic wing flexure of the aircraft. An adaptive Kalman filter uses innovation sequences based on the maximum likelihood estimated criterion to adapt the system noise covariance matrix and the measurement noise covariance matrix on line, which is used to estimate the misalignment if the model of wing flexure of the aircraft is unknown. From a number of simulations, it is shown that the accuracy of the adaptive Kalman filter is better than the conventional Kalman filter, and the erroneous misalignment models of the wing flexure of aircraft will cause bad estimation results of Kalman filter using attitude match method.展开更多
Submarine landslides can cause severe damage to marine engineering structures. Their sliding velocity and runout distance are two major parameters for quantifying and analyzing the risk of submarine landslides.Current...Submarine landslides can cause severe damage to marine engineering structures. Their sliding velocity and runout distance are two major parameters for quantifying and analyzing the risk of submarine landslides.Currently, commercial calculation programs such as BING have limitations in simulating underwater soil movements. All of these processes can be consistently simulated through a smoothed particle hydrodynamics（SPH） depth integrated model. The basis of the model is a control equation that was developed to take into account the effects of soil consolidation and erosion. In this work, the frictional rheological mode has been used to perform a simulation study of submarine landslides. Time-history curves of the sliding body＇s velocity, height,and length under various conditions of water depth, slope gradient, contact friction coefficient, and erosion rate are compared; the maximum sliding distance and velocity are calculated; and patterns of variation are discussed.The findings of this study can provide a reference for disaster warnings and pipeline route selection.展开更多
The effect of diffraction on the ultrasonic velocity measured by the pulse interference method has been investigated in VHF range theoretically and experimentally. Two silicate glasses are taken as the specimens, thei...The effect of diffraction on the ultrasonic velocity measured by the pulse interference method has been investigated in VHF range theoretically and experimentally. Two silicate glasses are taken as the specimens, their frequency dependences of longitudinal velocities are measured in the frequency range of 50-350 MHz, and the phase advances of ultrasonic signals caused by diffraction effect are calculated using A. O. Williams' theoretical expression. For the velocity error due to diffraction effect, the experimental results are in good agreement with the theoretical prediction. It has been shown that the velocity error due to diffraction effect is directly proportional to dθ21 (f)/df, whereθ21 (f) is the phase advances difference between the two partial reflection signals used in velocity measurement and f is the ultrasonic frequency.展开更多
The risks and uncertainties related to the storage of high-level radioactive waste （HLRW） can be reducedthanks to focused studies and investigations. HLRWs are going to be placed in deep geological repositories,enve...The risks and uncertainties related to the storage of high-level radioactive waste （HLRW） can be reducedthanks to focused studies and investigations. HLRWs are going to be placed in deep geological repositories,enveloped in an engineered bentonite barrier, whose physical conditions are subjected tochange throughout the lifespan of the infrastructure. Seismic tomography can be employed to monitor itsphysical state and integrity. The design of the seismic monitoring system can be optimized via conductingand analyzing numerical simulations of wave propagation in representative repository geometry.However, the quality of the numerical results relies on their initial calibration. The main aim of this paperis to provide a workflow to calibrate numerical tools employing laboratory ultrasonic datasets. The finitedifference code SOFI2D was employed to model ultrasonic waves propagating through a laboratorysample. Specifically, the input velocity model was calibrated to achieve a best match between experimentaland numerical ultrasonic traces. Likely due to the imperfections of the contact surfaces, theresultant velocities of P- and S-wave propagation tend to be noticeably lower than those a prioriassigned. Then, the calibrated model was employed to estimate the attenuation in a montmorillonitesample. The obtained low quality factors （Q） suggest that pronounced inelastic behavior of the clay has tobe taken into account in geophysical modeling and analysis. Consequently, this contribution should beconsidered as a first step towards the creation of a numerical tool to evaluate wave propagation innuclear waste repositories. 2014 Institute of Rock and Soil Mechanics, Chinese Academy of Sciences. Production and hosting byElsevier B.V. All rights reserved.展开更多
Predicting the failure time of unstable slopes is one of the most pivotal issues. In this paper, the inverse square root acceleration(INSRA) method was proposed to estimate the time-of-failure(TOF) of landslides. Four...Predicting the failure time of unstable slopes is one of the most pivotal issues. In this paper, the inverse square root acceleration(INSRA) method was proposed to estimate the time-of-failure(TOF) of landslides. Four collapsed slopes were presented in the three open-pit mines, two of them were probed by ground-based radar, and two of them were obtained from previous scientific papers. The inverse velocity(INV) method and INSRA method were adopted to analyze these four landslides and one slope which had great deformation but did not reach failure. Compared with the traditional INV method, the INSRA method can promote the forecasting effectiveness and has the advantage of higher accuracy.展开更多
The solution of complex rarefied flows with the BGK equation and the Discrete Velocity Method(DVM)requires a large number of velocity grid points leading to significant computational costs.We propose an adaptive veloc...The solution of complex rarefied flows with the BGK equation and the Discrete Velocity Method(DVM)requires a large number of velocity grid points leading to significant computational costs.We propose an adaptive velocity grid approach exploiting the fact that locally in space,the distribution function is supported only by a sub-set of the global velocity grid.The velocity grid is adapted thanks to criteria based on local temperature,velocity and on the enforcement of mass conservation.Simulations in 1D and 2D are presented for different Knudsen numbers and compared to a global velocity grid BGK solution,showing the computational gain of the proposed approach.展开更多
In most of partially averaged Navier-Stokes（PANS）methods,the Reynolds stress is solved by a linear hypothesis isotropic model.They could not capture all kinds of vortexes in tubomachineries.In this paper,a PANS mode...In most of partially averaged Navier-Stokes（PANS）methods,the Reynolds stress is solved by a linear hypothesis isotropic model.They could not capture all kinds of vortexes in tubomachineries.In this paper,a PANS model is modified from the RNG k-？turbulence model and is used to investigate the influence of the nonlinear shear stress on the simulation of the high pressure gradient flows and the large curvature flows.Comparisons are made between the result obtained by using the PANS model modified from the RNG k-？model and that obtained by using the nonlinear PANS methods.The flow past a curved rectangular duct is calculated by using the PANS methods.The obtained nonlinear shear stress agrees well with the experimental results,especially in the high pressure gradient region.The calculation results show that the nonlinear PANS methods are more reliable than the linear PANS methods for the high pressure gradient flows,the large curvature flows,and they can be used to capture complex vortexes in a turbomachinary.展开更多
基金This work was supported by The National Basic Research Program of China(Grant No.2007CB714600)Funding of Jiangsu Innovation Program for Graduate Education(Grant No.CXLX110170).
文摘Based on the field velocity method,a novel approach for simulating unsteady pitching and plunging motion of an airfoil is presented in this paper.Responses to pitching and plunging motions of the airfoil are simulated under different conditions.The obtained results are compared with those of moving grid method and good agreement is achieved.In the conventional field velocity method,the Euler solver is usually used to simulate the movement of the airfoil.However,when viscous effect is considered,unsteady Navier-Stokes equations have to be solved and the viscous flux correction must be taken into account.In this work,the viscous flux correction is introduced into the conventional field velocity method when non-uniform grid speed distribution is occurred.Numerical experiments for the flow around NACA0012 airfoil showed that the proposed approach can well simulate viscous moving boundary flow problems.
文摘There is a common difficulty in elastic-plastic impact codes such as EPIC[2,3] NONSAP[4], etc.. Most of these codes use the simple linear functions usually taken from static problem to represent the displacement components. In such finite element formulation, the stress components are constant in each element and they are discontinuous in any two neighboring elements. Therefore, the bases of using the virtual work principle in such elements are unreliable. In this paper, we introduce a new method, namely, the compatible stress iterative method, to eliminate the above-said difficulty. The calculated examples show that the calculation using the new method in dynamic finite element analysis of high velocity impact is valid and stable, and the element stiffness can be somewhat reduced.
文摘We apply the method of guidance by a required velocity for solving the optimal control problem over spacecraft’s reorientation from known initial attitude into a required final attitude.We suppose that attitude control is carried out by impulse jet engines.For optimization of fuel consumption,the controlling moments are calculated and formed according to the method of free trajectories together with principle of iterative control using the quaternions for generating commands to actuators.Optimal solution corresponds to the principle“acceleration-free rotation-separate corrections-free rotation-braking”.Rotation along a hitting trajectory is supported by insignificant correction of the uncontrolled motion at discrete instants between segments of acceleration and braking.Various strategies of forming the correction impulses during stage of free motion are suggested.Improving accuracy of achievement of spacecraft's final position is reached by terminal control using information about current attitude and angular velocity measurements for determining an instant of beginning of braking(condition for start of braking based on actual motion parameters is formulated in analytical form).The described method is universal and invariant relative to moments of inertia.Developed laws of attitude control concern the algorithms with prognostic model,the synthesized control modes are invariant with respect to both external perturbations and parametric errors.Results of mathematical modeling are presented that demonstrate practical feasibility and high efficiency of designed algorithms.
基金substantially supported by the Shuguang Program from Shanghai Education Development FoundationShanghai Municipal Education Commission, China (Grant No. 19SG19)+1 种基金National Natural Science Foundation of China (Grant No. 42072302)Fundamental Research Funds for the Central Universities, China
文摘The data-driven phenomenological models based on deformation measurements have been widely utilized to predict the slope failure time(SFT).The observational and model uncertainties could lead the predicted SFT calculated from the phenomenological models to deviate from the actual SFT.Currently,very limited study has been conducted on how to evaluate the effect of such uncertainties on SFT prediction.In this paper,a comprehensive slope failure database was compiled.A Bayesian machine learning(BML)-based method was developed to learn the model and observational uncertainties involved in SFT prediction,through which the probabilistic distribution of the SFT can be obtained.This method was illustrated in detail with an example.Verification studies show that the BML-based method is superior to the traditional inverse velocity method(INVM)and the maximum likelihood method for predicting SFT.The proposed method in this study provides an effective tool for SFT prediction.
文摘Cylindrical Cam Mechanism which is one of the best eq uipments to accomplish an accurate motion transmission is widely used in the fie lds of industries, such as machine tool exchangers, textile machinery and automa tic transfer equipments. This paper proposes a new approach for the shape design and manufacturing of the cylindrical cam. The design approach uses the relative velocity concept and the manufacturing approach uses the inverse kinematics concept. For the shape desig n, the contact points between the cam and the follower roller are calculated bas ed on relative velocity of which the direction is on the common tangential line, and then the whole shape of cam is determined from transformation of the coordi nate system. For the manufacturing procedures, the location and the orientation of cutter path can be allocated corresponding to the designed shape data. The in tegral NC code for multi-axis CNC machining center is created using the inverse kinematics concept from the data of the location and the orientation of cutter path. As the advantages of the proposed approach, the machine tool is designed t o having an alternative size in fabricating the general cam, while the tool must be fitted to diameter size of the follower in the conventional approach. Finally, CAD/CAM program, "Cylindrical DAM", is developed on C++ lan guage. This program can perform shape design, manufacturing and kinematics simul ation, which can make integral NC code for multi-axis CNC machining center. The proposed method can be applied easily on fields of industries.
基金Scientific Research Fund of the Institute of Engineering Mechanics,CEA under Grant No.2016B09,2017A02 and 2016A06the National Natural Science Foundation of China under Grant No,51378478,51408565,51678538 and 51161120360the National ScienceTechnology Support Plan Projects(2016YFC0701106)
文摘This paper employs a velocity plus displacement（V＋D）-based equivalent force control（EFC） method to solve the velocity/displacement difference equation in a real-time substructure test. This method uses type 2 feedback control loops to replace mathematical iteration to solve the nonlinear dynamic equation. A spectral radius analysis of the amplification matrix shows that the type 2 EFC-explicit, Newmark-β method has beneficial numerical characteristics for this method. Its stability limit of Ω = 2 remains unchanged regardless of the system damping because the velocity is achieved with very high accuracy during simulation. In contrast, the stability limits of the central difference method using direct velocity prediction and the EFC-average acceleration method with linear interpolation are shown to decrease with an increase in system damping. In fact, the EFC-average acceleration method is shown to change from unconditionally stable to conditionally stable. We also show that if an over-damped system with a damping ratio of 1.05 is considered, the stability limit is reduced to Ω =1.45. Finally, the results from an experiment with a single-degree-of-freedom structure installed with a magneto-rheological（MR） damper are presented. The results demonstrate that the proposed method is able to follow both displacement and velocity commands with moderate accuracy, resulting in improved test performance and accuracy for structures that are sensitive to both velocity and displacement inputs. Although the findings of the study are promising, additional test data and several further improvements will be required to draw general conclusions.
基金The National Science Foundation of China under contract Nos 40906017 and 41376038the National "863" Project of China under contract No.2013AA09A506+1 种基金the National Key Scientific Research Projects of China under contract No.2012CB955601the Special Projects on Public Sector under contract Nos 200905024 and 201409089
文摘A time splitting technique is common to many free surface ocean models. The different truncation errors in the equations of the internal and external modes require a numerical adjustment to make sure that algorithms correctly satisfy continuity equations and conserve tracers quantities. The princeton ocean model （POM） has applied a simple method of adjusting the vertical mean of internal velocities to external velocities at each internal time step. However, due to the Asselin time filter method adopted to prevent the numerical instability, the method of velocity adjustment used in POM can no longer guarantee the satisfaction of the continuity equation in the internal mode, though a special treatment is used to relate the surface elevation of the internal mode with that of the external mode. The error is proved to be a second-order term of the coefficient in the Asselin filter. One influence of this error in the numerical model is the failure of the kinetic boundary condition at the sea floor. By a regional experiment and a quasi-global experiment, the magni- tudes of this error are evaluated, and several sensitivity tests of this error are performed. The characteristic of this error is analyzed and two alternative algorithms are suggested to reduce the error.
文摘An improved algorithm for velocity field of general configurations ispresentd for low-order panel method based on the internal Dirichlet boundary condi-tion. A direct calculating method for the velocity distribution by means of a limit pro-cess combining with analytic evaluation of higher-order singular integrals instead of theconventional method of doublet strength gradient is devised in order to avoid the diffi-culty of edge extrapolation of doublet strength. The problem of substantialunderpredictions of the induced drag coefficient obtained from the VSAERO analysisdisappears for the present improved algorithm. Illustrative calculations for several testcases such as swept back wing, swept forward wing and wing-body combination showthat the accuracy of results may be improved and is competitive with high-order panelmethod. In addition, the present direct integral method can be used to evaluate the ve-locity distribution for external flow field correctly, where the method of gradient cannot be used at all.
基金funded by the Earthquake Science Fund of Xinjiang (201404)Seismic Tracing Oriented Task of China Earthquake Administration(2015010122)
文摘Based on the seismic observation report data provided by the Xinjiang Digital Seismic Network from 2009 to 2014,we calculate the wave velocity ratio and its background value for medium and small earthquakes by using the multi-station method in Tianshan,Xinjiang.This paper analyzes the variation of the wave velocity ratio disturbance value to highlight the abnormal,and also back-traces 7 moderate earthquakes at the research area.The results show that:(1)the background value of the wave velocity ratio is almost 1.70,the wave velocity ratio obviously decreases in the middle-eastern part of Tianshan and the region near the Puchang fault;(2)the wave velocity ratio disturbance value is mostly low in the epicenter before four earthquakes of M≥5.0 from 2011 to 2013 in the study area;(3)before 7 moderate strong earthquakes,the earthquake events with low value of the wave velocity ratio account for over 60% of corresponding total events near the epicenters,and the low value of the wave velocity ratio is relatively obvious before moderate earthquakes.
文摘Great interest has been aroused on deeply-situated Osaka clay since Kobe Earthquake in 1994. In this paper is presented the analysis on the results of a series of lab tests on Osaka clay situated from 100 m to 1500 m under the ground. The wave velocity method, bender element method, LDT and the formula derived by the authors are used, focus is put on the pre-failure mechanical behavior of the clay. The analysis shows that, (i) pore-pressure coefficient B is less than 1.0, (ii) the relationship between shear modulus and Poisson's ratio is not in agreement with that reported before, (iii) the modulus measured with LDT is still less than that measured with bender element method, and (iv) there are two threshold strains, within which the clay can be considered as elastic, and both of them are larger than that reported before.
基金funded by the National Natural Science Foundation of China(Grant NO.41772324)the Open Foundation of Chengdu Center of China Geological Survey。
文摘Predicting the failure time of a landslide is considered as challenging work in the field of landslide research,and inverse velocity is proved to be an effective and convenient method.The onset of acceleration(OOA)has a crucial effect on the prediction failure time from the inverse velocity method.However,a simple method to identify OOA points is lacked,and most of the identifications rely on expert experience.Therefore,this study presents an application of a simple framework developed to identify the OOA by analyzing monitoring velocity data in three steps,including selection of the absolute value of velocity,reliable area identification and OOA identification.A new parameter based on exponential moving average(EMA)is developed to identify the landslide OOA.The framework is applied to three historical case studies to test its practicability and effectiveness.The forecasting results show a good correspondence between the accuracy rate and the coefficient of determination(R2).The predicted failure time according to the linear extrapolation starting from the identified OOA points is acceptable with a high R2 and high accuracy.
文摘This paper focus on the accuracy enhancement of parallel kinematics machine through kinematics calibration. In the calibration processing, well-structured identification Jacobian matrix construction and end-effector position and orientation measurement are two main difficulties. In this paper, the identification Jacobian matrix is constructed easily by numerical calculation utilizing the unit virtual velocity method. The generalized distance errors model is presented for avoiding measuring the position and orientation directly which is difficult to be measured. At last, a measurement tool is given for acquiring the data points in the calibration processing. Experimental studies confirmed the effectiveness of method. It is also shown in the paper that the proposed approach can be applied to other typed parallel manipulators.
基金supported by the National Natural Science Foundation of China (Grant Nos. 11078009,U1331114,61107059 and 61307076)the 111 project (B13015) to Harbin Engineering University+1 种基金the Fundamental Research Funds for the Central Universitiesthe Opening Project of Key Laboratory of Astronomical Optics and Technology,Nanjing Institute of Astronomical Optics and Technology,Chinese Academy of Sciences
文摘A photonic lantem is studied which is formed by seven multimode fibers inserted into a pure silica capillary tube. The core of the tapered end has a uniform refractive index because the polymer claddings are removed before the fibers are inserted. Consequently, the light distribution is also uniform.Two theories describing a slowly varying waveguide and multimode coupling are used to analyze the photonic lantern. The transmission loss decreases as the length of the tapered part increases. For a device with a taper length of 3.4 cm, the loss is about 1.06 dB on average for light propagating through the taper from an inserted fiber to the tapered end and 0.99 dB in the reverse direction. For a device with a taper length of 0.7 cm, the two loss values are 2.63 dB and 2.53 dB, respectively. The results show that it is possible to achieve a uniform light distribution with the tapered end and a low-loss transmission in the device if parameters related to the lantern are reasonably defined.
文摘The alignment accuracy of the strap-down inertial navigation system （SINS） of airborne weapon is greatly degraded by the dynamic wing flexure of the aircraft. An adaptive Kalman filter uses innovation sequences based on the maximum likelihood estimated criterion to adapt the system noise covariance matrix and the measurement noise covariance matrix on line, which is used to estimate the misalignment if the model of wing flexure of the aircraft is unknown. From a number of simulations, it is shown that the accuracy of the adaptive Kalman filter is better than the conventional Kalman filter, and the erroneous misalignment models of the wing flexure of aircraft will cause bad estimation results of Kalman filter using attitude match method.
基金The Specialized Research Fund for the Doctoral Program of Higher Education under contract No.20120041130002the National Key Project of Science and Technology under contract No.2011ZX 05056-001-02the Fundamental Research Funds for the Central Universities under contract No.DUT14ZD220
文摘Submarine landslides can cause severe damage to marine engineering structures. Their sliding velocity and runout distance are two major parameters for quantifying and analyzing the risk of submarine landslides.Currently, commercial calculation programs such as BING have limitations in simulating underwater soil movements. All of these processes can be consistently simulated through a smoothed particle hydrodynamics（SPH） depth integrated model. The basis of the model is a control equation that was developed to take into account the effects of soil consolidation and erosion. In this work, the frictional rheological mode has been used to perform a simulation study of submarine landslides. Time-history curves of the sliding body＇s velocity, height,and length under various conditions of water depth, slope gradient, contact friction coefficient, and erosion rate are compared; the maximum sliding distance and velocity are calculated; and patterns of variation are discussed.The findings of this study can provide a reference for disaster warnings and pipeline route selection.
文摘The effect of diffraction on the ultrasonic velocity measured by the pulse interference method has been investigated in VHF range theoretically and experimentally. Two silicate glasses are taken as the specimens, their frequency dependences of longitudinal velocities are measured in the frequency range of 50-350 MHz, and the phase advances of ultrasonic signals caused by diffraction effect are calculated using A. O. Williams' theoretical expression. For the velocity error due to diffraction effect, the experimental results are in good agreement with the theoretical prediction. It has been shown that the velocity error due to diffraction effect is directly proportional to dθ21 (f)/df, whereθ21 (f) is the phase advances difference between the two partial reflection signals used in velocity measurement and f is the ultrasonic frequency.
文摘The risks and uncertainties related to the storage of high-level radioactive waste （HLRW） can be reducedthanks to focused studies and investigations. HLRWs are going to be placed in deep geological repositories,enveloped in an engineered bentonite barrier, whose physical conditions are subjected tochange throughout the lifespan of the infrastructure. Seismic tomography can be employed to monitor itsphysical state and integrity. The design of the seismic monitoring system can be optimized via conductingand analyzing numerical simulations of wave propagation in representative repository geometry.However, the quality of the numerical results relies on their initial calibration. The main aim of this paperis to provide a workflow to calibrate numerical tools employing laboratory ultrasonic datasets. The finitedifference code SOFI2D was employed to model ultrasonic waves propagating through a laboratorysample. Specifically, the input velocity model was calibrated to achieve a best match between experimentaland numerical ultrasonic traces. Likely due to the imperfections of the contact surfaces, theresultant velocities of P- and S-wave propagation tend to be noticeably lower than those a prioriassigned. Then, the calibrated model was employed to estimate the attenuation in a montmorillonitesample. The obtained low quality factors （Q） suggest that pronounced inelastic behavior of the clay has tobe taken into account in geophysical modeling and analysis. Consequently, this contribution should beconsidered as a first step towards the creation of a numerical tool to evaluate wave propagation innuclear waste repositories. 2014 Institute of Rock and Soil Mechanics, Chinese Academy of Sciences. Production and hosting byElsevier B.V. All rights reserved.
基金supported by the National Natural Science Foundation of China (Grant Nos.51839009 and 51679017)。
文摘Predicting the failure time of unstable slopes is one of the most pivotal issues. In this paper, the inverse square root acceleration(INSRA) method was proposed to estimate the time-of-failure(TOF) of landslides. Four collapsed slopes were presented in the three open-pit mines, two of them were probed by ground-based radar, and two of them were obtained from previous scientific papers. The inverse velocity(INV) method and INSRA method were adopted to analyze these four landslides and one slope which had great deformation but did not reach failure. Compared with the traditional INV method, the INSRA method can promote the forecasting effectiveness and has the advantage of higher accuracy.
基金the French National Research Agency(ANR)in the frame of the“Investments for the future”Programme IdEx Bordeaux(ANR-10-IDEX-03-02),Cluster of excellence CPU.“National Group for Scientific Computation(GNCS-INDAM).
文摘The solution of complex rarefied flows with the BGK equation and the Discrete Velocity Method(DVM)requires a large number of velocity grid points leading to significant computational costs.We propose an adaptive velocity grid approach exploiting the fact that locally in space,the distribution function is supported only by a sub-set of the global velocity grid.The velocity grid is adapted thanks to criteria based on local temperature,velocity and on the enforcement of mass conservation.Simulations in 1D and 2D are presented for different Knudsen numbers and compared to a global velocity grid BGK solution,showing the computational gain of the proposed approach.
基金supported by the National Natural Science Foundation of China(Grant Nos.51406010,51479166)
文摘In most of partially averaged Navier-Stokes（PANS）methods,the Reynolds stress is solved by a linear hypothesis isotropic model.They could not capture all kinds of vortexes in tubomachineries.In this paper,a PANS model is modified from the RNG k-？turbulence model and is used to investigate the influence of the nonlinear shear stress on the simulation of the high pressure gradient flows and the large curvature flows.Comparisons are made between the result obtained by using the PANS model modified from the RNG k-？model and that obtained by using the nonlinear PANS methods.The flow past a curved rectangular duct is calculated by using the PANS methods.The obtained nonlinear shear stress agrees well with the experimental results,especially in the high pressure gradient region.The calculation results show that the nonlinear PANS methods are more reliable than the linear PANS methods for the high pressure gradient flows,the large curvature flows,and they can be used to capture complex vortexes in a turbomachinary.
