【技术实现步骤摘要】
本专利技术涉及材料疲劳损伤分析技术，具体涉及一种基于时间双尺度分解的疲劳损伤演化分析方法。
技术介绍
1、对于工程中大量使用的混凝土以及岩石等准脆性材料，在受到反复荷载作用时易产生疲劳损伤，严重影响其结构安全性和可靠度。传统的疲劳分析方法多基于经验模型，通常在唯象学的框架下进行，给出的损伤变量的经验演化表达式缺乏理论背景和逻辑完善性，无法准确描述材料的微观损伤机制，可能导致实际应用中对材料疲劳损伤的预测结果存在较大误差，并且分析过程采用逐步有限元分析，存在分析耗时长的问题。
技术实现思路
1、本专利技术所要解决的技术问题是：提出一种基于时间双尺度分解的疲劳损伤演化分析方法，解决传统的疲劳分析方法不能准确反映疲劳过程中损伤演化的发展、分析耗时长的问题，为准脆性材料的疲劳寿命预测和损伤评估提供科学依据。
2、本专利技术解决上述技术问题采用的技术方案是：
3、一种基于时间双尺度分解的疲劳损伤演化分析方法，包括以下步骤：
4、s1、通过对准脆性材料微观层次微裂纹扩展过程的分析，确定纳米尺度的微裂纹的扩展速率；
5、s2、根据纳米尺度的微裂纹的扩展速率，结合损伤累积效应，确定疲劳过程能量耗散表达式；
6、s3、根据疲劳过程能量耗散表达式，结合细观随机断裂模型，建立准脆性材料疲劳损伤演化方程；
7、s4、引入时域摄动参数和时间平均算符，通过对准脆性材料损伤演化方程中的状态变量进行摄动展开和时间平均，分别建立微观时间尺度的控
8、s5、对所述微观时间尺度的控制方程和宏观时间尺度的控制方程进行交错迭代求解，直至达到准脆性材料完全损伤破坏的判定条件，获得该准脆性材料随时间变化的疲劳损伤分布结果。
9、进一步的，步骤s1中，确定的纳米尺度的微裂纹的扩展速率的公式表达为：
10、
11、其中，为两个相邻亚稳态之间的平均迁移速率；为材料的纳米微颗粒间距； e为自然常数；为特征频率，，其中 k、 t和 h分别是boltzmann常数、绝对温度和plank常数；
12、q0为能量势垒；为微裂纹扩展一个纳米微颗粒间距时释放的能量均值；
13、为微裂纹扩展时宏观损伤的微小增加距离； v为损伤影响的体积； y为损伤能释放率；
14、为均匀化表面能，，其中，为考虑疲劳扰动影响下的表面能；
15、为裂纹尖端张开位移的变化的度量，定义为在一定周期内微裂纹尖端所经历的平均应变幅：
16、
17、其中， t为当前时间；t0为疲劳加载周期；为应变率。
18、进一步的，步骤s2中，确定的疲劳过程能量耗散表达式为：
19、
20、其中， n为微裂纹的总个数，为材料体积中所有微裂纹扩展时的表面能， p为尺度常数； η为几何常数，为材料体积内微裂纹的平均长度；
21、为与温度有关的常数，。
22、进一步的，步骤s3中，建立的准脆性材料疲劳损伤演化方程表示为：
23、
24、其中，d为材料损伤变量；h(.)为heaviside函数；为疲劳累积能量耗散；为特征能量； x为所考虑材料空间坐标；为断裂应变随机场的空间微分；为对时间的微分； ψ为材料常数，宏观损伤对微观裂纹群的抑制；为弹性模量；为断裂应变随机场； β为一阶常系数。
25、进一步的，步骤s4中，将宏观时间尺度记为 t，将微观时间尺度记为 τ，；
26、其中，为预设的远小于1的尺度参数；
27、建立的宏观时间尺度的控制方程包括阶的宏观时间尺度上的平衡方程和本构方程；
28、其中，阶的宏观时间尺度上的平衡方程为：
29、
30、阶的宏观时间尺度上的本构方程为：
31、
32、约束边界条件为：
33、
34、其中，各类变量右下标“+”和“-”分别代表受拉和受压分量，右上标“0”和“1”分别代表对应摄动展开后对应阶和阶分量；
35、(x),t表示变量x对宏观时间t的偏导；为拉普拉斯算子；
36、为宏观尺度下的有效应力张量；为宏观尺度下的位移矢量；为时间平均算符，，为微观时间尺度上的疲劳加载周期；
37、和为阶宏观损伤变量；和为阶宏观有效应力张量；和为阶微观损伤变量；为材料所在空间区域；为外部力作用边界；为外部位移作用边界；n为所考虑边界的法向向量；f和分别为面力和预定位移；in表示作用于空间区域内；on表示作用于空间区域的边界上；
38、建立的微观时间尺度的控制方程包括阶的微观时间尺度上的平衡方程和本构方程；
39、其中，阶的微观时间尺度上的平衡方程为：
40、
41、 阶的微观时间尺度上的本构方程为：
42、
43、约束边界条件为：
44、
45、其中，为微观尺度下的有效应力张量；为微观尺度下的位移矢量；和为阶宏观有效应力，式中右下标“+”和“-”分别代表受拉和受压分量；为预定位移，u表示位移，f表示面力；( x),τ表示变量 x对微观时间 τ的偏导。
46、进一步的，步骤s5中，所述交错迭代求解具体为：
47、先求解微观时间尺度的控制方程，将求解结果用于宏观时间尺度的控制方程的求解，将宏观时间尺度的控制方程的求解结果再反馈到微观时间尺度的控制方程的求解中，以此类推进行循环求解。
48、进一步的，步骤s5中，对所述微观时间尺度的控制方程和宏观时间尺度的控制方程进行交错迭代求解，包括：
49、s51、微观时间尺度的控制方程求解：
50、已知条件为： τ时刻微观时间尺度应力和有效应力张量和，宏观时间尺度有效应力张量和损伤变量，有限元软件自动提供微观时间尺度应变增量；
51、计算更新：计算有效应力，并计算相应的损伤能释放率，再根据式
52、
53、计算新的微观时间尺度损伤变量，式中为微观时间尺本文档来自技高网...
【技术保护点】
1.一种基于时间双尺度分解的疲劳损伤演化分析方法，其特征在于，包括以下步骤：
2.如权利要求1所述的一种基于时间双尺度分解的疲劳损伤演化分析方法，其特征在于，
3.如权利要求2所述的一种基于时间双尺度分解的疲劳损伤演化分析方法，其特征在于，
4.如权利要求3所述的一种基于时间双尺度分解的疲劳损伤演化分析方法，其特征在于，
5.如权利要求4所述的一种基于时间双尺度分解的疲劳损伤演化分析方法，其特征在于，
6.如权利要求5所述的一种基于时间双尺度分解的疲劳损伤演化分析方法，其特征在于，
7.如权利要求6所述的一种基于时间双尺度分解的疲劳损伤演化分析方法，其特征在于，
8.如权利要求7所述的一种基于时间双尺度分解的疲劳损伤演化分析方法，其特征在于，
【技术特征摘要】
1.一种基于时间双尺度分解的疲劳损伤演化分析方法，其特征在于，包括以下步骤：
2.如权利要求1所述的一种基于时间双尺度分解的疲劳损伤演化分析方法，其特征在于，
3.如权利要求2所述的一种基于时间双尺度分解的疲劳损伤演化分析方法，其特征在于，
4.如权利要求3所述的一种基于时间双尺度分解的疲劳损伤演化分析方法，其特征在于，
<...【专利技术属性】
技术研发人员：梁俊松，周泽林，丁兆东，郭珊珊，李鑫，王建，
申请(专利权)人：中国十九冶集团有限公司，
类型：发明
国别省市：
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。