Proposal of new models of the bottom-side B0 and B1 parameters for IRI

The time series of hourly electron density profiles N(h) obtained from 27 ionosonde stations distributed world-wide have been used to obtain N(h) average profiles on a monthly basis and to extract the expected bottom-side parameters that define the IRI profile under quiet conditions. The time series embrace the time interval from 1998 to 2006, which practically contains the entire solar cycle 23. The Spherical Harmonic Analysis (SHA) has been used as an analytical technique for modeling globally the B0 and B1 parameters as general functions on a spherical surface. Due to the irregular longitudinal distribution of the stations over the globe, it has been assumed that the ionosphere remains approximately constant in form for a given day under quiet conditions for a particular coordinate system. Since the Earth rotates under a Sun-fixed system, the time differences have been considered to be equivalent to longitude differences. The time dependence has been represented by a two-degree Fourier expansion to model the annual and semiannual variations and the year-by-year analyses of the B0 and B1 have furnished nine sets of spherical harmonic coefficients for each parameter. The spatial–temporal yearly coefficients have been further expressed as linear functions of Rz12 to model the solar cycle dependence. The resultant analytical model provides a tool to predict B0 and B1 at any location distributed among the used range of latitudes (70°N–50°S) and at any time that improves the fit to the observed data with respect to IRI prediction.     

Behavior of the equivalent slab thickness over three European stations

The total electron content (TEC) derived from the global positioning system (GPS) and the F2-layer peak electron density obtained from Digisonde data have been used to study the diurnal, seasonal and solar activity variations of the ionospheric equivalent slab thickness (τ) over three European stations located at Pruhonice (50.0°N, 15.0°E), Ebro (40.8°N, 0.5°E) and El Arenosillo (37.1°N, 353.3°E). The diurnal variation of the τ is characterized by daytime values lower than nighttime ones for all seasons at low solar activity while daytime values larger than nighttime characterizes the diurnal variation for summer at high solar activity. A double peak is noticeable at dusk and at dawn, better expressed for winter at low solar activity. The seasonal variations of τ depend on local time and solar activity, the daytime values of τ increases from winter to summer whereas nighttime values of τ show the opposite. The effect of the solar activity on τ depends on local time and season, there being very sensitive for winter nighttime values of τ. The results of this study are compared with those presented by other authors.     

稠密气固两相流颗粒质量流量测量方法研究

为解决粉末火箭发动机和粉末冲压发动机供粉系统冷流标定过程中颗粒质量流量的测量问题,提出了一种稠密气固两相流中颗粒质量流量的测量方法,并基于旋风分离器和电子天平运用该方法设计了一套气固两相流中固体质量流量的测量装置,同时运用粉末火箭发动机的供粉系统对该测量方法和装置的测试性能进行了实验验证研究。结果表明:测量装置的粉末收集效率可高达98.5%以上,为固体颗粒的质量流量测量提供了保障;测量方法具有较好的可靠性和测量精度,其测量结果相对误差低于3%;测量方法具有较好的测量稳定性,重复性实验中测量结果相对误差低于2%;通过改变旋风分离器结构设计等参数,可拓宽测量方法的适用范围。     

液体火箭发动机超高转速泵浮动密封环研究

浮动密封环作为非接触式液体间隙密封,具有密封压差大、泄漏量较小、能自动调心和对中、理论上无磨损、对转子有辅助支撑作用等优点。但它在工作时受力状态很复杂,液膜形成和最小液膜厚度的大小受到许多因素的影响,泄漏量的计算与普通工作状态有较大的差别。本文通过对一些重点因素的分析研究,结合产品在不同工作转速和介质中的试验情况,提出了浮动环的工作模式,确定了最小液膜厚度和泄漏量的计算方法和在结构设计和产品生产过程中应重点考虑的问题。     

航天器介质盘环结构内带电特性三维仿真分析

针对航天器内特定结构的内带电（IDC）问题,以航天器典型复杂介质结构——太阳帆板驱动机构介质盘环为研究对象,开展了地球同步轨道恶劣充电环境（Flumic3）下介质内带电三维仿真分析。通过Geant4实现电荷输运模拟,根据电荷守恒定律数值计算得到电场,在此基础上,研究了屏蔽厚度对介质内带电的影响规律,提出了根据屏蔽厚度调整入射电子能谱能量下限以提高计算效率的方法。仿真结果表明,盘环结构内带电最严重部位是盘环最外圈上层介质与金属导电环接触的上边沿;增加屏蔽厚度可以减缓充电风险,但是随着温度降低,屏蔽效果会随之减弱。在地球同步轨道恶劣充电环境（Flumic3）下,当温度低至183 K时,由于辐射诱导电导率成为总电导率的主导部分,从而增大屏蔽的同时也会降低介质电导率,导致即使3 mm铝屏蔽下仍可能出现接近10⁷ V/m的场强峰值。     

基于CFD结果的轴流压气机附面层堵塞计算方法

为研究轴流压气机附面层堵塞计算方法,提出一种基于计算流体力学结果自动提取附面层厚度的多参数组合方法,该方法整合了多种附面层判断准则,可有效避免单参数判断准则的局限性。通过多个算例,包括二维平板、多级轴流压气机端壁、进口导流叶片以及带有分离区的静子叶片等,对该多参数判断准则进行了检验。结果表明,此附面层提取方法具有很好的有效性和普适性。该研究建立了多级轴流压气机中端壁堵塞从高维向低维模型的反馈,为压气机设计体系改进提供了更多可能。     

基于零厚度内聚力单元单向碳纤维增强树脂基复合材料微观切削机理研究

为探究碳纤维复合材料(CFRP)微观切削机理,通过有限元法,采用零厚度内聚力单元模拟界面相,碳纤维建模呈圆柱状并随机分布于基体中,以此来真实反应CFRP的微观结构。通过对各组成相设置不同的材料本构、材料失效和演化准则,对4种典型角度(0°、45°、90°、135°)进行直角切削仿真,探究不同纤维角度下单向碳纤维增强树脂基复合材料(UD-CFRP)在切削过程中的微观切削机理。结果表明:不同纤维角度下CFRP的微观破坏形式不同,切削0°CFRP时破坏主要以界面开裂和纤维折断为主,切削45°和90°CFRP时主要是刀具的侵入破坏,切削135°CFRP时则发生纤维的断裂和沿纤维方向的裂纹,纤维断裂点在刀刃下方。最后,通过实验验证了微观模型的准确性。     

P／M—SiC_p增强LY12的研究

 采用粉末冶金方法（Ｐ／Ｍ）完成了ＳｉＣ颗粒增强ＬＹ１２的成形工艺和不同重量百分比条件下的力学性能实验。测试表明,含ＳｉＣ＿ｐ２０％（ｗｔ）的复合材料在室温下平均抗拉强度可达５９４ＭＰａ,延伸率为６％。其抗拉强度比ＬＹ１２ＣＺ提高３０％；２００℃拉伸时σ＿ｂ＝４９３ＭＰａ；３１５℃时σ＿ｂ＝２９９ＭＰａ。用扫描电镜及光学显微镜观察了ＬＹ１２粉末和试样断口的组织形貌。     

团聚硼粉与HTPB混合物流变特性

利用HAAKE流变仪研究了硼粉团聚方法、粒度、填充比例对B/HTPB混合物流变性能影响。结果表明,对于LiF包覆的硼粉(B/F),小颗粒的LiF增强了与HTPB预聚物的物理吸附作用,使体系的表观黏度反而大于未包覆的硼粉;采用HTPB包覆的硼粉(B/H,B/F/H)与HTPB混合物的表观黏度和屈服值远小于未处理硼粉,且随混合时间增加而保持不变。在相同填充质量下,小粒度团聚硼粉的黏度系数K值较大,n值较小,严重偏离牛顿流体。当团聚硼粉/HTPB质量比为55/45时,混合物的流动方程转变为Herschel-Bulkley方程,大粒度团聚硼粉的屈服值较小。因此,与小粒度硼粉相比,大粒度团聚硼粉与HTPB的混合物具有较好的流变性能。含大粒度团聚硼粉富燃料推进剂药浆具有良好的工艺性能,6 h内表观黏度小于1500 Pa.s,屈服值小于150 Pa。     

复合材料变厚度加筋板后屈曲耐久性/损伤容限一体化设计研究

在先进的复合材料飞机结构上大量采用复合材料加筋板这种结构形式。因此,本文着重研究了复合材料变厚度加筋板后屈曲、冲击损伤与冲击损伤对承载能力的影响,以及复合材料变厚度加筋板冲击损伤、后屈曲、耐久性/损伤容限设计一体化综合试验方法。最后,作者根据多年从事飞机型号结构设计经验,并结合本文的研究结果,总结出15项复合材料变厚度加筋板后屈曲耐久性/损伤容限一体化设计技术,以期对我国预研新机的复合材料飞机结构设计,对已研飞机的复合材料飞机结构改进、评定有所启示。