《Results in Materials》：Probabilistic Modeling and Numerical Simulation of the Three-Point Bending Behavior of Wooden Beams
编辑推荐：
本文提出一种简化但具有现实性的数值模型，该模型基于均匀概率分布律（uniform probabilistic law）与Monte Carlo模拟（Monte Carlo simulations, MCS），用于研究木材在静态三点弯曲（three-point
本文提出一种简化但具有现实性的数值模型，该模型基于均匀概率分布律（uniform probabilistic law）与Monte Carlo模拟（Monte Carlo simulations, MCS），用于研究木材在静态三点弯曲（three-point bending, TPB）条件下的失效行为。所提方法将真实非均质梁替换为几何理想化梁，该梁由多根几何特性明确且完全黏结的薄层（lamellae）组成，其力学性能根据均匀概率分布（uniform probabilistic distribution）随机生成。通过将静力学基本方程与由均匀概率分布导出的数学表达式相结合，建立模拟算法以研究三点弯曲下梁的失效。数值结果与文献报道结果在定性上一致。所得的力-位移（force-displacement）与损伤-位移（damage-displacement）曲线，结合严谨的统计分析，揭示了木梁在弯曲下力学行为的离散性，从而有助于更好地控制其固有生物变异性。此外，应力-尺寸（stress-size）关系证实了Ba?ant尺寸效应律（Ba?ant’s size effect law, SEL）在表征木材弯曲响应方面的相关性。
研究背景与问题：木材作为本地可获得建筑材料，具有可再生、易加工和建筑结构应用广等特点，有利于降低建筑成本、减少进口依赖并促进可持续建设与减排。然而，木材本征生物变异性强，属非均质、各向异性和多孔材料，其力学性能跨尺度变化；纤维尺寸与取向、节子及纤维方向偏差等会影响抗弯承载和失效机制。常规有限元方法结合Tsai–Hill与Tsai–Wu等弹性失效准则可描述木材力学响应，但峰后行为和性能统计变异性表征仍困难；随机有限元与连续损伤力学提供了替代途径，但同时考虑缺陷、各向异性、渐进失效和试件尺寸效应仍是重要建模难题。复合材料与架构材料研究表明，拓扑结构、内部构造和加载条件会影响力学性能，为理解材料变异与结构响应关系提供背景。因此，需要简化且具物理可解释性的概率模型，以描述木梁弯曲行为并考虑材料变异、渐进失效和试件尺寸效应。研究人员据此构建基于均匀概率分布与Monte Carlo模拟（MCS）的框架，在薄层（lamella，本文中为计算离散层）尺度表征性能变异，模拟渐进薄层断裂及应力重分布，并将试件尺寸纳入概率分析。研究检验两个假设：概率模型可再现文献报道的弯曲响应离散性；结构抗力随尺寸变化可与Ba?ant尺寸效应律（Ba?ant’s size effect law, SEL）联系。
关键技术方法（不超过250字）：本文为数值模拟研究，无实验样本队列。研究人员将木梁理想化为高h、跨L、宽b、由N根完全黏结薄层组成的矩形梁；各薄层采用线弹性本构及统一杨氏模量（Young’s modulus, E），失效阈值按有符号均匀分布随机赋值，负值为拉伸阈值、正值为压缩阈值。基于Euler–Bernoulli弯曲应力关系与Huygens定理计算有效截面二次矩，按位移增量加载；薄层应力超过阈值即失效并从有效截面刚度中移除，重算中性轴与有效截面二次矩，引起应力重分布。以Monte Carlo模拟（MCS）重复生成梁样本并统计失效、损伤和尺寸效应；参数包括E=7260 MPa、位移增量Δ=0.1 mm、N=80、L=80 cm，h、b和模拟梁数k随算例变化，框架在Fortran中实现以自动化大规模模拟。
研究结果：1）荷载随挠度的演化：通过单梁及300根名义相同梁的三点弯曲模拟，响应先呈线性至约峰值荷载40–50%，随后因最弱薄层（主要在拉伸区）失效而偏离线性，峰后逐渐软化而非突塌；300根梁响应显著离散，反映薄层力学性能随机分布对损伤萌生和扩展的影响。虽局部为线弹性，但全局非线性源于渐进损伤导致的刚度退化与随机失效，而非材料非线性。2）损伤随挠度的演化：通过损伤-挠度曲线，损伤非线性增加并趋于饱和；损伤由最脆弱薄层失效启动，随有效刚度降低和剩余完整区域应力集中而累积加速；几何相同梁的损伤轨迹差异明显，说明损伤演化对弹性极限概率分布高度敏感。3）试件尺寸对结构强度的影响：通过改变试件尺寸的模拟，名义结构强度随试件尺寸增大而降低，与木材弯曲强度研究、尺寸效应文献及弱链概念（weakest-link concept）一致；但模型未显式表示节子、孔洞或纤维偏差，尺寸依赖性来自失效阈值分布与渐进薄层失效的相互作用，而非物理缺陷统计分布。h/L>1的梁在模拟中承载能力更高，说明几何构型影响结构响应。4）关键参数的统计研究：通过300次Monte Carlo模拟的统计分析，损伤萌生平均应力约1.36 MPa，失效受拉伸区控制；峰值应力主要在压缩区达到，压缩区对极限承载贡献更显著；损伤萌生应力分布较宽，说明失效不在单一阈值发生；薄层数增加时离散性降低。损伤峰值平均损伤（damage variable, D）为0.401、标准差0.109，失效薄层数为10–24，损伤峰值平均挠度约为首次失效时的3.8倍。5）最顶层薄层压缩失效概率：通过Monte Carlo框架按首次记录失效位置和应力符号统计，300次模拟中17次首次失效发生在最顶层薄层且为压缩应力，经验概率为5.67%，95% Wilson置信区间（Wilson confidence interval）约3.57%–8.89%；该结果仅针对首次记录失效的位置和模式，不代表整个加载过程最顶层薄层压缩失效的总概率。6）讨论：模型预测的力学响应由渐进损伤与局部失效阈值概率变异性共同控制；弹性向非线性转变由薄层失效、应力重分布和有效截面刚度降低驱动；响应与损伤演化离散性反映失效阈值随机性。名义强度随尺寸增大而降低说明木梁数值建模需考察尺寸效应；但模型未显式表示物理缺陷及其空间分布，模拟的尺寸依赖趋势不能直接归因于遭遇临界缺陷概率增加，应解释为所采用概率框架内的尺寸相关响应，而非缺陷控制弱链机制的直接证据；仍需进一步实验验证和显式断裂过程表示以确定其物理来源。
研究结论（翻译）：本研究建立了概率数值框架，用于分析木梁在三点静态弯曲下的弯曲行为。通过将均匀概率分布与Monte Carlo模拟及结构力学基本原理耦合，所提模型有效捕捉木材固有非均质性和各向异性。该框架在Fortran中实现，可自动化大规模模拟，提供整体结构响应并揭示局部损伤机制。结果表明，与文献实验发现具有良好定性和定量一致性。模拟响应显示初始线弹性阶段约至峰值荷载40–50%，随后为渐进薄层失效控制的非线性阶段和延性行为。损伤萌生平均应力约1.36 MPa，具有显著离散性，反映材料随机性。对大量模拟梁的统计分析证实损伤演化和峰值强度均有较大变异性。此外，研究显示明显尺寸效应，结构强度随试件尺寸增大而降低，与Ba?ant尺寸效应律一致。几何参数影响亦得到体现，较高长细比（h/l>1）梁具有更高承载能力。除再现观测力学行为外，模型将分析能力扩展至难以实验研究的大尺度结构。总体而言，所提方法为将材料变异性纳入木结构分析与设计提供高效稳健工具，在可靠性评估和概率工程应用中具潜力。未来工作将扩展至单轴拉伸等其他加载条件，并整合更多概率分布以提高预测精度和通用性。
