受限亚音速气流中倒置悬臂壁板静气弹稳定性的理论及实验研究

板壳结构在航空航天、高速列车、能量采集等诸多工程领域已经得到了广泛应用. 将悬臂壁板倒置于轴向气流中并在壁板周围流场中设置刚性壁面可有效地调控壁板的失稳速度, 是俘能器优化设计的重要措施之一. 但针对刚性壁面作用下亚音速气流中倒置悬臂壁板的失稳机制仍需要开展深入研究. 本文以受限亚音速气流中倒置的二维悬臂壁板为对象, 以理论分析及风洞实验为手段, 研究了单侧刚性壁面效应对倒置悬臂壁板静态失稳特性的影响规律. 在理论分析中, 首先应用镜像函数法来处理壁面约束条件, 基于算子理论研究获得了以Possio积分方程为表征的壁板气动力, 壁面效应实际表征为一包含移位Tricomi算子的复合算子; 然后将壁板失稳方程的求解问题转化为定区间上的函数逼近问题; 最后, 依据Wererstrass定理并利用最小二乘法求解该最优函数, 以获得系统的失稳临界参数. 在试验研究中依据压杆稳定原理设计了壁板静态失稳的测试方法并完成了风洞实验. 理论分析结果表明, 壁板会发生发散(静气动弹性)失稳, 临界动压随壁板与壁面间距的增加而增大并最终趋于稳定(无壁面情况); 通过理论与风洞实验结果的对比分析, 验证了本文气动力及理论分析的适用性及准确性. 针对倒置悬臂壁板结构的气动弹性失稳问题, 本文提出的方法不涉及系统方程的离散及特征值求解问题, 而是将其转化为了定区间上的函数逼近问题进行求解, 这为弹性结构静气动弹性失稳问题的研究提供了一个可行的新思路.      

三维位势问题的梯度边界积分方程的新解法

三维位势问题的边界元分析中,关于坐标变量的边界位势梯度的计算是一个困难的问题. 已有一些方法着手解决这个问题,然而,这些方法需要复杂的理论推导和大量的数值计算. 本文提出求解一般边界位势梯度边界积分方程的辅助边值问题法. 该方法构造了与原边界值问题具有相同解域的辅助边值问题,该辅助边值问题具有已知解,因此通过求解此辅助边值问题,可获得梯度边界积分方程对应的系统矩阵,然后将此系统矩阵应用于求解原边值问题,求解过程非常简单,只需求解一个线性系统即可获得原边值问题的解. 值得注意的是,在求解原边值问题时,不再需要重新计算系统矩阵,因此辅助边值问题法的效率并不很差. 辅助边值问题法避免了强奇异积分的计算,具有数学理论简单、程序设计容易、计算精度高等优点,为坐标变量梯度边界积分方程的求解提供了一个新的途径. 3个标准的数值算例验证了方法的有效性.      

板中热弹波传播: 一种改进的勒让德多项式方法

近年来, 超声导波因其衰减小, 传播距离远和信号覆盖范围广, 成为无损检测领域快速发展的方向之一. 然而, 基于超声导波的高温在线检测和激光超声技术却发展缓慢, 其关键在于热弹耦合波动方程求解难度大、传播与衰减特性研究困难. 作为一种有效的求解方法, 勒让德正交多项式方法已广泛应用于导波传播问题, 但该方法在求解热弹导波传播时存在两个不足, 限制其进一步的发展和应用. 这两个缺陷是: (1)求解过程中大量积分的存在, 致使计算效率低下; (2)仅能处理等热边界条件的热弹导波传播. 针对两项不足之处, 提出一种改进的勒让德正交多项式方法, 以求解分数阶热弹板中的导波传播. 推导求解方法中积分的解析表达式, 以提高计算效率; 引入温度梯度展开式, 发展适合勒让德多项式级数的绝热边界条件处理方法. 与已有文献结果对比表明改进方法的正确性; 与已有方法的计算时间对比说明改进方法的高效性. 最后将改进的方法用于求解分数阶热弹板中的导波传播, 研究分数阶次对频散、衰减曲线和应力、位移、温度分布等的影响.      

双曲冷却塔塔筒水平地震响应分析

为明确冷却塔在水平地震下的内力环向分布特征及内在原因,同时探究不同地震波时程与规范反应谱所得内力差异的原因,以某大型双曲冷却塔为例,在动力特性分析的基础上,通过反应谱方法和时程方法的水平地震响应计算及对比分析,对上述两个问题进行了研究。研究表明:由于仅侧弯振型对水平地震有贡献,而塔筒的侧弯振型和实际响应均表现为整体侧倾并伴随微弱的截面“流动”变形,这也使塔筒各内力的环向分布分别呈现正弦、余弦分布特征;其整体侧倾可类比于悬臂杆结构,塔筒子午向轴力FY、子午向弯矩MY、剪力FXY和扭矩MXY的环向分布可借助悬臂杆侧倾时截面正应力和剪应力分布来解释;而截面“流动”变形则决定了环向轴力FX和环向弯矩MX的环向分布;整体侧移显著而截面变形极小也使FY和FXY的幅值在塔筒中下部明显大于FX;由于冷却塔第1阶侧弯振型在水平地震响应中往往起绝对主导作用,因此可先对所选地震波计算得到其反应谱,对比第1阶侧弯振型周期对应的水平地震影响系数α值,即可初步推断不同时程及规范反应谱方法所得结果的大小关系。     

非齐次动力学方程的一种精细积分单步方法

针对非齐次动力学方程■,结合精细积分法和微分求积法,利用同阶的显式龙格-库塔法对计算过程中待求的v_(k+i/s)(i=1,2,…,s)进行预估,提出了一种避免状态矩阵求逆的高效精细积分单步方法。该方法采用精细积分法计算e~(Ht),而Duhamel积分项采用s级s阶的时域微分求积法,计算格式统一且易于编程,可灵活实现变阶变步长。仿真结果表明,与其他单步法及预估校正-辛时间子域法进行数值比较,该方法具有高精度、高效率及良好的稳定性,在求解大规模动力系统时间响应问题中具有较大的优势。     

泰勒定理的妙用

论述泰勒定理在不等式证明、行列式计算、定积分计算及金融数学债券定价中的应用.     

形心在积分学中的应用

介绍了计算和证明积分问题的一种方法—形心法     

二维正交各向异性位势问题的高阶单元快速多极边界元法

研制了一种适用于二维正交各向异性位势问题的高阶单元(线性单元和二次单元)快速多极边界元法. 在快速多极边界元法中, 源点对于远场区域的积分采用快速多极展开式计算, 而对于近场区域的积分则直接进行计算. 高阶单元的使用使得近场积分, 尤其是奇异积分和几乎奇异积分的计算更加复杂. 通过引入复数表达对其进行简化, 若边界采用线性单元插值, 近场积分可直接解析计算; 若采用二次单元插值, 则给出一个半解析算法计算近场积分. 高阶单元奇异积分和几乎奇异积分计算难题的解决, 使得高阶单元快速多极边界元法不仅能够计算一般结构, 也能被应用于超薄体结构, 拓宽了高阶单元快速多极边界元法的适用范围. 数值算例表明, 若计算精度一定, 高阶单元快速多极边界元法较常值单元快速多极边界元法使用的单元数量显著减少, 且高阶单元快速多极边界元法计算时间与自由度数量成线性关系, 其计算效率仍处于$O(N)$量级, 因此高阶单元快速多极边界元法可更加高效求解大规模问题.      

山东大学(威海)开展的核物理研究

综述了山东大学威海校区原子核物理研究团队在原子核精细谱学、核天体物理、探测器研制和高能核物理等方向开展的研究工作及最新进展;尤其重点介绍了$A\sim 110$核区原子核的形状共存和带交叉延迟,“订书机”和“雨伞”模式转动带,碳氮氧循环过程中关键核反应的测量进展,中子星参数化的状态方程及双中子星并合引力波研究,带电粒子探测器的设计与制作,相对论重离子碰撞物理中量子输运理论和高阶反常输运等研究工作,并展望了下一步的工作重点。     

共价有机框架催化材料:二氧化碳转化与利用

多相催化二氧化碳(CO₂)高效转化为高附加值的精细有机化学品和化工燃料,具有非常重要的研究价值和工业应用潜力。共价有机框架材料(Covalent Organic Frameworks, COFs)由于其高比表面积、有序的孔道结构、高化学和热稳定性、可控的催化位点等特性,在CO₂吸附和转化方面都展现出突出的优势。通过合理的策略从分子水平对其孔道或表面进行改性引入功能性物种或催化活性位点,能够有效实现其对特定反应的选择性调控并提供物质传输的有利微环境,因此COFs材料近年来在CO₂催化转化领域中得到了迅速的发展,也具有良好的应用前景。本文针对近几年COFs材料在CO₂转化为重要化学品的研究领域进行了简要的综述,总结当前存在的问题,并对COFs催化材料今后在CO₂转化领域的发展进行了展望。