Shedding vortex simulation method based on viscous compensation technology research

Owing to the influence of the viscosity of the flow field, the strength of the shedding vortex decreases gradually in the process of backward propagation. Large-scale vortexes constantly break up, forming smaller vortexes. In engineering, when numerical simulation of vortex evolution process is carried out, a large grid is needed to be arranged in the area of outflow field far from the boundary layer in order to ensure the calculation efficiency. As a result, small scale vortexes at the far end of the flow field cannot be captured by the sparse grid in this region, resulting in the dissipation or even disappearance of vortexes. In this paper, the effect of grid scale is quantified and compared with the viscous effect through theoretical derivation. The theoretical relationship between the mesh viscosity and the original viscosity of the flow field is established, and the viscosity term in the turbulence model is modified. This method proves to be able to effectively improve the intensity of small-scale shedding vortexes at the far end of the flow field under the condition of sparse grid. The error between the simulation results and the results obtained by using fine mesh is greatly reduced, the calculation time is shortened, and the high-precision and efficient simulation of the flow field is realized.     

钝头体中的广义雷诺比拟关系

钝头体壁面的摩阻和热流分布规律不同,平板流动中的雷诺比拟关系在钝头体壁面失效. 文章在前期高超声速广义雷诺比拟理论研究工作的基础上,利用数值仿真的方法对不同外形和来流参数条件下的钝头体广义雷诺比拟关系开展进一步研究. 通过建立钝头体绕流边界层的理论分析模型,得到了钝头体壁面雷诺比拟系数的线性分布预示公式. 采用数值求解 N-S 方程的方法,计算了圆柱和幂次体壁面的摩阻和热流以及二者之间的比拟系数. 通过与前期数值和理论结果对比,以及计算收敛性和网格无关性检验,对数值方法进行了验证. 通过在不同雷诺数 ($Re_\infty = 3.98\times 10^2 \sim 1.59\times 10^6$) 和马赫数 ($M_\infty = 3\sim 12$) 条件下的计算结果对比分析雷诺比拟系数的分布,总结了钝头体中广义雷诺比拟关系受外形和来流条件的影响,评估了广义雷诺比拟理论的适用性. 研究发现,在较高雷诺数条件下,离驻点较远的下游 ($\theta > 60^\circ$) 部位,雷诺比拟系数的分布不同程度地偏离理论预示的线性规律. 相比于圆柱外形,幂次体壁面的雷诺比拟系数分布的线性规律相对较好,其分布斜率略低于圆柱壁面的结果. 研究表明,如果针对实际外形和雷诺数进行适当修正,可以提高广义雷诺比拟关系的预示精度.      

缝合式C/SiC复合材料非线性本构关系及断裂行为研究

C/SiC复合材料具有高比强度、高比模量和优良的热稳定性能等一系列优点, 广泛应用于航空航天领域中. 裂纹扩展进而引起的脆性断裂是其主要失效形式之一, 因而材料的断裂性能分析对材料的结构设计和应用有重要的指导意义. 本文开展了缝合式C/SiC复合材料简单力学试验和断裂试验, 研究了材料在不同载荷下的力学响应及断裂特征. 基于缝合式C/SiC复合材料简单力学试验, 建立了材料宏观非线性损伤本构方程, 并模拟了缝合式C/SiC复合材料单边切口梁和双悬臂梁的断裂行为. 本构方程采用简单函数描述了材料在复杂应力状态下的非线性应力-应变曲线, 并考虑了反向加载过程中造成的裂纹闭合. 基于商业有限元软件ABAQUS, 通过编写UMAT子程序实现非线性损伤本构方程, 采用单个单元验证了建立的本构方程的有效性. 在此基础上, 采用线弹性损伤本构和非线性损伤本构分别模拟了缝合式C/SiC复合材料单边切口梁和双悬臂梁的断裂行为. 采用非线性损伤本构方程模拟的力-位移曲线结果与试验结果更为吻合, 非线性损伤本构预测的失效载荷与试验失效载荷更为接近, 验证了所建立的非线性损伤本构方程的准确性, 为C/SiC复合材料断裂行为的研究提供了借鉴, 为缝合式C/SiC复合材料结构的设计和应用提供了理论基础.      

介质与静压对激波管校准压阻式绝压传感器动态特性的影响

激波管通常用于动态压力传感器的校准，压阻式绝压传感器在激波管校准过程当中，会出现谐振频率等动态性能指标随着激波管静态压力环境、气体介质变化而改变的情况，影响传感器动态特性的校准。基于压阻式传感器的工作原理，对传感器的敏感膜片结构进行了机理分析，建立了膜片结构与校准环境中介质和静压关系的动态模型；通过ANSYS与SIMULINK软件开展了数值模拟验证工作，模拟结果与理论推导一致。通过激波管校准实验验证了气体介质与静压的影响关系，结果表明:传感器的谐振频率与静压间存在非线性关系，并且随着敏感膜片径厚比的增大而显著增大；系统阻尼比大小与气体介质有关，随着气体密度的降低而升高；传感器的灵敏度与气体介质和静压无太大直接关系。在使用激波管校准压阻式绝压传感器时，应当考虑介质与静压参数对校准结果的影响。     

基于相场描述的拉压不对称强度结构拓扑优化

运用了基于相场描述的拓扑优化方法，来寻找在拉伸和压缩中表现出不对称强度行为的连续体结构的最优布局。依据Drucker-Prager屈服准则和幂率插值方案，优化问题可以描述为在局部应力约束下的最小化结构的体积。用qp放松法来解决应力约束的奇异性，并采用基于P-norm函数的聚合方法对应力约束进行凝聚，该方法实现了约束个数的降低，同时引入了稳定转化法来处理大量的局部应力约束和高度非线性的应力行为，以修正应力，提高优化收敛的稳定性。在优化问题求解时，使用拉格朗日乘子法对目标函数和应力约束进行处理。利用伴随变量法进行灵敏度分析，并通过求解Allen-Cahn方程更新相场函数设计变量。数值算例证明了该优化模型和相应数值技术的有效性，相关算例还揭示了考虑拉压不同强度和考虑同拉压强度约束时得到的结构优化拓扑构型具有显著的差异。     

基于遗传算法的升力体外形优化设计

升力体由于低热流率再入物理特性和高效的内部容积利用率，是高超声速飞行器气动外形的一种典型布局.文章对升力体飞行器进行参数化数值建模，并提取其表征外形的参数作为设计变量，综合考虑飞行器再入过程中的气动力、气动热、容积利用率及稳定性等性能指标.运用多目标混合遗传算法对升力体进行了多变量、多约束下的气动外形优化设计，获得了再入飞行器外形的最优Pareto解.数值模拟结果表明，典型状态下最优Pareto解与CFD结果相差12%，验证了优化结果的准确性.      

两个dpa类配体的铜(Ⅱ)配合物的合成、结构、核酸酶活性及细胞毒性

以dpa衍生配体4-methyl-N,N-bis（pyridin-2-ylmethyl）aniline（L）合成2个单核铜配合物[CuL（NO₃）₂]（1）和[CuL（OAc）（H₂O）]ClO₄（2）,并对其进行了表征。单晶结构显示,配合物1中的Cu中心可以描述为畸变的五角双锥构型,而2的Cu中心为畸变的八面体构型。运用电子吸收和发射光谱法研究了配合物与CT-DNA的键合作用,结果表明2个配合物与DNA的相互作用均为部分插入模式。通过改变浓度、时间进一步检测配合物切割DNA的能力,以及验证其切割机理,结果表明在外界诱导剂存在下,2个配合物均表现出强的切割DNA的能力,其作用机理为氧化切割机理,其中活性氧可能为·OH和¹O₂。利用MTT法测定了配合物对体外HeLa、HepG-2和SGC-7901肿瘤细胞增殖的抑制能力。     

Phenylethynyl-terminated Imide Oligomers Modified by Reactive Diluent for Resin Transfer Molding Application

To meet the processing requirements of resin transfer moulding(RTM)technology,reactive diluent containing m-phenylene moiety was synthesized to physically mixed with phenylethynyl terminated cooligoimides with well-designed molecular weights of 1500-2500 g/mol derived from 4,4’-(hexafluoroisopropylidene)diphthalic anhydride(6 FDA),3,4’-oxydianiline(3,4’-ODA)and m-phenylenediamine(m-PDA).This blend shows low minimum melting viscosity(<1 Pa·s)and enlarged processing temperature window(260–361℃).FPI-R-1 stays below 1 Pa·s for2 h at 270℃.The relationship between the molecular weight of the blend and its melting stability was first explored.Blending oligoimides with lower molecular weights exhibit better melting stability.Upon curing at 380℃for 2 h,the thermosetting polyimide resin demonstrates superior heat resistance(T_g=420-426℃).     

ANALYSIS OF FRICTION INDUCED VIBRATION UNDER THE ACTIVE DISTURBANCE REJECTION CONTROL FRAMEWORK 1)

自抗扰控制(active disturbance rejection control, ADRC)是一种具有两自由度控制结构的工程化方法, 由于其能够直观有效地处理多种扰动, 近些年来在许多机电系统上得到了成功应用. 当采用ADRC对带有摩擦力的机电系统进行调节时, 可能会产生极限环振动. 目前, 还没有ADRC框架下摩擦力振动精确分析的相关工作. 因此, 本文采用非线性动力学系统的分析工具对这一问题进行研究. 首先, 考虑两种典型摩擦力模型, 静态切换模型和动态LuGre 模型, 对一类二阶运动系统设计不同阶次的ADRC, 得到控制器的等效形式, 并揭示出与比例积分微分(proportional-integral-derivative, PID)控制之间的联系. 然后, 采用打靶法结合拟弧长延拓方法求解系统中的极限环, 并根据Floquet理论判断极限环的稳定性、可能出现的分岔以及分岔类型. 此外, 通过雅克比矩阵和近似数值方法对系统平衡点集的局部稳定性进行了分析. 最后, 通过数值计算研究了摩擦力模型和参数、ADRC阶次和参数对极限环和平衡点集的影响. 计算结果表明, 决定摩擦力Stribeck效应负斜率的参数$\beta$作用较大. 当$\beta>1$时, 两种摩擦力模型下的闭环系统呈现出相同的特性, 极限环会出现环面折叠分岔(cyclic fold bifurcation, CFB)且平衡点集是局部稳定的. 然而当$\beta<1$时, 两种闭环系统呈现出完全不同的特性. 此外, 不同阶次的ADRC在极限环的存在性和稳定性、平衡点集的稳定性上面的结论是相同的, 而低阶次的ADRC能够更好地解决摩擦力补偿和稳定鲁棒性之间的矛盾问题. 这些结论对实际现象的理解、ADRC阶次的选择以及参数整定提供了一定指导.     

STABILITY OF CIRCULAR RING OF CIRCULAR CROSS-SECTION UNDER UNIFORMLY DISTRIBUTED TORSION 1)

圆环在均布扭力矩作用下的变形特殊且容易发生跳跃,不考虑稳定性而得到的一些结论将可能不存在。本文利用最小势能原理导出圆截面圆环在均布扭力矩作用下的平衡路径,并根据系统稳定性的能量判据,对圆环的平衡稳定性、变形与运动过程进行分析,最后得到稳定平衡状态下圆环横截面上的内力最大值。