求解矩形平板问题的通用应力函数

本文用变量分离法推导了求解矩形域平面问题的通用应力函数和用该通用应力函数表达的应力分量、应变分量以及位移分量公式。该通用应力函数可用来求解受任意载荷作用的矩形域平面问题。且所得解答不是放松边界条件的圣维南意义下的解答，而为精确级数解答。     

SH波对结构表面上含有多个半圆形介质夹杂的散射

利用复变函数和多极坐标系的方法,研究了在结构表面上含有多个半圆形介质夹杂的散射问题,首先给出 SH波对多个半圆形介质夹杂的散射波函数,并利用多极坐标系的方法,来满足介质界面上的连续条件和表面上的边界条件,归结为对一组无穷代数方程组的求解,简化了问题的求解。作为算例,对结构水平面上含有 2个半圆形夹杂介质进行分析,并讨论了夹杂介质参数的变化及孔距对位移的影响。     

受三点,四点弯曲半圆形单边裂纹SIF的研究

本文提出用以确定就力强度因子的一个新的由光弹性与权函数两种方法相结合的混合法，用此方法经济有效地求解了含单边裂纹半圆形板受三点及四点弯曲时的应力强度因子，并与和了比较，二者符合良好。同时还分别给出了应力强度因子随载荷位置和裂纹长度变化的情况。本方法为断裂力学的研究和工程应用提供了一个有效的手段。     

表面裂纹疲劳扩展和寿命计算的高效高精度数值分析方法

基于弹性力学理论,建立了三维裂纹应力强度因子计算的混合边界元法基本理论和数值求解技术;针对表面裂纹疲劳扩展过程中,需要计算每个裂纹扩展步下的应力强度因子,从而需要重复计算大型非对称系数矩阵问题,提出了仅在初始裂纹状态下一次计算主控矩阵,对于随后的疲劳裂纹扩展,只需做非常小规模矩阵的计算,且以显式形式给出应力强度因子解而无需求解大型线性代数方程组的方法,大大提高了计算效率;针对疲劳裂纹扩展过程中,单元需要不断重新划分问题,由于混合边界元法中的主控矩阵与裂纹无关,故只需对裂纹表面单元进行重新划分,对半椭圆表面裂纹,由于将其映射到单位半圆上划分单元,而单位半圆上的单元在疲劳扩展过程中不变,从而通过映射关系自动重新划分裂纹表面单元。最后,通过若干算例和试验,考核了本文方法的精度和可靠性。本文的研究为工程结构表面裂纹疲劳扩展和寿命计算的高效高精度数值分析建立了理论基础和实现方法。     

层合厚板混合状态Hamiltonian元的半解析解

给出任意角铺设或斜交铺设层合厚板的混合状态Ｈａｍｉｌｔｏｎ正则方程及其有效的半解析法。该方法在板厚方向未作任何有关应力或位移的人为假设,采用状态空间法给出解析解答,且通过状态转移矩阵的建立,保证了层间位移和应力的连续,减少了计算中未知量数目。     

具有任意边界的弹性力学梁固有特性的解法

 给出了可以用于求解具有任意边界条件的弹性力学梁的固有振动特性的双向展开方法和半解析方法。这两种方法均把梁看成为平面应力问题,不对纵向位移和横向位移作任何假设,涉及到的结构刚度和质量矩阵的元素均为显式。以简支梁固有特性求解问题为例,通过对本文结果和解析解进行比较,说明了给出的半解析方法和双向展开方法的精度和效率。     

在尾流激振情况下叶片振动应力预估技术

研究了前排静子叶片的尾流对后排转子叶片振动的影响.采用参数多项式方法和振荡流体力学理论, 求解静子叶片后的尾流场及尾流场作用下转子叶片通道内的非定常流场.把非定常气动分析给出的压力场, 转化为结构动力分析中的载荷压力场形式.根据试验结果得出响应分析中的阻尼, 先求解出各阶谐波作用下转子叶片的响应结果, 再把各阶响应结果进行叠加, 得到总响应.从流场求解到响应求解, 为工程应用初步建立起了一个尾流激振情况下叶片振动应力预估的半经验方法.      

基于均匀化方法的缝纫层压板拉伸弹性模量有限元预测

利用均匀化方法研究缝纫复合材料层压板的弹性性能，给出了相应的数学分析模型。建立了有限元分析的单胞模型，并进行求解。计算结果表明，利用这种半解析模型可以得到缝纫层压板的弹性模量的上下限，其平均值与实验结果平均值吻合得很好。     

纤维增强复合材料币状I型裂纹桥联问题

利用平面等效方法，将轴对称币状裂纹问题转化平面裂纹问题，由此给出复合材料内币状Ⅰ型裂纹的桥联积分控制方程。对权函数法应用于求解币状裂纹问题的可行性提出质疑，并给出修正的建议。然后对币状Ⅰ型裂纹桥联问题给予了断裂力学分析。利用一种新的分段插值方法求解了桥联积分控制方程。针对该积分程提出了一种近似分析解法，其结果与数值结果符合得相当好。     

Hamilton体系层合板动力响应半解析法

对时间方向引进正则变换函数, 给出了层合板动力学的Hamilton正则方程及其半解析法。通过分离变量, 就可以在板平面内采用通常的有限元离散, 而在板厚方向采用状态空间法给出解析解, 且通过传递矩阵, 保证了层间位移和应力的连续, 建立了层合板上下表面相变量的关系式, 利用打靶法进行求解。数值算例表明:本文方法未知量少, 精度高。      

非紧致气动噪声半空间传播基本解的边界积分方法

为了获得低马赫数流动诱发的非紧致气动噪声在半空间内传播的基本解,结合复等效源方法和边界元方法建立了半空间精确格林函数的边界积分方程,当半空间边界为刚度型阻抗边界时可避免奇异性积分。同时基于等效源方法提出一种半空间二维非紧致圆柱声散射模型,推导了静止介质中声散射基本解的理论表达式。对静止介质中的二维圆柱声散射,数值解在研究的频率与观察点处与理论解一致。采用数值方法计算二维圆柱绕流诱发的半空间声场基本解,结果显示半空间边界强化了声源的声辐射能力,小于马赫数0.2的介质运动对声传播的影响可以忽略。     

含孔边裂纹板的弯曲断裂计算

 采用复变函数理论和边界配置方法 ,对含孔边裂纹板的弯曲断裂进行了分析计算。首先假设挠度的复变函数式 ,进而可以求板的内力。它们能满足一系列的基本方程和支配条件 ,仅板的边界条件需要考虑 ,并且可用边界配置法和最小二乘法近似满足。对孔双边裂纹问题进行了应力强度因子计算。数值算例表明 ,本文方法精度较高 ,计算量小 ,是一种有效的半解析、半数值计算方法。     

考虑地面影响的半空间气动噪声数值计算方法

含非紧致边界的半空间气动噪声问题普遍存在于航空领域,在数值计算过程中不仅要考虑非紧致边界的散射效应,还需要考虑地面对声波的反射。基于Lighthill声模拟理论和广义格林函数方法,引入边界声压积分项和可压缩修正项可以考虑非紧致边界以及流体可压缩性对声波散射的影响。然而在半空间声场的计算过程中,包含沿无限大地面的积分项,为了克服沿地面积分所带来的巨大工作量,根据无穷大平板对声波的反射特性,推导得出频域下基于半空间格林函数的声学积分方程。同时开展了半空间中应力源噪声与圆柱绕流噪声的数值预测与分析,验证了该方法的准确性与高效性。     

Bump进气道中鼓包诱导的激波/边界层干扰特性

为了探索Bump进气道中鼓包诱导的锥形激波和机身发展而来的湍流边界层干扰问题,分析其气动优势,首先选取了半锥和半棱锥这两种与鼓包的流场结构具有一定相似性的构型作为参照,采用数值仿真方法,分别对这三类典型的三维激波/湍流边界层干扰问题进行了流场分析。在此基础之上,设计了三个不同马赫数的鼓包,并研究了设计马赫数对鼓包流场特性的影响。结果表明:当三类构型的无黏激波强度相等时,半锥诱导产生的旋涡强度最强,鼓包次之,半棱锥最弱。尽管鼓包诱导的流场非常复杂,其干扰流场却呈现出准锥形相似的特性。虽然半锥对边界层的排移能力最强,但是综合考虑边界层排移能力及进气道出口流场畸变下,鼓包最具优势,这也是其被选为超声速进气道前缘压缩面的重要原因之一。此外,在设计状态下,适当增加设计马赫数能改善鼓包排移边界层的能力,但设计马赫数太高,边界层排移能力基本不变,反而使得进气道总压损失急剧增加。      

任意翼型有小分离气泡跨音速绕流的迭代算法

 本文给出一种带小分离气泡的任意翼型粘性跨音速绕流的计算方法,采用有粘-无粘干扰迭代的概念。无粘流的全速势方程用AF差分格式在保角变换法生成的计算网格中求解,粘流附面层方程用C-S盒式法求解,用逆算法消除分离点处的奇性。本文对Ma_∞=0.8,Re_∞=2×10~6,迎角α=3.5°和4°的NACA64A010翼型粘性绕流进行了计算,结果与实验相比较,吻合良好。     

Predictive & Adaptive MPPT Perturb and Observe Method

The perturb and observe (P&O) best operation conditions are investigated in order to identify the edge efficiency performances of this most popular maximum power point tracking (MPPT) technique for photovoltaic (PV) applications. It is shown that P&O may guarantee top-level efficiency, provided that a proper predictive (by means of a parabolic interpolation of the last three operating points) and adaptive (based on the measure of the actual power) hill climbing strategy is adopted. The approach proposed is aimed at realizing, in addition to absolute best tracking performances, high robustness and promptness both in sunny and cloudy weather conditions. The power gain with respect to standard P&O technique is proved by means of simulation results and experimental measurements performed on a low power system. Besides the performance improvements, it is shown that the proposed approach allows possible reduction of hardware costs of analog-to-digital (A/D) converters used in the MPPT control circuitry.     

分流机匣对风扇/压气机性能影响的计算比较

采用流线曲率法、端部势流场分析加附面层修正等数值模拟手段,就分流机匣对风扇/压气机性能的影响展开研究,主要针对分流机匣端部中弧线形状、端弧形状,以及其轴向位置所引起的风扇/压气机特性变化规律进行了分析。研究结果表明:分流机匣几何参数及其轴向位置对整个风扇/压气机的性能有相当大的影响,对其几何参数的合理选取可以使风扇/压气机在较大的涵道分流比范围内,有效控制机匣内外表面附面层流动状况,在抑制紊流附面层分离,改善堵塞裕度方面有独特作用。      

Semi-analytic Solution of Steady Temperature Fields During Thin Plate Welding

Usually, it is very difficult to find out an analytical solution to thermal conduction problems during high temperature welding. Therefore, as an important numerical approach, the method of lines (MOLs) is introduced to solve the temperature field characterized by high gradients. The basic idea of the method is to semi-discretize the governing equation of the problem into a system of ordinary dif-ferential equations (ODEs) defined on discrete lines by means of the finite difference method, by which the thermal boundary condition with high gradients are directly embodied in formulation. Thus the temperature field can be obtained by solving the ODEs. As a numeri-cal example, the variation of an axisymmetrical temperature field along the plate thickness can be obtained.     

Numerical simulation of viscous flows with adaptively refined Cartesian grid

Quadtree-based Cartesian grid was automatically generated from specified geometry.Adaptive refinements were performed according to geometric parameters and solution of flow field.An altered CCST (curvature corrected symmetry technique) approach was proposed to apply solid wall boundary conditions.Driven flows in a square cavity and flows around NACA0012 airfoil were simulated and compared with the result of published structured grid and stretched Cartesian grid.The results show that solid wall boundary condition are accurately applied by current altered CCST approach,while incompressible/compressible subsonic,transonic and supersonic viscous flows are adequately simulated with adaptively refined Cartesian grid.