褶皱岩层初始挠曲的分叉特性

基于拖带坐标法建立了岩层大变形的一维非线性微分方程，然后用解析与数值方法研究了变刚度、可伸缩岩层在侧向压力及自重作用下的分叉特性，初步揭示了褶皱岩层初始挠曲的形成机理。     

重力插值方法研究

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

关于非球形粒子光散射的T-矩阵数值计算方法的研究

首先介绍了无规取向轴对称非球形粒子光散射的T－矩阵数值方法。然后，进一步研究了T－矩阵计算的收敛问题，提出了一种新的关于无规取向无吸收非线形散射物体的物理收敛方法，我们研究了椭球粒子的收敛问题，结果表明我们的数值方法和T－矩阵计算程序是有效的， 收敛精度与粒子的尺度和形状有很强的依赖性，在一定的条件下，我们的物理收敛速度优于NASA Mishchenko的数学收敛速序。     

三类密近双星自转同步性参量的理论计算和统计分析

利用双星自转同步性理论给出了６９个三类密近双星系统中９３个子星的临界同步自转参量和临界自转周期．并把利用临界自转同步参量所计算的临界自转周期与由气体星自转不稳定理论所计算的临界自转周期做了比较，其结果是两者均属同一量级．     

在重力分层的太阳大气中磁流浮现的数值解

为了解释太阳观测中磁流浮现形成活动区的主要特征——AFS弧顶的上升运动和弧脚的下沉运动,谱斑,Ellerman炸弹等,我们采用了MHD-Lagrangian数值模拟法计算了重力分层(压力密度指数递减),太阳大气中同基底磁场反向和同向的磁偶极子从光球下浮出时产生的动力学演化过程.计算证明:反向的偶极子会导致流体产生会聚和下降运动,为了能让反向偶极子浮出必须适当增压,演化的结果形成电流片,因而有可能引起如Ellerman炸弹等小型爆发现象.还证明了:在浮现处一定产生加热增温,即导致谱斑发亮现象.在弧顶气体是上升运动,速度小;弧脚处气体作下降运动,速度较大,这主要是重力分层气体压力密度非均匀引起的,和Brunt-Vaisala振动有同一的起因.我们还模拟了通常活动区形成时的初期现象:1500高斯强磁偶极子浮进100高斯的弱光球色球基底场时产生的动力学现象.计算表明:为了不致产生激波必须减压浮出,即类似于黑子降温式的浮现,浮现时间应是小时的量级,浮现后确实在色球里导致强的电流片,当这电流片升至低日冕后,会成为Heyvaerts,Priest等人期望的耀斑发生源,因而浮现过程有可能是耀斑储能和触发机理之一.另外我们利用计算结果试解释了北京天文台怀柔站的磁流浮现资料,即文[1]和文[2]的一些分析结果。     

介质横向不均匀对震后地表位移场的影响

利用三维有限元模拟两层地球模型(包括弹性层和粘弹性层)介质的横向不均匀和断层错动的复杂性对震后地表位移场的影响，其数值结果包括弹性解和粘弹性解。结果表明，震后位移场的分布特征不仅与粘弹性层的流变特性有关，还与弹性层介质的剪切模量以及断层的同震错动有关。     

完善征地争议裁决制度的立法思考

如何使土地争议救济渠道顺畅，切实保障被征地农民合法权益，可以说是当前我国土地争议解决中面临的一个大难题。现行法律法规中虽然规定了行政复议、征地补偿安置争议裁决、诉讼、信访等争议解决途径和相应处理部门，但从当前征地争议案件的处理情况看。每一条争议解决渠道都无法有效解决问题。     

液核地球自转速率变化的带谐潮效应Ⅱ——数值计算与比较

地球自转速率的潮汐变化可由无量纲参数k／Cm)(k和Cm分别为壳幔的有效Love数和有效板转动惯量)来表征。对于一个具有弹性地幔、平衡海潮和核幔不耦合的地球k／Cm=0．944，且与潮波频率无关。海潮的非平衡扰动使k／Cm为复数，且与频率有关。大气对自转速率有效勒夫数的贡献约为△kat=0．0075。同时地幔滞弹性对勒夫数也产生扰动。利用本文得到的理论公式和最新的潮汐数据计算了地球自转速率的潮汐变化，及其有关地球物理机制的影响。     

首都圈地区现今断层活动的数学模拟及地震危险性讨论

以首都圈地区现今活动断层上近２０年的位移测量资料为依据，用二维线弹性有限元对该区断层的活动特征进行了拟合，并结合有关资料讨论了该区的地震危险性。结果表明：１９７７－１９８６年间该区主压应力优势方位与华北较一致，约为Ｎ４５°－８０°Ｅ；１９８６－１９９０年其主压应力方向向北偏转，约为Ｎ５°－６０°Ｅ；１９９０年至现在其主压应力又向东偏转，角度大于第一阶段，约为Ｎ８０°－９５°Ｅ。该区的张家口－延庆一带近期有发生中强震的可能；丰镇－阳高－大同地区和凉城－古营盘地区应力较高，也应引起注意。