广义有限元法计算飞机开口裂纹应力强度因子

机身开口圆角损伤可能逐渐形成裂纹从而严重影响飞机使用安全,为此对圆角裂纹进行剩余强度估算显得十分重要。使用Williams位移场广义有限元法进行裂尖单元加强,并对裂尖区域应用放射形网格划分技术,编写广义有限元法的MATLAB程序,研究并计算飞机舱门开口圆角裂纹的应力强度因子。结果表明：在裂纹分析上,广义有限元法较常规有限元法具有更好的适应性和更高的精度。     

基于高阶间断伽辽金数值方法的浸没边界法

为了降低复杂边界条件下贴体网格生成难度，近年来基于笛卡尔网格的浸没边界法逐渐成为研究此类问题的主要数值方法之一，然而计算精度和计算效率仍然是此类方法目前面临的挑战。与传统的空间2阶精度的有限体积格式相比，空间精度为3阶或3阶以上的高精度方法具有空间精度高、数值分辨率高、数值耗散小的优点，而作为高精度数值方法之一的间断伽辽金方法在浸没边界方面的应用仍较少。本文将高阶间断伽辽金方法的高精度优势与浸没边界法无需贴体网格的优势结合起来，提出了适用于可压缩流动的高精度浸没边界法。其中边界条件采用体积惩罚方式实现，同时采用牛顿法迭代以及MPI并行提高计算效率，物面的数据重建采用插值点处的逆距离权重（IDW-IP）方法代替高阶格式下多项式插值方法。本文基于笛卡尔网格测试验证了二维定常和非定常情况下的数值模拟效果，并与传统贴体网格的计算结果进行了对比。     

用于空间探测的静电分析器轻量化设计与验证

针对空间探测载荷的轻量化需求，以空间等离子体探测的静电分析器为研究对象，进行了轻量化设计和验证。在材料及加工方式选择方面，充分利用增材制造这一新兴加工手段，打破传统轻量化设计时追求低密度材料的思路，对比多种材料及加工方式下设备重量及力学表现，确定铝合金为主要材料，3D打印为主要加工方式。在结构设计方面，基于3D打印加工的优势，在减小设备各部分结构厚度的同时，适当设置加强筋以解决薄壁在后处理时易产生形变的问题。采取上述方案设计的设备质量减少至1.2 kg，比使用镁合金及传统加工方式的设备质量（2.2 kg）下降45%。以典型航天任务的鉴定级力学试验条件作为输入，开展了设计模型的有限元仿真，完成实物加工装配以及力学试验，验证了该设计的抗力学性能。     

基于折纸的大型空间薄膜遮阳结构设计与分析（英文）

提出了一种轻质、展开顺畅、收纳比高且不需要额外驱动动力的新型自展开薄膜遮阳结构。利用4根等长的带状弹簧驱动叶内环绕折叠的聚酰亚胺薄膜，并计算了收纳比。为了防止薄膜在展开过程中出现钩挂和破损的问题，对薄膜展开的轨迹、薄膜的应力分布以及带状弹簧骨架的模态分布进行了有限元仿真，并制作了测试样机以验证设计的可行性。结果显示，利用带状弹簧驱动薄膜展开，展开过程顺利，带状弹簧展开时的冲击不会造成薄膜损伤，薄膜折痕处的最大应力约为42.3 MPa，低于材料的许用应力，且收纳比达到了69.4。     

星载双反射抛物面天线热变形分析

星载天线在轨运行时受到空间外热流的影响，会经历周期性的高低温交变，导致其反射面产生热变形。为了保证星载天线稳定运行，选取某型号星载双反射抛物面天线作为研究对象，采用有限元法对天线高温工况下的在轨温度场进行分析，进而将天线上分布的温度载荷作为边界条件映射到结构场中进行热变形分析；同时详细分析了材料属性、铝蜂窝芯厚度、碳纤维贴层厚度、反射面支撑约束位置等因素对双反射抛物面天线热变形的影响，以期为星载天线结构优化设计提供理论参考。     

涤纶增强橡胶基复合材料各向异性超弹性本构研究

基于一种适用于平纹涤纶增强橡胶复合材料的各向异性超弹性本构模型，将应变能分解为橡胶基体应变能、织物纤维拉伸应变能与织物增强橡胶剪切应变能3部分，并根据单轴拉伸试验数据确定了本构模型参数。编写了有限元材料子程序进行仿真分析，并与试验数据对比验证了本构模型的合理性。该模型从宏观出发，能更好地表征复合材料编织物在拉伸过程中由于大变形所引起的非线性各向异性力学行为，具有结果准确、简单实用等优点，为织物增强橡胶复合材料的设计应用提供了理论依据。     

航空磁探中潜艇磁场建模方法分析

为提高航空磁探中潜艇磁场模型精度,分析现有潜艇磁场建模方法,对边界积分法、有限元法、积分方程法、磁体模拟法等主要潜艇磁场建模方法的基本原理和推导过程进行分析。通过数值分析比较各方法的建模精度,最后归纳潜艇低空磁场建模方法的基本要点。     

非结构网格通量重构算法下三种紧致WENO限制器对比研究

在通量重构算法框架下对比研究了3种最新提出的用于间断伽辽金算法的紧致型WENO(weighted essentially non-oscillatory)限制器，即简单WENO限制器、p阶加权WENO限制器和多精度WENO限制器。这3种WENO限制器能够较高精度地模拟流场并捕捉激波，其紧致性在于问题网格单元中的重构只涉及问题网格单元及其相邻网格单元。同时，使用保精度的保正限制器用以避免可能出现的非物理的密度和压力负值。最后，采用非结构四边形网格在双马赫反射、激波与涡相互作用、激波反射、黏性弓形激波和激波与层流边界层相互作用等多个二维算例中，对这3种WENO限制器进行分析比较。结果显示：多精度WENO限制器与p阶加权WENO限制器能够高精度模拟流场并捕捉激波，同时一定程度抑制通量重构方法本身在激波附近的数值伪振荡；p阶加权WENO限制器与多精度WENO限制器相比，其稳态收敛性相对更好；简单WENO限制器则性能较差。根据以上研究，提出亟须发展能够使高阶WENO限制器在稳态问题中收敛的间断探测器。     

大型空间可展开结构热致振动研究

搭载高精度光学载荷的空间飞行器在进出地球阴影区时,飞行器上处于展开状态的太阳帆板因剧烈的温度变化会诱发振动问题.本文采用有限元法,首先对太阳帆板进行在轨进出阴影区瞬态温度求解;然后基于卫星-太阳帆板整体动力学模型,将时变温度场等效为时变热载荷加载到整星系统上,对整星系统的热致振动动力学响应进行了数值模拟.结果表明:处于...     

电阻应变计敏感栅结构参数对其测量精度的影响

为了提高电阻应变计测量精度，建立了“简支梁-过渡层-基底-栅丝”结构的3维有限元模型，利用数值模拟方法研究了 敏感栅结构参数对测量精度的影响；提出横向效应长度比参数，统一了栅丝长度、间距、弯数3种结构参数，建立了测量误差与横 向效应长度比的函数关系，基于此函数关系推导出结构参数优化公式，并利用优化公式完成了敏感栅结构参数组合设计。结果表 明：栅丝直径越小，测量误差越小；栅丝长度、间距、弯数、长宽比存在最优值；测量误差与横向效应长度比成分段线性函数关系，且 当横向效应长度比为6.15%时测量误差最小；利用函数关系式可以快速估算出给定敏感栅结构参数的应变计测量误差，估算结果 与有限元计算结果的相对误差不超过4%；利用优化公式设计出的参数组合的测量误差不超过0.06%，且较单独优化的最优值参 数组合的减小0.1%左右。