非齐次Poisson过程的若干新性质

本文在一些文献的基础上,进一步讨论非齐次Poisson过程的某些性质,给出了若干新性质。 设{N(t),t≥0}是累积强度为的非齐次Poisson过程。迄今的文献(如[1][2]等)指出,(?)n≥1,前n个到达时刻τ_1,τ_2,…,τ_n的联合概率密度为 本文定理1指出,上式不仅是必要的,而且是充分的,并给出了充分性的证明。从而,得到了描写过程统计规律的一个刻画。然后,在N(t_i)=n_i(i=1,…,k,0相似文献   

延性金属的断裂强度

由Griffith脆性断裂基础理论引伸,导出了延性断裂理论,求得含有穿透裂纹或表面裂纹非加劲平板结构断裂强度新的表达式。与常用的线弹性断裂力学使用一个材料参数不同,在表达式中使用两个材料参数。本理论独特之处在于两个参数可以由单向拉伸的应力一应变曲线求出;并且,对常用的结构金属,在很宽的裂纹尺寸范围内,应力超过或者低于金属屈服应力下,理论结果和试验数据相当符合。 A—半椭园表面裂纹临界面积,(πac)/2,in~2。(吋~2) Au—在σ=σ_U下半椭园表面裂纹临界面积,in~2。(吋~2) A—埃,0.394×10~(-8)in。(吋) a—半椭园表面裂纹的深度,in。(吋) a_U—在σ=σ_U下半椭园表面裂纹的深度,in。(吋) 2C—穿透裂纹或表面裂纹的长度,in。(吋) 2C_U—在σ=σ_U下穿透裂纹或表面裂纹的长度,in。(吋) 2C_L—在σ=σ_L下穿透裂纹或表面裂纹的长度,in。(吋) E—拉伸时的杨氏模量,Psi(磅/吋~2) h—滑移带的有效高度,in。(吋) h_F—裂纹前缘变形区城的有效高度,in,(吋) h_U—裂纹前缘附近变形区域的有效高度,in。(吋) K_O—线弹性平面应力或混合型的断裂韧性,Psi in~(1/2)。(磅/吋~(3/2)) K_(1C)—线弹性平面应变断裂韧性,Psi in~(1/2)。(磅/吋~(3/2)) K_(TC)—具有中心穿透裂纹的薄板或平板的断裂靱性,Psi(in)~(1/(2 ω)(磅/吋~((3 2ω)/(2 ω)) K_(pC)—具有中心表面裂纹的薄板或平板的断裂靱性,Psi(in.)~(1/(2 ω)(磅/吋~((3 2ω)/(2 ω))) K—厚度参数 L_G—单向拉伸试验中所用的应变片长度,in。(吋) n—ε_(TP)之Ramberg—Osgood关系的指数 P—单位厚度塑性能吸收率,L bs/in。(磅/吋) T—产生单位面积新裂纹表面所消耗的能量,Lbs/in。(磅/吋) t—断裂试件厚度,in。(吋) t—单向拉伸试件厚度,in。(吋) t_o—平面应力断裂的最大厚度,in。(吋) U_E—可用于产生新裂纹表面的单位厚度弹性能,Lbs(磅) U_S—产生新裂纹表面时单位厚度所消耗的能量,Lbs(磅) U_P—塑性变形时单位厚度所消耗的能量,Lbs(磅) U_F—裂纹前缘塑性变形时单位厚度所消耗的能量,Lbs(磅) U_(F1)—在σ=σ_U下,裂纹前缘塑性变形时单位厚度所消耗的能量,Lbs(磅) U_(F2)—在σ=σ_L下,裂纹前缘塑性变形时单位厚度所消耗的能量,Lbs(磅) U_U—裂纹前缘附近塑性变形时单位厚度所消耗的能量,Lbs(磅) U_(U1)—在σ=σ_U下,裂纹前缘附近塑性变形时单位厚度所消耗的能量,Lbs(磅) U_(U2)—在σ=σ_L下,裂纹前缘附近塑性变形时单位厚度所消耗的能量,Lbs(磅) W—试件宽度,in。(吋) W_F—在应力—应变曲线下面,从颈缩开始时的应变到σ_F的应变之间的塑性能密度, Psi(磅/吋~2) W_U—在应力—应变曲线下面,从σ_L的应变到颈缩开始时的应变之同的塑性能密度, Psi(磅/吋~2) β—厚度参数ε_L—在σ=σ_L下的单向拉伸应变ε_N—修正后的颈缩单向拉伸应变ε_U—颈缩开始(σ=0.995σ_U)时的单向拉伸应变ε_F—在σ=σ_F下的修正后的单向拉伸应变ε_F—在σ=σ_F下的平均单向拉伸应变(应变片长度内平均) ε_Y—在σ=σ_Y下的单向拉伸应变ε_(PL)—在σ=σ_L下的单向塑性应变ε_(PU)—在颈缩开始时的应力下的单向塑性应变ε_(PF)—断裂应力下的单向塑性应变ε_(TL)—在σ=σ_L下的单向真正拉伸应变ε_(TY)—在σ=σ_Y下的单向真正拉伸应变ε__(TU)—颈缩开始时的单向真正拉伸应变ε_(TF)—在σ=σ_F下的单向真正拉伸应变ε_(TP)—单向真正塑性拉伸应变ε_(TPU)—在σ=σ_L下的单向真正塑性拉伸应变ε_(TPY)—在σ=σ_Y下的单向真正塑性拉伸应变ε_(TPU)—颈缩开始时的单向真正塑性拉伸应变ε_(TPF)—在σ=σ_F下的单向真正塑性拉伸应变λ—裂纹形状因子μ—厚度参数ν—波松比σ—垂直于裂纹平面的总(毛)面积应力(单向拉伸应力),Psi(磅/吋~2) σ_L—相当于0.0005单向塑性应变的弹性极限拉仲应力,Psi(磅/吋~2) σ_Y—单向屈服拉伸应力,Psi(磅/吋~2) σ_U—单向极限拉伸应力,Psi(磅/吋~2) σ_(UF)—从σ_U至σ_F的平均单向拉伸应力,Psi(磅/吋~2) σ_F—单向断裂拉伸应力,Psi(磅/吋~2) σ_T—单向真正拉伸应力,Psi(磅/吋~2) σ_(TY)—单向真正屈服拉伸应力,Psi(磅/吋~2) σ_(TU)—单向真正极限拉伸应力,Psi(磅/吋~2) σ_(TUF)—从σ_(T_U)至σ(TF)的平均真正单向拉伸应力,Psi(磅/吋~2) σ_(TL)—单向真正极限拉伸应力,Psi(磅/吋~2) σ_(TF)—单向真正断裂拉伸应力,Psi(磅/吋~2) φ—裂纹形状参数ω—断裂靱性参数     

用消去法求逆矩阵的标准程序设计

在本文中,我们将给出一个求逆矩阵的方法和它的程序设计。如果 Γ_(nin)…Γ_(lil)AΓ_(lil)P_1…Γ_(njn)P_nQ=E 则我们有 A~(-1)=Γ_(ljl)…I_(njn)Γ_(njn)…Γ_(ljl)EΓ_(ljl)P_1…Γ_(nfn)P_nQΓ_(nin)…Γ_(lil) 此处 detΓ_(if)≠0,detP_m≠0,(m=1,2,…,n). 这样,我们得到了一个求逆阵的消去法。     

单向分类随机模型方差分量的经验Bayes估计

对于在质量控制、生物育种以及计量经济等方面有着广泛应用的单向分类随机模型 其中y_(ij)为观测数据,ε_i～N(0,σ_1~2)为随机效应,ε_i～N(0,σ_2~2)是随机误差,通常要求考虑对其参数μ,σ_1~2,σ_2~2作出估计,特别是对方差分量σ_1~2,σ_2~2的估计。而在目前常见的估计方法中,大都有一个共同的缺点,即不能保证方差分量估计的非负性。本文采用经验Bayes方法,在平方损失下,利用密度函数及其偏导数的核估计构造出方差分量的一个经验Bayes估计,这个估计是非负的。并证明了这种估计是渐近最优的,在一定条件下,其收敛速度可任意地接近于0(n~(-1/2)。     

关于Davidson—Lanczos方法的收敛率

本文对文[1]提出的求解大型对称矩阵A的极端(几个最大或最小)特征值及相应特征向量的Davidson—Lanczos方法,用Rayleigh—Ritz逼近理论,研究了该方法的收敛率。证明了由该方法产生的规范正交向量{v_i}_i~m=1是Krylov子空间K_m≡Span(v_1,Av_1,…,A~(m-1)v_1)的一组基。设A的k个最大特征值为又,λ_1>λ_2>…>λ_k,相应的近似特征值为λ_i~(m)(i=1,…,k),得到 这里γ_i(γ_i>1),W_i和W_i~(m)是常数。     

关于二阶马尔可夫过程的样本序列是否存在无后效性的讨论

本文指出了工程界关于高阶马尔可夫过程的一个错误定义,证明了(p=2)满足这个定义的平稳高斯过程是不存在的。 本文还指出由二阶微分方程 x″(t) a_1x′(t) a_2x(t)=ε(t) (其中ε(t)是白高斯过程)描写的随机过程x(t)的任意均匀采样序列都不能是AR(2)序列,而由下面微分方程 x″(t) a_1x′(t) a_2x(t)=ε′(t) βε(t)描写的随机过程x(t),当β~2>[max(c_1~2,c_2~2)]时(c_1、c_2是特征方程z~2 a_1z a_2=0的根),至少存在一个采样间隔Δ_1,使相应的样本序列是AR(2)模型,因此是一个二阶广义马尔可夫序列。     

对某些奇异积分的逐次分半加速法

对于f(x)∈C~(2p 1)[0,1],尤拉-麦克洛林求和公式可写为:(h_k=1/2~k) (?)(1)对此,在逐次减半加速法中计算序列: T_(mi)=(γ_mT_(m-1,i 1)-T_(m-1),:)/(γ_m-1)(m=1,2,3,……) (2)其中γ_m=4~m,形成γ序列(4,4~2,4~3,……)。 对于f(x)∈C~(2p 1)(0,1],我们建立公式: (?)(3)其中(?) (4)对于f(x)=x~βln~nx(β>-1,n=0,1,2,…)则有(?) (5)(3),(5)右边的前半部是由于f(x)在x=0上的奇异性而引起的误差。如将α=2~(β 1)以(n 1)次加入γ序列,就可怕除此项误差。选出十个计算实例作为说明。     

面向纤维增强复合材料低速冲击损伤的非线性混合模型

基于连续介质损伤力学提出了一种纤维增强复合材料(Fiber reinforced polyrner/plastic,FRP)结构低速冲击损伤预测的渐进损伤模型,包含非线性剪应力应变关系和归一化的混合模式基体损伤演化,用来预测复合材料层合板低速冲击损伤。模型区分了纤维拉伸/压缩、纤维间拉伸/压缩4种层内损伤以及层间分层损伤;纤维间损伤起始由Puck失效准则预测,损伤演化由断裂面上的等效应变控制,失效判定时考虑了就位效应对强度的影响;模型中加入单元特征长度以消除计算结果对网格密度的依赖性。以[45_4/-45_8/45_4],[0_3/45/-45]_S和[45/-45/0_2/90/45/0_2/-45/45]_3三种铺层的复合材料层合板为例,预测了不同冲击能量下复合材料层合板的低速冲击损伤响应参数,试验结果证明了本文模型的有效性。     

导电薄圆筒上的电荷分布——矩量法的改进处理

设薄圆筒以X轴为轴,二端面各为x=-1,x=1。当静电平衡时面电荷密度σ(x)满足积分方程: integral from -1 to 1 G(x∣x′)σ(x′)dx′=常数 (1)设:U(x):=integral from 0 to x σ(x)dx (2)并令:g(U,U′)=G(x∣x′) (3)(1)可表为:ingegral g(U,U′)dU′=常数 (4)对于二维(即圆筒半径为无穷大)情形,(4)的解为 U=2/πsin~(-1)x (5)现以此作为一般情形的尝试解: (ⅰ)把这U区间(-1,1)作2n等分,在与U=-1,-1 2/n,……,1-2/n,1相应的n 1圆环上分布线电荷,其密度各为q_1,q_1 q_2,……,q_(n-1) q_n,q_n,使它们在与U=-1 1/n,……,1-1/n相应的n个圆环上产生相同的电位,对应于(4)可列出n阶线性方程组。 (ⅱ)解出q_1,q_2,…q_n。对于二维情形,可证: q_1=q_2=……q_(n-1) (6)对于一般情形并不如此,但可由此构成新的x-U曲线。 反复(ⅰ),(ⅱ),直到(6)近似满足而使x-U曲线稳定为止。 本法对粗圆筒特别适用,沿圆筒长度不取等分点,而是电荷越密集,取点越密,因而节省计算量,但仍提高了精密度.     

关于A.M.Ostrowski的切线—平行线程序的一类收条件

本文给出非线性泛函方程P(x)=0用A. M. Ostrowski的切线——平行线程序: (?)_n=X_n-Γ_nP(Xn),X_n+1=X_n+Γ_nP(Xn)(n=0,1,2…)此处Γ_n=[P'(Xn)]~(+1),可以求解的一类收饮条件。它推广了M. A. [3]和徐利治[5]的结果。     

广义2一D系统Roess er模型的可接受输入序列研究

通过将有限矩形上的广义２－Ｄ Ｓｏｅｓｓｅｒ模型转化为等价的代数方程，给出了有限矩形上所有输入均为该系统的可接受输入的充要条件；利用２－Ｄ Ｚ－变换以及多项式矩阵的性质，给出了无限矩形上所有输入均可接受的充要条件。证明了当系统的输入矩阵行满秩时，只要系统在某个有限矩形上的所有输入均为可接受输入，则在无限矩形上的所有输入也均为可接受输入。基于矩阵的初等行变换，提出了确定广义２－Ｄ系统Ｒｏｅｓｓｅｒ模     

求多边形最小包容矩形的遗传算法

建立了求任意多边形包容矩形的数学模型 ,将求最小包容矩形问题转化为函数优化问题 ,并用遗传算法求得函数的最优解     

超声速进气道与弹体一体化外形研究

通过数值求解N S方程和风洞测力试验,研究了超声速进气道的布局、剖面形状及几何参数对全弹气动力和进气道内流的影响。结果表明:细长旋成体弹身背风面和侧边不宜安放进气道;矩形进气道能够产生涡升力,升阻比明显高于半圆形进气道,而且内部流动接近二维流动,流场畸变及流动损失情况均比半圆形进气道要好;矩形进气道横截面积增加,对升力系数影响不大,但阻力系数增加明显;其高宽比增加,升力系数增大,阻力系数减小。     

自动几何校正中摄像机的自动校正方法

为了解决多投影仪拼接显示系统自动几何校正中摄像机的自动校正问题,提出了一种基于图像识别的自动校正方法.首先对采集的图像进行分析,识别出目标区域,控制云台对摄像机进行定向,然后计算目标区域角点集的凸包,利用凸包的最小外接矩形判断出目标区域的姿态,并控制云台对摄像机姿态做出实时调整,实现了自动几何校正工作中摄像机的自动校正功能.     

椭圆喷管设计与数值模拟

提出了一种椭圆喷管的设计思路,并对椭圆喷管流场与矩形喷管流场进行了对比分析。通过数值计算表明,在相同的驻室参数、相同的长度和出口面积条件下,椭圆喷管出口马赫数略高于矩形喷管。湍流和高温真实气体效应均降低了喷管的有效面积比,改变了喷管流场的膨胀波系,通过选择合适的喷管出口位置可以获得较好的试验均匀区。椭圆喷管作为高超声速风洞特种试验装置,可以有效利用加热器的能量,提高设备的参数模拟能力,可适用于大尺寸扁平状前缘、舵、翼等模型的防热试验和大宽高比的冲压发动机试验研究。     

旋转直管流场的数值模拟

根据旋转坐标系下的 Favre平均 NS方程以及旋转条件下的流动特征 ,修改了 B/ L湍流模型 ,并提出了一种旋转条件下出口边界压力分布的给定方法 ,对已有的 NAPA通用湍流计算软件进行了改进和发展 ,并计算了旋转直管中的流动。通过与实验数据的比较 ,文中的计算结果与实验结果吻合较好 ,正确反映了流动特征和物理现象 ,验证了修改后的 NAPA软件的可靠性。     

微小矩形通道内流动与换热特性

采用试验和数值模拟方法研究了当量直径为0.633～1.079 mm的微小通道内的流动与换热特性.工质为 液体水,雷诺数Re为100~3 000.试验结果表明,该尺寸范围的通道内层流向紊流转变的临界雷诺数约为1700,层流表观阻力系数和努塞尔数的试验值与数值计算结果吻合良好.最后通过数值模拟对恒定泵功及恒定流 量条件下不同结构参数通道的换热进行了计算,给出了最佳换热结构.     

BTi6431S钛合金盒形件热拉深成形工艺数值模拟

为研究BTi6431S钛合金热拉深成形工艺,本文以某型号飞行器天线盒盖为例,通过设计正交实验方案,先后分别对无压边情况下拉深成形进行了有限元数值模拟,用最薄壁厚减薄率作为描述零件成形结果的评价标准,得出该钛合金件在成形温度为650°C,首次拉深凹模圆角半径11mm、凸模圆角半径7mm、转角半径12mm、拉深高度25mm参数下,最薄壁厚减薄率最小(20.3%)。     

地形跟随系统的仿真

本文探讨了低空突防地形跟随系统的仿真问题。分为如下几个方面:①低空突防飞行任务的仿真目的与基本要求;②地形跟随系统的主要组成部分及其仿真实现方案;③地形输入的实现方法和地形跟随雷达的仿真;④一个实用的多功能地形跟随系统仿真软件包;⑤一种以系统计算机实现的地形跟随系统实时仿真方案,及其实时控制技术原理;⑥地形跟随数字仿真和半物理实时仿真的仿真结果。     

圆柱-球体三自由度超声电机及其控制

研制了一种一阶纵振和二阶弯振模态组合的圆柱-球体三自由度超声电机。研究了其中两个自由度的运动轨迹跟踪控制策略,并将其应用于电机样机的控制。其中,运转速度的控制采用激励电压幅值调节法,针对位置控制提出了两种运动轨迹控制策略一种是逐点比较控制,另一种是矢量分解控制。本文提出了调节激励电压幅值的3种脉冲宽度调制(PWM)方法,对这些方法作了比较,并通过一系列实验得到验证。结果表明当改变激励电压幅值时,保持3个振动模态的相位差不变和激励电压信号不失真,是简化控制过程、提高控制性能的关键;矢量控制法效果较好,并有轨迹跟踪误差小、运动平稳和噪声小等优点。