固体火箭发动机药柱热老化结构分析方法

为研究推进剂热老化过程中的发动机结构分析方法,建立了一种含老化效应的粘弹性本构模型,基于增量有限元法对该本构进行了数值离散,导出了本构方程的增量形式,并提出了应力更新方法;基于Abaqus软件二次开发技术,编写了相应材料子程序并开展了算例分析。结果表明,老化粘弹性本构可表征分析过程的材料老化影响;本文有限元解与对应解析解吻合良好,分析方法及算法程序有效,可用于三维条件下的固体发动机药柱老化粘弹性有限元分析。     

三维粘弹性大变形动力响应分析

将大变形增量问题的Ｔｏｔａｌ Ｌａｇｒａｎｇｉａｎ法应用于粘弹性结构的动力响应分析，利用Ｐｒｏｎｇ级数导的Ｋｉｒｃｈｈｏｆｆ应力和Ｇｒｅｅｎ应变表示的粘弹性积分型增量本构方程及其递推公式，从而给出一种计算三维粘弹性大变形动力响应的有限元增量迭加方法，经算例验证，所述方法及其相应程序具有机器内存少、计算量小和效率高的优点。     

不可压缩粘弹性有限元法及其对推进剂药柱应力分析的应用

由于固体推进剂具有近似不可压缩的特点 ,利用位移法有限元求解这类材料的力学问题不可能得到满意的结果。本文根据 Herrmann变分原理导出一种适用于不可压缩和近似不可压缩粘弹性物质的新的本构关系式 ,介绍了相应有限元公式 ,并开发了用于平面问题分析的计算机程序 VFAPINP。计算结果表明 ,本方法和程序可以用于所有泊松比的粘弹性问题的计算 ,尤其适用于固体推进剂药柱的应力分析。     

粘弹性阻尼结构的动特性分析

本文研究了粘弹性阻尼结构动力特性的分析方法。根据分数导数模型表示的粘弹性材料本构方程,对一种阻尼材料用拟合法构造了它的力学性能解析表达式。建立了一种八自由度夹层梁有限元模型,从而给出了粘弹性阻尼结构的有限元动力学方程。为分析粘弹性阻尼结构的动力特性,提出了一种非线性特征值问题的迭代解法。应用上述方法对悬臂约束阻尼层梁的复模态参数进行了分析,并进行了实验验证。结果表明,分析与实验的结果一致,说明本文所提出的方法是可靠的且有较好的精度和实用性。     

固体推进剂的热粘弹性有限元分析

本文首先从弹性——粘弹性的相似性出发,导出张量形式的三维线性热粘弹本构方程,然后根据虚功原理建立相应的有限元公式,最后编制了平面应变条件下的有限元分析程序TVE,并计算了二个实例。计算结果表明所介绍的分析方法与计算程序具有相当高的精度,可用于固体推进剂药柱在固化冷却过程巾的应力分析。     

粘弹性结构振动特性分析的摄动求解技术

 本文根据粘弹性结构的特点,对其振动特性进行简化分析,构造了这类结构振动特性分析的摄动求解方法,给出了低阶摄动求解公式。还讨论了零阶摄动有重特征值时的求解。     

腐蚀电场平面问题边界元法研究

 研究腐蚀电场平面问题的边界元理论和分析方法。首先,详细论述腐蚀电场平面问题的边界元理论。其次,针对电极极化曲线的非线性特性,提出一个求解含非线性边界条件的边界元问题的迭代解法。最后,按本文所述理论和方法,具体计算了某些电偶腐蚀电场的电位与电流密度分布,并与有关实驻结果比较,发现二者相符甚好。     

时变泊松比HTPB固体推进剂体积和剪切松弛模量反演

为表征HTPB复合固体推进剂泊松比的时间相关性，以经典积分型粘弹性本构模型为基础，开展了从单轴受力响应到体积和剪切松弛模量反演问题研究。运用径向基配置法和带近似奇异核的高斯求积公式求解反演问题的第二类Volterra积分方程。将获得的体积和剪切松弛模量导入Abaqus有限元软件并开展材料力学行为仿真。将有限元结果与单轴松弛、蠕变和等速率加载实验数据进行对比，验证反演方法的正确性。分析了粘弹性泊松比Prony级数参数对相对松弛模量的影响规律。结果表明，当泊松比与时间相关时，相对体积和剪切松弛模量曲线出现分离，其分离程度受长期与瞬时泊松比之差的影响，而分离快慢受泊松比迟滞时间影响。采用正确的体积和剪切松弛模量，经典粘弹性本构模型不仅能合理预测HTPB固体推进剂松弛和蠕变等粘弹性行为，还能正确反映泊松比的时间相关性。     

一种新的玫瑰扫描目标跟踪边界重心法

针对自动寻的导弹玫瑰扫描目标跟踪方法中由于扫描过程非线性引起的目标采样非线性,导致求取目标重心算法的运算量和计算误差大等问题,提出了一种基于格林公式的、新的求取目标重心的边界重心法,该方法仅使用目标边界上的采样点。通过对运算量和计算误差的理论分析,与现有的权重重心法相比,本方法可以有效解决采样不均匀的问题,同时运算量减少了90%,误差仅为权重重心法的下界。另外,本方法也可用于基于图像信息的凝视目标跟踪等类似算法中来提高跟踪准确性及运算速度。     

近距格斗空空导弹发射边界判定研究

简要介绍中空空导弹发射边界的拟合方法，在此基础上提出一种基于工程设计经验的“定参数查表拟合法”，该方法选取对导弹发射边界影响较大的参数作为拟合公式的自变量，通过两次回归算法的构造出拟事公式，并对拟合公式进行了验证计算。计算结果表明，该方法满足空战中对空空导弹发射包线计算的实时性要求，并具有较高的拟合精度。     

时域推进与边界元混合瞬态热传导问题

对时间变量在小区间内应用多项式插值,对空间变量采用边界元求解,提出一种新的分析瞬态热传导问题的计算格式。与现行的利用含时间变量的基本解所得到的边界元计算格式相比较,具有基本解简单易得、系数矩阵容易计算以及无重复积分等优点;与向前差分消去时间因子的边界元计算格式相比较,则在算法稳定性、计算精度以及时间步长的选取等方面具有明显的优势,并且很容易推广到求解非线性问题。算例表明,采用时域推进与边界元混合方法分析瞬态热传导问题有效可行。     

表面裂纹的分区混合元分析

 本文对三维体表面裂纹问题进行分析,应用应力和位移的渐近解在表面裂纹区域构造一个奇异应力单元,采用分区混合有限元法计算裂纹顶端各个位置的应力强度因子。     

梁壳组合结构的一种非线性有限元计算方法

提出了梁壳垂直组合结构在有限元计算中存在的问题,并以几何大变形问题为依托提供了解决该问题的计算原理和有限元方法。数值算例证明了该理论和方法的正确性。该法不仅适用于梁壳组合结构的大变形非线性计算,也适用于静动力小变形问题的线弹性分析计算。     

复合材料胶接搭接接头应力分析方法研究

胶接是复合材料结构主要连接方法之一,对胶接接头进行应力分析是保证复合材料安全性、耐久性的关键。在初步设计阶段,一般采用解析方法对胶接接头进行应力分析及参数研究。针对复合材料双搭接和单搭接胶接接头,在Tsai等人的理论分析方法（TOM方法）基础上,提出了一种改进的搭接接头剪应力分析方法,该方法考虑了被胶接件的剪切变形,认为被胶接件只有在靠近胶层的半个厚度上产生剪切变形,剪应力沿该半厚度呈线性分布。算例分析结果表明：本文方法比现有的分析方法更接近于有限元模拟结果,可用于估算复合材料胶接接头剪应力分布。     

基于Kriging模型梯度解析解的改进一次二阶矩方法

改进一次二阶矩（AFOSM）法是一种基于功能函数梯度的结构可靠性分析方法，鉴于其对隐式函数的梯度较难求解，提出了一种基于Kriging模型梯度解析解的AFOSM方法，利用Kriging代理模型的解析表达式推导得到功能函数对输入变量的梯度解析解，为AFOSM中设计点的确定提供高精度的梯度信息。通过Kriging与AFOSM的结合，很好地解决了基于有限元模型的隐式情况下梯度计算量相当大、可靠性分析难的问题。数值与工程算例验证了所提Kriging梯度解析解的较高精确性，同时验证了所提基于Kriging解析解的AFOSM结构可靠性分析方法的正确性与较高精度。     

机翼和翼身融合体跨音速绕流有限元算法的研究

从全速位方程出发研究了机翼和翼身融合体跨音速绕流的有限元算法。给出适用于有限元法和复杂外形的块结构网格生成技术和使单元系数矩阵积分数目减少2/3的方法;并发展了一种内存量少的总体系数矩阵存贮方式。为了提高有限元离散方程组解的精度和加速收敛,提出α—GMRES算法代替传统全速位方程有限元解法中所习用收敛慢的线松弛迭代解法。通过机翼和翼身融合体实例考核,效果是好的,为飞行器先进气动布局的数值计算提供了一种新的有效算法。     

有限元的一种并行算法及实现

当前,由于有限元素法大都使用串行算法,处理时间较长。本文根据有限元素法的特点和采用的并行计算机结构,提出了一种并行有限元算法,能大大提高其处理速度,加快处理周期。且算法简单,实现方便。     

基于导热性能的复合材料微结构拓扑优化设计

基于均匀化理论和拓扑优化技术,提出了复合材料稳态导热性微结构构型的设计方法。根据均匀化理论给出的周期性材料微结构(单胞)等效导热系数的计算公式,结合有限元分析手段对单胞的等效导热系数求解;定义设计区域内微结构每个单元的导热系数为伪密度设计变量,约束材料用量,分别以材料微结构某个方向导热性能最好为目标和宏观结构的散热性能最佳为目标;采用实体各向同性惩罚函数法(SIMP)构造单胞拓扑结构优化数学模型,导出材料等效导热系数的灵敏度计算公式;采用凸规划对偶求解与周长控制约束相结合进行拓扑优化并获得了优化的微结构构型,数值算例结果验证了所提出的拓扑优化方法可以有效地进行复合材料微结构的导热性能设计。     

用应力杂交法计算Reissner型板复合型弯曲应力强度因子

本文在Reissner型平板裂纹尖端应力应变场一般解~[2,3]的基础上,采用应力杂交法构造高阶奇异元。由于采用了奇异元并解决了单元协调问题,使得Ⅱ型、Ⅲ型应力强度因子均获得较好的计算精度。本文计算了有限尺寸板复合型弯曲应力强度因子,并对Reissner型板若干力学特性进行了讨论分析。     

基于三维裂纹尖端应力场的应力强度因子计算方法

提出一种基于无限大体裂纹尖端弹性应力场理论解的前几项多项式函数,对实际裂纹体弹性应力场有限元解进行拟合来计算应力强度因子的方法.该方法在计算应力强度因子时不需要预先假设裂纹尖端的应力应变状态,应力强度因子计算结果更符合三维裂纹体裂纹尖端实际的应力应变状态.首先基于二维无限大板中心穿透裂纹应力场理论解验证了方法的有效性,探讨了拟合确定应力强度因子需要的多项式应力函数的项数.然后分别以二维大板中心穿透裂纹、三维大体内埋圆裂纹和三维有限厚板中心穿透裂纹的应力强度因子计算为例,通过与无限大板和无限大体应力强度因子理论解以及基于位移外推法和1/4节点张开位移法的应力强度因子有限元解的对比分析,验证了该方法的有效性和合理性.研究表明该方法能够合理反映三维裂纹体裂纹尖端的实际应力应变状态,计算得到的应力强度因子数值更合理.