逐步回归方法及其计算程序

(一)计算方法概述1.给出n一1个自变量X:,X:,…Xn一1和因变量Xn的m组观察值:xit(i二1,2,…n,t二1,2,一m)求线性回归方程Xn二bo bxX一 … bn一IXn一l(1)使与m组观察值在最小二乘意义下最佳拟合。由数理统计知识知道,求bi的正规方程为:n.1 艺101{m艺(Xit一又i)(Xj,一又j m_一、,一、=艺(Xjt一Xj)(Xnt一xn) t,l (j=1,n,!)‘2) .O、.,了」、.了t.l其中X 1一Mm艺Xit.1(Xit一又i)(Xjt一又j)=1,n) S‘石a‘j=百落了则(2)式可化为标准正规方程:n一l 乏aij ai== anj(j=l,n一l)(8)i一1在求得ai以后,则回归系数(盛) _n一1-bo=Xn一艺biX(5)解方…     

用有限元法计算锥壳基础的应力

前言实践表明,在一些大型烟囱的设计中,采用空腔锥壳基础新技术,在节约钢筋,水泥等建筑材料方面,具有良好效果。在包头二冶建研所模型实验的基础上,我们用弹性力学的有限元法对他们的Q—2,Q—3两实验壳进行了计算。计算结果与试验实测数值基本一致。     

用有限元法解板架结构屈曲问题

本文采用多项式位移函数建立压杆的刚度矩阵，采用逆迭代法求解板架结构单向受压时面外屈曲的欧拉应力和失稳波形。探讨了该位移函数的误差，并提出考虑材料物理非线性的两种处理方法。     

弹性体有限元法的平面问题分析

文章基于微分分析法，对限于二维范围内的弹性体所存在的平面应力、应变问题，得用内力平衡的条件，给出有限元法平面应力问题的基本方程的分析过程，得到由位移求应变的几何方程和由应变求应力的弹性方程。     

用有限元法对水坝的位移和应力分析

应用有限元法对一三角形截面水坝在静水压力作用下的位移和应力进行了计算．本方法也适用于其它类似的平面应变问题．     

用动态有限元法解一维气体动力学方程组

动态有限元(MFE:Moving Finite Elements)法是一种解波动问题的有效方法,它的特点是将结点视作动点,在计算中使它们自动地向大梯度的地方集中,从而提高了计算的精确程度,特别是在求解带有强间断面的波动问题中,显出了其特有的优点,激波面的厚度可以处理得十分薄。美国学者K.Miller,J.Gelinas等人近年来发表了不少的文章,介绍     

用有限元法解重调和方程边值问题的收敛性

本文对于三角形剖分分别应用六次、七次多项式作为试验函数。对于矩形剖分应用双三次扦值函数作为试验函数,来寻求重调和方程的数值解,得出离散化的误差估计,并指出其收敛的价数。     

基于有限元法的轮毂应力及变形分析

为研究轮毂应力大小及变形情况,针对安徽车桥自主研发的AQ5020127244410型轮毂进行建模,对静载和紧急制动两种工况分别添加载荷以及约束条件,分析轮毂受力及变形。结果表明,静止状态下,轮毂承受最大应力为101. 73MPa,产生最大变形0. 023mm。紧急制动状态下,轮毂承受最大应力122. 07MPa,产生最大变形0. 028625mm。紧急制动状态下,轮毂承受应力变形更大,轮毂承受应力小于材料许用应力,满足使用要求。     

用解加数法解运输问题的通用程序

物资的合理调拨问题是属于线性规划中的一种特殊类型——运输问题。这类问题运算复杂规模一般都是很大的。因此,在实际工作中,充分使用电子计算机是准确、迅速解决这类问题的好方法。解运输问题最常用的有位势法,单纯形法等。然而通过实际计算表明,采用在我国尚未使用的解加数法比位势法有更多的优点。它既能节省机器的内存贮量,又能大大缩短机器的运算时间。例如用DJS-8计算机计算m×n=19×419(m为发站,n为收站)的化肥运输合理方案时,用解加数法仅需机时为5分钟,而用位势法需机时约25分钟。虽然电子计算机的运算速度很快,但是在方法的选择和程序的编制中,如何达到既能节省机器的内存单元,又能提高机器的使用效率,仍然是一个很重要的问题。目前,在我国铁路运输工作中,已开始使用解加数法编制通用程序,计算物资的合理调拨方案,如计算铁矿石、化肥、汽油等的最优调拨方案,并取得较好的效果。本文将介绍用解加数法解运输问题的方法和程序的编制。     

用k_ε紊流模型及有限元法解二维异重流方程

本文用水沙混合水流的运动方程和质量守恒方程作为二维泥沙异重流的控制方程,用k-ε紊流模型闭合雷诺应力,用有限元法解上述方程组,计算了泥沙异重流的流速、含沙量、紊动动能和紊动粘滞系数的分布,并与试验的流速和含沙量资料进行了对比,最后分析了含沙量与紊动动能、流速的关系,阐述了在含沙水流中紊动动能有可能转化为水流的势能.     

粘弹性体接触时接触应力的有限元法分析

以粘弹性Maxwell模型材料的积分松弛型本构关系为基础，进行了本构关系递推，并建立了Maxwell模型粘弹性体的有限元分析模型，再利用节点对接触模型理论编制了Maxwell模型二维接触问题接触应力计算有限元程度，并验证了其正确性。最后应用该程序，研究了简单Maxwell模型轴对称粘弹性体接触时接触面上轴向应力分布及其变化规律。     

用非线性空间有限元法分析砌石拱坝应力的初步探讨

一、前言 近十多年来,砌石拱坝全国各地修建很多,目前坝型向着较高较薄方向发展。但是,目前坝体应力分析方法仍然沿用砼拱坝的方法,与实际砌石坝体的应力状况出入较大,特别是砌石拱坝梁向刚度较砼拱坝小,为了反映砌石拱坝这一特点,笔者试提出如下计算模式,如图1所示将砌石拱坝沿垂直方向分成二种不同材料的单元体相间构成,即一层六面体单元 (八结点)一层任意四边形薄板单元,也就是将几层砌体的砌缝砂浆填料集中在一层几十公分厚的薄板单元来考虑,这种单元的材料只能承受很小的拉力,而当压应力等超过一定值时将产生塑性变形。以此来考虑砌石…     

大烘缸的有限元法应力分析

用有限元法对Ф2500mm大烘缸(杨克式)进行了全面的应力分析和研究,求出了大烘缸各部位的应力分布规律,其中包括:缸体应力分析;端盖的应力分析;缸体与端盖连接处法兰口附近的局部应力分析和端盖人孔附近应力集中分析。通过分析研究,提出了大烘缸强度计算中应注意的几个问题。     

用有限元法计算大坝底孔的非稳定温度场和温度应力

本文介绍了用有限元法计算某坝底孔温度场和温度应力所得的结果,并根据温度应力的计算结果,分析了其对工程结构的影响。计算结果和原型观测成果进行了比较,两者甚为符合。文后附录介绍了计算方法,并讨论了计算格式的稳定性。     

鲁棒过程控制系统的设计方法及其CAD程序

过程控制中常用的二阶加纯滞后模型往往具有多线性相关的参数摄动形式，本文找到了计算这类系统的闭环鲁棒稳定裕度的简便方法。以此为基础，形成了满足指定鲁棒稳定裕度的鲁棒过程控制系统的设计方法，并给出了ＣＡＤ程序框图。按此框图编程容易，勿须借助任何软件包和工具箱。     

用有限元法解五层对称平面光波导的传播特性

将有限元法应用于五层对称平板波导传播特性的精确计算中 ,采用Matlab编程 ,计算了凸型和W型五层对称平板波导导模的传播常数和模场分布 ,与传统报道的结果比较表明 ,有限元方法计算简便易行 ,计算精度高     

用有限元法分析非轴对称载荷作用下波纹管的应力与位移

提出了非轴对称载荷作用下计算波纹管的位移和应力的有限单元法，将波纹管离散为锥壳单元，单元之间用节点圆相宫结，载荷和位移沿环向展开成傅立叶级数，这样降低了计算维数，减少了计算工作量。分析和算例表明，该方法能方便有效地用于计算波纹管的位移和应力，可为工程中波纹管的设计提供理论依据。     

用有限元法分析非轴对称载荷作用下波纹管的应力与位移

提出了非轴对称载荷作用下计算波纹管的位移和应力的有限单元法。将波纹管离散为锥壳单元，单元之间用节点圆相连结，载荷和位移沿环向展开成傅立叶级数，这样降低了计算维数，减少了计算工作量。分析和算例表明，该方法能方便有效地用于计算波纹管的位移和应力，可为工程中波纹管的设计提供理论依据     

基于依赖图的程序克隆分析及近似解求解方法

大型软件系统中的代码复制和修改现象可能导致程序缺陷的扩大以及无用代码的存在,克隆代码检测分析则有助于抽取可复用的软件组件和模式,在软件重构和软件演化中起重要的作用。在给出了一种基于子图同构进行克隆代码检测的方法基础上,加入依赖边类型约束求近似解的算法,以改善算法的时间复杂性。     

用原始对偶算法解运输问题及其程序

本文讨论用原始对偶算法解运输问题的计算步骤。其比点是计算工作量小,与一般解运输问题的算法相比约少50%。文中列出用“719”算法语言编写的程序及其说明。程序采用存储正变量的个数和它们的位置的方法,因此大大缩短了运算时间。在“719”电子计算机上,求解大约有100个变量的运输问題,包括打印计算结果的时间在内,约需15秒钟;机上允许出现的变量可达6000个以上。如果采用单纯形算法程序,则解同样的问题所需要的时间将超过2分钟;允许出现的变量也只有几百个。对于有大量变量的问題,原始对偶算法将更为有效。