[论文解读] Brittle Crack Roughness in Three-Dimensional Beam Lattices

B. Skjetne, Torbjørn Helle+1 more|arXiv (Cornell University)|May 25, 2005

Numerical methods in engineering被引用 4

一句话总结

本研究通过数值模拟，研究了在强无序条件下三维弹性杆格点中的裂纹粗糙度，以模拟脆性断裂。结果发现，裂纹粗糙度指数 ζ = 0.59 ± 0.01 具有普适性，与小尺度实验观测结果一致，表明在无序介质中缓慢、准静态裂纹扩展存在一个具有特征长度尺度的特定标度行为 regime。

ABSTRACT

The roughness exponent is reported in numerical simulations with a three-dimensional elastic beam lattice. Two different types of disorder have been used to generate the breaking thresholds, i.e., distributions with a tail towards either strong or weak beams. Beyond the weak disorder regime a universal exponent of 0.59(1) is obtained. This is within the range 0.4-0.6 reported experimentally for small scale quasi-static fracture, as would be expected for media with a characteristic length scale.

研究动机与目标

• 研究在不同无序度下三维弹性杆格点中裂纹粗糙度的标度行为。
• 确定在三维弹性杆网络中，粗糙度指数 ζ 是否在强无序条件下表现出普适性。
• 将模拟结果与实验数据进行比较，评估杆模型对真实脆性断裂过程的相关性。
• 探讨是否存在明显不同的标度行为——特别是小尺度下 ζ 值较低的情况——并与实验观测结果一致。
• 评估特征长度尺度在断裂形貌中的作用，尤其是在具有离散弹性元件（如杆格点）的系统中。

提出的方法

• 模拟采用具有固定连通性的立方杆格点，每个杆可抵抗轴向力、剪切力、弯曲力和扭转力。
• 通过共轭梯度法迭代求解机械平衡，以最小化弹性能量，确保每个节点处的力和力矩平衡。
• 通过幂律无序模型分配断裂阈值：r^D，其中 r ∈ [0,1]，D 控制无序强度及弱杆或强杆的分布方向。
• 通过在顶面施加均匀位移来启动断裂，杆按其轴向强度最弱的顺序依次断裂，每次断裂后进行应力重分布。
• 通过计算垂直于平均断裂平面的裂纹表面高度的均方根波动，确定粗糙度指数 ζ，使用 W ∼ L^ζ 的关系。
• 系统尺寸 L 范围为 3 至 40，通过不同 D 值的 W 与 L 的对数图上的线性拟合提取标度指数。

实验结果

研究问题

• RQ1 在强无序条件下，三维弹性杆格点中的粗糙度指数 ζ 是否收敛到一个普适值？
• RQ2 粗糙度指数 ζ 如何依赖于无序类型，特别是当阈值分布偏向弱杆或强杆时？
• RQ3 所观测到的粗糙度指数是否与小尺度、准静态断裂实验报告的值一致？
• RQ4 杆格点模型是否表现出如先前研究建议的小尺度下明显的标度行为？
• RQ5 与随机熔断模型或 Born 模型等标量模型相比，杆模型在预测裂纹粗糙度方面表现如何？

主要发现

• 在强无序条件下（D = 4），粗糙度指数收敛至 ζ = 0.59(1)，表明三维杆格点中存在普适的标度行为。
• 该指数 ζ = 0.59(1) 低于大尺度断裂实验报告的 ζ = 0.8，但处于小尺度、准静态实验中观测到的 0.4 ≤ ζ ≤ 0.6 范围内。
• 当 D = 2 时，指数为 ζ = 0.62(2)；当 D = 1 时，指数为 ζ = 0.77(1)，表明系统仅在强无序条件下进入普适标度区。
• 当 D = -4 时，结果为 ζ = 0.65(2)，略高但仍与普适值一致，表明在三维系统中，负 D 值的瞬态标度区可能延伸至更高的 |D|。
• 杆模型能够保持旋转对称性并真实模拟弹性变形，其粗糙度指数与小尺度实验数据一致，支持存在一个独特标度行为的存在。
• 研究结果支持如下假设：具有特征长度尺度的系统（如杆格点）在小尺度下表现出较低的粗糙度指数，这与 Cosserat 弹性理论和小尺度断裂行为一致。