第一作者:常鹏(1979—),男,陕西泾阳人,副教授,博士.研究方向为古建结构安全保护.email:pchang@bjtu.edu.cn.
古建筑木结构在长期竖向荷载作用下会产生累积变形,针对累积变形古建筑木柱在竖向荷载作用下平面应力问题,运用弹性力学方法求解了应力与位移分量及屈服荷载解析解.结合半逆解求解方式和等效荷载法,将累积变形柱进行了等效简化,对于古建筑端部半刚性特性,运用位移和应力边界条件进行处理.研究结果表明:柱跨中侧向挠曲变形对截面压应力影响大,而端部侧移影响较小;半刚性节点的非线性对截面压应力影响显著.通过Ansys数值计算与理论公式进行了对比,讨论了累积变形对木柱屈服荷载的影响,解析解与数值解结果吻合良好.研究结果对长期竖向荷载作用下有累计残损变形的古建木柱力学状态评估有一定的借鉴意义.
The ancient wood structure subject to accumulative deformation under long-term vertical load. Aiming to solve the problem of plane stress for the ancient building’s wooden column under vertical load, analytical solutions of the stress, displacement component and the yield load are proposed by using elasticity mechanics. The semi-inverse solution is combined with the equivalent load method to get the equivalent simplification of the accumulative deformation column. For the semi-rigid characteristic of the ancient building’s terminal part, the displacement and stress boundary are adjusted to solve it. The results show that the lateral deflection of the column cross section has a great influence on the sectional compressive stress and a smaller influence on the lateral displacement at terminal part. At the same time, the semi-rigid joints’ nonlinearity characteristic affects the sectional compressive stress significantly. The influence of accumulative deformation on the yield load of the wooden column is discussed through comparing the numerical simulation by Ansys with the theoretical formula, which are mutually consistent. This paper has referential significance on the state evaluation of wooden column’s accumulative deformation under long-term vertical load.
木柱是古建筑中主要的受力构件, “ 立木顶千斤” 足以说明其在古建筑中的重要角色.木柱的健康与否关系到结构的寿命及安全, 因此, 木柱受力状态的理论研究具有重要意义.梁和柱理论是力学分析中的经典课题, 在一定程度上, 梁和柱除了受荷类型不同, 理论分析具有相似性.戴瑛等[1]通过简化后的固定边界条件, 求出了均布载荷下两端固定短梁的平面应力解析解.黄德进等[2]通过Airy应力函数法, 研究了均布载荷作用下正交各向异性固支梁的平面应力问题, 并与数值解进行对比.榫卯连接是古建筑木结构特色, 连接特性介于固接和铰接之间, 为半刚性连接.姚侃等[3, 4]通过试验和数值相结合的方法研究了古建梁柱榫卯的连接特性.张雷等[5]推导出榫卯连接古建木梁在均布荷载作用下的应力、位移分量解析解, 并通过实例进行了讨论.国外学者Aristizabal等[6, 7, 8]在2008年推导了基于铁木辛柯理论和二阶效应, 考虑弯曲变形和剪切变形时半刚性连接的对称截面构件在轴力作用下的转角位移方程, 该公式考虑了在构件发生横向变形时轴力的剪切分量的影响, 并运用4个经典例子说明了所推公式的有效性; 在2010年推导了基于欧拉-伯努利理论和二阶效应, 考虑弯曲变形并忽略剪切变形时半刚性连接且具有对称截面的构件在偏心荷载作用下的转角位移方程; 在2011年运用弯矩分配法对静定梁和框架结构进行了稳定性和非线性二阶分析, 该分析过程均考虑了弯曲和剪切变形、轴向荷载和半刚性连接的影响, 并运用3个典型的例子证明了该研究的有效性.
综上, 国内外相关学者对于木梁构件的研究分为在各向同性的情况下引入变形或在各向异性的情况下不考虑变形影响, 而对于柱截面应力和承载力方面的研究未同时考虑正交各向异性和累积变形.
本文作者结合古建木结构特性, 选取西藏地区宫殿内古建木柱为研究对象, 考虑其细长形态及竖向荷载大的特点, 将线性阶段推进到非线性阶段, 考虑端部节点半刚性、端部侧移、跨中挠曲变形作为非线性因素, 研究了端部轴向荷载作用下累积变形木柱的解析解, 给出了累积变形对屈服荷载的影响.此处累积变形是指柱高方向的侧向挠曲变形和端部侧移.
古建木柱在柱端约束比较弱, 弯矩作用较小, 所以假定木柱挠曲形状为半波余弦曲线, 即y1(x)=δ 1cos(π x/L).累积变形木柱计算简图如图1(a)所示.
对于受竖向荷载发生挠曲变形的木柱, 可通过附加一个等效虚荷载来近似考虑挠曲变形的影响, 即将挠曲变形在跨中截面产生的弯矩和近似等效成横向均布荷载单独作用时产生的弯矩, 计算简图如图1(b)所示; 其中:δ 1为跨中挠曲变形值; L为柱高; Δ 0为柱端侧移; P为端部轴向力; Ky为柱端部半刚性刚度值; h为截面长度.
为了求出均布荷载q, 需先求半刚性理想构件在q下的解析解.本文研究的柱为方柱, 截面为矩形, 荷载均匀分布, 且横截面尺寸不随x变化, 所以该问题为平面应力问题, 木材本构方程为
式中:ε x为x向正应变; ε y为y向正应变; γ xy为xy向切应变; Ex为x向弹性模量; Ey为y向弹性模量; μ xy为xy向泊松比; μ yx为yx向泊松比; Gxy为xy面剪切模量; σ x为x向正应力; σ y为y向正应力; τ xy为xy向切应力.
几何方程表达为
应变协调方程为
采用半逆解法求解, 设应力函数为φ (x, y), 联立式(2)和式(3)得到
挤压应力σ y主要是由均布荷载q引起的, 且q不随x的变化而变化, 所以先假设σ y=f(y), 对σ y进行2次积分得
应力边界条件为
结合正应力σ x、σ y、τ xy的对称性, 得
式中:H、K为待定参数; a为代表各向异性的参数, a=-s11/(2s12+s66).
对式(2)两边积分得到位移分量
式中:u为x向位移分量; v为y向位移分量; u0、v0代表刚体运动; ω 为待定参数; J(x)=s11
古建梁柱连接为半刚性连接, 设端部弯矩为M, 其中M=Krotz· θ z, Krotz为端部绕z轴的转动刚度, 根据位移边界得
端部为小边界问题, 运用圣维南定理
式中:N为文献[5]中两端固支时的轴向力.结合式(8)、式(9)和式(10)得
将式(11)参数带到式(7)和式(8)中即可得到横向均布荷载作用下构件应力和位移分量解析解.
用等效荷载法将侧向挠曲在跨中截面产生的弯矩近似等效横向等效均布荷载单独作用时产生的弯矩, 通过施加横向均布荷载近似考虑累积变形影响.
式中b为截面宽度.
端部看成小边界, 在x=L/2处运用圣维南定理得
端部位移边界条件为
结合式(13)和式(14)得到未知参数
将式(15)参数带到式(7)和式(8)中即可得到累积变形影响下木柱应力和位移分量近似解析解.
木材主要是顺纹方向的力学性质起作用, 在累计变形下跨中截面不仅承受力P, 同时跨中还承受最大的截面弯矩, 所以当跨中左侧截面边缘应力首次达到屈服压应力
式中:R1=4LKyh5(
根据上述得到的解析解, 计算一根实例构件.按照《清式营造则例》[9], 木柱截面高度为4.5斗口, 取360 mm; 宽为3.6斗口, 取288 mm; 构件长度取3 600 mm.此时P=600 kN, Ky=105 N/m[10], Krotz=1.8× 105(N· m)/rad[11], Δ 0、δ 1挠度限值取《古建筑木结构维护与加固技术规范》[12]的L/250.红松材的弹性常数参照文献[13]取值, 见表1.
| 表1 红松材的弹性常数 Tab.1 Elastic constants of Korean Pine |
木柱中间截面左侧边缘压应力最大, 得到柱端侧移Δ 0、跨中挠曲变形值δ 1对截面最大压应力σ x 的影响如图2所示.由图2可见, 柱端侧移Δ 0、跨中挠曲变形值δ 1的相互作用均对截面最大压应力产生影响; 当Δ 0为0, δ 1达到限值L/250时, 使σ x为7.28× 106N/m2, 增大了25.52%; 当δ 1为0, Δ 0增大到挠度限值L/250时, 使σ x为6.16× 106N/m2, 仅增加了6.29%, 相比而言, Δ 0影响较δ 1小.
梁柱连接处易受环境影响, 其半刚性程度不断变化, 假设Δ 0=δ 1=L/1000, Ky与P同上, 得到半刚性程度对截面最大压应力σ x的影响, 见图3.
由图3可知, 半刚性对压应力影响为非线性, 且影响显著, 不能忽略.说明半刚性节点对结构节点受力状态评估影响很大, 需在性能评估时进行考虑.
采用Ansys建立模型, 柱的尺寸选用长度h× 宽度b=450 mm× 300 mm, 高度L=4 500 mm.单元选用Solid 45实体单元, 该单元是八节点六面体实体单元, 可定义各向异性材料.端部半刚性连接采用Matrix 27单元, 该单元可定义6个方向的刚度, 通过调整单元刚度矩阵中的参数模拟端部侧移和转动.
木材弹性参数参照文献[14]选取, 藏青杨的弹性常数如表2所示.
| 表2 藏青杨的弹性常数 Tab.2 Elastic constants of Tibetan cathay poplar |
木材具有正交各向异性, 这种特性体现在单元坐标系中.为了正确模拟这种特性, 选用半刚性连接的构件进行建模[16], 施加均布面荷q=1 kN/m2, 验证模型的计算简图如图4所示, 图4中Ka为A端节点刚度值, Kb为B端节点刚度值.理论与数值结果对比如表3所示.
| 表3 理论与数值结果对比 Tab.3 Comparison of theoretical and numerical models |
对比发现, 理论计算结果和数值模拟计算结果非常接近, 误差在1%以内, 证明了模型的正确性.
建立跨中挠曲变形值δ 1为0, L/1000, 2L/1000, 4L/1000的实体有限元模型.结合式(16)得到跨中节点解析解与数值解的荷载-位移曲线.荷载-位移曲线对比图如图5所示.
由图5可知, 弹性阶段侧向挠曲变形对轴向荷载影响较大; 当跨中挠曲变形δ 1为零时, 理想柱跨中位移v较小; 且轴向荷载p一定时, 随着跨中挠曲变形值δ 1的增大, 跨中位移v增大; 解析解与数值解吻合良好.
为了进一步研究累积变形对屈服荷载的影响, 结合式(15)分别得到累积变形对屈服荷载的影响, 跨中侧向挠曲变形δ 1影响下及柱端侧移Δ 0影响下柱屈服荷载pe及荷载曲线分别如表4和表5及图6和图7所示.
| 表4 跨中侧向挠曲变形影响下柱屈服荷载 Tab.4 Column yield load under midspan lateral deflection |
| 表5 柱端侧移影响下柱屈服荷载 Tab.5 Column yield load under lateral displacement at the end of column |
| 图7 柱端侧移影响下柱屈服荷载曲线Fig.7 Curves of column yield load under lateral displacementΔ 0 at the end of column |
由图6、图7可知, 跨中侧向挠曲变形δ 1和柱端侧移Δ 0对木柱屈服荷载pe的影响均呈线性, 但使其下降变化的梯度不同.当跨中侧向挠曲变形δ 1累积变形值达到4L/1000时, 屈服荷载pe降低39%; 当柱端侧移Δ 0累积变形值达到4L/1000时, 屈服荷载pe仅降低2%.因此, 跨中侧向挠曲变形δ 1对木柱屈服荷载pe的影响明显大于端部侧移Δ 0的影响.说明在残损调查中, 相较于端部侧移Δ 0, 更应关注跨中挠曲变形δ 1, 力学评估中也更应关注这一因素, 其对构件承载力也有较大影响, 古建筑的维修保护中应重点关注.
运用弹性力学方法结合半逆解求解方式和等效荷载法, 对累积变形下古建木柱的应力与位移分量及屈服荷载进行了解析求解并验证, 研究可为轴向荷载作用下存在累计残损变形的古建木柱力学状态评估提供依据.得出的主要结论如下.
1)综合考虑柱端侧移值、跨中侧向挠曲变形值以及半刚性节点的影响, 推导出了累积变形古建筑木柱在轴向荷载作用下的解析解.
2)通过算例对比分析验证了解析解的正确性.结果表明跨中侧向挠曲变形对截面压应力影响大, 柱端侧移影响较小, 同时, 半刚性节点的非线性对截面压应力影响明显.
3)建立了有限元模型, 进行数值解与解析解的比较, 得出弹性阶段跨中侧向挠曲变形对轴向荷载影响大; 当荷载一定时, 理想柱跨中位移很小, 而随着跨中侧向挠曲变形的增大, 跨中位移很大的结论.
4)跨中侧向挠曲变形对屈服荷载的影响明显大于端部侧移对屈服荷载的影响.当跨中侧向挠曲变形累积值达到L/250时, 屈服荷载降低39%; 而当柱端侧移变形累积值达到L/250时, 屈服荷载仅降低2%.
| [1] |
|
| [2] |
|
| [3] |
|
| [4] |
|
| [5] |
|
| [6] |
|
| [7] |
|
| [8] |
|
| [9] |
|
| [10] |
|
| [11] |
|
| [12] |
|
| [13] |
|
| [14] |
|
| [15] |
|
| [16] |
|

