一种采用解析积分的直线单元的边界元法

本文提出一种采用解析的直线单元计算受有离心力的二维问题的边界元法,推导了相应直线单元的积分解析式。算例表明,本文方法能获得较高的数值精度,尤其在计算靠近边界的内点应力时,效果更好。     

复合材料弹翼－弹身组合结构强度计算

本文利用自行开发的复合材料结构应力分析软件（CMSSA）对某型导弹的复合材料弹冀－弹身组合结构进行了静力分析。本文的计算软件采用动态链表数据结构存贮总体刚度矩阵，具有很高的存贮和计算效率。文中给出了任意四边形壳单元和三角形壳单元的形函数及几何矩阵，这两种单元能适应结构中的复杂曲线边界及单元过渡区域。     

应用Plemelj积分生成二维通道内的正交网格

首先选用的解析函数由速度模以及流向角组成，在已知通道边界上的流向角时，求得通道边界上的速度。由势函数和流函数的定义，求出壁面上的势函数及流函数。其次将解析函数用势函数表示，运用Ｐｌｅｍｅｌｊ边界积分关系的内点公式，求得通道内任意点的势函数和流函数。由此，从等势线和等流函数值线，可形成流体力学数值计算用的正交网格。文中给出两个算例，该方法无需迭代。     

基于比例边界有限元的结构强度分析方法研究

航天结构经历复杂严酷的载荷环境，且诸多大传载关键结构存在小尺寸倒角。为有效处理航天结构小尺寸倒角强度评估问题，提高结构应力应变求解精度和计算效率。本文基于比例边界有限元方法对典型小尺寸倒角的航天结构边界进行离散，构建强度计算模型，针对典型结构不同倒角尺寸开展强度分析，得到不同尺寸倒角对局部应力的影响规律，给出局部区域的应力分布，并与光测试验和常规有限元分析结果进行对比。结果表明，与常规有限元方法相比，比例边界有限元方法可以有效减少计算自由度，提高计算效率和计算精度，由于其半解析的特点，可以方便的得到内部域的应力场，为小尺寸倒角航天结构强度评估、寿命预测及优化设计提供有效的分析手段。     

有进动的高速旋转体的应力分析

本文讨论了采用有限元法对具有进动的高速旋转的轴对称体进行应力分析的方法，首先，推导了轴对称体中哥氏惯性力的分布公式，在柱坐标系下，哥氏力中含有沿径向的零阶对称不载荷和沿轴向的一阶对称及一载荷，首次，采用了具有计算精度，对边界适应性强的八节点四边等参数环形单元，并将受哥氏惯性力的轴对称体看作非轴对称问题，对有限元分析中位移模式的建立和哥氏惯性力的处理作了讨论。最后，简要地说明了所研制的计算机程序，并     

含孔复合材料层合板轴向拉伸三维有限元渐进失效分析

基于渐进失效分析方法,建立了带孔的三维有限元模型.该模型能有效计算应力集中区和自由边界区的应力以及层间应力,为渐进失效分析提供了可靠的基础.本文采用了修正的三维Hashin准则作为单元失效的判断依据,使用交互式退化模型对失效单元进行材料性能退化;同时研究了材料常数υ23值和刚度缩减系数(SCR)对失效载荷预测结果的影响;实现了失效扩展过程的仿真分析;计算了不同层合板的最终失效载荷,计算值与试验值误差小于10%.     

平面图形三角形网格的一种自动生成方案

本文对由若干个标准块串联组成的复杂平面图形提供了一个网格全自动生成方案。标准块分四边形和三角形两种。标准块的串联边是直线段,非串联边可以是直线段也可以是曲线段,因此它适用于常见的复杂平面图形,特别适用于再入飞行器端头温度场及热应力计算、发动机喉衬的温度场及热应力计算、大型压力容器颈口区应力分析等,在标准块内采用二次等参变换的方法来确定各节点的坐标值。本文提供的方案的主要优点是:第一、简单方便、只需把平面图形划分为几个标准块并输入与各标准块有关的少数信息就够了。不需人工预先划分单元和对节点作编号,因此它是全自动的。第二、当需要加密单元时,所有标准块的有关输入信息全部不变,只需改变单元控制边长一个数据就够了。因此对研究解的收敛情况特别方便。以上两点是许多通用程序(如SAP-6、7)采用的方案所不能做到的。第三、由于对标准块采用二次等参变换,因此边界模拟精度较高,能准确模拟任意直线段和二次抛物线段,对圆弧和椭圆弧段也能做到相当准确的模拟。第四、单元形状能较好地满足有限元法的要求,节点编号比较合理,即形成的总刚度矩阵带宽最小,要求计算机内存少和计算工作量小。本文提供的全自动网格生成方案与作者于今年6月完成的外载矢量信息的自动生成(在标准块的非串联边界上采用一维二次插值)和节点温度的自动插取(采用二维线性插值)结合在一起使得平面问题和轴对称物体有限元的算前准备工作变得十分简单。     

切削过程优化问题的边界极值解法

本文首先讨论切削过程优化问题目标函数的极值问题,并完成切削过程优化问题目标函数无极值的证明,给出了在刀具耐用度满足泰勒公式条件下的特例证明。在此基础上提出了该类问题的边界极值解法,把两维空间内的求极值问题转化为一维问题求解,大大减小了计算量。本文还讨论了边界极值求解方法在计算机实现过程中若干问题,包括约束分类方法、约束边界的组合方法、约束边界在计算机内的表达方式等问题。本文最后给出了一个应用边界相比求解的实际例子,并与采用传统的数学规划求解方法进行了比较。结果表明,两者计算结果一致,采用边界极值法的计算速度比数学规划方法提高了一个数量级。     

FL-9低速增压风洞主体结构有限元分析与气压试验

FL-9低速增压风洞是"十五"国家批准建设的大型航空基础设施,笔者简要介绍了该风洞基于有限元方法的应力、变形和模态分析,并与风洞气压试验结果进行了对比.结果表明:有限元计算中所采用的单元类型、模型简化、边界约束等处理方法合理可行,应力计算值与实测值比较一致,最大薄膜应力为147.6MPa,低于许用值42.1%,风洞结构具有较大的安全裕度.风洞的前六阶模态分析为风洞安全运行提供了参考依据.     

含孔边裂纹各向异性有限板的二维问题

应用边界元法，并结合一种特殊的复变函数基本解，研究了含孔边裂纹的各向异性有限板的应力强度因子修正系数的计算。这种特解边界元法避免了在裂纹附近划分单元，可以获得较好的数值计算精度。通过算例分析，结果表明本方法精度高，运算量小，具有一定的工程应用价值。     

复合材料叠层结构孔口附近应力分析方法研究

本文在总结前人对复合材料结构开口的应力集中研究的基础上,根据孔口附近应力场的特点,提出了一种适于工程应用的复合材料叠层结构开口附近应力分析的有限元改进算法。由于应力计算精度的改进仅在孔边、拐角等高应力梯度区域局部进行,且改进后的单元能进一步精确描述应力集中区的应力场。因而,使用该方法计算孔边应力,用较少的单元就能得到满意的精度,大大减少了计算费用和时间。 文中选择了一些有代表性的考题验证了该方法的有效性。并着重研究了带有圆角矩形孔的复合材料层板结构的孔边应力分布,给出应力集中系数随圆角半径的变化规律、孔边加强对应力集中的改善以及±45°纤维铺层比例对孔边应力分布的影响,得到了一些对工程部件开口设计有用的数据资料。     

三向应力状态分析的图解法

本文根据斜截面上应力分量的计算式,提出了“辅助应力圆”的概念,以“辅助应力圆”和莫尔圆的组合探讨了三向应力状态下任意斜截面上应力分量的图解法,以及获得主应力和主方位角的方法。     

用应力函数法求功能梯度梁的弹性理论解

假设弹性模量为厚度的指数函数,泊松比为常数,用应力函数法求解FGM梁的弹性理论解。得到了FGM简支梁弯曲问题的解析解。考察了不同弹性模量变化规律对FGM深梁位移和应力分量的影响。为校验FGM梁、板近似理论和数值方法的有效性提供了依据。     

飞机变速恒频电源系统可靠性的Monte-Carlo分析

随着电力电子技术的发展,用静止变换器取代机械式的恒速传动装置成为可能,因而出现了变速恒频电源系统。在变速恒频电源系统中,发电机是直接由飞机发动机来传动的。由于发动机的转速不是恒定的,因而发电机输出的交流电的频率也是变化的,故再用变换器把变频交流电变换为恒频交流电。飞机变速恒频电源系统的主要部件有无刷交流变频发电机、变换器和控制器。 一般情况下,变换器和控制器主要由电子元器件组成,其失效参数服从指数分布。但是发电机的失效参数为正态分布,因此,其可靠性的分析研究,用一般的解析方法是困难的。本文分析了飞机变速恒频电源系统的可靠性模型。应用Monte-Carlo方法对某型号飞机变速恒频电源系统方案的可靠性指标(R(t),MTBF)进行了数字仿真,得到了各主要部件的工作模式重要度,为该系统的研制设计提供了理论基础。     

振动时效消除铝合金试件残余应力的研究

对振动时效消除铝合金构件残余应力进行了探索，选择矩形截面的铝合金试件作为研究对象，通过加载的方式对其预制残余应力，然后在一定的支承条件下对其激振，实验结果表明，激振之后试件的残余应力数值得到了降低利均化。本文的工作对开展有色金属构件的振动时效研究提供了试验依据。     

Method for Measuring Residual Stresses Induced by Boring in Internal Surface of Tube and Its Validation with XRD Method

Residual stresses can have a strong effect on the usability of machined parts,and the X-ray diffraction(XRD)measuring equipment,which is commonly used to measure residual stresses,is very expensive.This paper presents a method of measuring the residual stresses induced by boring in the internal surface of a tube with much cheaper equipment.The method,called the strain-based method is mainly based on the strains measured on the external surface of the tube.It is proposed on the basis of the very long tube assumption.The finite element method(FEM)analysis is thus used to validate the length of the tube.Guided by the FEM results,an appropriate length of the tube is chosen,and the residual stresses are obtained from both the strain-based method and the XRD method.Stress profiles obtained from both two methods are compared.The comparison result indicates that the profiles of the two methods agree well with each other.Therefore,it can be concluded that the accuracy of the strain-based method is high enough,and it can be applied to residual stress measurement in practice.     

基于SFEM的结构元件疲劳可靠性分析

将结构元件的弹性模量、泊松比、几何尺寸定义为随机场,并将随机场离散成随机变量,将抗拉强度定义为随机变量,采用Monte Carlo随机有限元法(MCSFEM),给出了元件危险点应力的随机分布。由危险点应力分布、材料的p-S-N曲线和Miner理论给出了概率名义应力法,计算常幅载荷下元件的寿命分布,并利用剩余寿命模型计算元件的可靠度。给出了两个算例,结果表明:计算结果与试验结果吻合良好。     

结构元件疲劳可靠性估算的剩余寿命模型

从安全余量的定义出发，提出了一个估算结构元件疲劳可靠性的剩余寿命模型，基于常幅载荷下，疲劳寿命可用对数正态函数分布描述，本文分析了常幅和变幅载荷作用下材料的疲劳剩余寿命的分布形式，认为可用带位置参数的对数正态函数近似描述。在此基础上，系统阐述了元件在常幅和多级载荷下的疲劳可靠性分析方法。3个算例表明：预测结果与试验结果符合很好，本文给出的剩余寿命模型是合理的。     

高加速寿命试验(HALT)与高加速应力筛选(HASS)

本文简要介绍了HALT／HASS试验技术的产生、发展及在国外的应用情况，在深入分析HALT／HASS试验技术基本原理和环境应力与诱发的故障模式之间关系的基础上，阐述了HALT／HASS试验技术的核心－试验剖面图的建立方法，并对HALT／HASS试验技术的特点、所需资源和其它的相关内容作了说明。     

Influence of Patch Side of Heat-Ray Absorbing Film on One-Dimensional Unsteady Thermal Stresses in Window Glass

Heat-ray absorbing film is used to be bonded on the existing sheet glasses of the windows.It is effective for air-conditioning energy saving against the global warming,because it absorbs heat-ray in the thin film and decreases the incoming heat-ray into the room.On the other hand,the sheet glasses increase the temperature at the surface which the sheet is bonded and sometimes yield heat cracks by thermal stresses.It is important to know the state of thermal stresses accurately in order to develop the heat-ray absorbing film with higher performance and without heat cracks.In this paper,the analysis model is treated as the two-layer plate of the conventional soda sheet glass and the heat-ray absorbing film with different absorptivities.The unsteady temperature and thermal stresses are analyzed and calculated numerically.The influence of the patch side,which the heat-ray absorbing film is bonded at the exterior side or the interior side,on the heat-ray absorbing performance and the thermal stresses is discussed.It is found that the alternative patch side has no effect on the heat-ray absorbing performance and that the patch side is recommended to be interior side from a view point of decreasing thermal stresses against the heat crack of glasses.