半空间饱和土内置点载荷作用下的热弹性波动

基于Biot波动理论及热弹性动力理论，利用已建立的饱和多孔弹性介质热流固耦合控制方 程，研究半无限地基在内置点热力源作用下的动力响应问题. 求解过程引用Hankel变换技 术，得到了热力源作用下土体中温度增量、应力、位移和孔隙水压力的积分形式解答. 利用Hankel数值逆变换得到计算结果，分析了热流固耦合条件下激振频率对竖向位移和孔 隙水压力响应的影响. 对热流固耦合、热弹性和多孔弹性模型计算结果进行了比较.     

含损伤演化的TM耦合数值模型及其应用研究

从岩石材料的细观结构层次出发，应用损伤力学和热弹性理论，对热力耦合作用下岩石破裂 过程中热－应力相互作用关系进行了分析. 初步建立了细观岩石热－力(TM)耦合数值模型， 探讨了TM耦合作用下岩石材料的细观结构损伤及其诱发的材料力学性能演化机制，把岩石 热固耦合问题的研究从应力状态分析深入到损伤、破坏过程分析之中. 运用该数值模型对某 硬岩实验室开展的原位尺度实验中的废料处理井间柱稳定性进行了模拟分析，其应力场、岩 石剥离破坏形态及破坏诱发的AE特性等均与实验监测结果表现出了较好的一致性，证 明了该数值模型的合理性和有效性.     

热-机载荷下不规则夹杂的奇异性应力场分析

基于有限元特征分析法得到的夹杂角部场数值特征解开发了一种超级奇异单元模型,并将其与普通四节点单元紧密结合,用于热-机载荷下夹杂角端部的应力场分析。在数值计算中,考察了热-机载荷下不同弹性比和不同夹杂尺寸的应力强度因子,并将所得结果与文献解和传统有限元方法解比对。结果表明,本文方法对热-机耦合条件下的不规则夹杂角端部的热弹性应力分析极为有效,可避免局部网格的高度加密,并提高计算效率。模型在复合材料夹杂的局部强度问题分析方面具有很好的实用性。     

弹塑性变形油藏中多相渗流的数值模拟

基于流固耦合力学理论,建立了弹性变形油藏中多相渗流的数学模型,该模型考虑了渗流与变形的耦合作用,以及注采交变载荷作下油藏多孔介质的弹性变形特征,给出了耦合数值模拟方法和算例。     

热结构瞬态响应的耦合灵敏度分析方法与优化设计

研究结构瞬态热变形和热应力的灵敏度分析方法及其优化设计,灵敏度计算给出了直接法和伴随法两种算法.考虑了温度场的耦合作用,在直接法中需要计算温度场对设计变量的导数,在伴随法中需要计算热载荷对温度场的导数.数值算例验证了该方法的精度.伴随法在应用程序中的实现,为大型结构优化提供了高效率的灵敏度计算方法.     

基于光滑有限元法的热-弹相场断裂研究

基于相场理论，论文构建了热力耦合作用下热-弹性材料的复杂断裂控制方程，采用光滑有限元法对该模型进行了数值离散；设计编写了相应的Matlab计算程序；开展了热-力载荷下的二维拉伸断裂模拟，将结果和有限元计算结果对比，吻合一致.同时采用该模型开展热冲击下陶瓷板的多裂纹扩展过程模拟，计算结果和实验结果吻合良好，验证了论文方法的有效性.     

磁-电-弹性半空间在轴对称热载荷作用下的三维问题研究

考虑力-电-磁-热等多场耦合作用, 基于线性理论给出了磁-电-弹性半空间在表面轴对称温度载荷作用下的热-磁-电-弹性分析, 并得到了问题的解析解. 利用Hankel 积分变换法求解了磁-电-弹性材料中的热传导及控制方程, 讨论了在磁-电-弹性半空间在边界表面上作用局部热载荷时的混合边值问题, 利用积分变换和积分方程技术, 通过在边界表面上施加应力自由及磁-电开路条件, 推导得到了磁-电-弹性半空间中位移、电势及磁势的积分形式的表达式. 获得了磁-电-弹性半空间中温度场的解析表达式并且给出了应力, 电位移和磁通量的解析解. 数值计算结果表明温度载荷对磁-电-弹性场的分布有显著影响. 当温度载荷作用的圆域半径增大时, 最大正应力发生的位置会远离半无限大体的边界; 反之当温度载荷作用的圆域半径减小时, 最大应力发生的位置会靠近半无限大体的边界. 电场和磁场在温度载荷作用的圆域内在边界表面附近有明显的强化, 而磁-电-弹性场强化区域的强化程度跟温度载荷的大小和作用区域大小相关. 本研究的相关结果对智能材料和结构在热载荷作用下的设计和制造具有指导意义.      

THREE-DIMENSIONAL ANALYSIS OF A MAGNETOELECTROELASTIC HALF-SPACE UNDER AXISYMMETRIC TEMPERATURE LOADING 1)

考虑力-电-磁-热等多场耦合作用, 基于线性理论给出了磁-电-弹性半空间在表面轴对称温度载荷作用下的热-磁-电-弹性分析, 并得到了问题的解析解. 利用Hankel 积分变换法求解了磁-电-弹性材料中的热传导及控制方程, 讨论了在磁-电-弹性半空间在边界表面上作用局部热载荷时的混合边值问题, 利用积分变换和积分方程技术, 通过在边界表面上施加应力自由及磁-电开路条件, 推导得到了磁-电-弹性半空间中位移、电势及磁势的积分形式的表达式. 获得了磁-电-弹性半空间中温度场的解析表达式并且给出了应力, 电位移和磁通量的解析解. 数值计算结果表明温度载荷对磁-电-弹性场的分布有显著影响. 当温度载荷作用的圆域半径增大时, 最大正应力发生的位置会远离半无限大体的边界; 反之当温度载荷作用的圆域半径减小时, 最大应力发生的位置会靠近半无限大体的边界. 电场和磁场在温度载荷作用的圆域内在边界表面附近有明显的强化, 而磁-电-弹性场强化区域的强化程度跟温度载荷的大小和作用区域大小相关. 本研究的相关结果对智能材料和结构在热载荷作用下的设计和制造具有指导意义.     

基于热固耦合问题的多尺度拓扑优化设计

基于均质材料的拓扑优化逐渐难以适应现代化生产对工业产品高品质、轻量化的需求,同时很多工业设备经常面临着高温和高负载的工作环境.为了提高结构在温度载荷和机械载荷共同作用下的力学性能,本文提出了一种在稳态热源作用下的热固耦合连续体结构并行化拓扑优化方法.以结构刚度作为目标函数,材料的体分比为约束,利用能量均匀化方法预测微结构的等效属性,建立热固耦合结构并行化拓扑优化数学模型.为便于数值计算,利用直接法进行灵敏度分析,同时采用Heaviside非线性密度过滤技术抑制数值不稳定性,将OC准则算法用于求解优化问题.数值算例表明,本文方法能够有效地进行优化设计且能显著地提高结构的刚度性能和散热性能.     

载流圆形薄板的磁弹性应力与变形分析

本文在建立非定常电磁场和机械载荷作用下载流圆形薄板非线性变形状态下的磁弹性二维关系方程和运动方程的基础上，对其在非静态电磁场和机械载荷联合作用下的位移及应力进行了分析，获得了耦合场中一些参量的变化规律，通过具体算例，得出了载流圆形薄板在机械场和电磁场作用下的位移及应力与通电电流强度之间的关系，并解决了圆板中心处的奇异性问题，给出了在轴对称条件下的数值解，计算结果表明，改变通电电流及电磁感应强度的大小，可以改变载流薄板的应力和变形状态，达到控制薄板的受力和变形的目的。     

含热冲击预损伤的陶瓷基复合材料损伤本构模型

陶瓷基复合材料结构在服役过程中不可避免地经受热冲击(较高的热应力梯度)而产生热机械损伤, 因此, 建立含循环热冲击预损伤材料的损伤本构模型, 以描述材料在热机械载荷作用下的力学行为, 对材料结构损伤容限设计与结构完整性评估非常重要. 本文首先对经历了循环热冲击的材料进行单调拉伸损伤实验, 发现对于含循环热冲击预损伤的材料, 其弹性模量的下降与所施加的应变直接相关. 然后在连续介质损伤力学的框架下, 基于平面应力假设, 建立了含循环热冲击预损伤材料的损伤演化模型, 该模型所涉及的参数可通过一个偏轴(45$^\circ$)以及两个正轴(平行于两个主方向)的单调拉伸试验获得. 最后, 采用经典塑性理论对由基体损伤引起的非弹性应变进行了描述. 本文所提出的应变损伤宏观模型可以描述陶瓷基复合材料在热机械载荷作用下的损伤演化, 同时弥补了含预损伤的陶瓷基复合材料在机械载荷下损伤本构模型在理论及实验研究方面的不足.      

基于概率的平面连续体结构拓扑优化

构建了随机载荷作用下具有随机参数的平面连续体结构拓扑优化设计模型,以结构总质量均值极小化为目标函数,以满足单元应力可靠性为约束条件;对应力可靠性约束进行了等价显示化处理;导出了随机参数结构的应力反应数字特征的计算表达式;提出了相应的求解策略和方法。算例的优化结果表明文中模型的合理性和方法的有效性。     

多工况下结构拓扑优化设计

考虑单元删除和增加对结构应力约束的影响，提出了一种新的多工况下结构的双方向渐进优化方法．首先，基于在结构孔洞或边界周围附加人工材料的思路，建立了结构优化模型和应力灵敏度公式．然后，结合结构应力和应力灵敏度，给出了多工况静力载荷下考虑静应力约束的结构优化准则和算法．开展了结构仿真设计，结果表明给出的方法是正确和有效的，并具有较好的工程应用价值．     

单晶有限变形的热力耦合计算模型

以弹性变形梯度作为基本变量,结合热力学理论构造了单晶有限变形的热、力耦合计算模型。该模型考虑了温度、变温速率以及塑性耗散等条件对单晶有限变形的影响,相对于传统的以弹性变形梯度为基本变量的晶体塑性模型,算法能够体现温度效应的影响。采用隐式的积分方法对建立的控制方程进行计算以保证求解过程的稳定。以1100Al单晶为例计算了不同升温、降温速率,以及不同应变率影响下的材料应力-应变的响应。结果表明,模型能较好地反映变温过程中,单晶各向异性性质的演化以及应力、应变之间关系的变化。     

A Revisit to the Frozen Stress Phenomena in Photoelasticity

Stress freezing in photoelasticity is regularly employed for 3-D elastic stress analysis of geometrically complex models, wherein the models under load are soaked at the stress freezing temperature and slowly cooled to room temperature. During this process, a model may distort or fail undergoing large deformations at the stress freezing temperature owing to thermo-mechanical interactions. The contribution of thermal deformation to the model distortions is neglected in the available literature of stress freezing. This aspect of stress freezing is investigated in this paper, wherein Digital Image Correlation (DIC) is used to study the strain behavior during a complete stress freezing cycle for an epoxy made of CY230 resin cured with HY951 hardener. The results show that the thermal contributions to the model distortions at the critical temperature must be taken care of to estimate the failure margins. The distortions and failures would mainly depend upon the thermal and mechanical response of the model material and the complexity of the model.     

基于有限变形单晶塑性滑移的微观热-力耦合模型及其有限元计算

基于率相关的晶体塑性滑移理论,论文考虑晶体内部塑性变形产生的热以及快速热冲击作用下温度急剧变化产生热应力的热-力双向耦合效应,建立起微观单晶的瞬态热-弹-塑性耦合模型,推导出与温度有关的剪应变率和弹塑性切模量公式.根据论文建立的模型,对ABAQUS软件进行二次开发[1],数值模拟出<001>/{100}单晶Cu在单轴拉伸状态下的应力、应变与温度的关系和弹性模量的变化,结果如下:轴向应力随温度升高先呈线性增加再呈非线性减小,轴向应变随温度增加而增加;弹性模量随塑性变形的增加而降低,与分子动力学模拟的趋势[2]是一致的.数值实验表明,论文建立的模型和算法是正确合理的,且计算量远远小于分子动力学模拟.     

正交各向异性材料热弹性问题的三维无网格法计算模型及应用

建立了正交各向异性材料热弹性问题的三维无网格伽辽金(Element Free Galerkin, EFG)法计算模型。利用该计算模型对三维复合材料汽轮机叶轮和轴承座进行了热弹性分析,对比了材料方向角及热导率因子、热膨胀系数因子和拉压弹性模量因子不同组合情况下轴承座的最大热变形总位移和当量应力值,讨论了材料方向角及上述正交各向异性因子对热变形和当量应力的影响规律,并与各向同性材料进行了对比。结果表明:三维EFG模型的热变形总位移和当量应力相对误差范数分别比有限元法小0.1215%和0.1359%;材料方向角同时影响热变形的大小和方向,但对当量应力方向影响不大;正交各向异性材料因子主要影响热变形和当量应力的大小。在考虑热-机械载荷作用下的三维复合材料零件结构设计中,当以刚度或强度为主要需求时,材料方向角、热导率因子、热膨胀系数因子、拉压弹性模量因子分别在(45°~60°,8:1:4~10:1:5,(1/6):(1/5):1~(1/5):(1/4):1,(7/5):1:(9/5)~(3/2):1:2)或(0°~10°,(1/10):1:(1/5)~(1/8):1:(1/4),(1/5):1:(1/6)~(1/4):1:(1/5),1:(1/5):(1/10)~1:(1/4):(1/8))范围内取值能有效降低轴承座等结构的热变形和当量应力。     

含损伤演化的TM耦合数值模型及其应用研究

基于生物力学中的Wolff法则，发展了一种连续体拓扑优化的新方法. 有待优化的结构 被看作是一块遵从Wolff法则生长的骨骼，把寻找其最优拓扑的过程比拟为骨骼的重建/生 长过程. 采用骨骼的重建/生长规律作为准则更新材料分布，直至达到一个平衡状态并由此 获得结构的最优拓扑. 算例表明了所提出方法的有效性.