We investigated correlation of the slope of the initial part of P-wave envelope, hypocenter depth and plate boundaries by B-Δ method, which is used to determine epicentral distances in the Japan Meteorological Agency...We investigated correlation of the slope of the initial part of P-wave envelope, hypocenter depth and plate boundaries by B-Δ method, which is used to determine epicentral distances in the Japan Meteorological Agency and Japan Railway Company earthquake early warning(EEW) systems. The Tohoku region was chosen as the study region. 19,899 strong motion data for 265 events with magnitudes in the range from 5.0 to 7.6 from KiK-net(Kiban Kyoshin network) had been collected. The coefficient c to investigate is obtained from the linear relation between log B and log Δ. Compared to the hypocenter depth, the coefficients c of events is more likely to decide by the spatial correlation of the plate boundaries. The differences are likely to be due to earthquake characteristics, since deeper events in the subducting slabs the structural effects are likely to be larger than or comparable to those for shallow crust events.展开更多
The recently built China Digital Seismic Network consists of the China National Digital Seismic Network (CNDSN), 31 regional seismic networks and several small aperture arrays with more than 1 000 stations including...The recently built China Digital Seismic Network consists of the China National Digital Seismic Network (CNDSN), 31 regional seismic networks and several small aperture arrays with more than 1 000 stations including 850+ broadband stations. It forms a gigantic seismic array that provides an unprecedented opportunity to study the Earth's deep interior besides its routine task of seismic monitoring. Many modern seismic studies rely on rotation of vertical and horizontal components in order to separate different types of seismic waves. Knowledge of the orientations of the two horizontal components thus is important to perform a correction rotation. We analyzed particle motions of teleseismic P waves recorded by the network and used them to estimate the northcomponent azimuth of each station. An SNR-weighted-multi-event method was introduced to obtain component azimuths that best explain the P-wave particle motions of all the events recorded at a station. The method provides robust estimates including a measurement error calculated from background noise levels. We found that about one third of the stations have some sort of problems, including misorientation of the two horizontal components, mislabeling and polarity reversal in one or more components. These problems need to be taken into account for any rotation based seismic studies.展开更多
In this note we consider some basic, yet unusual, issues pertaining to the accuracy and stability of numerical integration methods to follow the solution of first order and second order initial value problems (IVP). I...In this note we consider some basic, yet unusual, issues pertaining to the accuracy and stability of numerical integration methods to follow the solution of first order and second order initial value problems (IVP). Included are remarks on multiple solutions, multi-step methods, effect of initial value perturbations, as well as slowing and advancing the computed motion in second order problems.展开更多
The present paper contributes in studying the phase velocities of P- and S-waves in a half space subjected to a compressive initial stress and gravity field. The density and acceleration due to gravity vary quadratica...The present paper contributes in studying the phase velocities of P- and S-waves in a half space subjected to a compressive initial stress and gravity field. The density and acceleration due to gravity vary quadratically along the depth. The dispersion equation is derived in a closed form. It is shown that the phase velocities depend not only on the initial stress, gravity, and direction of propagation but also on the inhomogeneity parameter associated with the density and acceleration due to gravity. Various particular cases are obtained, and the results match with the classical results. Numerical investigations on the phase velocities of P- and S-waves against the wave number are made for various sets of values of the material parameters, and the results are illustrated graphically. The graphical user interface model is developed to generalize the effect.展开更多
In this paper,we present a novel initial costates solver for initializing time-optimal trajectory problems in relative motion with continuous low thrust.The proposed solver consists of two primary components:training ...In this paper,we present a novel initial costates solver for initializing time-optimal trajectory problems in relative motion with continuous low thrust.The proposed solver consists of two primary components:training a Multilayer Perceptron(MLP)for generating reference sequence and Time of Flight(TOF)to the target,and deriving a system of linear algebraic equations for obtaining the initial costates.To overcome the challenge of generating training samples for the MLP,the backward generation method is proposed to obtain five different training databases.The training database and sample form are determined by analyzing the input and output correlation using the Pearson correlation coefficient.The best-performing MLP is obtained by analyzing the training results with various hyper-parameter combinations.A reference sequence starting from the initial states is obtained by integrating forward with the near-optimal control vector from the output of MLP.Finally,a system of linear algebraic equations for estimating the initial costates is derived using the reference sequence and the necessary conditions for optimality.Simulation results demonstrate that the proposed initial costates solver improves the convergence ratio and reduce the function calls of the shooting function.Furthermore,Monte-Carlo simulation illustrates that the initial costates solver is applicable to different initial velocities,demonstrating excellent generalization ability.展开更多
Summarized in the paper are the author's studies on incipient motion of sediment in recent 40 years. In addition to the forces of gravity, drag and lift, the cohesive force and the additional static pressure are i...Summarized in the paper are the author's studies on incipient motion of sediment in recent 40 years. In addition to the forces of gravity, drag and lift, the cohesive force and the additional static pressure are important for fine particles. The relations between three stages of incipient motion are defined by use of instantaneous velocity. Formulas for initial velocity and critical shear stress are given and overall verified by the author's and others' experimental data.展开更多
Constrained modeling and state estimation have attracted much attention in recent years. This paper focuses on target motion modeling and tracking in road coordinates. An improved initialization method,which uses the ...Constrained modeling and state estimation have attracted much attention in recent years. This paper focuses on target motion modeling and tracking in road coordinates. An improved initialization method,which uses the optimal fusion of the position measurements in different directions,is presented for the constraint coordinate Kalman filter(CCKF). The CCKF is evaluated with a comprehensive comparison to the state-of-art linear equality constraint estimation methods. Numerical simulation results demonstrate the better performance of the CCKF. Then the interacting multiple model CCKF(IMM-CCKF) is proposed to manifest the advantages of the CCKF in complex motion modeling and state estimations. The effectiveness of the IMM-CCKF in maneuvering target tracking with spatial equality constraints is demonstrated by numerical experiments.展开更多
Based on the kinematic relationship of the drill, drill point and workpiece and the coordinate systems founded,the mathematical models of the drill point wandering motion for the drill initial penetration, i. e. the f...Based on the kinematic relationship of the drill, drill point and workpiece and the coordinate systems founded,the mathematical models of the drill point wandering motion for the drill initial penetration, i. e. the formulae for the drill point center locus, hole shape and cutting edge locus, are established. According to these formulae the hole shape to be made is simulated and analyzed. It is verified that the bending chatter of the drill during the drill penetrating is the main cause to form a polygon like hole.展开更多
Harmonic suppression, non-periodic and non-closing in straightness profile error that will bring about harmonic component distortion in measurement result are analyzed. The countermeasure-a novel accurate two-probe me...Harmonic suppression, non-periodic and non-closing in straightness profile error that will bring about harmonic component distortion in measurement result are analyzed. The countermeasure-a novel accurate two-probe method in time domain is put forward to measure straight-going component motion error in machine tools based on the frequency domain 3-point method after symmetrical continuation of probes' primitive signal. Both straight-going component motion error in machine tools and the profile error in workpiece that is manufactured on this machine can be measured at the same time. The information is available to diagnose the fault origin of machine tools. The analysis result is proved to be correct by the experiment.展开更多
针对铁道车辆蛇行运动分岔图中不稳定极限环通过整车动力学仿真求解困难的问题,提出直接积分求解不稳定极限环的方法——初态激扰法。该方法利用多体系统动力学软件Simpack建立高速列车动力学模型,采用Matlab软件对拟周期解进行动态加...针对铁道车辆蛇行运动分岔图中不稳定极限环通过整车动力学仿真求解困难的问题,提出直接积分求解不稳定极限环的方法——初态激扰法。该方法利用多体系统动力学软件Simpack建立高速列车动力学模型,采用Matlab软件对拟周期解进行动态加权调整后作为车辆系统主要部件的初始状态,在光滑轨道上时域积分获取刚体运动状态,在轮对横移幅值随运行速度变化的分岔图中绘制不稳定极限环及平衡点和稳定极限环,从而得到完整的车辆蛇行运动分岔图。以某高速列车为例,基于初态激扰法求解不同轮轨接触工况和抗蛇行减振器故障工况下蛇行运动分岔曲线。结果表明:新轮和磨耗轮工况的车辆蛇行运动分别对应Hopf亚临界分岔和超临界分岔行为,且磨耗轮工况下蛇行运动由亚临界分岔变为超临界分岔;不改变抗蛇行减振器阻尼仅减小卸荷力,对车辆Hopf分岔临界速度没有影响,但会降低车辆LPC (Limit Point Bifurcation of Circles)分岔临界速度并减小不稳定极限环的幅值,从而降低车辆横向稳定性。展开更多
基金supported by the Director Foundation of the Institute of Seismology,China Earthquake Administration,China(IS200756046)
文摘We investigated correlation of the slope of the initial part of P-wave envelope, hypocenter depth and plate boundaries by B-Δ method, which is used to determine epicentral distances in the Japan Meteorological Agency and Japan Railway Company earthquake early warning(EEW) systems. The Tohoku region was chosen as the study region. 19,899 strong motion data for 265 events with magnitudes in the range from 5.0 to 7.6 from KiK-net(Kiban Kyoshin network) had been collected. The coefficient c to investigate is obtained from the linear relation between log B and log Δ. Compared to the hypocenter depth, the coefficients c of events is more likely to decide by the spatial correlation of the plate boundaries. The differences are likely to be due to earthquake characteristics, since deeper events in the subducting slabs the structural effects are likely to be larger than or comparable to those for shallow crust events.
基金supported by NSF grant EAR-063566(F.N.)National Natural Science Foundation of China grant 40774042(J.L.)
文摘The recently built China Digital Seismic Network consists of the China National Digital Seismic Network (CNDSN), 31 regional seismic networks and several small aperture arrays with more than 1 000 stations including 850+ broadband stations. It forms a gigantic seismic array that provides an unprecedented opportunity to study the Earth's deep interior besides its routine task of seismic monitoring. Many modern seismic studies rely on rotation of vertical and horizontal components in order to separate different types of seismic waves. Knowledge of the orientations of the two horizontal components thus is important to perform a correction rotation. We analyzed particle motions of teleseismic P waves recorded by the network and used them to estimate the northcomponent azimuth of each station. An SNR-weighted-multi-event method was introduced to obtain component azimuths that best explain the P-wave particle motions of all the events recorded at a station. The method provides robust estimates including a measurement error calculated from background noise levels. We found that about one third of the stations have some sort of problems, including misorientation of the two horizontal components, mislabeling and polarity reversal in one or more components. These problems need to be taken into account for any rotation based seismic studies.
文摘In this note we consider some basic, yet unusual, issues pertaining to the accuracy and stability of numerical integration methods to follow the solution of first order and second order initial value problems (IVP). Included are remarks on multiple solutions, multi-step methods, effect of initial value perturbations, as well as slowing and advancing the computed motion in second order problems.
基金supported by the Research Fellow of Indian School of Mines in Dhanbad (No. 2010DR0016)
文摘The present paper contributes in studying the phase velocities of P- and S-waves in a half space subjected to a compressive initial stress and gravity field. The density and acceleration due to gravity vary quadratically along the depth. The dispersion equation is derived in a closed form. It is shown that the phase velocities depend not only on the initial stress, gravity, and direction of propagation but also on the inhomogeneity parameter associated with the density and acceleration due to gravity. Various particular cases are obtained, and the results match with the classical results. Numerical investigations on the phase velocities of P- and S-waves against the wave number are made for various sets of values of the material parameters, and the results are illustrated graphically. The graphical user interface model is developed to generalize the effect.
基金This study was funded by the National Natural Science Foundation of China(Nos.11972077 and 12272039).
文摘In this paper,we present a novel initial costates solver for initializing time-optimal trajectory problems in relative motion with continuous low thrust.The proposed solver consists of two primary components:training a Multilayer Perceptron(MLP)for generating reference sequence and Time of Flight(TOF)to the target,and deriving a system of linear algebraic equations for obtaining the initial costates.To overcome the challenge of generating training samples for the MLP,the backward generation method is proposed to obtain five different training databases.The training database and sample form are determined by analyzing the input and output correlation using the Pearson correlation coefficient.The best-performing MLP is obtained by analyzing the training results with various hyper-parameter combinations.A reference sequence starting from the initial states is obtained by integrating forward with the near-optimal control vector from the output of MLP.Finally,a system of linear algebraic equations for estimating the initial costates is derived using the reference sequence and the necessary conditions for optimality.Simulation results demonstrate that the proposed initial costates solver improves the convergence ratio and reduce the function calls of the shooting function.Furthermore,Monte-Carlo simulation illustrates that the initial costates solver is applicable to different initial velocities,demonstrating excellent generalization ability.
文摘Summarized in the paper are the author's studies on incipient motion of sediment in recent 40 years. In addition to the forces of gravity, drag and lift, the cohesive force and the additional static pressure are important for fine particles. The relations between three stages of incipient motion are defined by use of instantaneous velocity. Formulas for initial velocity and critical shear stress are given and overall verified by the author's and others' experimental data.
基金Sponsored by the National Natural Science Foundation of China(Grant No.61201311)
文摘Constrained modeling and state estimation have attracted much attention in recent years. This paper focuses on target motion modeling and tracking in road coordinates. An improved initialization method,which uses the optimal fusion of the position measurements in different directions,is presented for the constraint coordinate Kalman filter(CCKF). The CCKF is evaluated with a comprehensive comparison to the state-of-art linear equality constraint estimation methods. Numerical simulation results demonstrate the better performance of the CCKF. Then the interacting multiple model CCKF(IMM-CCKF) is proposed to manifest the advantages of the CCKF in complex motion modeling and state estimations. The effectiveness of the IMM-CCKF in maneuvering target tracking with spatial equality constraints is demonstrated by numerical experiments.
文摘Based on the kinematic relationship of the drill, drill point and workpiece and the coordinate systems founded,the mathematical models of the drill point wandering motion for the drill initial penetration, i. e. the formulae for the drill point center locus, hole shape and cutting edge locus, are established. According to these formulae the hole shape to be made is simulated and analyzed. It is verified that the bending chatter of the drill during the drill penetrating is the main cause to form a polygon like hole.
基金National Nature Science Foundation of China.No.50075056
文摘Harmonic suppression, non-periodic and non-closing in straightness profile error that will bring about harmonic component distortion in measurement result are analyzed. The countermeasure-a novel accurate two-probe method in time domain is put forward to measure straight-going component motion error in machine tools based on the frequency domain 3-point method after symmetrical continuation of probes' primitive signal. Both straight-going component motion error in machine tools and the profile error in workpiece that is manufactured on this machine can be measured at the same time. The information is available to diagnose the fault origin of machine tools. The analysis result is proved to be correct by the experiment.
文摘针对铁道车辆蛇行运动分岔图中不稳定极限环通过整车动力学仿真求解困难的问题,提出直接积分求解不稳定极限环的方法——初态激扰法。该方法利用多体系统动力学软件Simpack建立高速列车动力学模型,采用Matlab软件对拟周期解进行动态加权调整后作为车辆系统主要部件的初始状态,在光滑轨道上时域积分获取刚体运动状态,在轮对横移幅值随运行速度变化的分岔图中绘制不稳定极限环及平衡点和稳定极限环,从而得到完整的车辆蛇行运动分岔图。以某高速列车为例,基于初态激扰法求解不同轮轨接触工况和抗蛇行减振器故障工况下蛇行运动分岔曲线。结果表明:新轮和磨耗轮工况的车辆蛇行运动分别对应Hopf亚临界分岔和超临界分岔行为,且磨耗轮工况下蛇行运动由亚临界分岔变为超临界分岔;不改变抗蛇行减振器阻尼仅减小卸荷力,对车辆Hopf分岔临界速度没有影响,但会降低车辆LPC (Limit Point Bifurcation of Circles)分岔临界速度并减小不稳定极限环的幅值,从而降低车辆横向稳定性。