热变形求解及其对高速点接触弹流润滑影响研究

提出了一种固体表面热变形求解新方法(ITD),由此研究了热变形对高速点接触弹流润滑行为的影响. 为此,基于计入流体惯性项的Reynolds方程获得了油膜压力,采用追赶法对润滑剂和接触固体的温度进行了求解,进而研究了不同工况下有无热变形的高速点接触非牛顿热弹流润滑性能. 采用有限元法和离散累加法对ITD法进行了验证,通过中心膜厚试验验证了考虑热变形的正确性. 结果表明:ITD法可准确快速地计算表面热变形;考虑热变形后,油膜厚度降低且向油膜出口倾斜,考虑热变形后的中心膜厚更接近试验结果.      

自适应无网格热弹塑性接触模型研究

提出一种自适应无网格热弹塑性接触求解模型,求解接触问题的线性规划-增量初应力法与基于应变能梯度的自适应无网格法相结合,给出了模型计算理论和算法实现.通过圆柱体与弹塑性平面热弹塑性接触算例对模型进行验证.对是否考虑材料应变硬化,是否考虑摩擦力和热输入,是否考虑材料屈服强度温度相关等情况的两种算例进行了讨论.结果表明,该模型能有效地求解考虑不同情况下的热弹塑性接触问题,在较真实地模拟接触状况的同时,具有较高的计算精度和计算效率.     

热-应力耦合结构灵敏度分析方法

研究稳态／瞬态热传导灵敏度分析、以及热与机械荷载同时作用的热结构应力灵敏度分析问题。考虑了温度场随设计变量的变化及其对应力的影响,提出温度场与结构热应力耦合问题的灵敏度计算方法。特别指出了热－应力耦合灵敏度分析中温度场导数的影响, 说明了在热－应力耦合结构灵敏度分析中必须考虑耦合灵敏度。在应用软件系统JIFEX中实现了所提出的方法,数值算例验证了灵敏度算法的精度。     

水流损失和热补偿共同作用对增强型地热系统(EGS)产能影响的研究

基于青海共和盆地-3705m地热田实测数据,结合流固耦合传热理论并运用Comsol软件,建立了离散型裂隙岩体流体传热模型。考虑水流损失和热补偿共同作用,模拟得到了开采过程中上、下岩层(盖层和垫层)为绝热不渗透、传热不渗透、渗透传热时,储层(上、下岩层和压裂层)温度场的变化特征,分析了产出流量、水流损失、产出温度、产热速率的变化规律。研究结果表明:采热过程中产出流量始终小于注入流量;产出流量增幅速率先增大后减小,最后趋于稳定,前3a产出流量增幅超过总增幅量的3/4;忽略水流损失,将高估产热速率,采热初期甚至达到考虑水流损失时产热速率的3倍以上;考虑水流损失,产热速率呈先快速上升再趋于稳定后逐渐下降的趋势,最优开采时间为3a^11a;研究上、下岩层对产出温度的影响,仅考虑传热,采热寿命延长5.43%,同时考虑渗流传热时,采热寿命延长2.71%;采热前9a,水流损失占主导作用,即流入上、下岩层水流损失对产热速率的影响高于热补偿效应,开采10a后,热补偿效应占主导作用;同时考虑水流损失和热补偿效应得到的产热速率变化规律与实际工程更为符合,建议选择低渗透能力的上、下岩层延长增强型地热系统(EGS)运行时间。     

考虑热渗效应的有限长圆柱热固结解

在作者已建立的饱和多孔介质耦合非线性热弹性理论基础上,考虑热渗效应,建立了饱和多孔介质耦合热弹性固结方程,并推导了有限长圆柱热固结问题的解析解,进而以温控三轴试验的试样为例进行了算例分析,同时利用COMSOL软件进行了数值模拟,并将解析结果和数值结果进行对比．结果表明：在不排水条件下,影响试样最终孔压大小的参数是：土的泊松比>弹性模量>水的体膨胀系数,渗透系数对孔压变化影响不大,考虑与不考虑水土压缩性不影响孔压的计算;相对线弹性情况,考虑非线弹性膨胀系数随温度变化时的孔压有所下降,轴向应变变化很小．     

考虑摩擦系数和微凸体相互作用的粗糙表面接触热导分形模型

基于三维分形理论,建立了考虑摩擦系数和微凸体相互作用的粗糙表面接触热导分形模型,并且考虑了微凸体的弹性变形、弹塑性变形和完全塑性变形. 通过该模型,分析了摩擦系数、分形维数、分形粗糙度和接触载荷对热接触热导的影响. 研究结果表明:接触热导随着摩擦系数和分形粗糙度的增大而减小,随着分形维数和接触载荷的增大而增大. 该研究为开展接合面的热传递提供了一定的理论基础.      

热变形下干气密封系统角向摆动稳定性分析

对考虑热变形下气膜的温度场分布进行分析,建立了气膜-密封环系统角向摆动计算模型,进而计算动环变形量。用Maple软件拟合得到了变形后的气膜厚度,继而对角向刚度和角向阻尼进行拟合,得到了一个含二次、三次项的非线性受迫振动方程。以Floquet指数为准则对系统进行稳定性分析,得出了系统稳定运行时的螺旋角范围,即当螺旋角α=75°9′55″时系统发生Hopf分岔。结果表明:在一定的精度下,考虑热变形对系统的稳定运行有一定影响,螺旋角的范围较不考虑热变形时略有减小。     

内压柱壳在激光束照射下热爆破数值模拟

为研究内压圆柱壳在激光束照射下的破坏现象提出了一个高效可行的数值分析模型 :采用了子结构方法把非线性区局限在较小的范围 ;采用了与应变速率有关的耦合的热传导分析 ;考虑了材料物理性质随温度的变化 ;考虑了热变形对结构刚度和荷载的几何非线性效应 ;考虑了内压是予先加载以及随后的热弹塑性过程。数值结果与文献 [1 ]的实验大致吻合     

航空发动机叶片-机匣碰摩热-结构耦合仿真

针对转静子碰摩过程中的摩擦热效应现象,建立了叶片-机匣的碰摩热-结构耦合模型。采用热-结构耦合单元,对航空发动机叶片-机匣进行了碰摩热效应特性研究。考虑摩擦因数及旋转过程中的离心力作用,对瞬态-热结构耦合场进行了有限元分析,研究了结构在温度场作用下的热-结构耦合应力和温度分布。与纯机械载荷下的应力分布对比,发现了碰摩热效应在叶片与机匣上的扩散规律。研究结果表明,考虑摩擦热效应时由于热应力的作用,导致结构的总应力水平升高进而产生热变形,从而使故障进一步恶化。由此可见,碰摩热效应的影响在实际航空发动机振动分析中不能忽视。     

脆性材料热-力耦合模型及热破裂数值分析方法

针对混凝土、岩石等脆性材料,利用热传导和热-力耦合的相关理论,并结合材料在细观尺度上的损伤演化规律,提出了一种考虑损伤的热-力耦合模型,并在原有材料破坏过程分析系统RFPA(Realistic Failure Process Analysis)模型的基础上建立了脆性材料热破裂过程分析的数值模拟方法.该方法考虑了脆性材料在细观层次上力学性质的非均匀性(包括强度、弹模、传导系数等),并通过统计分布函数建立了宏、细观力学性能之间的联系.对不同均匀程度材料的数值模拟结果表明:材料的非均匀性对热传导规律、热应力分布以及热破坏模式有较大的影响.材料热力学性质的非均匀性加剧了材料内部热应力分布的非均匀性,这是致使非均匀材料热破裂的一个重要因素.对稳态和瞬态热传导两种条件下的脆性介质破裂过程模拟分析表明,考虑瞬态热传导计算所得到的破裂区小于相同条件下稳态热传导所得到的结果,表明在热破裂过程分析中,应注重考虑瞬态热传导对破裂过程的影响.     

半无限大体广义电磁热弹耦合的二维问题

基于Lord和Shulman广义热弹性理论,研究了热、电可导的半无限大体电磁热弹耦合的二维问题。半无限大体受热和外加恒定磁场的作用,文中建立了电磁热弹性耦合的控制方程,利用正则模态法求解得到了所考虑物理量的解析解,并用图形反映了各物理量的分布规律,从分布图上可以看出,介质中出现了电磁热弹耦合效应,各物理量的非零值仅在一个有限的区域内。     

功能梯度材料矩形中厚板的受压/热致屈曲

考虑材料组份沿板厚度方向按幂律变化的情形,研究了温度均匀变化时固支功能梯度材料(FGM)矩形中厚板的受压屈曲、热致屈曲和考虑热／机械预应力时的屈曲问题,给出了基于Reddy高阶剪切理论研究板屈曲荷载和屈曲临界温度的半解析数值方法．并以Si3N4／SUS304板为例考虑了材料组份、预加应力、横向剪切变形及面内位移约束条件等对FGM板屈曲承载能力的影响．     

压电层合板热屈曲问题的精确分析

本文在作者们以前工作的基础上导出了热压电介弹性稳定性问题的控制微分方程，提出了热压电层合板热屈曲问题的数学和学解法。由于最后所得特征值问题的复杂性，热屈曲临界必须数值求解。文中对不同压电材非压电材料层构成的层合板作了具体计算，考虑了压电性对屈曲温度的影响。     

ANALYTICAL SOLUTION OF THERMAL CONSOLIDATION OF SATURATED SOILS WITH TWO PARALLEL CYLINDRICAL HEAT SOURCES1）

基于应力平衡条件、渗流连续方程、能量守恒方程,考虑土颗粒和孔隙水热膨胀系数的不同,建立考虑热水力耦合的饱和土体三维热固结控制方程。利用傅里叶变换和拉普拉斯变换导出变换域上的控制方程,解得点热源在变换域和实数域上的解析解,再利用区域积分给出两平行圆柱形热源热固结土体温度、孔压、位移的解析解,并对其进行分析,发现径距比增大会导致两热源温度相互影响程度减弱,热固结系数减小会导致孔压和位移的峰值增大。     

两平行圆柱热源作用下饱和土体热固结解析解

基于应力平衡条件、渗流连续方程、能量守恒方程,考虑土颗粒和孔隙水热膨胀系数的不同,建立考虑热水力耦合的饱和土体三维热固结控制方程。利用傅里叶变换和拉普拉斯变换导出变换域上的控制方程,解得点热源在变换域和实数域上的解析解,再利用区域积分给出两平行圆柱形热源热固结土体温度、孔压、位移的解析解,并对其进行分析,发现径距比增大会导致两热源温度相互影响程度减弱,热固结系数减小会导致孔压和位移的峰值增大。     

一种着装人体动态热湿传递模拟方法

针对服装热湿功能设计问题,结合改进的人体热调节模型和微元织物热湿耦合模型,实现了着装人体动态热湿传递过程的模拟。水分的蒸发凝结、纤维对湿的吸附解吸、液态水的毛细传递、人体汗水在体表积聚等现象,及其对服装热性能的作用在模型中都给予了考虑。并给出了计算流程。算例表明了模拟方法的有效性。     

高强度钢板热成形数值模拟-静力显式

在已建立的高强度钢热成形热、力、相变耦合本构方程的基础上,基于虚功率方程及持续平衡方程建立了热成形静力显式多场耦合有限元列武。将热成形过程中的相变潜热引入温度场,并进行了有限元分析;在自主开发的商业化金属成形CAE软件KMAS(King-Mesh Analysis System)基础上,开发了考虑多场耦合的非线性、大变...     

考虑横法向热应变的C~0型整体-局部高阶层合梁理论

通过考虑横法向热变形,本文建立了预先满足层间应力连续的C0型整体-局部高阶层合梁理论,并用于分析复合材料层合梁热膨胀和热弯曲问题。虽然考虑了横法向应变,不增加额外的位移变量。此理论位移场不含有横向位移一阶导数,便于构造多节点高阶单元。基于虚功原理推导了复合材料层合梁平衡方程,并分析了简支多层复合材料梁热膨胀和热弯曲问题。数值结果表明,建立的模型能准确分析复合材料层合梁热膨胀和热弯曲问题,忽略横法向应变的理论分析热膨胀问题误差较大。     

石墨烯增强功能梯度材料板条的热弹性响应

研究了石墨烯增强功能梯度复合材料板条的热弹性问题。基于推广后的Mian和Spencer功能梯度板理论,在板厚度方向上考虑热传导引起的稳态温度场,沿板厚方向考虑了三种石墨烯纳米片的分布形式,最终给出了不同边界条件下石墨烯增强功能梯度板条在温度场作用下的解析解。通过算例分析验证了本文方法的有效性,并讨论了边界条件和石墨烯的分布形式等因素对板条热弹性响应的影响。计算结果表明:石墨烯增强功能梯度材料板条在温度场作用下的热响应与在机械力作用下的弹性响应有较大不同。     

考虑气动热对粘性系数μt影响的压气机流场计算

根据气动热影响叶片壁面附面层粘性系数tμ的理论,把考虑气动热对粘性系数tμ影响的粘性体积力计算方法运用于压气机叶栅三维粘性流场计算。计算时使用有限体积显式时间推进法解算Navier-Stokes方程,用Baldwin-Lomax模型模拟湍流流动。所得结果比不考虑气动热影响的三维计算结果更接近实验结果,能更真实地反映压气机叶栅流场的流动。