斜坡应力分布的系统测试与反演分析

采用简易、快速、经济的洞壁应力恢复法对某工程边坡勘探平硐洞壁应力场作了系统测试, 并采用FLAC2D软件对测试成果作了斜坡初始应力的反演分析。结果表明, 斜坡应力场随深度变化具有明显分带性, 应力量值随硐深曲线具有典型的驼峰式特征, 与理论分析结果相符, 表明该方法能较好地应用于斜坡应力场的系统测试和研究。成果对认识斜坡岩体的变形破裂及卸荷特征等具有重要意义。     

Time-dependent stress change during failure of rocks

The energy balance of a solid subject to fracture has been explored using heat and mass transfer equations with regard to the volumetric and superficial components. In the suggested model, brittle fracture of a cracked solid considered as a heterogeneous two-phase medium is described by an equation analogous to the Griffith’s criterion for propagation of a single crack. The derived equation is used, together with estimates of relative change in specific interface area, to study the respective change of free strain energy and pressure in rocks associated with failure.     

西藏札达盆地控盆断裂有限元数值模拟

札达盆地是喜马拉雅构造带中的一个山间断陷盆地,其演化过程与盆地两侧的控盆构造密切相关。对控盆断裂的构造应力场进行模拟计算,将有助于进一步深化对本区构造控盆的认识。因此,在对盆地构造地质进行详细调查的基础上,结合本区的深部地质与地球物理资料,对札达盆地控盆断裂的构造应力场进行了模拟。计算结果表明,札达盆地的演化明显受盆地两侧边界断裂的控制,札达盆地是在整体南北向挤压应力的作用下,不同块体差异隆升作用的结果。其南侧的控盆断裂为北倾的正断层,北侧的控盆断裂为南倾的逆断层,二者共同形成了南降北升的翘板式断陷盆地运动过程,是喜马拉雅地块在陆内汇聚挤压构造环境中构造应力场调整的一种方式。     

On solving Kepler's equation for nearly parabolic orbits

We deal here with the efficient starting points for Kepler's equation in the special case of nearly parabolic orbits. Our approach provides with very simple formulas that allow calculating these points on a scientific vest-pocket calculator. Moreover, srtarting with these points in the Newton's method we can calculate a root of Kepler's equation with an accuracy greater than 0.001 in 0–2 iterations. This accuracy holds for the true anomaly || 135° and |e – 1| 0.01. We explain the reason for this effect also.Dedicated to the memory of Professor G.N. Duboshin (1903–1986).     

Extension of the solution of Kepler's equation to high eccentricities

The classic Lagrange's expansion of the solutionE(e, M) of Kepler's equation in powers of eccentricity is extended to highly eccentric orbits, 0.6627 ... <e<1. The solutionE(e, M) is developed in powers of (e–e*), wheree* is a fixed value of the eccentricity. The coefficients of the expansion are given in terms of the derivatives of the Bessel functionsJ _n (ne). The expansion is convergent for values of the eccentricity such that |e–e*|<(e*), where the radius of convergence (e*) is a positive real number, which is calculated numerically.     

On the physical model of earthquake precursor fields and the mechanism of precursors＇timespace distribution──origin and evidences of the strong body earthquake──generating model

Ｏｎｔｈｅｐｈｙｓｉｃａｌｍｏｄｅｌｏｆｅａｒｔｈｑｕａｋｅｐｒｅｃｕｒｓｏｒｆｉｅｌｄｓａｎｄｔｈｅｍｅｃｈａｎｉｓｍｏｆｐｒｅｃｕｒｓｏｒｓ＇ｔｉｍｅｓｐａｃｅｄｉｓｔｒｉｂｕｔｉｏｎ──ｏｒｉｇｉｎａｎｄｅｖｉｄｅｎｃｅｓｏｆ　ｔｈｅｓｔｒｏｎ...     

On the periodic solutions of a particular Lame equation

We consider the Hill's equation: % MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGGipm0dc9vqaqpepu0xbbG8F4rqqrFfpeea0xe9Lq-Jc9% vqaqpepm0xbba9pwe9Q8fs0-yqaqpepae9pg0FirpepeKkFr0xfr-x% fr-xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaWaaSaaaeaaca% WGKbWaaWbaaSqabeaacaaIYaaaaOGaeqOVdGhabaGaamizaiaadsha% daahaaWcbeqaaiaaikdaaaaaaOGaey4kaSYaaSaaaeaacaWGTbGaai% ikaiaad2gacqGHRaWkcaaIXaGaaiykaaqaaiaaikdaaaGaam4qamaa% CaaaleqabaGaaGOmaaaakiaacIcacaWG0bGaaiykaiabe67a4jabg2% da9iaaicdaaaa!4973!\[\frac{{d^2 \xi }}{{dt^2 }} + \frac{{m(m + 1)}}{2}C^2 (t)\xi = 0\]Where C(t) = Cn (t, {frbuilt|1/2}) is the elliptic function of Jacobi and m a given real number. It is a particular case of theame equation. By the change of variable from t to defined by: % MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGGipm0dc9vqaqpepu0xbbG8F4rqqrFfpeea0xe9Lq-Jc9% vqaqpepm0xbba9pwe9Q8fs0-yqaqpepae9pg0FirpepeKkFr0xfr-x% fr-xb9adbaqaaeGaciGaaiaabeqaamaabaabaaqcaawaaOWaaiqaaq% aabeqaamaalaaajaaybaGaamizaGGaaiab-z6agbqaaiaadsgacaWG% 0baaaiabg2da9OWaaOaaaKaaGfaacaGGOaqcKbaG-laaigdajaaycq% GHsislkmaaleaajeaybaGaaGymaaqaaiaaikdaaaqcaaMaaeiiaiaa% bohacaqGPbGaaeOBaOWaaWbaaKqaGfqabaGaaeOmaaaajaaycqWFMo% GrcqWFPaqkaKqaGfqaaaqcaawaaiab-z6agjab-HcaOiab-bdaWiab% -LcaPiab-1da9iab-bdaWaaakiaawUhaaaaa!51F5!\[\left\{ \begin{array}{l}\frac{{d\Phi }}{{dt}} = \sqrt {(1 - {\textstyle{1 \over 2}}{\rm{ sin}}^{\rm{2}} \Phi )} \\\Phi (0) = 0 \\\end{array} \right.\]it is transformed to the Ince equation: (1 + · cos(2)) y + b · sin(2) · y + (c + d · cos(2)) y = 0 where % MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGGipm0dc9vqaqpepu0xbbG8F4rqqrFfpeea0xe9Lq-Jc9% vqaqpepm0xbba9pwe9Q8fs0-yqaqpepae9pg0FirpepeKkFr0xfr-x% fr-xb9adbaqaaeGaciGaaiaabeqaamaabaabaaqcaawaaiaadggacq% GH9aqpcqGHsislcaWGIbGaeyypa0JcdaWcgaqaaiaaigdaaeaacaaI% ZaGaaiilaiaabccacaWGJbGaeyypa0Jaamizaiabg2da9aaacaqGGa% WaaSaaaKaaGfaacaWGTbGaaiikaiaad2gacqGHRaWkcaaIXaGaaiyk% aaqaaiaaiodaaaaaaa!4777!\[a = - b = {1 \mathord{\left/{\vphantom {1 {3,{\rm{ }}c = d = }}} \right.\kern-\nulldelimiterspace} {3,{\rm{ }}c = d = }}{\rm{ }}\frac{{m(m + 1)}}{3}\]In the neighbourhood of the poles, we give the expression of the solutions.The periodic solutions of the Equation (1) correspond to the periodic solutions of the Equation (3). Magnus and Winkler give us a theory of their existence. By comparing these results to those of our study in the case of the Hill's equation, we can find the development in Fourier series of periodic solutions in function of the variable and deduce the development of solutions of (1) in function of C(t).     

中国三个大地震的震源参数及讨论

大地震的震源机制及地震矩、应力降等参数,不(又对认识地震的破裂过程,而且对预报强地面运动都是非常重要的。一些文章对板内地震应力降的特点、地震矩与高频源谱的关系及区域特征的讨论,使我们可直接由长周期体波或面波得出的地震矩M₀估计出不同地区中大地震的高频源谱,以供强地面运动预报参考。 由于对发生在中国板内的1976年11月6日盐源地震（M_s=6.5）和1976年11月15日宁河地震（M_s=6.3）的震源参数还没有详尽的报道;另外,对1973年7月14日     

Large eddy simulation of hydrocyclone—prediction of air-core diameter and shape

Hydrocyclones are widely used in the mining and chemical industries. An attempt has been made in this study, to develop a CFD (computational fluid dynamics) model, which is capable of predicting the flow patterns inside the hydrocyclone, including accurate prediction of flow split as well as the size of the air-core. The flow velocities and air-core diameters are predicted by DRSM (differential Reynolds stress model) and LES (large eddy simulations) models were compared to experimental results. The predicted water splits and air-core diameter with LES and RSM turbulence models along with VOF (volume of fluid) model for the air phase, through the outlets for various inlet pressures were also analyzed. The LES turbulence model led to an improved turbulence field prediction and thereby to more accurate prediction of pressure and velocity fields. This improvement was distinctive for the axial profile of pressure, indicating that air-core development is principally a transport effect rather than a pressure effect.     

Increasing critical sensitivity of the Load/Unload Response Ratio before large earthquakes with identified stress accumulation pattern

The Load/Unload Response Ratio (LURR) method is proposed for short-to-intermediate-term earthquake prediction [Yin, X.C., Chen, X.Z., Song, Z.P., Yin, C., 1995. A New Approach to Earthquake Prediction — The Load/Unload Response Ratio (LURR) Theory, Pure Appl. Geophys., 145, 701–715]. This method is based on measuring the ratio between Benioff strains released during the time periods of loading and unloading, corresponding to the Coulomb Failure Stress change induced by Earth tides on optimally oriented faults. According to the method, the LURR time series usually climb to an anomalously high peak prior to occurrence of a large earthquake. Previous studies have indicated that the size of critical seismogenic region selected for LURR measurements has great influence on the evaluation of LURR. In this study, we replace the circular region usually adopted in LURR practice with an area within which the tectonic stress change would mostly affect the Coulomb stress on a potential seismogenic fault of a future event. The Coulomb stress change before a hypothetical earthquake is calculated based on a simple back-slip dislocation model of the event. This new algorithm, by combining the LURR method with our choice of identified area with increased Coulomb stress, is devised to improve the sensitivity of LURR to measure criticality of stress accumulation before a large earthquake. Retrospective tests of this algorithm on four large earthquakes occurred in California over the last two decades show remarkable enhancement of the LURR precursory anomalies. For some strong events of lesser magnitudes occurred in the same neighborhoods and during the same time periods, significant anomalies are found if circular areas are used, and are not found if increased Coulomb stress areas are used for LURR data selection. The unique feature of this algorithm may provide stronger constraints on forecasts of the size and location of future large events.