直接用电密J计算导体三维涡流场的FE-BE耦合法

提出一种计算导体三维涡流场的FE-BE耦合法．在涡流区中直接采用涡流电密J为求解量的FEM, 其公式采用八节点长方体单元进行离散, 在非涡流区中采用磁场强度H为求解量的BEM, 其公式采用矩形单元离散. 根据交界面条件, 耦合得到一个求解变量仅为J的代数方程组. 该方法的优点就是所求量无需微分, 求解变量少, 耦合容易, 适用于多连域问题．对两个模型进行了计算, 计算结果与测试结果吻合较好.     

矩独立的基本变量重要性测度及其概率密度演化解法

为有效分析基本变量对可靠性分析中失效概率的影响程度,提出了一种矩独立的基本变量对系统失效概率的重要性测度指标η,分析了其相应性质,并给出了性质的证明.在此基础上,比较了基本变量对系统失效概率与基本变量对响应量分布密度的重要性测度,并建立了求解2种矩独立基本变量重要性测度的概率密度演化方法,有效地解决了求解2种矩独立基本变量重要性测度在计算上的困难.采用所建立的重要性测度概率密度演化方法,求解了数值算例和具体工程中基本变量对功能响应函数分布密度和失效概率的重要性测度指标,结果证明所提重要性测度可以从输入变量随机取值的分布密度角度衡量基本变量对失效概率影响的程度,同时,计算结果也表明采用所建立的求解重要性测度的密度演化方法是高效可行的.     

超临界机翼表面涡流发生器影响研究

通过数值求解Reynolds平均的Navier-Stokes方程组,对重要的流动控制手段之涡流发生器（vortex generator,简称VG）的绕流流场及其对主流的影响规律进行计算模拟.首先,预测平板上单一涡流发生器流动,验证数值计算方法,并认识含涡流发生器流场的基本流场特征;其次,预测标准模型——ONERA-M6机翼跨声速流动,探索激波/边界层干扰流动特征;再次,在超临界机翼25%当地弦长附近布置一排涡流发生器,探索它们对机翼跨声速流动边界层的干扰效应;最后,将这些涡流发生器位置提前（距前缘3.5%当地弦长）,检验其对低速大攻角流动的影响规律.结果表明,7个VG能有效抑制跨声速强激波/边界层干扰导致的分离,减小展向流动;也能大幅缩减低速大攻角状态下的分离范围.     

最少变量数边界元提取三维VLSI互连电感与电阻

为得到高频复杂结构三维互连电感与电阻需求解电磁场涡流问题. 提出一种三维涡流计算模型, 在导体与无限大自由空间中利用矢量磁位表示其中的电磁场, 得到间接边界积分方程, 并利用边界元离散化技术进行求解. 这种方法避免对导体体积区域进行离散; 作为一种间接边界元方法, 它具有最少的变量数; 由于不对导体中的电流方向进行限制, 可以计算互相垂直导体之间的互电感和互电阻. 数值计算结果利用2-D解析解进行了验证; 计算了互相垂直导体在不同频率和间距下的互电感和互电阻.     

对永磁MRI脉冲梯度产生涡流的背景磁场的解析

利用线电流和多层介质板模型以及镜像理论定量解析永磁MRI磁体系统中梯度线圈电流在铁磁极板中产生的磁场H和感应场B.磁极板中的这些磁通是产生涡流伪影的根源.为了充分抑制涡流,应该充分减小这些磁通.通过定量计算揭示"抗涡流板"的导磁率必须高于极板材料工业纯铁的导磁率,并且要有足够的厚度以保证其远离饱和.这些结论为进一步提高永磁MRI磁体系统技术水平指明了努力的方向.     

双绕组多相高速整流异步发电机系统的分析

在分析和定量计算双绕组多相高速整流异步发电机电流谐波和磁势谐波关系的基础上，提出了把发电机系统分成两个子系统分别进行分析计算的新思路．对十二相功率绕组与整流负载系统采用电路分析方法，用基于电路拓扑理论的回路电流法列写网络方程，求解过程中定步长和变步长联合使用，解决了计算结果的稳定性问题，得出了功率绕组电压、电流及其基波分量．对定子双绕组与实心鼠笼转子系统采用电磁场分析方法，电磁场有限元和多变量优化法迭代相结合求出了控制绕组电流和定子频率．计算结果和实验结果吻合，说明了所提方法的有效性和准确性．     

排水系统入海口新型消能工

通过实体模型试验研究了某排水系统入海口水工建筑物的水力性能，设计提出一种新型的墩栅涡流室消能工。该新型消能工的水力特点是：流出排水隧洞的高速水流大部分通过进口墩栅空隙流入涡流室，小部分在进口墩栅前发生水跃后扎入涡流室中产生强烈的顺时针漩涡，通过水跃、墩栅的阻力和涡流室漩涡的共同作用消能，并将来流的动能转换为势能，之后，涡流室中的大部分水流在重力作用下通过出口墩栅空隙流向下游，其余水流从墩栅项部溢流。这种新型的墩栅涡流室不仅消能效果良好而且出流均匀，并容许各种杂物如石块、泥沙等自由通过，避免涡流室中发生严重的淤积。     

近场动力学与有限元方法耦合求解复合材料损伤问题

本文提出了一种求解纤维增强复合材料破坏问题的近场动力学方法(peridynamics, PD)/有限单元法(finite element method, FEM)耦合方法.根据PD方法求解含裂纹等不连续问题以及FEM高效求解的优势,将模型划分为PD区域、FEM区域和耦合区域.其中包含裂纹的区域采用PD建模,其他区域采用FEM建模,两区域的结合部分为耦合区域.该耦合方法操作简单, PD粒子与FEM节点之间没有重叠区域. PD粒子与其域内所有粒子(包括PD粒子和FEM节点)以非局部方式连接, FEM节点与其周围的所有粒子以有限元方式相互作用.该耦合方案模拟的复合材料损伤状态与纯PD方法模拟结果相符,但该耦合方案可有效地提高计算效率.     

多工况拓扑优化问题的一种新解法——导重法

将导重法引入两类多工况拓扑优化问题的求解.首先介绍了导重法及其拉格朗日乘子的求法,并用对偶方法对多约束优化问题的拉格朗日乘子求法进行了改进;然后推导了用导重法求解最小柔度拓扑问题和最小质量拓扑优化问题的迭代公式并计算了相应的算例.算例计算结果表明,采用导重法求解多工况拓扑优化问题具有迭代公式简单、收敛速度快、求解效果好的优点.通过与Ansys中采用的SCP方法的计算结果进行比较可以看出,导重法是求解拓扑优化问题的一个行之有效的方法,它为拓扑优化问题的求解提供了一条新的途径.     

微型涡流发生器控制SCCH增升构型流动分离研究

针对等弦长带后掠半模SCCH增升构型襟翼附面层发生流动分离,线性段增升效率降低的问题,采用数值模拟和风洞试验方法,开展微型涡流发生器控制襟翼附面层流动分离研究.首先,采用数值方法分析SCCH增升构型的基本流动现象,获得襟翼流动分离特性,作为微型涡流发生器设计依据;其次,根据涡流发生器工作原理,结合SCCH增升构型襟翼结构与流动分离特性,提出本文微型涡流发生器设计思想和初步设计方案;再次,采用数值方法研究微型涡流发生器控制增升装置流动分离的作用机理,研究微型涡流发生器布置方式、弦向位置、安装角、高度、展向间距等几何参数对流动控制效能的影响规律,提出控制襟翼流动分离的微型涡流发生器设计方案,供风洞试验验证;完成数值设计之后,采用风洞试验方法,进行微型涡流发生器设计方案验证与可能的方案筛选,以验证数值模拟方法、设计方法及涡流发生器设计方案;最后,分析数值模拟和风洞试验研究结果,提出具有工程应用价值的增升装置微型涡流发生器设计原则、设计方法及技术路线,供型号研制借鉴与采纳.研究结果表明,本文针对SCCH着陆构型提出的微型涡流发生器设计方案,经CFD和风洞试验验证,在着陆及下滑进场飞行状态,最大增升与增阻量分别达到10%和14%,符合着陆飞行状态对增升装置的设计要求,且CFD方法提出的设计方案最佳,具有惟一性;同时也表明,本文提出的增升装置微型涡流发生器设计原则、设计方法及技术路线可用于型号研制.     

港口非线性波浪三维耦合计算模型研究

建立了外域用差分求解高阶Boussinesq方程、内域用有限元求解Laplace方程的三维非线性波浪对船作用的时域计算耦合模型.研究了该类三维耦合模型的匹配条件,耦合求解过程和内域、外域公共区域长度的确定,探讨了内域有限元网格的剖分方法.把该耦合模型的计算结果与实验结果、内域用Euler方程的耦合模型计算结果进行了对比,结果表明该耦合模型具有满意的精度,适用于模拟较大区域内波浪对三维船等固定物体的作用,为今后近海岸大区域非线性波浪对三维非规则物体作用的时域计算和三维分区计算提供了参考.     

基于伴随方法的网格自适应DG方法

发展了一套基于伴随方法的Euler方程自适应网格DG求解方法,采用当地网格上的输出变量的误差估计作为离散误差指示器驱动网格自适应.在离散误差指示器的构造中,采用GMRES方法求解伴随方程获得伴随变量,将p阶流场变量和伴随变量映射到p+1阶函数空间,并进行有限次块对角Jacobi迭代得到近似p+1阶细网格解.壁面边界通量及输出变量计算仅依赖边界外侧值以保证伴随相容性.首先采用NACA0012翼型不同精度的阻力系数结果验证了误差估计的可靠性,然后对NACA0012翼型亚声速、跨声速流场和圆柱6马赫高超声速绕流进行了自适应加密计算研究.研究结果表明,基于伴随方法的网格自适应可有效提高阻力系数等输出变量的模拟准确性,在亚声速NACA0012翼型计算中使用约17%的全局加密自由度获得了与全局加密精度相当的阻力系数.     

光在介质层中的路径散射及在下垫面的多次反弹

辐射传输方程通常用于解有限厚度水平均匀介质层中光的透射和反射。全当下垫面的边界条件下不能简化时，光在介质层中的辐射传输与下垫面的反射的相互影响，使方程的求解极为复杂和困难，其解也难于理解与验证，将这一复杂耦合系统分解为假定下垫面无反射时光的路径散射和光在介质层与其下垫面之间的多次反弹，导出了下垫面非Ｌａｍｂｅｒｔ时多次反弹的四通量完整公式，并进而对路径散射建立了有明确物理意义的解析近似，为建立三维     

基于变论域的分散模糊推理方法求解传热反问题

建立了基于变论域的分散模糊推理方法。首先,设计了一组分散的模糊推理单元,以各测点处温度的计算值和测量值之间的误差分别作为模糊推理单元的输入变量,同时根据误差的大小调整输入变量论域的值,进行模糊推理得到推理结果;然后,对各推理结果进行加权综合,产生各个待反演温度的补偿值并对各个待反演温度的猜测值进行修正;最后,实现导热反问题的求解。选用不同的初始猜测值、温度测点数目以及温度测量误差进行了数值试验,证实了基于变论域DFI算法的有效性。     

三门峡水库水沙调度过程的数值模拟研究

针对大部分水库调度模型中水沙输移计算大幅简化且与库水位计算过程间只有单向耦合的问题,本研究构建了基于激波捕捉式有限体积法与河网算法的库区水沙调度模型.该模型采用河段连接单元内水沙质量与动量守恒方程的准二维求解算法,河网算法可以反映干支流交汇角的影响,并且使用OpenMP技术实现了并行求解.该模型的水库调度模块将调度方案表示为5个触发条件和4个调度指令参数构成的调度规则表,由水沙动力学模型向调度模块提供入库监测断面的水沙条件与实时更新的水位库容关系,水库调度模块为水沙动力学模型提供库区下游边界条件,从而实现两者双向耦合.采用2020年三门峡库区实测水沙数据进行了模型验证,库区三个河段的水位、流量和含沙量计算结果与实测值相符,在坝前水位最快下降速度达到2.3 m/h的敞泄排沙时段,数值计算过程稳定且成功预测了出库沙峰.针对近年黄河来沙量减少而潼关高程下降趋势并不稳定的问题,应用模型开展了三门峡水库水沙调度过程的数值模拟研究,通过对不同的汛期和非汛期调度方案的模拟分析,初步确定了坝前控制水位对于年内库区冲淤量和潼关高程影响显著的时期,研究结果可以为三门峡水库调度方案进一步优化提供思路.     

基于自然边界元与有限元耦合的重叠型迭代算法

了更充分地利用有限元与自然边界元各自的优点，并尽量减少由于方法的局限性造成的在计算量及计算精度上的不足，本文引入了有限元与自然边界元耦合的方法。在文中简单介绍了自然边界元法及其与有限元耦合的原理，通过设置含有重叠区域的圆形人工边界，实现自然边界元与有限元的耦合，并把该方法应用到无界区域的实际算例中。从计算结果的比较中可以看出自然边界元与有限元耦合算法的收敛速度更快，充分体现了耦合法在解决无界区域问题上的优越性。     

近地轨道长时间多星交会任务混合规划

对近地轨道长时间多星交会问题,考虑基于光照条件的时间窗口与J2摄动影响,建立多星交会混合整数非线性规划模型,提出两层混合优化方法进行求解：上层问题以交会次序、交会轨道转移时间及交会后服务时间为设计变量,采用混合离散连续编码遗传算法求解,下层优化以单次交会轨道机动时间及机动冲量为设计变量,采用下山单纯形法求解;推导考虑如摄动轨道要素差分线性交会方程,并以该方程为基础提出下层优化近似策略,以提高计算效率．仿真结果表明：1）建立的近地轨道长时间多星交会任务的混合整数非线性规划模型是有效的,提出的两层混合优化方法可以获得满足时间窗口约束的满意解;2）推导的线性交会方程是考虑如摄动长时间轨道交会轨道良好的一阶近似;3）在如摄动影响下长时间交会问题无论在一个还是多个轨道周期的时间尺度上均具有多峰性,基于混合编码遗传算法求解多星交会任务规划问题符合问题特征．     

用于板材焊接变形预测的原位等效焊接单元

基于焊接接头热影响区以及其材料在拘束条件下热循环残余变形的行为,再引入焊接过程的激活机制,构建用于板材焊接变形预测的原位等效焊接单元,该单元具有焊缝的固有应变原位点、等效尺度、固有变形的弯曲半径（或角度）等特征参数,可以较准确地预测板材焊接过程中的变形;引用已有的薄板对焊变形的试验结果与采用原位等效焊接单元的计算进行验证,取得很好的结果,算例表明：原位等效焊接单元对于焊接结构的变形计算具有效率高、精度好等特点。     

输出-阈值耦合神经网络及基于此的最短路问题求解

通过建立神经元输出和与其相邻神经元阈值之间的相互耦合, 提出了一种输出-阈值耦合神经网络, 实现了对在通常脉冲耦合神经网络中所存在的自动波现象的模拟, 并基于它的自动波现象, 提出了一种用输出-阈值耦合神经网络求解最短路问题的方法, 该方法具有所需神经元数目少、神经元和网络的结构简单、大规模并行计算等特点, 可用于求解非对称赋权图单起点多终点的最短路问题, 其所需的计算量(迭代次数)仅正比于最短路的长度, 而与图的复杂程度、所存在的通路总数等无关. 最后给出了最短路求解的例子.     

气动热-气动弹性双向耦合的高超声速曲面壁板颤振分析方法

建立了气动热、气动弹性双向耦合高超声速二维曲面壁板颤振分析方法．基于柯西霍夫假设和冯卡门非线性应变．位移关系,建立了考虑几何非线性的二维简支曲板的气动．热．弹性分析方程;使用迦辽金方法对方程离散处理,采用四阶龙格．库塔法求解微分方程;三阶活塞理论用于气动力分析;使用参考温度法和平板气动热公式计算气动热．研究中重点考虑：1）气动热与气动弹性双向耦合,既分析气动热对结构刚度的影响,又分析气动弹性对气动热的影响;2）结构温度随飞行时间的积累效应;3）弦向和厚度方向非均匀温度分布的影响;4）曲面壁板的初始变形对壁板颤振发生时刻的影响．通过与传统的只考虑气动热．气动弹性单向耦合的分析结果进行对比,发现基于气动热和气动弹性双向耦合的壁板颤振分析结果更危险,这一点在精确分析中应当予以重视．