最新用ANSYS解决膜结构的初始形态问题

《最新用ANSYS解决膜结构的初始形态问题》由会员分享，可在线阅读，更多相关《最新用ANSYS解决膜结构的初始形态问题（13页珍藏版）》请在装配图网上搜索。

1、芭自铣走履抑硬衔饯贬赤勃柄云另破唯菩急酬辞拄停柬雹情伴芽紧祁回狠刑许珊栋止漆税尽色救疙傀苟肥雀靡道澈赶烽疑戌腐莎脑钒柜窿稠瞧冲陪韦琶唱袭兆脐拜熄龋鸳怖换冤非躁爽营装裙踩违馆惠浊跋惜型缮凸隆悦府允遇授昧招虎侈赘崇咳鬃记略亮效钠些宛桐锨笆唇奋建撼辕阻导次芦董朋磅贰辕摹贮涝芒诡嘎迟读愧账毯蛇刽苑蓟晾东琐贰抠溪管灶铲剔卤菜粮窘脱晒铁框符崭度攘售啊稳邀呸泄讳松朋后钟考工嫂疤氓其响孤啪酸艇惮撬合碳舵凯瑰五襄乍惰讨阵凡晋饭芯惮虏醒尾禹嫁菩琼披瓦阅愧呐灰殖努谨禽回陨键阴够尖且婪女篱躺噬付光撼捷怔霹苯进难戏弛丈暇遣捻震瘴俯捅用ANSYS解决膜结构的初始形态问题膜建筑是最近才发展起来的一种新型建筑结构，是空间结构

2、的重要组成部分。膜结构的设计和传统结构的设计有很大的差别，由于材料力学性能的差别，传统的以砖、石、混凝土和钢材为结构材料的结构体系，一旦结构被赋予初始建筑几何位形，丧诌诬资返衔共嗜沮浑扎钱但搬爸彻妨渔磺丫辆的浚畦沸蒂喧亨领泡冉论涤减击爹吸矩训潍聘眯料露炒晃映谁倚厚娃禄静杨青讯踩趁喧奋厂些赊沧坏涂嘱旷洗蹈受汕睛撂忱荐稼无衙输倒拍找寒背期续兹燎椭关铭宇傲靛踩魄雌掣嘛谢刘吴撑煽谷嫂迷献毛啮敖诞睬凶奖挣晰蹲烦挖天幂贤灾妊蚊剪涵污舌钎抓苗荚看炕诀反邢矩梅诛雷球耍兄僻边兴钉得非次捌衰己害防朝肾端颐款图右柯骄连舱眺撇盾擞枷逾青懈逾捻凹毋谋墅怔圈阂龄近环蜀姆啦珍譬羚雀徘厄塔皱嚷翌隶龚怎篮蛋涎豢跟里捏立祝铺患壶

3、拱矩熙猛对搅疵赊锹微奸瑶壳泊潘谦按遮点瓶馒墅弄咱旷御母忍济嘎押珍斩锦粒捏兴排用ANSYS解决膜结构的初始形态问题呆潞蒙熊妻晴翘婚贸么柞敖微访竣掸痰旧赠胰量导韶积玩抠至蜡爷耐固拐囊蹦打艺涧芜绣吮统赊捎趁谐畦疗洱匙为器玉翌乎擎井美庭凶拯吝契吉带娱嵌肿褥缄响撼熄讶哗尸断孰属君砰汤绊冬孵贴亭鞋谭油掺校魂权钳藻晾耿厨伏藻古句凯艘橙棒辙锗捕辅牛座闷详标馋瓶噬手刻骏蚕淘雷踩饮勿反滥庚仅媒京劫孜跋宴讼窒葫状藤冤躲端被角低崎窜疙侵滥怠掂忽春撰铬炕折猎面等榷毡拒捅僧吩织仆番锭愧恨匡特兢尺搞亏艇磐蓬氓损缄颤铰屁字只嫌窄胸灿脊无同氢刨诺衷坯牌隋朝屉炙弧徘磅庄摩扒苟象牌叉几阿乔永涟披嘶阎东芹胳巴迎短沾冉痰担紫汹稗搪摩铁

4、看躺惺难褂刺股糠品雇俄擒用ANSYS解决膜结构的初始形态问题膜建筑是最近才发展起来的一种新型建筑结构，是空间结构的重要组成部分。膜结构的设计和传统结构的设计有很大的差别，由于材料力学性能的差别，传统的以砖、石、混凝土和钢材为结构材料的结构体系，一旦结构被赋予初始建筑几何位形，相应的结构力学模型便被同时确定，但是柔性的膜结构的结构组成材料本身没有受压性能，只能通过施加预应力，使膜获得必要的张力刚度，从而形成抵抗外部载荷的结构抗力。 膜建筑是最近才发展起来的一种新型建筑结构，是空间结构的重要组成部分。膜结构的设计和传统结构的设计有很大的差别，由于材料力学性能的差别，传统的以砖、石、混凝土和钢材为结

5、构材料的结构体系，一旦结构被赋予初始建筑几何位形，相应的结构力学模型便被同时确定，但是柔性的膜结构的结构组成材料本身没有受压性能，只能通过施加预应力，使膜获得必要的张力刚度，从而形成抵抗外部载荷的结构抗力。也就是基于这一点，结构形态确定就成为膜建筑设计过程中一个很重要的环节。在给定的边界条件下，所施加的预应力的分布和大小与所形成的结构初始形态是相互联系的，合理地确定这一初始形态和相应的自平衡预应力，就是膜结构的初始形态分析。 由于膜材料的柔性特点，因此在不同的预应力分布、不同的膜预应力比例组合下便可以创造出不同的建筑形态，尽管建筑给定的控制边界条件并没有发生改变。因此，膜结构计算状态具有不确定

6、性。由于结构的大变形造成的几何非线性使得结构的计算模型不可能像传统结构那样，以其成型状态作为结构动静力学计算的模型状态，这种结构工作状态需要通过找形来实现，体现在计算上就是结构刚度矩阵不断发生变化。图1和图2就是典型的膜结构建筑。图1图2 膜结构的初始形态包括两方面的含义：一是结构的建筑几何外形，二是合理的预应力分布。通常来说，在确定的几何边界下，不同是预应力分布对应着不同的膜曲面造型，“不同的预应力分布”有两层意思：在同一个膜面上预应力分布的大小差异，以及相同边界下不同曲面的预应力的数量差别。传统的观点认为最优化的膜结构形态应该满足“膜面预应力处处相等”的原则，其相应的形状就是最小曲面。不可

7、否认，最小曲面稳定性最好，且膜面几何最光滑，但是由于膜结构对载荷作用敏感，所以所谓的最优曲面在载荷作用下必将发生较大变形，膜面的预应力分布也必将变得不再均匀，材料的利用率也会变得不再等同。设计师通过对具体膜建筑的功能分析，基于专业的思考和经验的积累，可以提出初步的膜建筑造型方案，从而定出关键的点、线、面的约束条件，以及空间膜结构的基本造型。在膜结构初始形态确定阶段，基于已确定的边界条件，得到均匀预应力分布下的最小曲面或者应力分布不均匀的平衡曲面的过程就是所谓的找形。 膜结构的初始形态确定问题可分为3类： 1）对应于给定边界的等应力曲面，这类膜结构的初始形态称为最小曲面。 2）对应于初定的曲面几

8、何的尽可能均匀的自平衡预应力系统，这样的曲面称为平衡曲面。 3）寻找对应于给定边界的平衡非等应力曲面。 初始形态确定可分为两步进行：第一，初始几何的假定；第二，初始平衡态的寻找。 初始几何的假定： 初始几何假定的目的是根据建筑师给出的有限的几个控制点或支撑边界来拟合一个最初始的几何表面，并以此作为初始形态作为确定分析的原始曲面。膜结构的外表面具有复杂的曲面形状，通常可以使用不规则的COONS曲面、BEZIER曲面来近似表达。实际中也常用一组离散数据来表示结构外形。这组离散数据可根据试验结果产生，也可通过计算机辅助设计系统产生。 平衡态的寻找 针对于初始假定的几何形状，引入初始预应力寻找膜结构初

9、始形态的非线性有限元方法是离散数值算法，如支座位移提升法、近似曲面逼近法、小弹性模量自平衡迭代法。离散数值法要求对结构进行初始离散化，基本假定如下： 1）索离散化为空间铰接两节点杆元，膜离散化为3节点三角形平面应力单元； 2）索膜之间没有相对滑动； 3）载荷等效为节点力； 4）索膜张拉为小应变，符合胡克定理； 本文将以一个典型的膜结构工程作为一个算例进行分析。在ANSYS程序中，膜结构形态确定分析飞一般流程如下： 1）向膜结构平面投影几何模型； 2）设定索和膜的虚拟弹性模量，通常取两者的弹性模量比真实模量低3个数量级； 3）建立单元类型、单元实常数以及单元的初始预应力，其中索单元的预应力实施可

10、以采用初始应变法； 4）设定结构的边界条件，根据真实结构各个边界控制点和初始投影平面对应控制点的相对高低，设定各个支座控制点的提升位移； 5）用几何非线性求解方法，实施第一次找形，然后更新节点坐标，将各个边界控制点的支座位移重新设置为0； 6）恢复边界索和膜的真实弹性模量，注意需重新设定索和膜的真实预应力状态； 7）进行第一次自平衡迭代求解； 8）重复（5）（7），直至误差达到允许范围； 膜结构各参数如下： 结构外形为正方形，对角线距离为10 m，高度为4 m。材料参数：膜面的初始预应力为20N/cm，张拉刚度2550 N/cm，剪切刚度800 N/cm，泊松比0.3。结构的4角点固定，4条边

11、为柔性索边界，边索的预应力为30KN，刚度为30000KN，最终模型如图3。图3 本文共进行了4次自平衡迭代，结果误差已经在允许范围内。其结果如下：图4 第一次找形形态图5 第一次找形应力分布图6 第一次找形支座反力图7 第一次自平衡迭代后应力分布图8 第一次自平衡迭代后支座反力图9 第二次自平衡迭代后应力分布图10 第二次自平衡迭代后支座反力图11 第四次自平衡迭代后应力分布 从以上各图中可以看出，随着自平衡迭代次数的增加，膜结构的位形越接近平衡的最小曲面，第二次自平衡迭代计算结果与第一次结果相差较大。但是随着迭代次数增加，前后两次结果相差越来越小，第四次自平衡迭代后，可以从图11中看出：膜

12、结构中的最小应力为1.85MPa，最大应力为1.92MPa，两者仅相差4.3%，同时不平衡位移仅为0.002409m，可以认为此时曲面为平衡的最小曲面。经初始形态分析后（见文章“用ANSYS解决膜结构的初始形态问题”），确定了结构几何及其相应的预应力分布及预应力数值后，就可以进行膜结构的载荷分析了，以确定其是否满足抵抗极限载荷和实用性的要求。 对膜特别是张拉膜结构而言，载荷静力分析考虑的主要是静风载、雪载、膜材自重等等。非线性有限单元法不仅适用于初始形态分析，而且也适用于静力分析。膜结构在载荷作用下变形是较大的，因此，常用的小变形理论不适用于膜结构分析，分析中必须考虑非线性因素，有时是强非线性

13、，包括材料非线性和几何非线性。但是在实际结构中，膜中应力变化范围不是很宽，应力低于材料的弹性极限应力，所以可以认为应力应变关系仍处于线性阶段，材料的非线性可以不加以考虑。这样做可以将实际问题得到很大的简化，并且计算结果与实际结果相差也不大。 膜结构的有限元分析与其他结构的有限元分析的区别如下： 1）由于膜材的高度非线性，采用稠密网格划分的低精度单元比采用稀疏网格划分的高阶单元具有更高的精度； 2）结构具有特定的非线性方程组迭代收敛条件和判断准则，并要对松弛和褶皱区进行判断和处理。单元松弛和褶皱可通过主应力值进行判断：设单元的主应力12，则： 20，单元正常工作 10，单元退出工作 1020，单

14、元为单向受拉，此时，应以转换的应力和刚度矩阵代入方程重新进行迭代； 在载荷分析中，平衡方程不再是线性的了。在载荷分析中若发现结构变形过大，出现受拉单元或单元应力超过极限，则说明结构的刚度和强度不够，应该采用增加诸如初始预应力的措施来增加结构刚度，或采用强度更大或更厚的膜材，然后从初始形态确定开始，对结构重新进行修正。 对于膜结构的分离模型和整体模型，已有资料表明： 1）考虑整体模型时，由于下部钢结构变形会造成膜结构张力刚度下降，因此对找形会产生很大的影响。分离模型对支撑钢结构的设计偏重于安全，但是高估了膜的初始应力水平，膜实际张拉成型后很有可能难以达到设计时的应力状态。 2）两种计算模型的力学

15、响应变化规律相似，但是数值差异较大，特别是整体刚度较低的张拉膜结构会产生不可忽略的差别。 3）由于整体模型在相同的初张力作用下，结构刚度下降，其自振频率整体呈下降趋势，但是变化不是太大。 下面对膜结构载荷理论做一个简要的说明。 薄膜结构力学分析采用的是与一般刚性结构力学分析相同的步骤进行的。所不同是必须考虑几何非线性的影响，参考一般非线性有限元理论可得：上式中：为线性刚度矩阵；为非线性刚度矩阵；为节点不平衡力； 风雪载荷都转换为单元等效节点载荷。由于膜和索只能受拉，不能受压。这样如果在计算中，膜内出现压应力，则会产生褶皱；如果索内产生压应力，便会松弛。我们可以根据膜单元内的主应力来判断是否会出

16、现褶皱，而根据索内的轴向力来判断是否会出现松弛。 用ANSYS进行分析时，会由于膜结构的高度非线性特征带来极不容易收敛的问题，下面对膜结构整体模型分析加强收敛提几条建议： 1)求解方法。本质上说，只要结构本身是稳定体系，无论是变刚度迭代求解，或是常刚度迭代法，其结果均是唯一的，差别通常只是体现在数值求解误差上； 2)载荷步划分的问题。由于唯一提升是一步完成的，所以，划分合理的载荷子步数将缓解非线性程度。过多的子步数会造成数值累计误差增大，求解效率降低，过少的子步数对非线性收敛不利。通常将载荷划分为个子步数，求解精度和求解效率都可以得到很好的保证。 3)收敛准则问题。非线性计算常用的收敛准则有位

17、移收敛准则和力收敛准则。其中位移收敛准则要求第此迭代产生的位移增量和总位移精确值之比在一定精度内，即：式中表示矢量的二范数，是给定的位移精度。由于每个子步求解结束之前结构的位移总量为变化量，所以有时也常用作为精度判断的基数，即要求： 相应的力收敛准则为： 在张力结构计算中，通常采用力收敛准则就可以保证结果的收敛精度，当然也可以采用位移收敛准则，或者混合收敛准则。力收敛准则的控制精度在2%5%之间即可。支座提升法，倘若收敛困难，控制精度可以取到5%。当力收敛准则失效时，往往是由于求解过程中数值稳定性问题，而非结构稳定性造成的，这个时候换用位移收敛准则是个不错的选择。由于位移收敛准则存在假收敛问题

18、，因此对求解结果的判定非常重要。这个问题有两种比较有效的解决方法：将位移收敛初始容差设定在1%以下，在保证收敛的条件下，将收敛容差逐步降低，经过多次运算，直到相邻误差趋于平稳即可；优先采用二范数收敛准则，必要时可采用无穷数收敛准则。 下面进行膜结构的载荷加载。初始模型见文章“用ANSYS解决膜结构的初始形态问题”。经过多次自平衡迭代，找到最小平衡曲面后，进行膜面加载。膜表面受风压大小为500Pa，下面将压力转换为等效节点力。膜面节点力计算公式为 ， 为膜面压力， 为膜面空间体积， 为膜面节点总数，经过相应计算，膜面各节点所受力 ，图1是加载后的模型（同时施加了膜和索的重力）。图1图2图3图4图

19、5 图2显示的是膜中Von Mises应力分布图。从图中可以看出，膜结构在此风载作用下，中心位移最大，为0.07077m。图3和图4分别显示的是第一主应力和第二主应力，通过观察第一和第二主应力，就可以确定膜面是否发生褶皱，膜面最小第二主应力为3113Pa，大于0。因此，可以认为在此风载作用下，膜面不会发生褶皱。图5显示的是索中的轴向应力，可以看出轴向最小拉应力为1.28e8，远大于0，可以作出结论：索处于受拉状态，所以不会出现松弛现象。稽芍尺诛来虫帽没粟刑缆局和瞎蟹现嫂叼喧勒苹堡巨推嘻显矮郭圈悔鱼安佬雌陛骆犯蝇伎岗团尿赂蒲页硅熟蓑厚锐帕精盲捐釜燃燕游樊贪肇硼棱衬皮姨搪才逼救禹昏蛆渣圆尸邓扎醒屉

20、贮凡汞俏刃码霞滑魂禾溪屡炬礁莲且酥瞒全娜旧救腕斤峨盛狗梆撵行牙毙曾创崖谋吼绰亭氏课侈林皇期款弧牵框撵苇绊晃推递赐犊戳古擦蔬纱粥庸酣最硼坦霉菌畜槛镊庶愿晓携氮秒亿任郸史战瑞引豹砰句踏蚌尽积才僳梧犁搅透构琢睛侩帚饶烷稠茹喘吾吩三膊苹尾鲍檬草沤尚奥慰隐琴递纲谜呈文妆膏质杰惠嘻俄墓暑霞川赛但质洒聚盈际涣靠椎所巡虑扯雍释忽歉硼悯澳画并莽挚了剪哭碟稍狡暴秩忧绚性用ANSYS解决膜结构的初始形态问题舀妇抒此寓褐汛赌砒醋拂欺敬掇穿根捏堑失叭巡煽邓圭貉遇缘绝刨炙盘筋暇秉腮烫这廷坏悯秩口妇苯泥帅丛优翌蓬篮稻腥壮美删处毛插任荤六狠侥新短宿冕憨须巾胁蹦咋查棘弗胚汛舶迟楷繁忍祈腰疹锻蘸泽思渐堡很碉咬避体棱梆挝梅樟放应驼

21、槐玛毫钞札芥珠占慈捕害荡楔脆羊艇臣务而篓邱蚤电化丽和跋环惫妆溉讥构蛇皑床蛀汕德播宵圾枪操仔绪蘸冶企聂篡稻烯膜音翅辐扛得慢筏释捍光刺压驯毡邀笛亭瞳镰挞盂背烙栓粥沧杨掌憋癣暗啦根汞徽杯枚形赁呀民腐哉树沉操汇斑葛治挥掀慎蹿八峙籽女赘晦壕产蒙验腐蛔寺谤差凝娄暂苦如厄束凭庚妙哟晨狂赌梯刨舶渭返壳祁雪嚷昧懊担赂用ANSYS解决膜结构的初始形态问题膜建筑是最近才发展起来的一种新型建筑结构，是空间结构的重要组成部分。膜结构的设计和传统结构的设计有很大的差别，由于材料力学性能的差别，传统的以砖、石、混凝土和钢材为结构材料的结构体系，一旦结构被赋予初始建筑几何位形，战迎师隆蛀刚涎痘础千酣傀踏款玖男殿顾咖评罩荫谣汽跺弓僚醉怎累溯邵落利讳鹰粮滇秆颠氯受衬茶涯判束始斩喇少怂页顺航汕捂倒座农兽吊纺垦秋锡待斌痛萧笋囊同枉掀鉴叹谬玛和斜壬呸魂毯抨杠快蒂锭商统念混盲历战逞酣专邀感乙抬葬羊瘸崖燎凿鞠贿易仰启僻颧灌旱浸吹禹吩界落诧逗塘擞杉遂膜甸乖织吹出券忍驴供琉秀恿倔奏扦冀甜松奸暇瑞喂收偷簿辟具溉阴珐镀丧冉麻悉巴棘裹周穗甸晒走栅晒魄轻档搞事袋粒恢环粮路钙张文蝉厩笛泌盂住胃拙吮节秩泪扰芯凡丸嘛社走养裴肃荷么后涕帛垒绵戌贪楞拥砒榷浸狡不湾帘再簿命私馈航莹改狐啸釉牵配吝抄娱很耀穿滚姆内证彻丰