一柱式藏式古建筑木构架结构静力分析.pdf
第3 4卷 第6期 2 0 1 5年1 2月 兰州交通大学学报 J o u r n a l o f L a n z h o u J i a o t o n g U n i v e r s i t y V o l . 3 4N o . 6 D e c . 2 0 1 5 收稿日期2 0 1 5-0 3-1 8 学报网址 h t t p / /x b.l z j t u. e d u. c n 作者简介 靳文强(1 9 7 0-) , 男, 山西襄汾人, 副教授, 主要研究方向为结构工程. E-m a i l j i n w q @m a i l .l z j t u. c n 文章编号1 0 0 1-4 3 7 3( 2 0 1 5)0 6-0 0 4 5-0 4D O I1 0. 3 9 6 9/ j . i s s n . 1 0 0 1-4 3 7 3. 2 0 1 5. 0 6. 0 0 9 一柱式藏式古建筑木构架结构静力分析 靳文强 ( 兰州交通大学 土木工程学院, 甘肃 兰州 7 3 0 0 7 0) 摘 要 为研究藏式古建筑中较为普遍的木梁柱构架的力学性能, 利用非线性有限元方法, 建立梁柱构架的实体模 型.在上层集中荷载和楼面均布荷载作用下, 利用工具软件得到构架应力、 应变及剪力图形.分析结果表明 在给定 荷载作用下, 梁端约束部位的应力最大, 柱底的反方向应力最大; 梁和弓木均发生受弯变形, 且方向相反; 组件之间 通过接触面挤压、 摩擦传递荷载.结果发现垫木作为横纹承压构件最先破坏, 同时这种连接方式极易发生侧移和扭 转. 关键词 藏式古建筑;M i d a s; 有限元; 木梁柱构架; 静力分析 中图分类号TU 3 6 6. 2 文献标志码 A T h e S t a t i c A n a l y s i s o f a P i l l a r o f A n c i e n t T i b e t a n A r c h i t e c t u r a l W o o d F r a m e S t r u c t u r e J I N W e n-q i a n g (S c h o o l o f C i v i l E n g i n e e r i n g,L a n z h o u J i a o t o n g U n i v e r s i t y,L a n z h o u 7 3 0 0 7 0,C h i n a) A b s t r a c tI n o r d e r t o s t u d y t h e m e c h a n i c a l p r o p e r t i e s o f a k i n d o f p o p u l a r l y u s e d w o o d e n b e a m- c o l u m n f r a m e o f T i b e t a n a n c i e n t a r c h i t e c t u r e,t h e s o l i d m o d e l o f b e a m-c o l u m n f r a m e i s b u i l t b y u- s i n g n o n l i n e a r f i n i t e e l e m e n t m e t h o d .T h e s t r e s s,s t r a i n a n d s h e a r d i a g r a m o f t h e f r a m e w o r k i s o b t a i n e d b y t h e t o o l s o f t w a r e u n d e r t h e c o n c e n t r a t e d l o a d a n d u n i f o r m l y d i s t r i b u t e d l o a d o n t h e u p p e r f l o o r . T h e a n a l y s i s r e s u l t s s h o w t h a t i n t h e r o l e o f a g i v e n l o a d,t h e m a x i m u m s t r e s s a p- p e a r s a t t h e b i n d i n g s i t e o f t h e b e a m a n d t h e m a x i m u m s t r e s s o f o p p o s i t e d i r e c t i o n a p p e a r s a t t h e b o t t o m o f t h e c o l u m n;t h e b e n d i n g d e f o r m a t i o n o c c u r s o n b o t h t h e b e a m a n d b o w w o o d,a n d t h e i r d i r e c t i o n s a r e o p p o s i t e;t h e l o a d t r a n s f e r s b y e x t r u s i o n a n d f r i c t i o n b e t w e e n t h e c o n t a c t s u r f a c e o f t h e c o m p o n e n t . T h e r e s u l t s s h o w t h a t t h e p e d e s t a l f i r s t l y f a i l s a s a c o m p r e s s i o n m e m b e r, a n d t h i s k i n d o f c o n n e c t i o n i s e x t r e m e l y p r o n e t o l a t e r a l d i s p l a c e m e n t a n d t o r s i o n. K e y w o r d sT i b e t a n a n c i e n t a r c h i t e c t u r e;M I D A S;f i n i t e e l e m e n t;w o o d e n b e a m-c o l u m n f r a m e; s t a t i c a n a l y s i s 藏式建筑破坏日益严重, 分析构架在荷载作用 下破坏机理是藏式古建木结构加固的关键.姜怀英 等[ 1]总结了布达拉宫维修前的破坏情况, 从受力的 角度对木结构连接的优点做了评价, 并且给出了典 型梁柱节点的尺寸详图. 李鹏等[ 2]引入空间弹簧单 元, 得到藏式古建筑梁柱节点连接的简化模型. 曲扎 江措等[ 3]指出藏式建筑的柱网最基本的单元是一柱 式, 其他布置形式均由一柱式发展而来.本文运用有 限元软件 M i d a s对藏式古建筑木构架中的一柱式 结构形式的梁柱框架进行研究, 在静力作用下, 对其 进行竖向位移、 应力、 应变进行分析, 为藏式古木结 构物的维修和保护提供一些可靠的科学依据. 兰州交通大学学报第 3 4卷 1 木构架相关尺寸 一柱式藏式古建筑木构架结构单元从下往上依 次由柱、 斗、 垫木、 弓木、 梁构成, 各构件垂直叠压, 采 用暗销连接( 见图1) , 梁端结合处采用燕尾榫卯增 强联系( 见图2). 参照某典型藏式古建节点尺寸 [4], 模型梁、 柱等的几何尺寸如表1所示. 图1 藏式古建木构架节点详图 F i g . 1 T h e d e t a i l e d d i a g r a m o f t h e n o d e o f T i b e t a n a n c i e n t w o o d e n f r a m e 图2 梁左、 梁右燕尾榫尺寸( 单位mm) F i g . 2 D i m e n s i o n o f d o v e t a i l(u n i tmm) 表1 各构件尺寸( 单位mm) T a b . 1 D i m e n s i o n s o f e v e r y p a r t o f t h e s t r u c t u r e(u n i tmm) 构件名称顶面(LB)底面(LB)高(H) 梁左 3 7 0 01 8 0 3 7 0 01 8 0 4 0 0 梁右3 7 0 01 8 0 3 7 0 01 8 0 4 0 0 弓木 2 2 6 01 5 0 1 7 0 01 5 0 2 8 0 垫木 4 6 01 5 0 2 2 01 5 0 1 3 0 斗 2 1 02 1 0 2 0 02 0 0 1 5 0 柱 2 0 02 0 0 2 5 02 5 0 2 1 1 0 暗销6 04 0 6 04 0 1 0 0 木构架中梁、 弓木、 垫木、 斗为顺纹横向放置, 而 柱是竖向顺纹放置. 根据 木结构设计手册 中规 定[ 5], 木材( 藏青杨) 为正交各向异性材料, 密度取ρ =0. 3 5 2g /c m 3. 弹性模量E、 剪切模量G 和泊松比 等材料常数如表2所示. 2 有限元模型的建立 在建立模型时, 为了简化分析过程, 对其进行荷 载及约束的简化. 两个梁端、 柱底的约束根据其搁置 方式, 分别简化为竖向铰接和三向铰接[ 6] ( 见图3 ). 梁端燕尾榫增强联系, 其模型的尺寸如图2所示, 其 它各构件垂直叠压, 结构间由暗销联系. 表2 木材材性参数( 单位MP a) T a b . 2 P a r a m e t e r s o f w o o d p r o p e r t y(U n i tMP a) Ex/MP a Ey/MP a Ez/MP a E x y /MP a E y z /MP a Ex z/MP a μ x y μ y z μx z 7 2 6 0 7 2 6 3 6 3 5 4 5 1 3 1 4 3 6 0. 2 0. 4 3 0. 4 7 图3 计算模型图示 F i g . 3 S c h e m a t i c d i a g r a m o f t h e a n a l y s i s m o d e l 根据藏式古建筑传力机理[ 6], 对楼板不考虑建模, 只 是将楼板和椽子的重力换算为均布荷载q作用在梁 上, 由于该结构为多层建筑, 则一层柱构架还会受到 二层柱的集中力F作用, 为防止应力集中, 此处将 该集中力换算成均布荷载作用在节点正上方的一定 区域, 其计算模型如图3所示. 图4 有限元模型 F i g . 4 S i m p l i f i e d m o d e l 藏式建筑木构件的连接较为简单, 构件垂直叠 压, 以暗销定位, 暗销有一定约束构件滑移的作用, 这里接触面决定其滑移程度, 在下面接触面的接触 摩擦系数选择较大的值. 根据表1的尺寸构件实体 简化模型如图4所示. 有限元模型中采用一般接触, 其接触面罚刚度取 1, 渗透度取0 . 1 . 木材之间的摩擦系数在0 . 2 ~ 0 . 5之 64 第6期靳文强 一柱式藏式古建筑木构架结构静力分析 间, 这里接触摩擦系数取0 . 5, 接触面如图5所示. 图5 构件接触对 F i g . 5 P i c t u r e o f c o n t a c t p a i r s 图6 网格划分图 F i g . 6 M e s h g e n e r a t i o n m a p 由于各个构件特性不一样, 为了耦合, 采用布尔 运算差集来对网格进行尺寸控制. 这里楼板和椽子 等效均布荷载取q=7 2k N/m, 上层柱传下来的集 中力F=6 0k N, 均换算成面荷载[ 7], 网格划分图约 束及其荷载如图6所示. 3 木构架静力作用下的弹塑性分析 3. 1 木结构的本构关系 木材顺纹静态压缩下的应力应变关系, 可简 化为图7所示曲线; 对于木材横纹方向静态压缩下 的应力应变关系, 可简化为如图8所示的曲线, 参 数的意义见文献[ 8]. 图7 木材顺纹静压的简化本构关系 F i g . 7 T h e s i m p l i f i e d c o n s t i t u t i v e r e l a t i o n o f t h e w o o d u n d e r l o n g i t u d i n a l s t a t i c p r e s s u r e 3. 2 木构架有限元弹塑性计算结果分析 通过对梁上的均布荷载逐级加载至q=7 2k N/ m, 集中荷载保持F=6 0k N, 得到藏式古建筑结构 静力作用下模型的位移、 应变、 应力图如图9~图1 2 所示, 模型的主应力分布图如图1 3~图1 5所示. 图8 简化木材横纹静压的本构关系 F i g . 8 T h e s i m p l i f i e d c o n s t i t u t i v e r e l a t i o n o f t h e w o o d u n d e r l a t e r a l s t a t i c p r e s s u r e 图9 模型的竖向位移图( Y方向) ( 单位mm) F i g . 9 D i s p l a c e m e n t m a p o f v e r t i c a l d i r e c t i o n ( Y d i r e c t i o n)(u n i tmm) 图1 0 模型的水平向位移图(X方向) ( 单位mm) F i g . 1 0 D i s p l a c e m e n t m a p o f h o r i z o n t a l d i r e c t i o n (Xd i r e c t i o n)(u n i tmm) 通过图9可以看出在荷载作用下, 结构的竖向最 大位移发生在梁的跨中. 梁的变形是上部受压, 下部 受拉, 而弓木则是上部受拉, 下部受压. 从图1 0~图 1 1可以看出, 构件靠摩擦和相互挤压协同工作, 同时 也说明这种连接方式容易发生侧移和扭转. 从图1 2 可以看出, 在垫木的右下方的部位剪力最大, 出现应 力集中现象. 垫木一般是横纹放置, 此处可知, 垫木是 横纹受压破坏, 也是整个结构中最先破坏的构件. 74 兰州交通大学学报第 3 4卷 图1 1 模型的Z方向位移图( 单位mm) F i g . 1 1 D i s p l a c e m e n t m a p o f Zd i r e c t i o n( u n i tmm) 图1 2 模型最大剪力图( 单位k N) F i g . 1 2 D i a g r a m o f m a x i m u m s h e a r f o r c e o f t h e m o d e l(u n i tk N) 图1 3 模型的第一主应力值( 单位MP a) F i g . 1 3 V a l u e o f t h e f i r s t p r i n c i p a l s t r e s s(u n i tMP a) 图1 4 模型的第二主应力值( 单位MP a) F i g . 1 4 V a l u e o f t h e s e c o n d p r i n c i p a l s t r e s s(u n i tMP a) 图1 5 模型的第三主应力值( 单位MP a) F i g . 1 5 V a l u e o f t h e t h i r d p r i n c i p a l s t r e s s(u n i tMP a) 由图1 3 ~ 图1 5显示模型中的3 D单元( 梁柱) 的 P1、P2、P3值的最大正应力、 最小负应力的绝对值, 可 以看出应力最大处为梁端约束部位和柱底. 栌斗、 弓 木以及垫木所在的部位, 中间部分的竖向应力相对较 大, 说明这些构件将上部荷载实现了向下传递. 4 小结 以一柱式藏式古建筑木构架结构为研究对象, 使用 M i d a s F E A 中的接触分析基本理论, 对模型的 实体尺寸、 材料属性、 接触方式、 网格划分、 荷载的加 载方式及边界条件的定义等问题进行了说明. 梁上施加均布荷载同时在柱上施加集中荷载的 情况下, 结构的竖向最大位移发生在梁的跨中. 荷载 足够大时, 梁和弓木均会发生受弯破坏, 而且两者受 弯方向相反. 垫木是所有构件中最先破坏的构件, 发 生的是横纹承压破坏. 两个梁端约束部位和柱底应 力最大. 同时发现这种连接方式极易发生侧移和扭 转, 在构件设计和加固时应充分考虑这些因素. 参考文献 [1] 姜怀英.西藏布达拉宫修缮工程报告[R]. 北京 文物 出版社,1 9 9 4 2 8-3 2. [2] 李鹏. 藏式古建筑木梁柱节点力学性能研究[J] . 土木 工程学报,2 0 1 0, 4 3(Z 1) 2 6 3-2 6 8. [3] 曲扎江措, 泽仁玉珍.藏式古建筑的结构特点及其防雷 措施[J].西藏研究, 2 0 1 2(3) 1 0 5-1 0 9. [4] 张岩.藏式结构有限元模型修正与变形关系分析[D]. 北京 北京交通大学,2 0 1 0. [5] 中华人民共和国建设部. 木结构设计手册[M]. 3版. 北京 中国建筑工业出版社,2 0 0 5. [6] 杨娜, 闫会春. 典型藏式古建木构梁柱连接传力机理分 析[J]. 北京交通大学学报, 2 0 1 2,3 6(4) 1 1 2-1 1 9. [7] 郭婷.藏式古建结构的人致荷载分析[D].北京 北京 交通大学,2 0 0 9. [8] 陶俊林.木材静压大变形本构关系研究[J].力学与实 践,2 0 0 0, 2 2(5) 2 5-2 7. ( 责任编辑 李爱军) 84