ROMBERG方法的改进

本文研究了全部根构成一个等比数列的多项式的性质，并将所得的结果用到数值积分的Ｒｏｍｂｅｒｇ方法＾「１」，得到了相应遥改进。     

两类求解y“=f（x,y）具有极小相位延迟的显式方法

本文提出了两类数值积分二阶周期性初值问题ｙ〃＝ｆ（ｘ，ｙ），ｙ（ｘ０）＝ｙ０，ｙ（ｘ０）＝ｙ０具有检小相位延迟的显式两步法。这些方法推广和改进了文献「１」１－「７」中的某些方法。数值试验表明本文中的某些方法优于「１」－「７」中的某些方法。     

空间情形下旋转体体积的近似计算

基于定积分的数值计算方法,解决了空间一般式曲线绕任意轴所形成的旋转体体积的近似计算问题,给出了相应的Matlab程序,并用两个实例进行说明.     

基于国产仪器研究钨丝原子吸收光谱

用标准的150 W幻灯灯泡的钨丝作原子化加热材料,设计制作了一种新颖的、结构简单的钨丝电热原子化器,同时还设计了该原子化器的电源电路和信号的数据采集电路.该原子化器能的最高原子化温度可达3100 K左右;编写了硬件的控制软件和信号处理的应用软件.用该原子化器将实验室一台火焰原子吸收光谱仪改装成钨丝原子吸收光谱仪,并以铜元素标准溶液和PerkinElmer公司多元素标准溶液对仪器性能作了研究.仪器对铜的检出限为0.0133 μg/mL,线性范围为0.10～4.0 μg/mL;对1.0 μg/mL的铜标准溶液测试的相对标准偏差为RSD=4.1%(测试次数n=10),每次分析所需样品量20 μL.     

多体系统动力学微分／代数方程组数值方法

多体系统动力学微分／代数混合方程组又称Ｅｕｌｅｒ－ｌａｇｒａｎｇｅ方程,是近十年来动力学和计算数学领域研究的热点之一．本文介绍这两个领域中引入的传统的数值积分方法与新的理论．      

多体系统动力学微分／代数方程组数值方法

多体系统动力学微分／代数混合方程组又称Ｅｕｌｅｒ－ｌａｇｒａｎｇｅ方程，是近十年来动力学和计算数学领域研究的热点之一．本文介绍这两个领域中引入的传统的数值积分方法与新的理论．     

无条件稳定数值积分方法在拟动力实验中的应用研究

用无条件稳定的隐式α－方法代替目前的条件稳定的显式方法，通过数值迭代控制方法实现了拟动力实验．在此基础上进一步应用子结构技术，解决了刚度大、自由度较多结构采用传统的条件稳定数值积分方法进行拟动力实验的问题．应用不同数值积分方法进行的三层模型结构的实验结果表明，无条件稳定的数值积分方法的实验精度优于条件稳定的数值积分方法的精度．文中还就实验误差的控制补偿问题进行了讨论，提出了采用瞬时割线刚度进行补偿的方法     

非线性动力学积分方程分块积分解法

对于非线性动力学方程组分块地应用精细积分算法,使其化成积分方程表达式,求解的表达式中具有相对低阶的转换矩阵,从而使精细积分更适用于多自由度、强非线性、变系数、非保守系统,针对积分方程提出了一个显示预测-校正的单步四阶精度自起步的精细积分算法。算例表明本方法是有效的。     

延性单晶非均匀变形的三维有限元模拟

对延性单晶在拉伸载荷作用下的应变局域化和颈缩等非均匀变形过程进行了三维有限元数值模拟。将相关晶体塑性本构模型及一种新的数值积分方法补充到ABAQUS6.1商用有限元软件中。该方法的特点是，利用晶体塑性的动力学方程，获得一个关于晶体弹性变形梯度的演化方程，采用半隐式积分方案进行求解。本文推导出一种新的应力变本构矩阵。按此方式更新本构矩阵，计算速度和计算稳定性大大提高。加载方式，边界条件和变形程度等因素影响着滑移系的启动状况，这是平面模型所不能预测的。本文利用三维有限元方法模拟了不同取向下滑移系的启动状况，全面地考虑了FCC单晶材料12个可能滑移系在变形过程中的启动状况，合理地模拟了FCC面心立方单晶沿不同取向加载时晶轴旋转导致的应变局域化和颈缩等非均匀变形过程。