有限元法与伪谱法混合求解弹性波动方程

在地震波场数值模拟中,有限差分法、有限元法和伪谱法都是常用的基本方法,但它们各有不同的适应性和优缺点,如有限差分法、有限元法都存在减弱网格频散和提高计算效率的矛盾,而伪谱法的网格频散小且计算效率高.有限差分法和伪谱法在处理地表结构复杂或地表剧烈起伏以及地下结构复杂的情况时存在较大的难度,而有限元法可较为理想地拟合起伏地表和任意弯曲界面,且可方便地处理自由边界条件和界面边界条件.尝试将有限元法和伪谱法相结合,形成地震波场数值模拟的一种混合方法,利用二者的优点,克服二者的缺点,达到既减弱网格频散又提高计算精度和效率的目的.并采用所谓的‘过度区域‘技术解决两种不同算法的衔接问题.模拟实例表明,给出的混合模拟方法不失为弹性波场数值模拟的一种有效方法.     

山区公路软基病害研究

随着高等级公路建设规模的扩大，山区公路所面临的路基病害特别是软土路基病害问题也逐渐增多。结合对陕南勉(县)一宁(强)高速公路工程实例分析，对山区软土成因及其特性进行了初步分析，认为山区软土属以坡洪积、湖积和冲积为主的软土，也有少量是由坡残积物堆积而形成。山区软土的特殊性表现在成分的复杂性、分布的不均匀性、隐蔽性和物理力学性质的特殊性。根据其特性以及中国山区公路软基普遍存在的病害问题，总结了山区公路软基的主要病害类型为剪切拉裂破坏、浸水沉陷破坏、剥蚀坍塌破坏、推挤滑动破坏。还提出了相应的软基处理方法及建议。     

Boundary curves for families of planar orbits

For monoparametric familiesf(x,y)=c of planar orbits, created by a planar potentialV(x,y), we introduce the notion of the family boundary curves (FBC). All members of the familyf(x,y)=c are traced in an allowable region of thexy plane, defined by the corresponding FBC, with total energyE=E(c) varying along the family. Family boundary curves are also found for two-parametric familiesf(x,y,b)=c. The relation of equilibrium points and asymptotic orbits, possibly possessed by the potentialV(x,y), to be FBC is studied.     

模糊划分矩阵在岩土参数概率分布中的应用

讨论如何在小样本条件下用已有的过程经验与试验资料确定岩土参数概率分布,用模糊划分矩阵与BAYES方法相结合,给出由小样本试验数据确定岩土参数的概率分布。     

量子化学方法研究氯代苯酚类化合物对水生物毒性的构效关系

对20个氯代苯酚类化合物进行了CNDO/2量子化学计算,讨论了化合物电子结构与其对水生物翻车鱼(Bluegill)、花鳉鱼(Guppy)、和虾(Shrimp)半致死量之间的关系,分别获得了表示其构效关系的三个线性方程,它们的显著性均远高于α=0.01水平。结果表明:氯代苯酚化合物苯环上碳原子的兀电荷之和(∑Qπ_R)越大苯环上碳原子的净电荷之和(∑Q_R)越大,或LUMO轨道能(E_(LUMO))越低,均使化合物对Bluegill和Guppy水生物的毒性增大;而化合物苯环上羟基邻位的碳原子的兀电荷布居(P_3)越大,同时HOMO轨道能E_(HOMO)越高,则化合物对Shrimp毒性也越大。据此,可预测氯代苯酚类化合物对上述三种水生物的毒性。     

山西省石炭系与二叠系界线

１８８２年李希霍芬在太原西山建立石炭系的太原系含煤建造以来，石炭系与二叠系的分界意见，归纳起来有三：①以Ｐｓｅｕｄｏｓｃｈｗａｇｅｒｉｎａ（ｓ．１）带之底为界，下面为石炭系，上面为二叠系；②以晋东南的灰白色“佳祥砂岩”与太原西山黄绿色“骆驼脖子砂岩”对比，砂岩底面为石炭系与二叠系的分界线，下面是石炭系山西组，上面是二叠系石盒子组，山西组的地层时代定为晚石炭世；③以晋东南灰白色第三砂岩（陵川砂岩）与太原西山“北岔沟砂岩”对比，砂岩底面为石炭系与二叠系的分界线，下面是石炭系太原组，上面是二叠系山西组，山西组的地层时代为早二叠世。笔者认为石炭系与二叠系的划分，以第③方案较适宜。     

重力插值方法研究

为了根据离散观测数据构制连续空间重力变化图像，分析和讨论了３种数值插值方法，计算结果表明多面函数方法插值精度最高。由于逐步回归分析筛选核函数中心点的计算繁琐，文中提出根据分形理论和Ｓｈａｎｎｏｎ取样定理来确定核函数中心点。对滇西试验场进行模拟试算，插值精度可达到４～５（１０－８ｍｓ－２）。     

A MESO-β SCALE SIMULATION OF THE EFFECTS OF BOREAL FOREST ECOSYSTEM ON THE LOWER ATMOSPHERE

Based on the Intensive Field Campaign(IFC-1)data of Boreal Ecosystem-Atmosphere Study(BOREAS).a three-dimensional meso-β scale model is used to simulate the effect of boreal forests onthe lower atmosphere.A fine horizontal resolution of 2 km×2 km is used in order to distinguish thevegetative heterogeneity in the boreal region.A total of 20×25 grid points cover the entire sub-modeling area in BOREAS' South Study Area(SSA).The ecosystem types and their coverage ineach grid square are extracted from the North American Land Cover Characteristics Data Base(NALCCD)generated by the U.S.Geographical Survey(USGS)and the University of Nebraska-Lincoln(UNL).The topography of the study area is taken from the Digital Elevation Map(DEM)of USGS.The model outputs include the components of the energy balance budget within the canopyand at the ground.the turbulence parameters in the atmospheric boundary layer and the wind.temperature and humidity profiles extending up to a height of 1500 m.In addition to the fine timeand spatial step,the unique feature of the present model is the incorporation of both dynamic andbiological effects of the Boreal forest into the model parameterization scheme.The model resultscompare favorably with BOREAS' IFC-1 data in 1994 when the forest was in the luxuriant growingperiod.     

Numerical calculations in the orbital determination of an artificial satellite for a long arc

Numerical methods are usually used for the computation of ephemerides with perturbations for the precise orbital determination of an artificial satellite. But their numerical stability will be encountered in a long arc. In this case the use the improved Encke special perturbation methods has been suggested. The results of this paper show that Encke's method does indeed have a certain effectiveness, but cannot yet completely resolve the numerical stability, and the more efficient method is to use the energy integral or its variational relation to control the growth of the along-track error in general numerical calculations so that the aim of stabilization can be achieved.