单层网壳结构随机缺陷模态法研究pdf(共7页)

1、第 23卷 第 3期河北 建筑科 技 学院 学报Vol 23 No 32006年 9 月Journalof HebeiInstituteofArchitecturalScienceand TechnologySep 2006文章编号: 1007-6743( 2006) 03- 0011- 04单层网壳结构(jigu)随机缺陷模态法研究范志飞, 曾国华(u hu), 董 聪( 清华大学土木工程系, 北京(bi jn) 100084)摘要: 随机缺陷模态法是分析网壳结构缺陷的常用方法, 传统的方法以节点位移为随机变量进 行缺陷结构的稳定性分析, 为了使理论上更加合理及能够考虑初应力的影响, 本文选

2、取结构的 基本参数( 杆件长度、钢管外径、弹性模量、钢管厚度) 为基本的随机变量。根据这种思想, 通过 数值计算发现了杆件长度对结构稳定承载力的影响最为显著, 拟合了网壳结构极限承载力的概 率分布, 确定了规范规定的可靠度指标下的临界荷载标准值。关键词: 网壳结构; 随机缺陷; 初应力中图分类号: TU356文献标识码: A单层网壳的稳定性验算已成为结构设计中的 关键问题, 由于单层网壳是缺陷敏感性结构, 初始 缺陷对网壳结构稳定性的影响仍然是目前学术界 和工程界关注的重要课题。实际的网壳结构不可避免地具有各种初始缺陷, 包括杆件的初弯曲、对节点的初偏心、节点的 位置偏差和残余应力等。缺陷可以

3、分为两类: 第 一类为节点缺陷, 是指在实际结构中, 节点约束是弹性约束而非理想的无摩擦铰, 从而引起力不通 过刚心而形成偏心荷载, 影响结构的刚度和最大承载力; 第二类为几何缺陷, 指部分( 或全部) 结构 与原设计有偏差, 将改变网壳的初始几何形状从 而大大影响网壳的屈曲行为。研究指出, 将模拟缺陷结构的屈曲荷载与原理想结构相比, 第二类 缺陷的变化更为敏感, 而第一类缺陷对其影响不 大。缺陷分析的典型方法有: 一致缺陷模态法和 随机缺陷模态法。本文将对随机缺陷模态法作进 一步研究。1 随机缺陷模态法的新思路1. 1 初应力对网壳结构稳定性的影响确定结构的临界荷载就是寻找使结构几何非线性方

4、程的切线刚度矩阵成为奇异时的荷载临界 值 pcr 。对于U . L. 列式, 即求解以下特征值问题。t K = 0( 1)其中 t K = t KL + t KNL( 2)t KL =tBLT tD t BLtdV( 3)etVt KNL =lBNLTSl BNLldV( 4)elV式中 ! 对应(duyng)于临界值p cr 的时刻(shk); t BL 和t BNL ! 线 性应变(yngbin)和非线性应变与位移的转换矩 阵; t S ! Cauchy 应力矩阵, 所有这些矩阵和向量的元素都是对应于 时刻, 并且参考于 t 位形。 从式( 4) 中可以看到结构的刚度矩阵和应力状态有关。因

5、此, 结构的稳定临界荷载和结构内部 的初应力密切相关, 在计算缺陷对网壳结构稳定 性的影响时不能忽略结构的初应力效应。1. 2 随机变量的确定随机缺陷模态法把节点的几何偏差作为基本 的随机变量对结构进行分析, 抓住了影响承载力 的主要缺陷因素, 获得了很好的效果。然而我们也 注意到: 随机缺陷模态法所采用的节点的几何偏 差严格说不应作为一种缺陷因素, 而是缺陷的反 映。缺陷的影响是初应力引起的, 最后量测到的几 何构形是初应力作用下已经自平衡后的几何构 形。初始缺陷影响的研究应从缺陷的实质上加以 分析: 首先要明确网壳结构缺陷是指网壳结构最 基本的几何参数的变异, 并且这些最基本的几何 参数要

6、保证彼此之间独立无关; 其次要明确缺陷 的产生主要是由于杆件的制造误差和装配误差引 起的, 制造误差包括杆件长度、钢管壁厚、钢管外 径和弹性模量等基本的几何参数产生的变异。这收稿日期: 2006- 05- 09作者简介: 范志飞( 1980- ) , 男, 江西丰城人, 硕士研究生, 从事结构工程研究。12河 北 建 筑 科 技 学 院 学 报2006年样, 我们(w men)就确定了以杆件长度 l 、钢管(gnggun)壁厚 t 、钢管(gnggun) 外径 D 和弹性模量E 为基本的随机变量生成随机 缺陷, 分析结构考虑随机缺陷影响的稳定极限承 载力。在这些基本变量中, 由于网壳结构是多次

7、超 静定结构, 杆件长度 l i 的扰动不仅改变结构的几 何形状而且还会使网壳结构产生初应力, 这种初 应力必然使结构的极限承载力有所降低。根据以 上分析, 下面的分析应着重考察杆件长度 li 的扰 动对结构承载力产生的影响。2 数值计算结果K 6 型网壳的网格布置如图 1 所示, 跨度为18m, 矢高为 1. 15m, 网壳的杆件为 180 5 的钢 管, A = 27. 49cm2, Ix = 21. 6cm4, Iy = 1053cm4, Iz =1053cm4, E = 210GPa。网壳承受的荷载转化为竖向 节点集中荷载, 周边节点为刚接支承。有限元计算模型中, 杆件均采用两节点的

8、Timoshenko 单元, 每根杆件划分为三个单元, 考虑 几何非线性, 采用大变形、小应变假定。非线性求 解采用基于Newton- Raphson 法的弧长法, 采用力 收敛准则。2. 1 初应力对网壳结构稳定性的影响为考察初应力对网壳结构稳定性的影响, 现针 对杆件长度 l 以均值为 !( !为设计值) , 变异系数为 = 0. 002 进行随机扰动, 分别比较不考虑初应力影 响和考虑初应力影响的稳定承载力的区别。可以看到针对本例, 考虑初应力影响的极限 荷载 F = 350. 9kN, 不考虑初应力影响的极限荷载 F = 364. 7kN, 前者比后者降低 3. 8% , 这和前面的

9、分析是一致的。因此, 网壳结构的稳定分析必须要考虑初应力(yngl)的影响。2. 2 网壳结构(jigu)基本参数扰动对稳定性的影响在缺乏各参数(cnsh)分布的统计资料前提下, 为了比 较各参数对网壳结构稳定性的影响, 假定各参数服 从均值为 !( !为设计值) , 变异系数为 = 0. 002 的 概率分布。利用ANSYS 的参数化设计语言APDL, 对 随机变量按照概率分布进行抽样, 实现变量的扰 动。然后对结构进行非线性屈曲分析, 得到网壳结 构的屈曲临界荷载。对于随机变量为杆件长度 l 时, 由于杆件长度的扰动使杆件产生初应力, 所以还须 把初应力读入杆件元素中, 重复这样的过程进行

10、 m ( m 为抽样次数) 次, 即可以求出屈曲临界荷载的 经验分布函数。这样, 我们就可以根据可靠性指标 确定网壳结构极限承载力的设计值。取抽样次数 m = 1000, 由于数据量较大, 仅出 示临界荷载数据的直方图( 横坐标为节点荷载的 临界值, 纵坐标为临界荷载出现在该区间的次 数) , 如图 2、图 3、图 4、图 5 所示。从承载力的变异系数可以得知, 基本参数对网 壳结构承载力影响的程度从大到小的顺序是: 杆件 长度 l 、钢管外径 D 、弹性模量 E 和钢管壁厚 t 。3 临界荷载分布规律的研究对于临界荷载的计算结果采用假设检验的方法来分析临界荷载的分布规律。基于子样分组的#2

11、检验法在工程界得到广泛的应用, 但检验结果 依赖于分组方式, 子样分组方式不同可能导致不 同的结论。近年来发展起来的以经验分布函数和 理论分布函数之差为基础构造的EDF 统计量比 #2 统计量更有效。3. 1 统计量的构造与算法设随机变量子量容量为的随机样本为 x 1, x 2,#, xn , 定义 x ( 1) , x ( 2) , #, x ( n) 为按升阶排列的顺 序统计量。进一步假设 X 的理论分布函数为F( x , ) , 其中 F ( x , ) 为连续函数, 为参数矢量。X 的 经验分布函数 Fn( x , ) 可由式( 4) 确定。表 1各参数对结构临界荷载的影响参数极限荷载

12、的均值 !Pcr ( kN)极限荷载的方差 %Pcr ( kN)极限荷载的变异系数Pcr = !Pcr / %Pcr杆件长度 l360. 218. 950. 025钢管外径 D360. 631. 910. 0053钢管壁厚 t360. 690. 690. 0019弹性模量 E368. 860. 74130. 002第 3 期范志飞等: 单层网壳结构随机(su j)缺陷模态法研究130( x x1)Fx ( x ) = i/ n( x i x x ( i+ 1) ) i = 1, 2, #, n - 11( x n x )( 5) 当 n % & 时, 可以(ky)证明| Fn( x , ) -

13、 F ( x , ) |以概率(gil)1 趋近于 0。也就是说, 经验分布函数 Fn ( x ,) 是理论分布函数 F( x , ) 的一致估计。 我们将度量 Fn( x , ) 和 F( x , ) 之间差异的一类统计量定义为 EDF 统计量。EDF 统计量通常 划分为两种类型: 上确界型统计量和平方差型统 计量。这里采用平方差型统计量进行假设检验。平方差型统计量目前称为 Cramer - von Mises族统计量, 它以经验分布函数 Fn( x , ) 和理论分 布函数F ( x , ) 之差的均方积分为基础。为提高尾 区形态的检验效率, 在构造统计量时, 通常采用加 权方法。Cram

14、er - Von Mises 族统计量中的 A 2 统计量 可用通式表示为+ &Qn = n Fn( x ) - F( x ) 2 & ( x , ) dF ( x , )- &( 6)其中, & ( x , ) 为权函数。当 & ( x , ) 时, 该统计量(jling)为Cramer - von Mises 统计(tngj)量( W2 统计(tngj)量) 。当& ( x , ) = F( x , ) 1- F( x , ) - 1 时, 该统计 量为 Anderson- Darling 统计量( A 2 统计量) 。可以 证明, 当z 按升阶排列z ( 1) z ( 2) # z (

15、n) , 可导 出 W2 统计量和 A 2 统计量的表达式n222i - 11Wn =i= 1z i -2n+12nn12(2i - 1) lnz i + ln(1- z( n+ 1- i) ) - nA n= -ni= 1( 7)对于参数矢量 未知的情况, 采用经样本容量修正后的统计量 A *和 W*为A*=20. 75+2. 25,A n 1+nn2W* = W2n1 +0. 5( 8)n我们用OSL表示观测显著性水平。OSL 的物理含义是, 对于给定数据, 当原假设分布为真时, 拒 绝原假设分布的最大风险率。也就是说, OSL 在数 值上相当于不拒绝原假设分布的显著性水平的上 限值。OS

16、L 的计算公式为OSL =11+ exp - 0. 079+ 1. 24ln( A * ) + 4. 48A * (9)14河 北建 筑科技学 院学 报2006年表 2EDF 统计量检验结果参数分量极大似然估计结果EDF 统计量不同显著性水平( ) 的检验结论分布模型1n1nA*OSL! =lnx i % =( lnx i- !) 20. 05 0. 10 0. 25W *n i= 1n - 1 i= 1X LN ( !, %)5. 88640. 024610.58430. 1332 接受 接受 拒绝0.07973. 2 拟合优度检验结果以上的数值实验的结果, 可以得到以下结论:1)钢管外径、

17、壁厚、弹性模量的扰动对网壳结由于钢管外径、钢管壁厚和弹性模量对承载构承载力的影响远小于杆件长度扰动对网壳结构力的影响相对于杆件长度 l 而言较小, 故只对杆件承载力的影响,在进行缺陷对网壳结构稳定承载长度扰动得到的网壳结构承载力数据进行母体分力计算中可以忽略不计。布拟合优度检验, 检验结果如下。2)杆件长度的扰动不仅使网壳结构位形发生只要显著性水平小于 0. 1332, 对数正态分布变化, 而且使结构产生初应力, 根据数值计算的结的模型是可以接受。根据规范, 假设检验的显著性果, 初应力对网壳结构稳定性产生较大的影响, 在水平可取为0. 05。这样, 我们可以认为选取 != 5.设计中需要予以

18、考虑。8864, % = 0. 02461 对数正态分布作为描述网壳结3)以杆件长度为基本的随机变量进行随机缺构稳定承载力的模型是可以接受的。陷模态法的分析, 从理论上来说更加合理, 得到的网壳结构失稳破坏为脆性破坏,网壳结构安结果更加可信。全等级通常为一级, 根据规范的规定, 结构承载力参考文献:极限状态的可靠度指标 (=4. 2, 对于本题, 根据1 CBORRI C, SPINELI P. Buckling and post - buckling be公式( 11) 可以算得网壳结构承载力的标准值为havior of single reticulated shells effected

19、by random impe re324. 74kN。fection J . Comput Struct, 1988, 30( 4) : 937- 943.Fk =F- 1 ) (- () ( 10) 2 沈世钊, 陈昕. 网壳结构稳定性 M . 北京: 科学出版式中 F !对数正态分布的累积分布函数。社, 1999. 3 唐敢, 尹凌峰, 马军, 等. 单层网壳结构稳定性分析的改4 结论进随机缺陷法 J . 空间结构, 2004,10( 4) : 44- 47.4 王勖成. 有限单元法 M . 北京:清华大学出版社,在整个计算过程中, 基本的随机变量杆件长2003. 5 董聪. 现代结构系统

20、可靠性理论及其应用 M . 北京:度 l 、钢管外径 D 、钢管壁厚 t 、弹性模量 E 变异系科学出版社,2001.数相同, 这样可以保证计算结果具有可比性。根据 6 GB 50068- 2001, 建筑结构可靠度设计统一标准 S .Study of stochastic imperfection method for stability analysis of single layer lattice shellFAN Zhi fei, ZENG Guo hua, DONG Cong( Department of Civil Engineering, Tsinghua University, Beijing 100084, China)Abstract: Stochastic imperfection method is a common method to analyse lattice domes.Custom stochastic imperfection method assumes node displacement as a stochastic variable to study a lattice shell,h