金属力学性质的物理力学计算Ⅰ——用Morse势计算弹性常数

讨论了根据内聚能和体积弹性模量实验值,用优选法中黄金分割法,计算立方金属Morse势的方法。利用Morse势,计算出多种立方金属元素单晶和多晶的弹性常数。给出用基本FORTRAN算法语言编写的计算程序。     

金属力学性质的物理力学计算Ⅲ——用Morse势计算拉伸与压缩理论强度

讨论了利用Morse势,按照Born稳定性准则,计算立方金属理论强度的方法。用此法计算了多种立方金属元素(Na,Al,Ca,V,Cr,Fe,Ni,Cu,Sr,Nb,Mo,Rh,Ag,Ba,W,Ir,Au,Th)在单轴应力作用下的拉伸与压缩理论强度,给出用基本FORTRAN算法语言编写的计算程序。     

金属铝弹性常数从头计算的一种新方法

本文报导用原子簇模型计算金属铝的弹性常数。具有面心立方结构的Al_(13)原子簇在不同晶格常数下总能由多重散射兄方法(MS—X_a)算出。采用Morse势函数对原子簇总能进行最小二乘法拟合,然后用拟合出的势参数算出金属铝的弹性常数。我们的计算结果与其他理论计算值及实验值符合甚好。     

惰性气体扩散系数的物理力学计算

根据本文方法求出的惰性气体原子间的相互作用势,计算了He-Ne、He-Ar、He-Kr、He-Xe等气体的扩散系数。其结果与实验值符合得很好,相对误差都在2%以下。     

考虑压缩性影响的超空泡流动数值计算

为了研究水下航行体航速接近水中声速时水的可压缩性对超空泡外形的影响,采用理想可压缩流体定常流动以及Riabouchinsky空化封闭模型,基于势流理论和有限体积方法,对绕圆盘空化器的超空泡流动进行了数值计算.计算时考虑到水的压缩性,水的状态方程采用Tait幂次率.提出了一种新的超空泡外形迭代方法并分析了流体压缩性对超空...     

舰船偏离航道中心运动时水压场的计算方法

本文给出了舰船以亚临界航速偏离航道中心运动时水压场的计算方法。在船体细长、水深较浅的情况下,通过渐近匹配展开法可以得到舰船运动在矩形航道中的扰动速度势控制方程。用纵向分布源涡代替舰船的存在,为满足边界条件需要使用映像原理。从文献知道,源强正比于舰船横截面面积的变化率,涡强需从一个积分方程求取。通过数值计算得到源涡强度分布后,即可求出航道中舰船的扰动压力分布。     

以eAT 为基础的数字仿真方法

本文所叙述的以e~(AT)为基础的数字仿真方法,是将被仿真系统中的各个连续线性定常部分看作子系统,利用计算机先求出每个子系统的离散状态方程,然后根据各子系统的相互联接、其间的非线性或变系数特性,再进行数字仿真。这种仿真方法计算速度快、精度高、有较强的适应性。文章中着重介绍了用解析法求解子系统离散状态方程的系数阵。解析法的优点是:对线性定常系统而言,仿真前的计算量及仿真精度与仿真步距T的选取无关,其缺点是程序较繁。     

磁性合金A1-xBx相图的理论计算

采用变分累积展开方法,研究了块体磁性合金的相图.利用随机键海森堡模型计算了面心立方、体心立方的居里点和耐尔点,给出了磁性合金临界点的计算公式,进而分别给出自旋1/2和3/2情况下Tc-x相图,并加以分析.     

低能电子散射中的定域交换势(Ⅰ)电子被He，Ne和Ar原子的散射

文中利用自由电子气体近似(HFEGE)定域交换势和修正的半经典近似(MSCE)定域交换势,加上正交化方法(Orthogonalization technique),对能量在0.1～20eV的低能电子被He,Ne和Ar 原子的弹性散射过程进行了计算。其中靶原子的电荷分布以及入射电子和靶原子的静电作用利用原子的Roothaan-Hartree-Fock 波函数求出;相关和极化作用则用一个无参数的相关极化势近似表示。通过本文计算结果和严格计及交换作用的结果及实验数据相比较,对HFEGE和MSCE 交换势的准确性以及正交化过程对计算结果的影响进行了分析。     

低能电子散射中的定域交换势(I)电子被He,Ne和Ar原子的散射

文中利用自由电子气体近似(HFEGE)定域交换势和修正的半经典近似(MSCE)定域交换势,加上正交化方法(Orthogonalization technique),对能量在0.1～20eV的低能电子被He,Ne和Ar原子的弹性散射过程进行了计算。其中靶原子的电荷分布以及入射电子和靶原子的静电作用利用原子的Roothaan-Hartree Fock渡函数求出;相关和极化作用则用一个无参数的相关极化势近似表示。通过本文计算结果和严格计及交换作用的结果及实验数据相比较,对HFEGEE和MSCE交换势的准确性以及正交化过程对计算结果的影响进行了分析。     

利用Roothaan-Hartree-Fock 原子波函数计算的低能电子被He和Ar原子弹性散射截面

利用原子的Roothaan-Hartree-Fock波函数计算了低能e~--He和e~--Ar弹性散射截面。计算采用模型势方法,其中包括静电势、交换势和极化势。计算结果和实验吻合较好。通过把我们的计算结果和其他作者利用不同原子波函数的计算结果进行比较,讨论了不同原子波函数之间的差别对计算结果的影响。     

轴对称物体空泡流动的数值计算

基于势流理论和边界元方法,对水下运动的回转体的空泡流进行了数值模拟.将半经验空泡封闭模型由平面流动推广至轴对称流动.对头部为圆锥的圆柱体空泡流动进行了数值计算.分析了空化数、半锥角、空泡长度间的关系.计算与实验结果符合较好.     

爆炸合成金刚石中温度的计算

本文讨论材料在强冲击波作用下,状态参数温度的计算。在数十万大气压的冲击波作用下,固体材料的强度及电子对比热的贡献都可忽略。这样,可采用流体力学模型,用压力代替应力,材料的本构方程即为状态方程,本文取格留乃森状态方程形式。下面结合爆炸合成金刚石课题,讨论石墨、金刚石在强击波作用下,对其温度的数值计算。当然,对于其它材料,在动高压状态下的温度也可用这个方法估算。本文最后给出了适用于爆炸合成金刚石中对温度作工程估算用的解析拟合式。     

利用酉几何构作non-Cartesion认证码

利用酉几何构作了一个non-Cartesion认证码,计算了它的参数;当编码规则按等概率分布选取时,求出了模仿攻击和替换攻击成功的概率。     

平面飞片动态高压装置的一维差分计算

本文讨论了炸药接触爆轰推动金属飞片,与多层材料组成的靶板高速碰撞,获得高压的一维差分数值计算方法,并对中国科学院物理所设计的爆炸法合成金刚石装置作了细致计算。此外,本文还给出了用D=C_0+S_0u及Gruneison状态方程导出的冲击波压力及冲击波恒压的时间宽度等计算公式,它们被用来直接计算平面飞片增压装置中的有关参量。     

微计算机解析计算电偶腐蚀速度

本文介绍了微计算机在金属腐蚀研究中的一种应用。用微计算机进行数据处理,运用电量法进行电偶腐蚀速度的测定。我们采用辛普生数值积分方法,按BASIC语言编制程序,用电子计算机进行处理,求出了电偶腐蚀的阳极溶解量,并进行了907/921、902/921、4/921,902/907四个偶对的电偶腐蚀速度的测定。实验表明,用微计算机进行电偶腐蚀速度的计算,速度快、准确性高,这种方法是可行的,给研究工作人员带来很大的方便。     

KCl 和KBr 晶体热膨胀系数的计算

本文利用Gordon-Kim 模型给出的离子间相互作用势,计算了KCl 和KBr 晶体在室温下的热膨胀系数。所得的结果与实验值进行了对照。     

基于BWRS方程的压缩空气压缩因子计算

使用压缩空气作为介质进行流体计算时,为保证计算精度须考虑压缩因子的影响。为了精确得到压缩空气的压缩因子,对压缩因子的计算方法进行了研究。将空气简化为氮气与氧气的混合气体,通过3个真实气体状态方程计算压缩因子,并与已有数据进行比较,验证了将空气模型简化,再计算其压缩因子的方法是可行的。BWRS方程计算出的压缩因子精度较高,但在不同压力下的计算误差有较大差距。为进一步提高计算精度,考虑了偏心因子的影响,通过计算分析、经验修正,得到了精度较高的压缩因子计算算法,并进行了试验验证。结果表明：不同压力下改变偏心因子的数值对计算结果有不同程度的影响;改进后的BWRS方程计算精度得到了提高,从而验证了BWRS方程改进算法计算压缩因子的有效性和可行性。     

航速对水面舰船极限风压值的影响

本文利用苏联计算水面舰艇抗风性能之极限风速法[1]中有关计算极限风压P_i之方法,对我国3000k吨级以下各型水面战斗舰艇在有、无航速时之P_i进行了计算分析,并与我国水面舰艇稳性规范的计算结果进行了比较。结果表明按极限风速法求出之P_i值在高速时与规范法大致相当,在零速时按极限风速法求出之P_i较规范法大20％-40％。进而论证了我国规范法把倾侧力臂取为受风面积中心至水线的距离(Z_f-T)这处理方法尽管在物理意义上是不明确的,但在应用中是可行的。     

冲压式翼伞开伞仿真计算

本文分析了冲压式翼伞开伞的特点,提出了一个展弦向二维开伞模型,建立了伞衣的径向运动方程。在推导系统运动方程和分离体运动方程的基础上,建立翼伞开伞仿真状态方程,用龙格-库塔法求解,给出开伞充满时间、开伞动载、翼伞速度和位置坐标等参数的变化值。本文提供的开伞仿真计算软件包,可用于各类冲压式翼伞的开伞仿真计算。