二阶迎风有限体积法方腔流数值模拟

介绍了一种基于非结构网格的二阶迎风有限体积离散格式，在对流项的离散过程中，为了达到二阶精度，在界面上对物理量咖作Taylor展开，在处理展开项中的梯度时经入NND格式的优点，克服了中心差分格式不稳定的缺点，并对方腔流动进行了系统的数值模拟．计算网格采用三角形网格，节点数为12960，单元数为25600．压强修正基于SIMPLEC方法．给出了雷诺数Re达到5000时的定常流动结果，与以往方腔流计算的标准解非常吻合，但所用的网格数要少．计算结果表明该二阶迎风有限体积算法有很好的收敛性和稳定性，且松弛因子对计算结果影响很小．     

二维非结构网格的一个TVD型有限体积方法

文章将有关文献提出的有限差分算法应用到有限体积方法中,在无结构任意三角形网格上建立了一种守恒型有限体积的高分辨率算法,该方法简单易行,应用改进后的算法到二维浅水波方程的求解上,进行了二维部分溃坝、圆坝两个算例的数值模拟,模拟结果表明:改进后的算法既有较高的激波捕捉能力,又能处理界面任意复杂的计算问题。     

非结构网格中子输运方程的离散纵标节块数值解法

基于离散节块方法（DNTM）的思想，研究了一种新的适用于二维三角形节块的离散坐标数值解法．使用坐标变换将任意三角形变换为正三角形，源项和通量的空间分布使用二次正交多项式展开，横向泄漏项的空间二次分布通过使用二元二次多项式近似正三角形节块的通量分布来获得、针对小尺寸易出现不收敛的情况，提出了一种改进的节块平衡有限差分方法．数值计算结果表明，该方法能满足非结构几何区域问题的求解需求，与细网SN方法相比，该方法在较粗的网格下具有高的精度，同时还适用于较细网格下的中子输运计算．     

NS方程的非结构化网格方法及其差分格式

采用同位非结构化网格上的有限体积法,对ＮａｖｉｅｒＳｔｏｋｅｓ 方程的ＳＩＭＰＬＥ 算法及差分格式进行了研究．提出了适用于非结构化网格的不用求解单元顶点变量值的二阶混合差分格式．该格式的优点在于：１） 减少了计算工作量;２） 避免了普通混合差分格式因使用简单的一阶迎风差分格式所引起的网格界面方向相关性的问题．最后,采用三角形网格,利用提出的方法及差分格式,对方腔内的驱动层流及绕翼型湍流进行了数值计算,计算结果与基准解或实验值的符合程度优于普通混合差分格式     

非结构化网格上N-S方程的离散

采用有限容积法离散控制方程,实现了二维非结构化同位网格上的SIMPLE算法,并编写了四边形网格求解不可压缩N-S方程组的程序.分别采用三角形网格与四边形网格,对计算传热学中有基准解的顶盖驱动空腔流问题进行了模拟.结果表明:所提算法既可用于三角形网格,又可用于四边形网格；三角形网格与四边形网格的模拟解均与基准解吻合,验证...     

TVD有限体积法在喷流计算中的应用

该文根据ＷａｎｇＪ．Ｃ．Ｔ．的二阶ＴＶＤ有限体积法求解守恒型Ｅｕｌｅｒ方程，对原文中的数值格式进行了部分修改。计算网格单元是非均匀的任意四边形。文中计算了楔形风洞中的斜激波反射问题和无粘喷流近场流动结构问题。数值计算结果表明，该文方法具有激波分辨率高、计算稳定、伪振荡小等优点。同时该文方法可以处理任意复杂边界的二维及轴对称流动问题的计算，也可以将该文的有限体积法思想推广到三维具有复杂边界的流动的计     

用改进的强隐式方法求解三维不可压粘性流动问题

提出了一种用改进的强隐式方法求解三维不可压缩粘性流动的交错网格有限体积法．该方法在ＳＩＭＰＥＲ算法的基础上，对Ｎ－Ｓ方程的对流项采用ＱＵＩＣＫ格式离散，利用ＭＳＩ方法求解代数方程．数值计算表明本方法具有计算稳定、收敛速度快等特点．     

Cauchy微分方程的一种新的数值解法

众所周知，求解微分方程（组）常用的数值方法有有限差分法，有限元素法等，这些方法都是将微分方程（组）分离散化后求解．若将网格划分得粗了，则求解精度不高，不能满足工程实际需要，若将网格划分得细了，则所需计算机内存量和计算量都太大．为解决上述问题，本文给出微分方程（组）的解的概率表达式的一种新的数值解法──概率数值解法．     

基于RNG k-ε湍流模型钝体绕流的数值模拟

采用RNG k-ε湍流模型对钝体绕流流场进行了数值模拟．计算采用非交错网格的有限体积法求解二维不可压N-S方程．计算结果表明，RNG k-ε湍流模型可以成功地模拟绕钝体的不稳定、非定常和剧烈分离流动。     

椭圆型方程问题的两套网格剖分的混合有限体积法

本文简要叙述了混合有限体积方法的相关结论,主要是三角形和四边形网格上的基于两套网格剖分的混合有限体积法,同时列出了相应格式的收敛性结果。     

防止高拱坝坝踵开裂和屈服的工程措施研究

对防止高拱坝坝踵开裂和屈服的两种常用工程措施周边缝和底缝进行研究．研究中坝体混凝土、缝和坝基材料采用不同的非线性材料模型．结合292m高的小湾拱坝对设缝前后的坝体位移场、应力场、开裂区和屈服区进行分析研究,重点讨论设缝对坝体屈服和开裂的影响．对小湾拱坝的研究结果表明：(a)未设缝时,在坝基附近的垫座内出现较大范围的开裂和屈服;设置周边缝时,垫座内的开裂和屈服明显减少;设置底缝时,坝踵附近的开裂和屈服明显减少,但在设缝两侧,尤其是两侧缝端附近,垫座内的屈服现象比无缝情况时更严重．(b)小湾拱坝设周边缝的效果要优于设底缝,设缝范围可限制在缝面开裂的范围内．     

碾压混凝土拱坝诱导缝断裂特性研究

通过数值计算分析了诱导缝的断裂特性;以实际工程为例,用断裂力学方法模拟了诱导缝的开裂过程,并对诱导缝的作用进行了探讨.结果表明:诱导缝开裂特性与诱导缝的位置、形式、两侧约束距离及相邻缝端的影响等因素有关;在诱导缝的开裂过程中,坝体应力有所降低,特别是诱导缝完全开裂时,应力下降非常明显;诱导缝设置可以首次开裂条件作为判别依据;为减少施工工序、提高施工效率,可在保证始终满足诱导缝开裂判别式的前提下适当加大中部诱导板的缝间距.     

底缝对拱坝工作性态的影响

结合一试验拱坝数值方法研究了底缝对拱坝的非线性影响,坝体、坝基和接缝均考虑了材料的弹塑性,坝体,地基处理用Ｄｒｕｃｋｅｒ－Ｐｒａｇｅｒ准则模拟,接缝材料用Ｍｏｈｒ－Ｃｏｕｌｏｍｂ准则模拟,模拟了接缝的滑移、开裂等非线性行为,计算结果和试验值一致性较好。     

含横缝碾压混凝土拱坝的变形和应力重分布

研究含横缝碾压混凝土拱坝工作状态。含横缝碾压混凝土拱坝未灌浆前若出现洪水强迫蓄水，横缝或诱导缝在施工末期（汛前）灌浆，由于蓄水后坝体进一步降温重新拉开灌浆缝都会破坏拱坝整体性。当缝宽不大时（＜０．５ｍｍ），在水压自重作用下梁向变位加大，可将部分拱重新压紧，形成应力重新分布。作者用体积向量提高简化计算中缝的变形场精度。坝踵水平裂缝出现后，拱作用逐渐加大，梁向大变位和接触形成新拱可使水平裂缝停止扩展。     

重力坝纵缝非连续接触地震反应分析

采用非线性动力时程分析方法,通过超载响应分析,初步研究了设纵缝重力坝的抗震能力.运用弥散裂缝模型模拟混凝土在强震下的开裂行为,采用基于罚函数法的接触模型模拟纵缝在动力荷载下的力学行为.对比不设纵缝的计算模型,结果表明纵缝改变了坝体应力分布,降低了坝体的抗震能力.通过对开裂比和抗滑安全系数的分析,判定该设纵缝重力坝段的极限抗震能力为0.6g.     

周边缝拱坝应力的弹性有限条元分析

本文应用壳体理论和有限条元法解周边缝拱坝的坝体变位与应力，导得坝体条元的刚度矩阵及等效结线荷载，编制了计算机程序和计算实例。     

沙牌碾压混凝土拱坝诱导缝开裂研究

根据碾压混凝土拱坝诱导缝构造特点,分别将矩形诱导缝简化为穿透形裂缝和深埋椭圆形裂缝,假定为I型开裂;考虑相邻裂缝的影响．采用虚拟裂缝模型和线弹性断裂力学相结合的方法,建立了诱导缝失稳扩展时的亚临界扩展量和等效应力强度因子的解析表达式．结合碾压混凝土的试验断裂参数,计算了沙牌碾压混凝土拱坝诱导缝的等效强度,并将此结果应用到三维非线性应力场仿真分析中,对沙牌拱坝进行了开裂计算分析．计算结果表明,等效强度可以作为诱导缝的开裂准则．     

横缝间隙对拱坝应力状态的影响分析

考虑拱坝分期施工过程,采用常刚度加速迭代算法研究了横缝间隙对小湾拱坝应力状态的影响．对灌浆后横缝没有间隙和存在0．5mm的间隙两种情况进行的对比计算和分析结果表明,横缝间隙的存在对坝体应力分布规律的影响不大,但对应力的大小有较大影响,即坝体上游面拉应力大小和拱冠梁底部拉应力范围均有明显增大．     

坝后背管与拱坝共同工作的有限元分析

本文给出用以分析拱坝及其背管共同工作的计算模型.在背管与坝体之间布置夹层单元以模拟施工的薄弱面,并根据施工程序分级加载,对结构进行非线性分析.它可以分析背管对坝体应力分布的影响,还能给出背管与坝体之间接缝处和背管本身的应力状态.通过对李家峡工程的计算,分析了这种结构的工作情况,并得出一些有意义的结论.     

拱坝横缝配筋动力响应的模型试验和计算分析

配置横向跨缝钢筋后,拱坝在动力荷载下的非线性响应是高拱坝领域的一个重要课题。该文在采用脆性材料制作的带横缝配筋拱坝模型的基础上,通过振动台模型试验研究了带横缝拱坝在谐波和地震波作用下的动力响应,并运用改进的拱坝非线性分析有限元程序ADAP-88(增加配筋数值模型)进行计算比较。试验结果验证了横向配筋数值分析模型的正确性,揭示出带横缝拱坝在地震荷载作用下的一些反应规律,包括配置跨缝钢筋可明显减小强震作用下横缝的张开度,防止止水遭到破坏,并明显提高拱坝的整体性。