浅论经济增长与物价水平的关系

从交易方程出发，初步分析了经济增长与物价水平的关系，得到了一些有益的结果。     

Room-temperature equation of state of Li<Subscript>2</Subscript>VOSiO<Subscript>4</Subscript> up to 8.5 GPa

High-pressure single-crystal X-ray diffraction measurements of lattice parameters of the compound Li₂VOSiO₄, which crystallises with a natisite-type structure, has been carried out to a pressure of 8.54(5) GPa at room temperature. Unit-cell volume data were fitted with a second-order Birch-Murnaghan EoS (BM-EoS), simultaneously refining V ₀ and K ₀ using the data weighted by the uncertainties in V. The bulk modulus is K ₀ = 99(1) GPa, with K′ fixed to 4. Refinements of third order equations-of-state yielded values of K′ that did not differ significantly from 4. The compressibility of the unit-cell is strongly anisotropic with the c axis (K ₀(c) = 49.7 ± 0.5 GPa) approximately four times more compressible than the a axis (K ₀(a) = 195 ± 3 GPa).     

复杂网络河的数值模拟

依据复杂网络河的特点及其内在联系，本文建立了其特征矩阵及由Newton-Raphson格式表述的隐式圣－维南差分方程组的系数矩阵，经输入边界条件及解方程组，可得到不同各支段上，下游端点上的水位和流量值。     

长江口海水的密度

为了核对1980年国际海水状态方程应用于中国近岸海水的可靠性,高精度的磁力浮沉子密度计在15—25℃,盐度5—35,一个大气压条件下,测量了长江口附近海水及其人工稀释海水的相对密度。结果表明,海水密度直接测量值系统地偏高于状态方程计算值,二者之间的偏差随水样盐度值降低而增大,且与海水[Ca~(2+))/s比值成明显正相关关系。经计算机拟合,导出了长江口附近海水及其稀释海水的状态方程,该方程计算值与实验值的平均标准偏差为±3.9×10~(-3)kg/m~3。     

基于麦克斯韦理论的海底大地电磁探测技术

麦克斯韦理论为海底大地电磁探测技术奠定了坚实基础。大地电磁场的测量通常测取Ex，Ey，Hx，Hy,Hz。大地电磁仪的设计要点包括：高灵敏度、智能化、同步采集数据、锚系硬连接、承压密封和合理的下沉速度所观测的数据经计算机处理后，可获得被测海域的地球物理信息。     

O(T2)精度下椭球界面Dirichlet边值问题的积分解

研究了边界是参考椭球面的Laplace方程Dirichlet边值问题的求解，在O(ε4·T)精 度下给出了参考椭球界面上扰动重力位Dirichlet外问题的积分解式. 该结果理论上优于目 前常用的球近似下的积分解式，从而为研究物理大地测量中边值问题的求解提供了新的依据     

Quantum-thermodynamic treatment of intrinsic anharmonicity; Wallace’s theorem revisited

Wallace (in Thermodynamics of crystals, 1972) developed a theorem, rooted in rigid lattice dynamics, which incorporates intrinsic anharmonic effects in solids. The practical application of this theorem in mineral physics is computationally involved and this is the main reason for the theorem not getting the attention it deserves. Because intrinsic anharmonicity is an important issue at the extreme conditions in planetary mantles, we derived a method which removes the computational obstacles in applying this theorem. We extended the theorem to incorporate details of the phonon spectrum and tested our algorithm on forsterite (Mg₂SiO₄). Using a least squares inversion technique applied to all available experimental data, we show that it results in an accurate representation of thermodynamic properties and sound wave velocities of Mg₂SiO₄ in its complete pressure–temperature stability range. We also show that the accuracy of our results is not significantly affected by the use of a different equation of state.     

Incremental elastoplastic FEM for simulating the deformation process of suction caissons subjected to cyclic loads in soft clays

This paper presents an incremental elastoplastic finite element method (FEM) to simulate the undrained deformation process of suction caisson foundations subjected to cyclic loads in soft clays. The method is developed by encoding the total-stress-based bounding surface model proposed by the authors in the ABAQUS software package. According to the model characteristics, elastoplastic stress states associated with the incremental strains of each iteration are determined using the sub-incremental explicit Euler algorithm, and the state parameters describing the cyclic accumulative rates of strains are updated by setting state variables during the calculations. The radial fallback method is also proposed to modify the stress states outside the bounding surface to the surface during determination of the elastoplastic stress states. The stress reversals of soil elements are judged by the angle between the incremental deviatoric stress and the exterior normal vector at the image stress point on the bounding surface to update the mapping centre and state variables during cyclic loading. To assess the general validity of the method, the reduced scale model tests and centrifuge tests of suction caissons subjected to cyclic loads are simulated using the method. Predictions are in relative good agreement with test results. Compared with the limit equilibrium and quasi-static methods, the method can not only determine the cyclic bearing capacity, but can also analyse the deformation process and the failure mechanisms of suction caisson under cyclic loads in soft clays.     

Sensitivity of penman-monteith reference crop evapotranspiration in Tao’er River Basin of northeastern China

A non-dimensional relative sensitivity coefficient was employed to predict the responses of reference crop evapotranspiration （ET0） to perturbation of four climate variables in Tao＇er River Basin of the northeastern China. Mean monthly ET0 and yearly ET0 from 1961 to 2005 were estimated with the FAO-56 Penman-Monteith Equation. A 45-year historical dataset of average monthly maximum/minimum air temperature, mean air temperature, wind speed, sunshine hours and relative humidity from 15 meteorological stations was used in the analysis. Results show that： 1） Sensitivity coefficients of wind speed, air temperature and sunshine hours were positive except for those of air temperature of Arxan Meteorological Station, while those of relative humidity were all negative. Relative humidity was the most sensitive variable in general for the Tao＇er River Basin, followed by sunshine hours, wind speed and air temperature. 2） Similar to climate variable, monthly sensitivity coefficients exhibit large annual fluctuations. 3） Sensitivity coefficients for four climate variables all showed significant trends in seasonal/yearly series. Also, sensitivity coefficients of air temperature, sunshine hours and wind speed all showed significant trends in spring. 4） Among all sensitivity coefficients, the average yearly sensitivity coefficient of relative humidity was highest throughout the basin and showed largest spatial variability. Longitudinal distribution of sensitivity coefficients for air temperature, relative humidity and sunshine hours was also found, which was similar to the distribution of the three climate variables.     

岩石圈挠曲-弹性薄板小挠度弯曲的新方程

岩石圈挠曲研究采用的是弹性薄板小挠度弯曲方程(即克希霍夫方程),克希霍夫方程基于薄板的前提,忽略并假设薄板内垂向应力为0。本文在无需垂向应力为0的这一与地质事实不相符的假设的情况下,由弹性体几何方程、物理方程和静力平衡方程推导出岩石圈挠曲-弹性薄板小挠度弯曲的新中面方程,具有同等的数学简洁性。取泊松比为0.25时,有DFF/D=1.125,即新挠曲方程中的挠曲刚度DFF要比经典的克希霍夫挠曲方程的D值大12.5%。本文推导的新方程不仅可以在岩石圈动力学,也可以在弹性力学中获得应用。