共查询到20条相似文献,搜索用时 78 毫秒
1.
2.
3.
基于精细积分思想,提出了一种有效的病态代数方程组求解方法。类似于稳态热传导方程可视为瞬态热传导方程的极限形式,将具有正定对称实系数矩阵的病态代数方程组归结为一个常微分方程组初值问题的极限形式,并在此基础上建立了病态代数方程组的精细积分解法。该方法不仅精度高,而且能以指数速度收敛,具有较高的效率。本文还讨论了病态代数方程... 相似文献
4.
采用双重互易边界元法结合精细积分法求解二维含热源的瞬态热传导问题。针对边界积分方程中热源项和温度关于时间导数项引起的域积分,采用双重互易法处理,将域积分转换为边界积分。采用边界元法将边界积分方程离散后,得到关于时间的微分方程组,并利用精细积分法处理其中的指数型矩阵;对于微分方程组中由边界条件和热源项引起的非齐次项,采用解析的方法计算。为了比较精细积分-双重互易边界元法的计算效果,同时使用有限差分法计算温度对时间的导数项。通过数值算例验证了本文方法的有效性和精确性。计算结果表明:时间步长对于精细积分-双重互易边界元法的结果影响较小,而有限差分法对时间步长比较敏感且只在时间步长选取较小时有效;当选取较大时间步长时,精细积分-双重互易边界元法依然具有良好的计算精度。 相似文献
5.
采用物理信息神经网络PINN(Physics-informed Neural Networks)求解稳态和瞬态功能梯度材料(FGMs)热传导问题。该方法利用控制方程、边界及初始条件的残差构造损失函数,在无任何响应数据的情况下得到了更具泛化能力的神经网络模型,同时避免了传统数值方法在求解计算力学问题时所需的微分、积分公式推导以及繁重的建模和划分网格等前处理工作。本文探究了PINN及其域分解的扩展物理信息神经网络XPINN(eXtended Physics-informed Neural Networks)在求解稳态和瞬态FGMs热传导问题时的适用性,讨论了网络结构对预测结果的影响。研究结果表明,PINN/XPINN在解决几何复杂的稳态和瞬态FGMs热传导问题时仍具有较高的可靠性和简洁的求解流程,同时,为极端环境下求解复杂多场耦合和夹杂等问题提供了新思路。 相似文献
6.
瞬态热传导方程精细积分方法中对称性的利用 总被引:3,自引:0,他引:3
采用精细积分法求解瞬态热传导方程时,对指数矩阵进行变换后使其具有对称性,利用这一特性可使存贮量和计算量降低一半。变换后指数矩阵的带宽特性不变,采用子域精细积分可进一步提高算法的计算与存储效率。 相似文献
7.
热传导问题灵敏度分析的伴随法 总被引:5,自引:1,他引:5
在热传导灵敏度分析的直接法的研究基础上,进一步探讨了稳态和瞬态热传导问题灵敏度分析的伴随法.推导了伴随法的计算列式,对于瞬态热传导问题,研究了瞬态约束处理的关键点方法,并提出伴随方程的精细积分解法。算例表明,稳态问题灵敏度计算,伴随法与直接法的结果是一致的;瞬态问题灵敏度计算,两种方法的精度相当。 相似文献
8.
9.
10.
提出一种基于奇异边界法结合双重互易法的数值模型来求解瞬态热传导问题。奇异边界法属于配点型边界无网格方法,相对于网格方法,其具有无需划分网格,只需边界配点的优势。运用差分格式来处理热传导方程中的时间变量,将原热传导方程化为非齐次修正Helmholtz方程。修正Helmholtz方程的解由齐次解和特解两部分组成,齐次解通过奇异边界法求出,特解由双重互易法求出,源项由径向基函数近似。通过数值算例检验了本文数值模型的精度及有效性;算例结果表明,该数值模型计算精度较高,误差基本都在1%以内,具有很好的稳定性,能有效地应用于求解多连通域的瞬态热传导问题。 相似文献
11.
获得热传导问题“拟解析解”的精细积分算法 总被引:8,自引:1,他引:8
研究并讨论了用于热传导分析有限元解的精细积分算法,算法很好地克服了传统方法求解时的单调性问题,且对空间离散后所获得方程的解是解析的,因而算法的解将具有“拟解析解”的意义,论文证明了算法单调性. 相似文献
12.
对线性定常结构动力系统提出的精细积分方法,在数值精度等方面表现出极大优越性,但是当矩阵尺度很大时在数值计算与存储中将产生困难,对此,本文对瞬态热传导方程,根据结构的概念,将结构分为若干个子结构,对各子结构分别进行指数矩阵运算并通过了结构间界面的物理量相联系,从而提高精细积分方法的计算效率。 相似文献
13.
基于精细积分技术的非线性动力学方程的同伦摄动法 总被引:2,自引:0,他引:2
将精细积分技术(PIM)和同伦摄动方法(HPM)相结合,给出了一种求解非线性动力学方程的新的渐近数值方法。采用精细积分法求解非线性问题时,需要将非线性项对时间参数按Taylor级数展开,在展开项少时,计算精度对时间步长敏感;随着展开项的增加,计算格式会变得越来越复杂。采用同伦摄动法,则具有相对筒单的计算格式,但计算精度较差,应用范围也限于低维非线性微分方程。将这两种方法相结合得到的新的渐近数值方法则同时具备了两者的优点,既使同伦摄动方法的应用范围推广到高维非线性动力学方程的求解,又使精细积分方法在求解非线性问题时具有较简单的计算格式。数值算例表明,该方法具有较高的数值精度和计算效率。 相似文献
14.
提出了一种基于特征正交分解(POD)和有限元法的瞬态非线性热传导问题的模型降阶快速分析方法, 建立了导热系数随温度变化的一类瞬态非线性热传导问题有限元格式的POD降阶模型. 在隐式时间推进方法的基础上有效结合单元预转换方法和多级线性化方法发展了一种加速求解瞬态非线性热传导降阶模型的新型计算方法,并通过二维和三维算例验证了该方法的准确性和高效性. 研究结果表明: (1)降阶模型解的均方根误差在经过初始时段轻微的脉动后稳定于0.01%以下, 而其计算效率比有限元全阶模型提高2$\sim $3个数量级, 并且自由度数量(DOFs)愈大提高的幅度也愈加显著; (2)新型算法解决了常规算法在计算非线性降阶模型时加速性能差的问题, 即使是在DOFs比较小的时候也能够明显提高计算效率; (3)常数边界条件下得到的POD模态可以用来建立相同求解域在各种复杂时变边界条件下的瞬态非线性热传导降阶模型, 并对其传热过程和温度场进行快速准确的分析与预测, 具有很好的工程应用价值. 相似文献
15.
提出了一种基于特征正交分解(POD)和有限元法的瞬态非线性热传导问题的模型降阶快速分析方法,建立了导热系数随温度变化的一类瞬态非线性热传导问题有限元格式的POD降阶模型.在隐式时间推进方法的基础上有效结合单元预转换方法和多级线性化方法发展了一种加速求解瞬态非线性热传导降阶模型的新型计算方法,并通过二维和三维算例验证了该方法的准确性和高效性.研究结果表明:(1)降阶模型解的均方根误差在经过初始时段轻微的脉动后稳定于0.01%以下,而其计算效率比有限元全阶模型提高2~3个数量级,并且自由度数量(DOFs)愈大提高的幅度也愈加显著;(2)新型算法解决了常规算法在计算非线性降阶模型时加速性能差的问题,即使是在DOFs比较小的时候也能够明显提高计算效率;(3)常数边界条件下得到的POD模态可以用来建立相同求解域在各种复杂时变边界条件下的瞬态非线性热传导降阶模型,并对其传热过程和温度场进行快速准确的分析与预测,具有很好的工程应用价值. 相似文献
16.
针对u-p形式的饱和两相介质波动方程,采用精细时程积分方法计算固相位移u,采用向后差分算法求解流体压力p,建立了饱和两相介质动力固结问题时域求解的精细时程积分方法。针对标准算例,对该方法的计算精度进行了校核。开展了该方法相关算法特性的研究,对采用不同数值积分方法计算非齐次波动方程特解项计算精度的差异进行了对比研究,并对采用不同积分点数目的高斯积分法计算特解项条件下计算精度的差异进行了对比研究。研究结果表明,(1)该方法具有良好的计算精度。(2)计算非齐次波动方程特解项的数值积分方法中,梯形积分法的计算精度最差,高斯积分法、辛普生积分法和科茨积分法都具有较好的计算精度。(3)增加高斯积分点数目对于提高计算精度的作用并不显著。 相似文献
17.
18.
瞬态热传导问题普遍存在于航空航天、土木和冶金等领域中, 对这类问题精确、高效的数值求解方法一直备受关注. 为此, 本文提出了一种无条件稳定的单步时间积分方法. 在所提出的方法中, 拉格朗日插值函数被用来近似真实的温度场及其一次导数场, 之后, 加权残量法被用来建立二者之间的关系. 通过对算法参数的巧妙设计, 本文提出的方法具有二阶精度和L型数值耗散, 其中, L型耗散使得本文方法能够快速过滤掉虚假的高频振荡. 多数现有时间积分方法仅对线性系统具有无条件稳定性, 对非线性系统则是条件稳定的. 为此, 本文改进了Hughes针对一阶非线性热传导问题提出的时间积分方法稳定性评估理论, 并将改进的理论用于方法的参数设计中. 理论分析的结果表明本文方法对线性和非线性热传导系统均是无条件稳定的. 即使对于著名的Crank-Nicolson方法失稳的非线性热传导问题, 本文方法仍能给出稳定且精确的预测. 数值测试结果显示, 所提出的方法相较于当前流行的方法具有明显的精度、耗散和稳定性优势. 相似文献
19.
相较于传统有限元,数值流形方法(numerical manifold method, NMM) 的一个显著优点是在处理裂纹问题时网格无需与裂纹重合,这就方便了岩体破坏过程的模拟. 基于包含裂尖增强函数的NMM,采用Newmark 隐式动力学算法进行时间积分,重点研究了动力载荷条件下裂纹动态应力强度因子(dynamic stress intensity factor,DSIF) 的求解方法. 针对典型的线弹性动力裂纹问题,给出了NMM 的数值算例. 结果表明NMM 能够准确计算动载荷条件下裂纹的DSIF,并且具有较好的收敛性. 相似文献
20.
相较于传统有限元,数值流形方法(numerical manifold method, NMM) 的一个显著优点是在处理裂纹问题时网格无需与裂纹重合,这就方便了岩体破坏过程的模拟. 基于包含裂尖增强函数的NMM,采用Newmark 隐式动力学算法进行时间积分,重点研究了动力载荷条件下裂纹动态应力强度因子(dynamic stress intensity factor,DSIF) 的求解方法. 针对典型的线弹性动力裂纹问题,给出了NMM 的数值算例. 结果表明NMM 能够准确计算动载荷条件下裂纹的DSIF,并且具有较好的收敛性. 相似文献