一种含多裂纹损伤曲梁自由振动的自适应求解方法技术

本申请公开了一种含多裂纹损伤曲梁自由振动的自适应求解方法，包括：建立含裂纹损伤的Timoshenko曲梁坐标系，建立曲梁中裂纹的截面惯性矩模型和截面面积模型；建立曲梁的自由振动控制方程，基于给定的有限元网格，得到网格下的有限元解；通过超收敛拼片恢复法，计算得到振型超收敛解，再通过Rayleigh商计算得到频率超收敛解；对振型有限元解进行误差估计处理，得到相对误差值，输出目标值，或重复上述步骤优化处理。本申请可以得到更加精确的振型解和频率解；适用于不同边界条件、不同裂纹数量、不同裂纹分布情况下的含裂纹曲梁自由振动问题。题。题。

【技术实现步骤摘要】
一种含多裂纹损伤曲梁自由振动的自适应求解方法


[0001]本申请属于曲梁损伤评价领域，具体涉及一种含多裂纹损伤曲梁自由振动的自适应求解方法。

技术介绍

[0002]Timoshenko梁在土木、机械、船舶、航空航天等工程中具有普遍的应用，研究曲线结构的抗震与直线结构不同，它的力学特性复杂，分析起来更加困难。裂纹是梁结构损伤的主要形式之一，这在很大程度上影响结构的可靠性，裂纹尖端的应力集中现象会使裂纹继续扩张并导致结构失效。除此之外，在裂纹到达临界值前对其进行准确识别并进行维护、替换，可以在极大程度上避免巨大安全隐患。因此，对含裂纹Timoshenko梁结构的自由振动特性进行研究具有很强的实际意义。目前，对于复杂结构形式和边界条件的中厚圆柱壳自由振动分析多采用有限元法进行求解。
[0003]常规有限单元法使用的误差估计一般通过事先提供的上下限来分析计算误差，但是理论分析的困难使得这种方法无法满足设定精度的误差分析，且计算的误差范围与实际误差限通常相差较大，生成网格的技术也较为复杂。目前的有限元方法往往要借助经验的判断和定性的分析，精心地设计网格尺寸和阶次，如需要调整网格划分，需要重新计算分析网格参数，较为繁杂。即使这样，对于有限元结果，也只能得到一个相对合理但精度难以满足的解答。

技术实现思路

[0004]本申请提出了一种含多裂纹损伤曲梁自由振动的自适应求解方法，基于裂纹损伤惯性矩、面积，建立自由振动控制方程，采用常规有限元法求得当前网格下的有限元解，进一步使用超收敛拼片恢复法得到超收敛解，通过能量模形下的误差估计，不断提高求解精度，最终得到满足误差要求的解答。
[0005]为实现上述目的，本申请提供了如下方案：
[0006]一种含多裂纹损伤曲梁自由振动的自适应求解方法，包括如下步骤：
[0007]S1.建立含裂纹损伤的曲梁坐标系，在所述曲梁坐标系下，建立曲梁中裂纹的截面惯性矩模型和截面面积模型，所述截面惯性矩模型和所述截面面积模型均采用裂纹深度表征；
[0008]S2.基于所述截面惯性矩模型和所述截面面积模型，建立所述曲梁的自由振动控制方程，基于所述自由振动控制方程和给定的有限元网格，得到网格下的有限元解，所述网格下的有限元解包括振型有限元解；
[0009]S3.基于所述振型有限元解，通过超收敛拼片恢复法，得到振型超收敛解，基于所述振型超收敛解，通过Rayleigh商计算得到频率超收敛解；
[0010]S4.基于所述振型超收敛解，对所述振型有限元解进行误差估计处理，得到所述振型有限元解的相对误差值，当所述相对误差值大于预设误差值时，调整所述自由振动控制
方程，重复所述S2
‑
S4，当所述相对误差值小于等于所述预设误差值时，所述振型超收敛解即为目标振型解，所述频率超收敛解即为目标频率解，完成所述曲梁的自由振动求解。
[0011]优选的，所述S1中，建立所述曲梁的截面惯性矩模型和所述截面面积模型的方法包括：
[0012]基于所述曲梁坐标系，采用裂纹截面损伤缺陷比拟法，建立裂纹截面损伤模型，所述裂纹截面损伤模型用于表征裂纹深度；
[0013]基于所述裂纹截面损伤模型，建立所述截面惯性矩模型和所述截面面积模型。
[0014]优选的，所述裂纹截面损伤模型为
[0015][0016]α为表征裂纹损伤深度的截面损伤率，α＝0表示梁截面完整无损伤，h
c
为裂纹绝对深度，h表示曲梁高度。
[0017]优选的，所述截面惯性矩模型为：
[0018][0019]I
c
表示考虑裂纹损伤的横截面惯性矩，b表示曲梁厚度；
[0020]所述截面面积模型为：
[0021]A
c
＝bh(1
‑
α)
[0022]A
c
表示考虑裂纹损伤的截面面积。
[0023]优选的，得到所述网格下的有限元解的方法包括：
[0024]基于所述自由振动控制方程，建立特征值方程，所述特征值方程采用频率和振型函数向量表示；
[0025]基于所述特征值方程和给定的有限元网格，建立线性矩阵特征值方程；
[0026]基于所述线性矩阵特征值方程，采用逆幂迭代法得到所述网格下的有限元解。
[0027]优选的，所述振型超收敛解w
*
为：
[0028]w
*
(x)＝Pa
[0029]其中，P为给定函数向量；a为待定系数向量。
[0030]优选的，得到所述振型有限元解的相对误差值的方法包括：基于所述振型超收敛解，对所述振型有限元解进行能量模形式下的误差估计处理，得到所述相对误差值。
[0031]优选的，当所述相对误差值大于等于所述预设误差值时，采用单元均匀细分加密的h型网格自适应方式调整所述自由振动控制方程。本申请的有益效果为：
[0032]本申请公开了一种含多裂纹损伤曲梁自由振动的自适应求解方法，在常规有限元解的基础上，引入超收敛拼片恢复法，得到了比常规有限元解更加精确的振型解和频率解；通过建立误差估计分析，评判求解精度，可以不断提高计算精度，进而得到符合精度要求的求解过程。本申请适用于不同边界条件、不同裂纹数量、不同裂纹分布情况下的含裂纹Timoshenko曲梁自由振动问题，具有广泛的适用性。
附图说明
[0033]为了更清楚地说明本申请的技术方案，下面对实施例中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本申请的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。
[0034]图1为本申请实施例一种含多裂纹损伤曲梁自由振动的自适应求解方法的流程示意图；
[0035]图2为本申请实施例中建立的曲梁坐标系示意图。
具体实施方式
[0036]下面将结合本申请实施例中的附图，对本申请实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本申请一部分实施例，而不是全部的实施例。基于本申请中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本申请保护的范围。
[0037]为使本申请的上述目的、特征和优点能够更加明显易懂，下面结合附图和具体实施方式对本申请作进一步详细的说明。
[0038]如图1所示，为本申请实施例一种含多裂纹损伤曲梁自由振动的自适应求解方法的流程示意图。
[0039]S1.在曲梁坐标系下，建立曲梁中裂纹的截面惯性矩模型和截面面积模型。
[0040]如图2所述，为本实施例建立的曲梁坐标系，其中，曲梁中性轴坐标为s，坐标系为xyz，其中x、y为曲梁平面内坐标，x沿轴线切向，y沿轴线法向，z垂直于轴线所在平面。面内振动的位移为：沿x轴位移振幅u、沿y轴位移振幅v和绕z轴的转角振幅Ψ
z
。记曲梁曲率半径为R(s)，截面剪切刚度修正系数本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种含多裂纹损伤曲梁自由振动的自适应求解方法，其特征在于，包括如下步骤：S1.建立含裂纹损伤的曲梁坐标系，在所述曲梁坐标系下，建立曲梁中裂纹的截面惯性矩模型和截面面积模型，所述截面惯性矩模型和所述截面面积模型均采用裂纹深度表征；S2.基于所述截面惯性矩模型和所述截面面积模型，建立所述曲梁的自由振动控制方程，基于所述自由振动控制方程和给定的有限元网格，得到网格下的有限元解，所述网格下的有限元解包括振型有限元解；S3.基于所述振型有限元解，通过超收敛拼片恢复法，得到振型超收敛解，基于所述振型超收敛解，通过Rayleigh商计算得到频率超收敛解；S4.基于所述振型超收敛解，对所述振型有限元解进行误差估计处理，得到所述振型有限元解的相对误差值，当所述相对误差值大于预设误差值时，调整所述自由振动控制方程，重复所述S2
‑
S4，当所述相对误差值小于等于所述预设误差值时，所述振型超收敛解即为目标振型解，所述频率超收敛解即为目标频率解，完成所述曲梁的自由振动求解。2.根据权利要求1所述的含多裂纹损伤曲梁自由振动的自适应求解方法，其特征在于，所述S1中，建立所述曲梁的截面惯性矩模型和所述截面面积模型的方法包括：基于所述曲梁坐标系，采用裂纹截面损伤缺陷比拟法，建立裂纹截面损伤模型，所述裂纹截面损伤模型用于表征裂纹深度；基于所述裂纹截面损伤模型，建立所述截面惯性矩模型和所述截面面积模型。3.根据权利要求2所述的含多裂纹损伤曲梁自由振动的自适应求解方法，其特征在于，所述裂纹截面损伤模型为α为表征裂纹损伤深度的截面损伤率，α＝0表示梁截面完整无损伤...

【专利技术属性】
技术研发人员：王永亮，王建辉，
申请(专利权)人：王永亮，
类型：发明
国别省市：