离心力场中的三维裂纹的J积分有限元计算

本文研究离心力场中的三维裂纹问题。在W.S.Blackburn给出的不计及离心力时的三维J积分公式的基础上,作者导出了离心力场中的三维裂纹的J积分公式为: J=integral from Γ(Wdy-T_iu_(i,1)ds)-integral from A((σ_(i3)u_(i,1),_3dA)-integral from A(F_iu_(i,1)dA) 本文将上式化为有限元表达式并研究了转换时的简化技巧;对中心含半圆表面裂纹的标准块在单向拉伸载荷作用下的J积分进行了计算和激光光弹实验;对某发动机涡轮盘的角隅裂纹进行了计及离心力的J积分计算。实践证明本文列出的方法是简便和切实可行的。     

断裂韧性K_(1C)和断裂门槛值ΔK_(th)可靠性测定方法

采用定量方程随机化方法 ,建立了断裂p JR Δa曲线方程 ,断裂p JR Δa曲线与钝化线的交点即为试验测定具有可靠度p的J积分临界值 (J1C) p值 ,根据断裂力学测试原理 ,由试验测定 (J1C) p 值可转化出断裂韧性 (K1C) p 值 ;还建立了断裂门槛ΔKth值的可靠性测定方法 ,并对试验数据进行了统计分析处理 ,得到了一些具有重要应用价值的结果和数据 .研究发现 :断裂门槛值ΔKth服从正态分布     

基于T-H方程的多脉冲最优交会方法

针对椭圆参考轨道交会问题,采用T-H方程描述相对运动,提出了一种时间固定燃料最省的多脉冲最优交会方法,优化参数为交会脉冲及其施加时刻.当考虑到J2摄动或航天器初始相对距离较大时,用T-H方程进行状态预测其线性化误差一般不容忽略,而若用轨道积分预测则耗时较多,进而导致优化时间过长.针对此问题,提出一种采用前一优化脉冲节点的状态导出的轨道根数预测当前节点状态的预测方法.该方法简单实用,有效地加快了优化收敛速度.最后基于多脉冲优化解进行了数值轨道积分以验证交会精度.仿真结果表明,即使加上J2摄动,在初始相对距离为1 000 km时,该方法的终端位置精度仍能达到75 m.     

用再结晶法研究三维弹塑性裂纹问题

本文用再结晶法测定了中心裂纹拉伸试件裂纹尖端高应变区的塑性应变分布,求得了高应变区参数与J积分及裂纹尖端张开位移COD的定量关系,并用三维有限元法进行了计算。实验和计算结果表明,用再结晶法研究三维弹塑性裂纹问题是一种有效的方法。     

北京航空学院学报1981年总目录(总第16期—总第19期)

(总第16期一总第19期)北京航空学院学报1981年第1期(总第16期)周盛沈孟育j 111 121 131 141 151 161 171 181 191201 21跨音叶栅流爆的相仿律伞翼气动力特性的实验研究二维广义J积分用复位表示的一般二维热弹性问题带冠空心冷却叶片的三维应力分析金属宏观断口形态的数学分析乙C4高强度铝合金应力腐蚀断裂机理的研究(一)关于非线性控制系统的近似方法相似理论及其在发动机调节系统中的实现现代滤波理论襟翼三予定位置瞬时螺旋运动的图解及解析计算一九八O年度我院获得上级嘉奖的科技成果简况张行熊昌炳朱永昌毛士艺林保具胡继忠吴国勋黄执…     

广义正交多项式在时变系统跟踪问题中的应用

采用广义正交多项式的展开技术,利用其积分、乘积运算阵,对于时变线性系统最优控制中的跟踪问题,直接从最优控制的性能泛函入手,将最优控制问题转化为代数极值问题,从而避免了求解非线性Riccati方程,得到了一个较为简便的算法,并结合实例说明了方法的可行性.      

航天飞行器的两大数学问题

本文考虑火箭燃烧室燃烧不稳定性,固体推进剂不稳定燃烧,液体火箭贮箱推进剂晃动引起火箭系统不稳定性以及航天飞行不稳定性等问题的偏微分方程组定解问题。并指明航天飞行器的断裂问题是断裂力学的时代问题之一,研究它的数学方法如奇异积分方程法,J积分方法,研究复合材料的线性断裂力学原理以及处理边缘断裂问题的边界单元法都是值得探讨的。     

金属混合断口的宏观定量分析

本文在分析拉伸试样发生缩颈时的屈服方式的基础上,给出由拉伸数据估算J积分值的近似公式:J=ak(a/b,m)ηI;并给出计算拉伸断口的纤维区尺寸、剪切唇尺寸的粗略的近似公式: 2a/d≈(n 3.87/n~2 7.74n 14.17)ψ 0.05和ψ≈((n′)~2 2n′ 0.3)/(n′ 1)(S/(d/2 0.77)     

一种微波及毫米波滤波器的精确设计方法

介绍了一种微波及毫米波滤波器精确设计方法:由滤波器的原型参数,导出滤波器的广义阻抗变换系数;采用模匹配法或耦合积分方程技术等数值算法精确计算各个耦合结的S参数,进而转换成广义阻抗变换系数;通过插值,获得满足原型要求的滤波器各个耦合结的结构尺寸以及谐振腔长度,将结构参数代入模匹配程序、HFSS或CST软件进行分析,分析结果与预期值一致。设计了26GHz矩形波导带通滤波器,模匹配法和HFSS以及CST软件的仿真结果、测试数据与设计要求一致,说明该种设计方法是有效的。     

1996-2008年期间活动区与耀斑活动的特征分析

在1996—2008年期间,可见日面共出现3067个活动区,其中活动区面积S满足10μh≤S≤90μh(μh为太阳半球面积的百万分之一)的活动区数量占活动区总数的58.52%,100μh≤S≤490μh的活动区数量占活动区总数的31.43%,而S≥500μh的活动区数量仅约占活动区总数的6%,S≥1000μh的活动区数量仅约占活动区总数的1%,面积为0的活动区约占活动区总数的4%.1996—2008年期间,除面积低于50μh的活动区外,不同面积活动区的数量占该年活动区总数比例的最大值通常都在太阳活动峰年之后,即相对太阳黑子数平滑月均值都有一定的延迟.在1997—2008年期间,至少能够产生一个A级耀斑的活动区数量为1408个,而在此期间,共产生了7706个C级耀斑、1141个M级耀斑和110个X级耀斑.C级耀斑在1~8A波段流量的积分值为22.711 J·m~(-2),M级耀斑在1~8 A波段流量的积分值为29.443 J·m~(-2),X级耀斑在1~8 A波段流量的积分值为36.178 J·m~(-2).第23太阳活动周期间,南半球的耀斑活动比北半球强.     

基于状态空间模型的无刷直流电机控制器设计

在无刷直流时机的位置伺服系统中，采用了基于状态空间模型的优化设计方法，综合出了次优控制器、建立了广义误差系统状态方程，并利用积分罚函数分段一经法处理伺服系统不等式约束的问题，得出了次优控制算法，经仿真证实，有杉该优化方法设计出的伺服系统的稳态跟踪误差为零；稳成输出不受阶跃干扰的影响；并具有期望的瞬态响应特性。     

基于矢量观测的陀螺故障诊断算法

卫星在轨运行期间, 采用反作用轮进行姿态控制及推力器进行位置保持, 星体惯量会发生改变并存在控制力矩. 提出了一种基于运动学方程的矢量观测算法应用于陀螺故障诊断, 采用广义罗德里格参数描述卫星姿态, 避免了四元数在计算过程中出现模值不为一的情况; 计算过程中采用四元数进行姿态运动学计算, 避免了广义罗格里德参数复杂的积分运算. 针对陀螺突变和缓变两种故障, 分别在角速度幅值变化很小和幅值缓慢减小两类情况下进行仿真. 仿真结果表明, 该算法可有效估计星体角速度, 在上述两类情况下均能准确地对突变和缓变故障进行诊断, 不受星体惯量变化和控制力矩的影响, 并且计算量较小, 对工程应用具有较好的参考价值.     

伽略金法分析双缝波导定向耦合器

本文利用伽略金法分析了由宽边重合的矩形波导,通过公共壁上两个纵向窄缝隙耦合而形成的定向耦合器特性。基函数系和检验函数系分别取为三角函数。根据等效原理和边界连续条件,分别在两个缝隙上下口径上建立四个耦合的积分方程,待求口径上的切向电场。通过矩量法把积分方程变成矩阵方程,对矩阵中的元素(或内积)进行处理之后,利用计算机进行数值求解,求出口径切向场、散射场、耦合系数,方向性系数、隔离度和电压驻波比等参数。给出了标准X波段矩形波导(B J—100或WG16)定向耦合器参数随频率的变化关系。     

铁木辛柯梁固有振动频率的边界元解法

根据傅里叶变换推导具有两个广义位移的铁木辛柯梁固有振动的基本解.利用加权残量方法,从控制微分方程出发建立边界积分方程,进而根据边界条件得到频率方程,采用代数特征值方法和影响系数方法求解频率,并分析了两种方法的特点.以杆为例证明了对于一维均匀结构,对不同的边界条件利用边界元方法(BEM,Boundary Element Method)都可以得到精确频率.将铁木辛柯梁的BEM结果与有限元结果和精确解进行了比较.     

由高频无线电波反射回波参数反演电离层运动和结构的高度剖面

本文从广义电波射线微分方程组出发,导出了时变、不均匀和各向异性电离层介质中高频无线电波反射回波的可测参数与电离层的运动和结构参数之间的积分方程组,并给出了利用现代数字测高仪中的多普勒图、定向图、频高图和测距图,反演电离层的结构和运动随高度的分布剖面的计算方法.根据这种反演方法,单站探测电离层运动和结构剖面时有多种可测参量的组合方式.文中提出了一种最佳的参量组合方案,为探测设备的改进和观测记录的分析提供了理论依据.     

广义非线性强度与广义滑动面

广义Mises准则满足八面体面上的剪应力和正应力之比为常数,SMP(Spatially Mobilized Plane)准则满足空间滑动面上的剪应力和正应力之比为常数,那么对于包含或逼近现有非线性强度准则的广义非线性强度,在物理空间中应存在一平面,且满足该平面上的剪应力和正应力之比为常数.为此,提出了介于八面体面和空间滑动面之间的广义滑动面,与广义非线性强度相对应,广义滑动面上的剪应力与正应力之比为常数,并给出了广义滑动面上应力比的表达式,及主应力空间中的强度面形状.通过深入分析和对比,表明了广义非线性强度和广义滑动面在材料强度描述上的一致性.      

J2轨道模型对月地激光通信预瞄准精度影响研究

基于J_2轨道预测模型,设计了不同倾角和轨道高度圆轨道月球卫星,通过将J_2轨道预测模型预瞄准仿真结果与直接积分RKF7(8)法仿真结果相对比,研究了不同类型月球卫星轨道对-Y面卫星舱板激光通信终端瞄准精度的影响。仿真结果表明,J_2轨道预测模型可以满足月球极地卫星月地激光通信要求,当轨道高度为1 000 km和2 000 km的时候,10 min方位角偏差不超过40μrad,而俯仰角偏差仅为7μrad时,在一定程度上J2轨道预测模型可以满足月地激光通信预瞄准要求。     

带包边的对称正交铺层复合材料层板分层问题的解析——广义变分解法

对于对称正交铺层复合材料层板单向受拉分层问题,本文根据穆什海里什维利的求解各向同性平面问题与列赫尼兹基的求解各向异性平面问题的复变函数方法,得到了满足所有基本方程、裂纹表面边界条件与层间连续条件的应力场、位移场的本征展开式。对于包边与不包边这两种情况,本文利用分区广义变分原理代替裂纹表面以外的边界条件,确定奇异应力场的控制量。由于所有基本方程预先得以满足,在变分方程中只有线积分而无面积分。计算表明,本方法前期准备工作简单,计算节省机时,结果收敛迅速。     

断裂韧性K1C和断裂门槛值 Δ Kth可靠性测定方法

采用定量方程随机化方法,建立了断裂p-J_R-Δa曲线方程,断裂p-J_R-Δa曲线与钝化线的交点即为试验测定具有可靠度p的J积分临界值(J_1C)_p值,根据断裂力学测试原理,由试验测定(J_1C)_p值可转化出断裂韧性(K_1C)_p值;还建立了断裂门槛ΔK_th值的可靠性测定方法,并对试验数据进行了统计分析处理,得到了一些具有重要应用价值的结果和数据.研究发现:断裂门槛值ΔK_th服从正态分布.