CN114112676B - 一种木材横纹抗压全时程本构关系的构建方法 - Google Patents
一种木材横纹抗压全时程本构关系的构建方法 Download PDFInfo
- Publication number
- CN114112676B CN114112676B CN202111463925.6A CN202111463925A CN114112676B CN 114112676 B CN114112676 B CN 114112676B CN 202111463925 A CN202111463925 A CN 202111463925A CN 114112676 B CN114112676 B CN 114112676B
- Authority
- CN
- China
- Prior art keywords
- load
- straight line
- wood
- point
- displacement
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000007906 compression Methods 0.000 title claims abstract description 142
- 230000006835 compression Effects 0.000 title claims abstract description 138
- 238000010276 construction Methods 0.000 title claims abstract description 18
- 239000002023 wood Substances 0.000 claims abstract description 161
- 238000000034 method Methods 0.000 claims abstract description 72
- 238000012360 testing method Methods 0.000 claims description 128
- 238000006073 displacement reaction Methods 0.000 claims description 125
- 238000012669 compression test Methods 0.000 claims description 16
- 238000012545 processing Methods 0.000 claims description 7
- 238000005056 compaction Methods 0.000 claims description 5
- XOFYZVNMUHMLCC-ZPOLXVRWSA-N prednisone Chemical compound O=C1C=C[C@]2(C)[C@H]3C(=O)C[C@](C)([C@@](CC4)(O)C(=O)CO)[C@@H]4[C@@H]3CCC2=C1 XOFYZVNMUHMLCC-ZPOLXVRWSA-N 0.000 claims description 4
- 238000011156 evaluation Methods 0.000 abstract description 8
- 238000004088 simulation Methods 0.000 abstract description 8
- 239000004566 building material Substances 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 230000003068 static effect Effects 0.000 description 3
- 241000219000 Populus Species 0.000 description 2
- 238000012795 verification Methods 0.000 description 2
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 238000005452 bending Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 229910052799 carbon Inorganic materials 0.000 description 1
- 238000004134 energy conservation Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000001727 in vivo Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000004806 packaging method and process Methods 0.000 description 1
- 230000000630 rising effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N3/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N3/08—Investigating strength properties of solid materials by application of mechanical stress by applying steady tensile or compressive forces
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2203/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N2203/0058—Kind of property studied
- G01N2203/0069—Fatigue, creep, strain-stress relations or elastic constants
- G01N2203/0075—Strain-stress relations or elastic constants
Landscapes
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
一种木材横纹抗压全时程本构关系的构建方法。本发明的目的是提供一种能够准确定量地评估木材横纹抗压的全时程本构关系，用于古建筑木结构和现代木结构的安全评估的木材横纹抗压全时程本构关系的构建方法。在该本构关系构建方法中，基于能量等效原理方法构建的应力－应变关系模型采用直线来拟合，包含了横纹抗压过程中第一阶段线弹性段、第二阶段平台段和第三阶段密实段，能够简洁准确地描述木材横纹抗压全时程本构关系，克服传统参数回归拟合方法缺乏物理含义且拟合参数多的问题；定量确定的屈服点和硬化点，能够克服传统方法取比例线性载荷为屈服点导致模拟第二阶段平台段起点过低的问题，也克服了传统方法仅能定性描述硬化点位置的问题。
Description
技术领域
本发明涉及绿色建材领域，特别是涉及一种用于对古建筑木结构和现代木结构的精准定量安全评估的木材横纹抗压全时程本构关系的构建方法。
背景技术
木材是一种绿色可持续的建材，促进其在建筑结构中的安全应用，可为我国建筑的节能减排、双碳目标的实现提供助力。在古建筑木结构的斗拱、榫卯节点处，木材横纹承压现象普遍存在，且对于逐渐兴起的现代木结构，在节点处也普遍存在木材横纹承压现象。由于木材材性的各向异性，横纹方向的模量和线性比例极限强度较弱，在经受较大荷载时，木材横纹抗压会由第一阶段线弹性段进入第二阶段平台段，对于遭受较大荷载时，例如罕遇地震时，木材横纹抗压将会由第二阶段平台段进入第三阶段密实段。因此，需要准确描述木材横纹抗压全时程的应力－应变关系，才可能对古建筑木结构和现代木结构实施精准定量的安全评估。
但是目前，由于缺乏木材横纹抗压全时程本构关系的构建方法，从而无法准确定量地评估木材横纹抗压全时程的应力－应变关系，导致目前古建筑木结构和现代木结构安全评估的精度较低，大多仍然停留在定性经验分析的水平。因此，通过提出木材横纹抗压全时程本构关系的构建方法，来准确定量地评估木材横纹抗压全时程的应力－应变关系，对古建筑木结构和现代木结构的精准定量安全评估来说，具有十分重要的意义。
经有关文献和资料检索，如中国标准GB/T 1939、美国标准ASTM D143、日本标准JIS Z2101、德国标准DIN52192等标准仅规定了线性比例极限强度的确定方法，未涉及应力－应变本构关系；已公开的专利“基于木材抗压强度预测木材抗剪强度的方法”(申请号：202011087516.6)、“基于ABAQUS的木材本构关系数值模拟方法”(申请号：201810039597.9)、“一种体内预应力胶合木梁抗弯承载能力的分析方法”(申请号：201711128330.9)，以及论文包装工程2019年第40卷第11期第86～93页的“杨木静动态压缩本构模型研究”、湖南大学学报(自然科学版)2018年第45卷第3期第86～93页的“木材径向反复抗压应力－应变模型研究”、土木工程学报2017年第50卷第4期第80～88页的“一种描述木材抗压的非线性本构模型及试验验证”、安徽农业大学学报2010年第37卷第4期第665～668页的“杨木压缩应力－应变关系研究”、力学与实践2000年第22卷第5期第25～27页的“木材静压大变形本构关系研究”等所涉及的横纹抗压，主要是规定了第一阶段线弹性段和第二阶段平台段的应力－应变本构关系，且其屈服强度是由线性比例极限载荷计算得到，导致了第二阶段平台段起点过低，并且第二阶段平台段也多是采用多参数回归拟合的方法得到的。
从上述文献公开的方法之中，可知目前对于横纹抗压过程中第一阶段线弹性段、第二阶段平台段和第三阶段密实段全时程的应力－应变关系尚缺乏相应的相应技术，尚不能满足本领域安全评估的需要。
发明内容
本发明的主要目的是提出一种木材横纹抗压全时程本构关系的构建方法，能精准定量地评估木材横纹抗压的本构关系，及其用于古建筑木结构和现代木结构的安全评估的应用。
为了解决上述技术问题，本申请提供了如下技术方案：
本发明提出的木材横纹抗压全时程本构关系的构建方法，包括以下步骤：
首先加工试件并进行横纹抗压力学性能试验，然后将木材横纹抗压力学性能试验全时程中的荷载－位移曲线定为三个阶段：包含了横纹抗压过程中第一阶段线弹性段、第二阶段平台段和第三阶段密实段，对荷载－位移曲线三个阶段对应的应力与应变关系采用直线来进行定义，第二阶段平台段对应的应力与应变关系直线是基于能量等效原理确定，然后根据上述关系确定木材横纹抗压全时程本构关系。
其中，较佳的是，本发明木材横纹抗压全时程本构关系的构建方法，其具体步骤可包含如下内容：
(1)将待测木材加工成为长方体试件或立方体试件，优选为试件的最小边长不小于5mm，试件的高度、宽度和厚度分别记为h、w、t，对试件进行横纹抗压力学性能试验，h方向为横纹抗压力学试验的加载方向，横纹抗压力学性能试验的应力记为σ、应变记为ε、加载位移记为Δ，加载位移Δ达到0.8倍的试件高度h：Δ＝0.8h，停止加载；
(2)将木材横纹抗压力学性能试验全时程中的荷载－位移曲线定为三个阶段，第一阶段线弹性段，第二阶段平台段，第三阶段密实段，对荷载－位移曲线三个阶段对应的应力与应变关系采用直线来进行定义，第二阶段平台段对应的应力与应变关系直线将基于能量等效原理确定；
(3)通过取试件表面应变值，确定木材横纹抗压全时程本构关系中第一阶段线弹性段的应力σ与应变ε关系，
其具体步骤可为如下所述：
(3-1)采用应变片直接接触法或数字散斑图像非接触方法，测试木材横纹抗压试验中试件高度中央位置沿加载方向的表面应变值；
(3-2)取木材横纹抗压试验全时程中第一阶段线弹性段的两个加载点，分别记为加载点A和点B，点A对应的荷载为FA、加载位移为ΔA、加载方向的表面应变值为εA，点B对应的荷载为FB、加载位移为ΔB、加载方向的表面应变值为εB；
(3-3)以上述点A作为第一点，以其对应荷载为FA作为纵坐标，加载位移为ΔA作为横坐标，点B作为第二点，以其对应荷载为FB作为纵坐标，加载位移为ΔB作为横坐标，作直线I：直线I与坐标轴横轴交点记为O，直线I与木材横纹抗压力学性能试验全时程中的荷载－位移曲线的上限交点记为C；
(3-4)根据上述步骤(3-2)得到的加载点A、点B的荷载FA和FB，以及上述步骤(1)试件的宽度w和厚度t，得到加载点A、点B对应的应力σA和σB：σA＝FA/(wt)、σB＝FB/(wt)；
(3-5)根据上述步骤(3-2)得到加载点A、点B的表面应变值εA和εB，以及上述步骤(3-4)得到加载点A、点B的应力σA和σB，确定木材横纹抗压全时程本构关系中第一阶段线弹性段的应力σ与应变ε关系定义为：
(4)通过取试件加载位移，确定木材横纹抗压全时程本构关系中第三阶段密实段的应力σ与应变ε关系，
该具体步骤可如下所述：
(4-1)采用位移计测量方法，测试木材横纹抗压试验全时程中试件加载位移；
(4-2)取木材横纹抗压试验全时程中第三阶段密实段的两个加载点，分别记为加载点K和点L，点K对应的荷载为FK、加载位移为ΔK，点L对应的荷载为FL、加载位移为ΔL；
(4-3)以上述点K作为第一点，以其对应荷载为FK作为纵坐标，加载位移为ΔK作为横坐标，点L作为第二点，以其对应荷载为FL作为纵坐标，加载位移为ΔL作为横坐标，作直线II：直线II与坐标轴横轴交点记为G，直线II与木材横纹抗压力学性能试验全时程中的荷载－位移曲线的下限交点记为J；
(4-4)根据步骤(4-2)得到的加载点K、点L的荷载FK和FL，以及上述步骤(1)试件的高度h、宽度w和厚度t，得到对应的应力σK和σL：σK＝FK/(wt)、σL＝FL/(wt)；
(4-5)根据上述步骤(1)得到的试件高度h，步骤(4-4)得到的加载点K、点L的应力σK和σL，以及上述步骤(4-2)得到的加载位移ΔK和ΔL，确定木材横纹抗压全时程本构关系中第三阶段密实段的应力σ与应变ε关系定义为：
(4-6)将木材横纹抗压力学性能试验全时程中的荷载－位移曲线与直线II、坐标轴横轴包络的面积记为S1；
(5)通过能量等效原理，确定木材横纹抗压全时程本构关系中第二阶段平台段的应力σ与应变ε关系，
该具体步骤可如下所述：
(5-1)根据上述步骤(1)试件的高度h和步骤(3)得到的直线I，将直线I沿坐标轴横轴方向向右平移5％h，得到直线III：直线III与坐标轴横轴交点记为M，直线III与木材横纹抗压力学性能试验全时程中的荷载－位移曲线的交点记为E，其对应荷载为FE、加载位移为ΔE；
(5-2)过上述交点E作一直线，直线IV的斜率记为SLi，i＝1,…k…,n，记为直线IV：F＝SLi(Δ-ΔE)+FE，直线IV斜率SLi的取值以0为起始点，以0.001递增，直线斜率SLi的取值范围为直线IV与直线I、直线II、坐标轴横轴包络的面积记为S2，直线IV与直线I、木材横纹抗压力学性能试验全时程中的荷载－位移曲线包络的面积记为S3，直线IV与木材横纹抗压力学性能试验全时程中的荷载－位移曲线包络的面积记为S4，直线IV与直线II、木材横纹抗压力学性能试验全时程中的荷载－位移曲线包络的面积记为S5，直线IV与直线I的交点记为D，其对应荷载为FD、加载位移为ΔD、应力σD＝FD/(wt)，直线IV与直线II的交点记为H，其对应荷载为FH、加载位移为ΔH、应力σH＝FH/(wt)；
(5-3)当直线IV的斜率取为SLk，直线IV与直线I、直线II、坐标轴横轴包络的面积S2满足即/>则斜率SLk为符合能量等效原理的最终确定值，因此木材横纹抗压全时程本构关系中第二阶段平台段的应力σ与应变ε关系为：
(6)根据上述步骤(1)～(5)，确定木材横纹抗压全时程本构关系，
综上所述该模型可为：
本发明还提供了一种木结构的定量安全评估方法，是将以上任一方法构建的木材横纹抗压全时程本构关系用于对古建筑木结构和现代木结构进行精准的定量安全评估。
与现有技术相比，本发明木材横纹抗压全时程本构关系的构建方法至少具有以下有益效果：
本发明木材横纹抗压全时程本构关系的构建方法中基于能量等效原理方法构建的应力－应变关系模型采用直线来拟合，包含了横纹抗压过程中第一阶段线弹性段、第二阶段平台段和第三阶段密实段，能够简洁准确地描述木材横纹抗压全时程本构关系，克服传统参数回归拟合方法缺乏物理含义且拟合参数多的问题；定量确定的屈服点和硬化点，能够克服传统方法取比例线性载荷为屈服点导致模拟第二阶段平台段起点过低的问题，也克服了传统方法仅能定性描述硬化点位置的问题。本发明方法能够准确定量地评估木材横纹抗压的本构关系，用于古建筑木结构和现代木结构的安全评估。
下面结合附图对本发明木材横纹抗压全时程本构关系的构建方法做进一步说明。
附图说明
图1为本发明方法中涉及的木材横纹抗压力学性能试验全时程荷载－位移曲线三个阶段的示意图；
图2为本发明方法中涉及的木材横纹抗压力学性能试验全时程本构关系确定方法的示意图；
图3为本发明方法中涉及的木材横纹抗压力学性能试验荷载－位移曲线所包络的面积S1的示意图；
图4为本发明方法中涉及的木材横纹抗压本构关系所包络的面积S2的示意图；
图5为本发明方法中涉及的木材横纹抗压本构关系构建的能量等效原理的示意图；
图6为本发明方法中涉及的木材横纹抗压本构关系构建的流程图；
图7为本发明方法中涉及的具体实施例木材横纹抗压的有限元模型图；
图8为本发明方法中涉及的具体实施例木材横纹抗压的有限元模拟计算结果。
具体实施方式
本发明提供的一种木材横纹抗压全时程本构关系的构建方法，是为了能够准确定量地评估木材横纹抗压的全时程本构关系，用于古建筑木结构和现代木结构的安全评估。
如图6所示，为本发明方法中涉及的木材横纹抗压本构关系构建的流程图。本发明的本构关系构建方法是基于能量等效原理方法构建的应力－应变关系模型采用直线来拟合，包含了横纹抗压过程中第一阶段线弹性段、第二阶段平台段和第三阶段密实段。
其能够简洁准确地描述木材横纹抗压全时程本构关系，克服传统参数回归拟合方法缺乏物理含义且拟合参数多的问题；定量确定的屈服点和硬化点，能够克服传统方法取比例线性载荷为屈服点导致模拟第二阶段平台段起点过低的问题，也克服了传统方法仅能定性描述硬化点位置的问题。
如图1～6所示，本发明木材横纹抗压全时程本构关系的构建方法，包括以下步骤：
(1)加工试件并进行横纹抗压力学性能试验，
将待测木材加工成为长方体试件或立方体试件，且通常试件的最小边长不应小于5mm，将试件的高度、宽度和厚度分别记为h、w、t，并对该试件进行横纹抗压力学性能试验，其中，h方向为横纹抗压力学试验的加载方向，横纹抗压力学性能试验的应力记为σ、应变记为ε、加载位移记为Δ，加载位移Δ达到0.8倍(为根据测试结果总结出来的经验值)的试件高度h：Δ＝0.8h，停止加载；
(2)将木材横纹抗压力学性能试验全时程中的荷载－位移曲线定为三个阶段，
如图1所示，第一阶段线弹性段，第二阶段平台段，第三阶段密实段，其中，对荷载－位移曲线三个阶段对应的应力与应变关系采用直线来进行定义，第二阶段平台段对应的应力与应变关系直线将基于能量等效原理确定；
(3)通过取试件表面应变值，确定木材横纹抗压全时程本构关系中第一阶段线弹性段的应力σ与应变ε关系，具体步骤如下：
(3-1)采用应变片直接接触法或数字散斑图像非接触方法，测试木材横纹抗压试验中试件高度中央位置沿加载方向的表面应变值；
(3-2)取木材横纹抗压试验全时程中第一阶段线弹性段的两个加载点，分别记为加载点A和点B，如图2所示，点A对应的荷载为FA、加载位移为ΔA、加载方向的表面应变值为εA，点B对应的荷载为FB、加载位移为ΔB、加载方向的表面应变值为εB；
(3-3)以上述点A作为第一点，以其对应荷载为FA作为纵坐标，加载位移为ΔA作为横坐标，点B作为第二点，以其对应荷载为FB作为纵坐标，加载位移为ΔB作为横坐标，作直线I：如图2所示，直线I与坐标轴横轴交点记为O，直线I与木材横纹抗压力学性能试验全时程中的荷载－位移曲线的上限交点记为C；
(3-4)根据上述步骤(3-2)得到的加载点A、点B的荷载FA和FB，以及上述步骤(1)试件的宽度w和厚度t，得到加载点A、点B对应的应力σA和σB：σA＝FA/(wt)、σB＝FB/(wt)；
(3-5)根据上述步骤(3-2)得到加载点A、点B的表面应变值εA和εB，以及上述步骤(3-4)得到加载点A、点B的应力σA和σB，确定木材横纹抗压全时程本构关系中第一阶段线弹性段的应力σ与应变ε关系定义为：
(4)通过取试件加载位移，确定木材横纹抗压全时程本构关系中第三阶段密实段的应力σ与应变ε关系，具体步骤如下：
(4-1)采用位移计测量方法，测试木材横纹抗压试验全时程中试件加载位移；
(4-2)取木材横纹抗压试验全时程中第三阶段密实段的两个加载点，分别记为加载点K和点L，如图2所示，点K对应的荷载为FK、加载位移为ΔK，点L对应的荷载为FL、加载位移为ΔL；
(4-3)以上述点K作为第一点，以其对应荷载为FK作为纵坐标，加载位移为ΔK作为横坐标，点L作为第二点，以其对应荷载为FL作为纵坐标，加载位移为ΔL作为横坐标，作直线II：如图2所示，直线II与坐标轴横轴交点记为G，直线II与木材横纹抗压力学性能试验全时程中的荷载－位移曲线的下限交点记为J；
(4-4)根据步骤(4-2)得到的加载点K、点L的荷载FK和FL，以及上述步骤(1)试件的高度h、宽度w和厚度t，得到对应的应力σK和σL：σK＝FK/(wt)、σL＝FL/(wt)；
(4-5)根据上述步骤(1)得到的试件高度h，步骤(4-4)得到的加载点K、点L的应力σK和σL，以及上述步骤(4-2)得到的加载位移ΔK和ΔL，确定木材横纹抗压全时程本构关系中第三阶段密实段的应力σ与应变ε关系定义为：
(4-6)将木材横纹抗压力学性能试验全时程中的荷载－位移曲线与直线II、坐标轴横轴包络的面积记为S1，如图3所示；
(5)通过能量等效原理，确定木材横纹抗压全时程本构关系中第二阶段平台段的应力σ与应变ε关系，具体步骤如下：
(5-1)根据上述步骤(1)试件的高度h和步骤(3)得到的直线I，将直线I沿坐标轴横轴方向向右平移5％h，得到直线III：如图2所示，直线III与坐标轴横轴交点记为M，直线III与木材横纹抗压力学性能试验全时程中的荷载－位移曲线的交点记为E，其对应荷载为FE、加载位移为ΔE；
(5-2)过上述交点E作一直线，记为直线IV，如图2所示，直线IV的斜率记为SLi，i＝1,…k…,n，记为直线IV：F＝SLi(Δ-ΔE)+FE，直线IV斜率SLi的取值以0为起始点，以0.001递增，直线斜率SLi的取值范围为直线IV与直线I、直线II、坐标轴横轴包络的面积记为S2，如图4所示，直线IV与直线I、木材横纹抗压力学性能试验全时程中的荷载－位移曲线包络的面积记为S3，直线IV与木材横纹抗压力学性能试验全时程中的荷载－位移曲线包络的面积记为S4，直线IV与直线II、木材横纹抗压力学性能试验全时程中的荷载－位移曲线包络的面积记为S5，如图5所示，直线IV与直线I的交点记为D，其对应荷载为FD、加载位移为ΔD、应力为σD＝FD/(wt)，直线IV与直线II的交点记为H，其对应荷载为FH、加载位移为ΔH、应力为σH＝FH/(wt)；
(5-3)当直线IV的斜率取为SLk，直线IV与直线I、直线II、坐标轴横轴包络的面积S2满足即/>则斜率SLk为符合能量等效原理的最终确定值，因此木材横纹抗压全时程本构关系中第二阶段平台段的应力σ与应变ε关系为：
(6)根据上述步骤(1)～(5)，确定木材横纹抗压全时程本构关系为：
下面通过一具体实施例，并结合附图，对本发明的技术方案做进一步具体的说明：
(1)试件加工及横纹抗压力学性能试验
将待测杉木木材加工成为高度h、宽度w和厚度t分别为30mm、20mm、20mm的长方体试件，对该试件进行横纹抗压力学性能试验，其中，h方向为横纹抗压力学试验的加载方向，横纹抗压力学性能试验的应力记为σ、单位为MPa，应变记为ε、单位为με，加载位移记为Δ、单位为mm，当加载位移Δ达到0.8倍的试件高度(h：Δ＝0.8h＝24mm)时，停止加载；
(2)将杉木木材横纹抗压力学性能试验全时程中的荷载－位移曲线定为三个阶段，如图1所示，其包括第一阶段线弹性段，第二阶段平台段，第三阶段密实段，对荷载－位移曲线三个阶段对应的应力与应变关系采用直线来进行定义，第二阶段平台段对应的应力与应变关系直线将基于能量等效原理确定；
(3)通过取试件表面应变值，确定杉木木材横纹抗压全时程本构关系中第一阶段线弹性段的应力σ与应变ε关系，具体步骤如下：
(3-1)采用应变片直接接触法或数字散斑图像非接触方法，测试木材横纹抗压试验中试件高度中央位置沿加载方向的表面应变值；
(3-2)取木材横纹抗压试验全时程中第一阶段线弹性段的两个加载点，分别记为加载点A和点B，如图2所示，点A对应的荷载FA为0.501kN、加载位移ΔA为0.212mm、加载方向的表面应变值εA为126με，点B对应的荷载FB为1.143kN、加载位移ΔB为0.471mm、加载方向的表面应变值εB为298με；
(3-3)以上述点A作为第一点，以其对应荷载为FA作为纵坐标，加载位移为ΔA作为横坐标，点B作为第二点，以其对应荷载为FB作为纵坐标，加载位移为ΔB作为横坐标，作直线I：如图2所示，直线I与坐标轴横轴交点记为O，点O对应的荷载FO为0kN、加载位移ΔO为0.010mm，直线I与木材横纹抗压力学性能试验全时程中的荷载－位移曲线的上限交点记为C，点C对应的荷载FC为1.182kN、加载位移ΔO为0.487mm；
(3-4)根据上述步骤(3-2)得到的加载点A、点B的荷载FA和FB，以及上述步骤(1)试件的宽度w和厚度t，得到加载点A、点B对应的应力σA和σB：
(3-5)根据上述步骤(3-2)得到加载点A、点B的表面应变值εA和εB，以及上述步骤(3-4)得到加载点A、点B的应力σA和σB，确定木材横纹抗压全时程本构关系中第一阶段线弹性段的应力σ与应变ε关系定义为：
(4)通过取试件加载位移，确定杉木木材横纹抗压全时程本构关系中第三阶段密实段的应力σ与应变ε关系，具体步骤如下：
(4-1)采用位移计测量方法，测试木材横纹抗压试验全时程中试件加载位移；
(4-2)取木材横纹抗压试验全时程中第三阶段密实段的两个加载点，分别记为加载点K和点L，如图2所示，点K对应的荷载FK为5.316kN、加载位移ΔK为为17.709mm，点L对应的荷载FL为6.973kN、加载位移ΔL为18.033mm；
(4-3)以上述点K作为第一点，以其对应荷载为FK作为纵坐标，加载位移为ΔK作为横坐标，点L作为第二点，以其对应荷载为FL作为纵坐标，加载位移为ΔL作为横坐标，作直线II：如图2所示，直线II与坐标轴横轴交点记为G，点G对应的荷载FG为0kN、加载位移ΔG为16.670mm，直线II与木材横纹抗压力学性能试验全时程中的荷载－位移曲线的下限交点记为J，点J对应的荷载FJ为4.558kN、加载位移ΔJ为17.561mm；
(4-4)根据步骤(4-2)得到的加载点K、点L的荷载FK和FL，以及上述步骤(1)试件的高度h、宽度w和厚度t，得到对应的应力σK和σL：
(4-5)根据上述步骤(1)得到的试件高度h，步骤(4-4)得到的加载点K、点L的应力σK和σL，以及上述步骤(4-2)得到的加载位移ΔK和ΔL，确定木材横纹抗压全时程本构关系中第三阶段密实段的应力σ与应变ε关系定义为：
(4-6)将木材横纹抗压力学性能试验全时程中的荷载－位移曲线与直线II、坐标轴横轴包络的面积记为S1，S1面积为38.554kN·mm，如图3所示；
(5)通过能量等效原理，确定杉木木材横纹抗压全时程本构关系中第二阶段平台段的应力σ与应变ε关系，具体步骤如下：
(5-1)根据上述步骤(1)试件的高度h和步骤(3)得到的直线I，将直线I沿坐标轴横轴方向向右平移5％h，得到直线III：如图2所示，直线III与坐标轴横轴交点记为M，点M对应的荷载FM为0kN、加载位移ΔM为1.510mm，直线III与木材横纹抗压力学性能试验全时程中的荷载－位移曲线的交点记为E，其对应荷载FE为1.589kN、加载位移ΔE为2.151mm；
(5-2)过上述交点E作一直线，记为直线IV，如图2所示，直线IV的斜率记为SLi，i＝1,…k…,n，记为直线IV：F＝SLi(Δ-ΔE)+FE，直线IV斜率SLi的取值以0kN/mm为起始点，以0.001kN/mm递增，直线斜率SLi的取值范围为 直线IV与直线I、直线II、坐标轴横轴包络的面积记为S2，如图4所示，直线IV与直线I、木材横纹抗压力学性能试验全时程中的荷载－位移曲线包络的面积记为S3，直线IV与木材横纹抗压力学性能试验全时程中的荷载－位移曲线包络的面积记为S4，直线IV与直线II、木材横纹抗压力学性能试验全时程中的荷载－位移曲线包络的面积记为S5，如图5所示，直线IV与直线I的交点记为D，其对应荷载为FD、加载位移为ΔD，直线IV与直线II的交点记为H，其对应荷载为FH、加载位移为ΔH；
(5-3)当直线IV的斜率取为SLk＝0.121kN/mm，直线IV：F＝SLi(Δ-ΔE)+FE＝0.121×(Δ-2.151)+1.589，直线IV与直线I的交点记为D，其对应荷载FD为1.398kN、加载位移ΔD为0.574mm、应力直线IV与直线II的交点记为H，其对应荷载FH为3.427kN、加载位移ΔH为17.340mm、应力/>直线IV与直线I、直线II、坐标轴横轴包络的面积S2为38.758kN.75，满足直线IV与直线I、木材横纹抗压力学性能试验全时程中的荷载－位移曲线包络的面积S3为0.586kN·mm，直线IV与木材横纹抗压力学性能试验全时程中的荷载－位移曲线包络的面积记S4为1.347kN·mm，直线IV与直线II、木材横纹抗压力学性能试验全时程中的荷载－位移曲线包络的面积记S5为1.923kN·mm，如图5所示，满足/>斜率SLk＝0.121kN/mm为符合能量等效原理的最终确定值，因此木材横纹抗压全时程本构关系中第二阶段平台段的应力σ与应变ε关系为：
(6)根据上述步骤(1)～(5)，确定杉木木材横纹抗压全时程本构关系为：
在验证实施例中，根据上述步骤(1)所描述的杉木木材试件的尺寸，采用有限元软件ABAQUS建立杉木试件的有限元几何模型，如图7所示，并将上述步骤(6)确定的木材横纹抗压全时程本构关系输入到有限元软件ABAQUS材料参数定义项中，并通过有限元软件ABAQUS模拟施加竖向位移，得到有限元模拟荷载-荷载位移曲线，与试验荷载-位移曲线吻合较好，如图8所示。
以上所述的实施例仅仅是对本发明的优选实施方式进行描述，并非对本发明的范围进行限定，在不脱离本发明设计精神的前提下，本领域普通技术人员对本发明的技术方案作出的各种变形和改进，均应落入本发明权利要求书确定的保护范围内。
Claims (4)
1.一种木材横纹抗压全时程本构关系的构建方法，其特征在于，首先加工试件并进行横纹抗压力学性能试验，然后将木材横纹抗压力学性能试验全时程中的荷载－位移曲线定为三个阶段：包含了横纹抗压过程中第一阶段线弹性段、第二阶段平台段和第三阶段密实段，其中，对荷载－位移曲线三个阶段对应的应力与应变关系采用直线来进行定义，第二阶段平台段对应的应力与应变关系直线是基于能量等效原理确定，即：
通过取试件表面应变值，确定木材横纹抗压全时程本构关系中第一阶段线弹性段的应力σ与应变ε关系的具体步骤如下：
(3-1)采用应变片直接接触法或数字散斑图像非接触方法，测试木材横纹抗压试验中试件高度中央位置沿加载方向的表面应变值；
(3-2)取木材横纹抗压试验全时程中第一阶段线弹性段的两个加载点，分别记为加载点A和点B，点A对应的荷载为FA、加载位移为ΔA、加载方向的表面应变值为εA，点B对应的荷载为FB、加载位移为ΔB、加载方向的表面应变值为εB；
(3-3)以上述点A作为第一点，以其对应荷载为FA作为纵坐标，加载位移为ΔA作为横坐标，点B作为第二点，以其对应荷载为FB作为纵坐标，加载位移为ΔB作为横坐标，作直线I：直线I与坐标轴横轴交点记为O，直线I与木材横纹抗压力学性能试验全时程中的荷载－位移曲线的上限交点记为C；
(3-4)根据上述步骤(3-2)得到的加载点A、点B的荷载FA和FB，以及上述步骤(1)试件的宽度w和厚度t，得到加载点A、点B对应的应力σA和σB：σA＝FA/(wt)、σB＝FB/(wt)；
(3-5)根据上述步骤(3-2)得到加载点A、点B的表面应变值εA和εB，以及上述步骤(3-4)得到加载点A、点B的应力σA和σB，确定木材横纹抗压全时程本构关系中第一阶段线弹性段的应力σ与应变ε关系定义为：
通过取试件加载位移，确定木材横纹抗压全时程本构关系中第三阶段密实段的应力σ与应变ε关系的具体步骤如下：
(4-1)采用位移计测量方法，测试木材横纹抗压试验全时程中试件加载位移；
(4-2)取木材横纹抗压试验全时程中第三阶段密实段的两个加载点，分别记为加载点K和点L，点K对应的荷载为FK、加载位移为ΔK，点L对应的荷载为FL、加载位移为ΔL；
(4-3)以上述点K作为第一点，以其对应荷载为FK作为纵坐标，加载位移为ΔK作为横坐标，点L作为第二点，以其对应荷载为FL作为纵坐标，加载位移为ΔL作为横坐标，作直线II：直线II与坐标轴横轴交点记为G，直线II与木材横纹抗压力学性能试验全时程中的荷载－位移曲线的下限交点记为J；
(4-4)根据步骤(4-2)得到的加载点K、点L的荷载FK和FL，以及上述步骤(1)试件的高度h、宽度w和厚度t，得到对应的应力σK和σL：σK＝FK/(wt)、σL＝FL/(wt)；
(4-5)根据上述步骤(1)得到的试件高度h，步骤(4-4)得到的加载点K、点L的应力σK和σL，以及上述步骤(4-2)得到的加载位移ΔK和ΔL，确定木材横纹抗压全时程本构关系中第三阶段密实段的应力σ与应变ε关系定义为：
(4-6)将木材横纹抗压力学性能试验全时程中的荷载－位移曲线与直线I I、坐标轴横轴包络的面积记为S1；
通过能量等效原理，确定木材横纹抗压全时程本构关系中第二阶段平台段的应力σ与应变ε关系的具体步骤如下：
(5-1)根据上述步骤(1)试件的高度h和步骤(3)得到的直线I，将直线I沿坐标轴横轴方向向右平移5％h，得到直线III：直线III与坐标轴横轴交点记为M，直线III与木材横纹抗压力学性能试验全时程中的荷载－位移曲线的交点记为E，其对应荷载为FE、加载位移为ΔE；
(5-2)过上述交点E作一直线，直线IV的斜率记为SLi，i＝1,…k…,n，记为直线IV：F＝SLi(Δ-ΔE)+FE，直线IV斜率SLi的取值以0为起始点，以0.001递增，直线斜率SLi的取值范围为直线IV与直线I、直线II、坐标轴横轴包络的面积记为S2，直线IV与直线I、木材横纹抗压力学性能试验全时程中的荷载－位移曲线包络的面积记为S3，直线IV与木材横纹抗压力学性能试验全时程中的荷载－位移曲线包络的面积记为S4，直线IV与直线II、木材横纹抗压力学性能试验全时程中的荷载－位移曲线包络的面积记为S5，直线IV与直线I的交点记为D，其对应荷载为FD、加载位移为ΔD、应力为σD＝FD/(wt)，直线IV与直线II的交点记为H，其对应荷载为FH、加载位移为ΔH、应力为σH＝FH/(wt)；
(5-3)当直线IV的斜率取为SLk，直线IV与直线I、直线II、坐标轴横轴包络的面积S2满足即/>则斜率SLk为符合能量等效原理的最终确定值，因此木材横纹抗压全时程本构关系中第二阶段平台段的应力σ与应变ε关系为：
根据上述关系确定木材横纹抗压全时程本构关系。
2.根据权利要求1所述的木材横纹抗压全时程本构关系的构建方法，其特征在于，该方法包括以下步骤：
(1)将待测木材加工成为长方体试件或立方体试件，所述试件的高度、宽度和厚度分别记为h、w、t，对试件进行横纹抗压力学性能试验，h方向为横纹抗压力学试验的加载方向，横纹抗压力学性能试验的应力记为σ、应变记为ε、加载位移记为Δ，加载位移Δ达到0.8倍的试件高度h：Δ＝0.8h，停止加载；
(2)将木材横纹抗压力学性能试验全时程中的荷载－位移曲线定为三个阶段，第一阶段线弹性段，第二阶段平台段，第三阶段密实段，对荷载－位移曲线三个阶段对应的应力与应变关系采用直线来进行定义，第二阶段平台段对应的应力与应变关系直线是基于能量等效原理确定；
(3)通过取试件表面应变值，确定木材横纹抗压全时程本构关系中第一阶段线弹性段的应力σ与应变ε关系；
(4)通过取试件加载位移，确定木材横纹抗压全时程本构关系中第三阶段密实段的应力σ与应变ε关系；
(5)通过能量等效原理，确定木材横纹抗压全时程本构关系中第二阶段平台段的应力σ与应变ε关系；
(6)根据上述步骤(1)～(5)，确定木材横纹抗压全时程本构关系。
3.根据权利要求2所述的木材横纹抗压全时程本构关系的构建方法，其特征在于，所述木材横纹抗压全时程本构关系为：
4.根据权利要求1-3任一项所述的木材横纹抗压全时程本构关系的构建方法，其特征在于，所述试件的最小边长大于等于5mm。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111463925.6A CN114112676B (zh) | 2021-12-03 | 2021-12-03 | 一种木材横纹抗压全时程本构关系的构建方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111463925.6A CN114112676B (zh) | 2021-12-03 | 2021-12-03 | 一种木材横纹抗压全时程本构关系的构建方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN114112676A CN114112676A (zh) | 2022-03-01 |
| CN114112676B true CN114112676B (zh) | 2023-12-19 |
Family
ID=80366824
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202111463925.6A Active CN114112676B (zh) | 2021-12-03 | 2021-12-03 | 一种木材横纹抗压全时程本构关系的构建方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114112676B (zh) |
Citations (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101303280A (zh) * | 2008-04-11 | 2008-11-12 | 北京联合大学 | 采用曲线平移法处理材料载荷变形曲线的方法 |
| WO2009124421A1 (zh) * | 2008-04-11 | 2009-10-15 | 北京联合大学 | 确定材料载荷变形曲线初始直线段和力学性能参数的方法 |
| CN102364489A (zh) * | 2011-10-25 | 2012-02-29 | 陈志勇 | 木材复杂各向异性本构关系模型的数值模拟方法 |
| CN103175735A (zh) * | 2012-12-06 | 2013-06-26 | 西南交通大学 | 材料拉伸真实本构曲线测试技术 |
| CN103967161A (zh) * | 2014-05-23 | 2014-08-06 | 香港华艺设计顾问(深圳)有限公司 | 一种对应构件不屈服性能状态所需钢筋面积的配筋方法 |
| CN104181040A (zh) * | 2014-07-30 | 2014-12-03 | 中国林业科学研究院木材工业研究所 | 一种古建筑木构件的剩余抗压承载力测定方法 |
| CN104931348A (zh) * | 2015-06-08 | 2015-09-23 | 西南交通大学 | 圆环径向压缩能量预测材料单轴本构关系测定方法 |
| CN108256211A (zh) * | 2018-01-08 | 2018-07-06 | 大连理工大学 | 基于abaqus的木材本构关系数值模拟方法 |
| CN108844824A (zh) * | 2018-06-27 | 2018-11-20 | 西南交通大学 | 一种基于圆锥形压头的已知材料残余应力测定方法 |
| WO2020143284A1 (zh) * | 2019-01-11 | 2020-07-16 | 华东理工大学 | 一种基于abaqus的多轴蠕变疲劳预测方法 |
| CN112432849A (zh) * | 2020-10-13 | 2021-03-02 | 北京交通大学 | 基于木材抗压强度预测木材抗剪强度的方法 |
| WO2021037036A1 (zh) * | 2019-08-28 | 2021-03-04 | 肖锋 | 一种测定材料真实应力应变曲线的试验与计算方法 |
| CN112541285A (zh) * | 2020-11-11 | 2021-03-23 | 北京交通大学 | 一种适于中国古建筑木结构用材本构关系的数值模拟方法 |
-
2021
- 2021-12-03 CN CN202111463925.6A patent/CN114112676B/zh active Active
Patent Citations (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101303280A (zh) * | 2008-04-11 | 2008-11-12 | 北京联合大学 | 采用曲线平移法处理材料载荷变形曲线的方法 |
| WO2009124421A1 (zh) * | 2008-04-11 | 2009-10-15 | 北京联合大学 | 确定材料载荷变形曲线初始直线段和力学性能参数的方法 |
| CN102364489A (zh) * | 2011-10-25 | 2012-02-29 | 陈志勇 | 木材复杂各向异性本构关系模型的数值模拟方法 |
| CN103175735A (zh) * | 2012-12-06 | 2013-06-26 | 西南交通大学 | 材料拉伸真实本构曲线测试技术 |
| CN103967161A (zh) * | 2014-05-23 | 2014-08-06 | 香港华艺设计顾问(深圳)有限公司 | 一种对应构件不屈服性能状态所需钢筋面积的配筋方法 |
| CN104181040A (zh) * | 2014-07-30 | 2014-12-03 | 中国林业科学研究院木材工业研究所 | 一种古建筑木构件的剩余抗压承载力测定方法 |
| CN104931348A (zh) * | 2015-06-08 | 2015-09-23 | 西南交通大学 | 圆环径向压缩能量预测材料单轴本构关系测定方法 |
| CN108256211A (zh) * | 2018-01-08 | 2018-07-06 | 大连理工大学 | 基于abaqus的木材本构关系数值模拟方法 |
| CN108844824A (zh) * | 2018-06-27 | 2018-11-20 | 西南交通大学 | 一种基于圆锥形压头的已知材料残余应力测定方法 |
| WO2020143284A1 (zh) * | 2019-01-11 | 2020-07-16 | 华东理工大学 | 一种基于abaqus的多轴蠕变疲劳预测方法 |
| WO2021037036A1 (zh) * | 2019-08-28 | 2021-03-04 | 肖锋 | 一种测定材料真实应力应变曲线的试验与计算方法 |
| CN112432849A (zh) * | 2020-10-13 | 2021-03-02 | 北京交通大学 | 基于木材抗压强度预测木材抗剪强度的方法 |
| CN112541285A (zh) * | 2020-11-11 | 2021-03-23 | 北京交通大学 | 一种适于中国古建筑木结构用材本构关系的数值模拟方法 |
Non-Patent Citations (6)
| Title |
|---|
| Effects of Grain Pattern on the Rolling Shear Properties of Wood in Cross-Laminated Timber;Guofang Wu等;Forests(第12期);全文 * |
| STRESS-STRAIN RESPONSE OF WOOD UNDER RADIAL COMPRESSION. PART I. TEST METHOD AND INFLUENCES OF CELLULAR PROPERTIES;Taghi Tabarsa等;Wood and Fiber Science;第32卷(第2期);全文 * |
| 复杂应力状态下木材力学性能的数值模拟;陈志勇;祝恩淳;潘景龙;;计算力学学报(第04期);全文 * |
| 木材径向反复受压应力 － 应变模型研究;谢启芳等;湖南大学学报（自然科学版）;第45卷(第3期);全文 * |
| 杨木静动态压缩本构模型研究;许威;《包装工程》;第40·卷(第11期);全文 * |
| 重组竹基本力学性能的试验分析与研究;吴文清等;武汉理工大学学报;第39卷(第4期);全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN114112676A (zh) | 2022-03-01 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Buchanan et al. | Testing, simulation and design of cold-formed stainless steel CHS columns | |
| Zhao et al. | Structural performance of stainless steel circular hollow sections under combined axial load and bending–Part 1: Experiments and numerical modelling | |
| Zhao et al. | Behaviour and design of stainless steel SHS and RHS beam-columns | |
| Netto et al. | The effect of corrosion defects on the burst pressure of pipelines | |
| Uy et al. | Behaviour of short and slender concrete-filled stainless steel tubular columns | |
| Ahmed et al. | Behavior of eccentrically loaded double circular steel tubular short columns filled with concrete | |
| Arrayago et al. | Experimental study on ferritic stainless steel RHS and SHS beam-columns | |
| Somodi et al. | Flexural buckling resistance of cold-formed HSS hollow section members | |
| CN106124313B (zh) | 混凝土及类似材料受压下综合变形性能的测试装置 | |
| Alashti et al. | Experimental and numerical investigation of ductile damage effect on load bearing capacity of a dented API XB pipe subjected to internal pressure | |
| CN102455263A (zh) | 一种基于载荷-深度曲线获得金属材料力学性能的方法 | |
| Shu et al. | Experimental and theoretical study on the behavior of cold-formed stainless steel stub columns | |
| CN110795879A (zh) | 一种连续的钢管混凝土构件压弯承载力计算方法 | |
| Peng et al. | A novel semi-analytical method based on equivalent energy principle to obtain J resistance curves of ductile materials | |
| Dobrić et al. | Design procedures for cold-formed stainless steel equal-leg angle columns | |
| Shu et al. | Experimental and numerical study of cold-drawn duplex stainless steel square tube columns | |
| CN111539142A (zh) | 一种管道裂纹扩展驱动力计算方法及系统 | |
| CN114112676B (zh) | 一种木材横纹抗压全时程本构关系的构建方法 | |
| Zhou et al. | Global buckling behaviour of bamboo scrimber box columns under axial compression: Experimental tests and numerical modelling | |
| CN108548720B (zh) | I型裂纹弹塑性理论公式获取延性材料j阻力曲线的方法 | |
| Zavalis et al. | Analysis of bed joint influence on masonry modulus of elasticity | |
| Zhang et al. | Experimental investigation into mechanical performance of novel inorganic-bonded bamboo composite columns under axial compression | |
| Howlader et al. | Cold-forming effect on stainless steel sections | |
| RU2599069C1 (ru) | Способ определения предела выносливости материала при растяжении-сжатии | |
| Jang et al. | Predictions of Tensile Strain Capacity for Strain-Based Pipelines With a Circumferential and Internal Surface Flaw |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |