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

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

重力插值方法研究

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

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

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

基于Windows消息机制的人机交互的研究与应用

结合实例着重介绍了如何通过串行通讯利用Windows的消息机制来解决目前工业设备控制中人机交互的模拟问题,并对这两种技术进行了详细的论述.     

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

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

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

为了解释太阳观测中磁流浮现形成活动区的主要特征——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。同时地幔滞弹性对勒夫数也产生扰动。利用本文得到的理论公式和最新的潮汐数据计算了地球自转速率的潮汐变化，及其有关地球物理机制的影响。     

原子钟噪声的蒙特卡络模拟方法

讨论了原子钟噪声的蒙特卡络模拟方法．具体过程是：先产生模拟所需的正态分布随机数，从实际钟的阿仑方差估算出噪声电平，然后用１组递点函数产生高质量的模拟钟时差的随机序列．在对１台商品铯钟进行具体模拟后，给出了真实数据与模拟值的阿仑方差之问的相符程度，最后对实际问题进行了分析．