M.E.Shvez迭代解法及其在力学中的应用

当微分方程中含有微量项时,可用M.E.Shvez迭代法求解。但当微量项出现奇性,或在某一区间内微量项并非微量时。用此法求解将遇到困难。本文针对这类问题,把原M.E.Shvez迭代解法稍加改变。算例表明,用改进了的M.E.Shvez法求解上述问题。其精度比原M.E.Shvez法的有所提高。     

KdV-Burgers-RLW方程的高精度差分格式

初值问题的差分解法,参数ε≥0,μ≥0. 这一方程当ε=μ=0时为KdV方程,δ=ε=0时为Burgers方程,而当δ=μ=0时为RLW方程.对于方程(1),已设计了许多计算格式.对于KdV方程,最早的格式当推Zabusky-Kruskal,后来有[2—6].对于RLW方程,也有许多工作.对于Burgers方程,格式就更多了.非线性波动方     

行列式(续)

四.行列式的基本性質2) 主对角元素都是1,其余的元素都是0的行列式,它的值等于1. 很容易可以用对阶数n作归納法来証明这个命題.当n=1时命題显然成立.对于任意n,所說行列式有形狀     

一类干扰序列性质的研究

设α是正实数且α()Z~+,Hegyvari研究干扰序列Q_α={[αn~2]|n∈Z~+}并证明了:1)α具有形如2k+1/2或2k+2/3的结构,其中k是非负整数,当且仅当:Q_α中元素均为偶数;2)α具有形如3k+3/4的结构,其中k是非负整数,当且仅当:Q_α中元素均为3的倍数.类似得到干扰序列αn(n+1)/2‖n∈Z~+的相关性质.     

三维快速多极边界元法分析地埋管群传热问题北大核心CSCD

基于三节点三角形线性单元,为克服单元跨叶子积分难题,将三维位势问题快速多极边界元法与几乎奇异积分的半解析算法相结合,实现了三维边界元法中几乎奇异积分的准确计算,该方法适用于U型地埋管薄体结构的换热分析.在制冷、制热两种工况下研究了U型地埋管壁厚对换热量的影响,并进一步分析了管群间的热相互作用.计算结果显示,当管壁导热系数一定时,管壁越厚,对管内流体和土壤之间的换热影响越大.当钻孔间距一定时,管群中埋管数量越多,热干扰现象越强烈,提高管群换热量的主要措施是降低管群间热干扰.因准确计算了几乎奇异积分,三维快速多极边界元法可以有效计算薄体和厚体耦合的三维热传导问题.该文方法和分析结果可为地埋管换热器系统的工程应用提供参考.     

中子活化示踪沙技术的研究

本文选取铱(Ir)作为信号元素,采用化学吸附标记的方法制作示踪沙,应用中子活化分析方法检测泥沙中信号元素的含量,依据信号元素的运动特征判断泥沙的运移规律.该方法克服了流明沙技术及放射性示踪沙技术的缺点,并能真实地模拟天然泥沙的运动特性,是一种先进的泥沙运动观测技术.在天津港进行的国内首次实际应用试验中取得了预期效果.     

基于倾向值匹配的干扰数据检测方法

数据中掺杂干扰数据的现象十分常见,对于随机出现的干扰数据处理,目前已有很多方法可以借鉴,但对于人为的干扰数据,若继续使用传统方法,则可能不会达到很好的效果.倾向值可以用一维数值来描述多维数据的特征,且当数据具有相近的倾向值时,其本身常常也很相似,并可能来自同一总体.因此,文章提出一种应用倾向值匹配检测干扰数据的新方法,...     

交错符号矩阵数的界

满足性质A1.所有元素都等于 1、-1或0;A2.每行、每列元素之和都等于 1;A3.每行、每列非零元素的符号交替变化的 n 阶方阵称为 n 阶交错符号矩阵,其全体记为 A_n,A_n=|A_n|为 n 阶交错符号矩阵的数目。通过计算已经验证:当1≤n≤10时,     

具有常量红利界的带干扰项的风险模型及最优红利界

对于保险公司来说,如何确定其红利策略,使得投保人利益最大化是一个需要研究的课题.研究了具有常量红利界的带干扰项的经典风险模型下,索赔量为混合指数分布情形时的最优红利界的计算方法.     

改进的分块模方法求解对角占优线性互补问题

为了高效求解中小型线性互补问题,本文提出了改进的分块模方法,并证明了关于严格对角占优(对角元素均为正数)线性互补问题的收敛性.对于广义对角占优线性互补问题,先将其转化为严格对角占优线性互补问题,再采用改进的分块模方法求解.数值结果表明,改进的分块模方法在求解广义对角占优线性互补问题时在内迭代次数和计算时间上均明显优于分...     

侧向加热腔体中的多圈型对流斑图

基于流体力学方程组的数值模拟,研究了倾角θ=90°时侧向加热的大高宽比腔体中的对流斑图.对于Prandtl数Pr=6.99的流体,在相对Rayleigh数2≤Ra r≤25的范围内,腔体中发生的是单圈型对流斑图.对于Pr=0.0272的流体,取Ra r=13.9,随着计算时间的发展,腔体中由最初的单圈型对流斑图过渡到多圈型对流斑图,这是出现在侧向加热大高宽比腔体中的新型对流斑图.对不同Ra r情况的计算结果表明,Ra r对对流斑图的形成存在明显的影响.当Ra r≤4.4时是单圈型对流滚动;当Ra r=8.9~11.1时是过渡状态;当Ra r≥13.9时是多圈型对流滚动.对流最大振幅和Nusselt数Nu随着相对Rayleigh数的增加而增加.该对流斑图与Pr=6.99时对流斑图的比较说明,对流斑图的形成依赖于Prandtl数.     

非奇异阵特征极值和条件数的近似计算

本文从共轭梯度法的公式推导出对称正定阵A与三对角阵B的相似关系,B的元素由共轭梯度法的迭代参数确定.因此,对称正定阵的条件数计算可以化成三对角阵条件数的计算,并且可以在共轭梯度法的计算中顺带完成.它只需增加O(s)次的计算量,s为迭代次数.这与共轭梯度法的计算量相比是可以忽略的.当A为非对称正定阵时,只要A非奇异,即可用共轭梯度法计算ATA的特征极值和条件数,从而得出A的条件数.对不同算例的计算表明,这是一种快速有效的简便方法.     

固体内任意元素的边界积分变分定理——任意裂纹开展时能量释放率的计算

根据[1～2]中提出的离散型固体力学及其变分原理的基础上,本文形成四种类型的元素的边界积分变分定理.当进行断裂分析时,可用它们计算沿裂纹边界法线方向的能量释解率;在有孔洞时,当在孔洞边界存在或不存在外力作用的情况下,可用它们计算沿孔洞边界法线方向的能量改变量;当进行离散分析时,也可以用来建立离散方程.以便求解待解函数值.并且由本文分析可知,在[3]中提出的J积分形式是不确切的.     

有限元广义伽略金方程,边界变分方程,边界积分方程

在[1]的基础上,我们进一步应用可动边界的变分原理于固体体系的离散分析,得到有限元广义伽略金方程,边界变分方程,边界积分方程.这些方程描述了待解函数在元素内部与元素的边界上应满足的方程.当对固体体系进行离散分析时,可以应用这些方程去建立不同情况下的求解待解函数的离散方程.亦可作为相应情况下的简化计算的依据.由本文得到的边界积分方程可知,在[2]中提出的J积分形式,应用于内部元素边界的围道积分计算是不适宜的.     

布尔代数与推理

对于只含有两个元素 0和 1的集合 ,规定三种运算 :+ (加法 ) ,·(乘法 )和′(补运算 ) ,用下表表示+ 0 10 0 11 1 1· 0 10 0 01 0 1ｘｘ′0 11 0对于加法和乘法 ,交换律、结合律、乘法对加法的分配律和加法对乘法的分配律都成立 ;对于补运算 ,德·摩根律成立 .这便是 (二值 )布尔代数 .我们约定某事用字母表示 ,某事为真时取值为 1 ,某事为假时取值为 0 .例如事件Ａ为真 ,事件Ｂ为假 ,则记为Ａ =1 ,Ｂ =0 .根据加法表 ,我们得到只要事件Ａ、Ｂ中有一个为真 ,则Ａ +Ｂ就为真 ,只有当Ａ和Ｂ都假时 ,Ａ +Ｂ才假 ,即只有当Ａ =0且Ｂ =0时 ,才有…     

正规稀疏幻方的存在性

设d＜n为两个正整数.一个密度为d的n阶正规稀疏幻方,记为Sn(d,o),是一个n×n的整数矩阵,其每行每列恰有d个非零元素、n-d个零元素,其非零元素集为1到nd的所有整数构成的集合,其每行每列每两条对角线上元素和都相等.正规稀疏幻方是幻方的推广且在图的标号中有很好的应用.本文证明存在一个Sn(d,o)当且仅当n为奇数时d≥3,n为偶数时d也为偶数且d≥4.     

轴对称弹性体的有限元分析

轴对称弹性力学问题的有限元分析长期以来都是采用三角圆环有限元和线性形状函数.由于积分困难,常用近似积分求得刚度矩阵,这种近似积分对于靠近旋转对称轴的元素,误差很大,所以,长期以来,被认为不满意的办法.也有用精确积分计算刚度矩阵的.但本文指出,这种积分只适用于有中孔的轴对称体.对于实心的轴对称体而言,这种刚度矩阵都不收敛,计算是无效的.本文提出了一种新的形状函数,当径向座标r接近于零时,这种形状函数的径向位移u自然地接近于零.如果用这种新的形状函数,则由此计算求得的刚度矩阵,不论三角圆环有限元的位置是否靠近轴线.都是存在的.这种有限元,就能用于计算实心的轴对称体的问题.     

特征2的有限域中一些特殊元素的存在性(II)

设q=2s.s,n为正整数,Fqn为qn元素的有限域.在本文中,我们考虑Fqn中一些特殊元素的存在性.主要结果是:当下面的条件之一成立时,在Fqn中存在ξ使得ξ和ξ+ξ-1都是本原元并且ξ+ξ-1还是一个正规元:1.当n|(q-1)时,n37,s>6,或者2.当n|■(q-1)时,n≥34,s>6.进一步,如果n是奇数,则当下列条件之一成立时,存在ξ∈Fqn使得ξ和ξ+ξ-1都是Fqn的本原正规元:1.当n|(q-1)时,n≥257,s>9,或者2.当n■(q-1)时,n≥43,s≥9.     

最高阶元素个数不同的有限群

本文首先讨论了当群 G 的某 r 阶元素的集合 M_r(G)只含两个元素时 G 的性质,然后给出了当群 G 的最高阶元素的集合 M(G)只含两个元素时群 G 的一个刻划,最后得到了当|M(G)|为奇数或不大于4时,群 G 为超可解;当|M(G)|=2p,p 为素数,G 为可解群.     

脑内各向异性扩散传质与吸附反应过程数值分析

对脑组织内传质过程的机理及其影响因素进行了分析,建立了综合考虑脑内物质各向异性扩散、吸附和反应过程的数学模型,模型方程采用隐式控制容积法进行数值求解.计算结果表明:组织迂曲度越大,物质的扩散越慢,当某一方向迂曲度较小时,物质浓度明显增大,物质扩散变快,由于脑组织的非均质性,脑内物质的扩散传递存在着竞争现象;吸附与反应作用会抑制脑内物质传递,吸附速率越大,抑制现象越明显,对于脑内非线性的米氏反应过程,当反应速率常数增大时,稳定浓度会显著减小,同时米氏常数的增大则会使得稳定浓度值增大.相较于吸附过程,米氏过程的抑制性作用更为明显.