k6型单层球面网壳非线性稳定性分析.docxVIP

k6型单层球面网壳非线性稳定性分析 1 考虑两种非线性方法的对比 网络外壳结构的稳定性是网络外壳结构设计的中心主题。对于仅考虑几何非线性的单层网壳结构的全过程分析方法,文献已经解决得比较好,并在大规模参数分析的基础上提出了可供一般设计部门采用的稳定承载力公式。但由于当时计算机条件的限制,这些分析没有同时考虑材料非线性的影响,其原因主要有两个方面:(1)如果采用考虑几何、物理两种非线性的全过程分析方法,对于当时的计算条件,所需计算时间需增加许多倍,实际上是很困难的;(2)网壳结构的正常工作状态是在弹性范围内,材料非线性对结构全过程曲线及其极限点的影响实际上是使结构稳定承载力的安全储备有所下降。对这种影响已有可能从定量上作出适当判断。当然我们一直认为最可靠和最有效的方法还应是结合不同类型的网壳结构形式,在基本参数(包括几何参数、构造参数、荷载参数等)的常用变化范围内,运用双重非线性分析方法,进行大规模的实际结构计算。现在随着计算机技术的发展以及各种大型通用有限元程序的出现,使得对大型网壳结构进行大规模弹塑性全过程分析成为可能。 基于以上认识,本文利用大型通用有限元软件ANSYS对K6型单层球面网壳结构(见图1)进行1200余例双重非线性全过程分析,得到了较为精确的网壳结构稳定极限承载力、并对所得计算结果进行了统计和归纳,考察了网壳弹塑性稳定的变化规律,比较真实地反映了此类结构的稳定性能,与以往仅考虑几何非线性的分析相比,更具有理论与实用价值。 2 beam189单元梁元理论 本文的有限元分析采用ANSYS软件中的BEAM189单元(见图2).BEAM189单元为3结点二次梁元,基于Timoshenko梁元理论,该梁元在非线性分析中能考虑大变形、大转角和大应变效应。材料的本构关系采用理想弹塑性模型,屈服准则采用Von Mises屈服准则。 2.1 跨国范围划分 从符合实际的角度出发,我们选择如下几种跨度和矢跨比进行分析。跨度:L=40m、50m、60m、70m;矢跨比f/L=1/5、1/6、1/7、1/8。杆件长度通常控制在3～5m范围内,因而K6型网格的分割频数(NF)随跨度的增大而增加,如表1所示。 2.2 网壳的截面尺寸 为了便于比较,杆件一律采用圆钢管。在参数分析中,对每一种跨度的网壳均选用四套不同的截面尺寸,按截面增强顺序排列为①、②、③、④号。每套截面中主肋和环杆采用较大截面,斜杆采用较小的截面。表2列出各种跨度的具体截面尺寸,它们都是工程中常用的规格。 2.3 恒载和活载组合 只考虑满跨均布和半跨均布两种荷载分布。恒载满跨均布,活载可满跨也可半跨均布,但活载(p)和恒载(g)可有几种不同比例:p/g=0、1/4、1/2。 2.4 网壳初始缺陷 在用一致缺陷模态法进行的分析中,以最大偏差等于(L/1000～L/300)的情形作为基本的初始缺陷对所有网壳进行分析,以便对计算结果进行统计;然后对两种跨度的网壳考虑由小到大共9种初始缺陷尺寸进行了专门的初始缺陷影响分析。 2.5 铰接承接情况 本文所研究的K6型网壳结构的支承条件均为固接支承,但对于实际工程中可能出现的铰接支承情形在后面也做了比较分析。 根据上述参数分析方案,我们对K6型单层球面网壳进行了1200余例结构的弹塑性全过程分析。 3 全装配曲线和限制负荷 3.1 网壳弹塑性分析 结构的稳定性能可以从其荷载—位移全过程曲线中得到完整的概念,我们首先在这里讨论一种基本情形,即网壳在满跨均布荷载作用下的全过程分析结果。对每例结构进行全过程分析之后,每个节点都可画出一条荷载-位移曲线。我们所关心的是网壳的极限承载力、相应的最大位移以及塑性发展程度,我们将每例结构计算迭代结束时刻位移最大节点的荷载-位移曲线作为典型代表。本文取第一个临界点处的荷载值作为结构的稳定极限荷载。 为节省篇幅,本文只给出了60m跨度、矢跨比1/7、四种截面尺寸4个网壳在均布荷载作用下的弹塑性全过程曲线,对每例结构在按无缺陷的完整结构进行分析的同时,还考虑了L/1000大小的缺陷进行分析,两组曲线放在同一图中,分别用实线和虚线表示,见图3。首先可以看到这些曲线非常有规律性,对应文献中的弹性网壳结构的荷载-位移曲线可以看出,考虑材料非线性网壳的极限承载力明显降低;考虑初始缺陷影响的荷载-位移曲线在极值点后变化较为平缓,同时临界荷载对应的位移值明显大于理想网壳。对于材料非线性对极限承载力的影响本文做了系统的统计,在本文研究参数范围内,K6型单层球面网壳考虑材料非线性与不考虑材料非线性的极限承载力平均比值Kp列于图4中,此比值最低可达38%。这进一步说明材料非线性对于网壳极限承载力的影响是至关重要的。 需要特别指出的是,与弹性网壳结构不同,K6型球面网壳的弹塑性荷载-位移曲线中没有出现过多的上下波动,也就是