空间网架有限元分析

《空间网架有限元分析》由会员分享，可在线阅读，更多相关《空间网架有限元分析（8页珍藏版）》请在金锄头文库上搜索。

1、论文摘要:空间网架结构是现代大跨度结构工程中最常用的结构形式。本文针对运用ANSYS 进行网空间网架结构的有限元分析，讨论了单元类型，以及材料模型的选择，并且定义材料 的极限应力，必要时甚至考虑压杆失稳状态（屈曲分析）。在各种荷载（永久荷载，风荷载， 地震荷载，自重）的作用下，利用强大有限元分析软件ANSYS对网架结构进行静力学分析， 之后利用ANSYS强大的优化设计功能，在结构安全的情况下，对网架杆件的截面进行优化 设计，然后根据现有的杆件的截面，合理的选择杆件截面，减少材料的消耗量，达到最经济 合理的设计，实现可持续发展，作出令人满意的设计。1前言空间网架结构由于具有轻便、通透的特点，适合

2、大面积、大跨度的使用，现已经成为土 木工程大跨度结构中采用的最多的结构形式。土木工程师在对网架结构进行受力分析时，往 往希望通过采用诸如ANSYS等大型通用有限元分析软件进行计算。考虑到大跨度网架结构 的材料用量，必然希望在结构安全的情况下，材料的实用量能够达到最少，ANSYS的后期 强大的优化设计功能很好地解决了这个问题。在运用ANSYS进行优化设计时应合理选用迭 代次数，在以尽量少的计算量的情况下，得出最优解。2单元类型的选取空间网架结构组成部分主要为杆单元，网架的受力特点是杆件均为铰接，不能承受弯矩 和扭矩，因此所有的杆件只受拉或受压。LINK8三维杆（或桁架）单元是有着广泛的工程应用的

3、杆单元，比如可以用来模拟： 桁架、缆索、连杆、弹簧等等。这种三维杆单元是杆轴方向的拉压单元，每个节点具有三个 自由度：沿节点坐标系X、Y、Z方向的平动。就像在铰接结构中的表现一样，本单元不承 受弯矩。本单元具有塑性、蠕变、膨胀、应力刚化、大变形、大应变等功能。对于LINK8 单元有如下假设：杆单元假定为一直杆，轴向荷载作用在末端，自杆的一端至另一段均为同 一属性。杆长应大于零，即节点I和J不重合。面积也必须比零要大。假定温度沿杆长线性 变化。位移函数暗含着在杆上具有相同的应力。初始应变也被用来计算应力刚度矩阵，即 便是对于第一次累计迭代。Link8单元示意图网架结构中的杆件的受力形式为二力杆，

4、使用Link 8来模拟是最合理的选择。3材料模型的选取网架结构中的材料模型的选取也对整个仿真工程的精确程度有着重大的影响，正确的选 择材料模型，将对整个虚拟工程的可信度起到重大的影响，再次我们选择低碳钢Q235作为 网架结构的材料，为了减少计算量，避免由于材料非线性引起的非线性问题，在此假设Q235 为线弹性材料，并且定义其弹性模量为201Gpa,泊松比为0.3，许用应力210Mpa。4荷载选取在荷载选择方面，在此我们选择一个最常规的荷载模式，就是在永久荷载以及本身自重 的作用下，其中忽略了一些次要的荷载。5 ANSYS优化设计介绍ANSYS优化设计大致可分为以下八个步骤：第一步：生成分析文件

5、分析文件的生成是ANSYS优化设计过程的关键部分。ANSYS程序运用分析文件构造 循环文件,进行循环分析。在分析文件中模型的建立必须参数化的（通常是优化变量为参数）, 结果也必须用参数来提取（用于状态变量和目标函数）。生生分析文件主要有以下几个步骤:（1）参数化建立模型。用设计变量作为参数建立模型的工作是在PREP7中完成的。（2）求解。求解用于定义分析类型和分析选项，施加荷载，指定荷载步，完成有限元 计算。分析中所用到的数据都要指出：凝聚法分析中的主自由度，非线性分析中的收敛准则， 谐波分析中的频率范围等。载荷和边界条件也可以作为设计变量。（3）参数化提取结果。在本步中，提取结果并赋值给相应

6、的参数。这些参数一般为状 态变量和目标函数。提取数据的操作用*GET命令（Utility MenuParametersGet Scalar Data） 实现。通常用POST1来完成本步操作，特别是涉及到数据的存储，加减或其他操作。（4）分析文件准备。第二步：建立优化工程中的参数。这里主要包括一些优化设计变量，优化状态变量以及目标函数的定义第三部：进入OPT，指定分析文件（OPT）第四步：申明优化变量这一步主要指定哪些参数是设计变量，哪些参数是状态变量，哪些参数是目标函数，允 许有不超过60个设计变量和不超过100个状态变量，但目标函数只能有一个。每个设计变 量和状态变量都可以定义最大值和最小值

7、，以及这些变量的公差，目标函数不需要给定范围。第五步：选择优化工具或优化方法。优化方法主要有：零阶方法、一阶方法、用户提供的优化方法、单步运行、随机摸索法、 等步长搜索法、乘子计算法、最优梯度法、用户提供的优化工具。第六步：指定优化循环控制方式每种优化方法和工具都有相应的循环控制参数，比如最大迭代次数等。第七步：进行优化分析在所有的控制选项设定好以后，就可以进行分析了，此时可能需要花去一段时间，如果 结果不收敛，程序可能会中断。在此类情况下，可增加优化模型的假设，调节失代次数等来 实现。第八步：查看优化设计序列结果主要可以看到各种满足条件的不同的解，以及优化设计变量，优化状态量和目标函数的 失

8、代分析曲线。6算例本文以一个受均匀荷载的四角锥平板网架的为例做优化设计，优化网架的截面。在网架 上定义七种不同的截面面积变量，分别为al,a2,a3,a4,a5,a6,a7（如图1所示），网架上平面结构 有9个节点，间距lm,下平面有16个节点间距lm,上下平面相距0.7m,优化设计前各杆 件的截面面积均为5e-5平方米，材料为Q235钢材。在荷载作用下的网架上平面的边缘节点 固定，下平面的每个节点承受垂直方向荷载10000N。并且限制网架在工况下的竖向最大位 移不超过lcm，各杆件的截面应力小于钢材的许用应力。图1网架平面示意图定义模拟网架结构的单元为Link8单元，杆件的截面(AREA)均

9、为5e-5平方米，以及 材料的弹性模量(Ex=2e11pa)和泊松比(Prxy=0.3)。之后就是建立网架有限元模型了，首先 建立平面上的节点，定义单元的属性，在直接由节点生成有限元模型。在网架的上平面边缘 施加约束，限制X，Y，Z三个方向的线位移。在下平面的9个节点上分别施加向下的集中 荷载，大小为10000N。(如图2所示)图2网架空间受力示意图进行完前处理之后，定义分析方式，为静力学分析，然后求解。使用通用后处理模块 （/P0ST1），显示整个模型的体积（如图3所示），体积约为0.35901e-2。A 55UM CommandFileIl F I %.Img e=Hr gA/ a/Ci

10、I | rSUM ALL THE fiCTIUE ENTRIES IN THE ELEMENT TABLETABLE LABEL TOTALEUOLUME St .359098E-02图3网架结构体积显示整个模型的轴应力云图，（如图4所示），其中最大应力为148.363Mpa。图4上述数据均为在优化设计前的一系列数据，现在进入ANSYS优化设计，定义优化设计 变量为AREA （桁架杆件的截面），其中1e-6 mvAREAvdmaxv0.003。然后定义目标函数， 这里体积为目标函数。指定优化方法（一阶优化方法），进行优化分析。 pavsmaxe显示优化状态变量最大应力（smaxe）迭代图形，（

11、如图5所示）。从图中可以网架杆 件的最大正应力有增大趋势，但是并未超过刚才的极限应力，这说明结构整体在不考虑屈曲 破坏的情况下还是安全的。8亍1/Vgg真 OPI-ERr 土丫 g-tI 肯1IL1VallieAPR 10 蚀DP-1,9图5显示最大位移（dmax）迭代图形，（如图6所示）。从图中可以清楚的看出，最大位 移有增大趋势，最终稳定在定义的区间里。u图66-bSet Ntirtibe显示杆件各截面(al,a2,a3,a4,a5,a6,a7)迭代图形，如图7所示，从图中可以看出，各杆 件的界面的值都有所下降，说明达到了预期的效果。CHA：-：LEk：!-APM 9 2BD9SO4.-B

12、上 g氐丁ifl,5Bet FfuraharV-s Lua4.33.&23:5:32prH lu 吨JPKj_u 1_b W HIKJ *I W gpifcHIMHM JU显示目标函数(volume)迭代图形，如图8所示。从图中可以明显的看出，整个桁架的用钢 量急剧下降，最终达到一个稳定值。jPTIHIZATICfJS e-t Nuto-b e-r本次优化分析中也也导出了最优解，整个网架总的体积约为0.29275e-2。比原优化前的设计 方案节约了大量的材料。7结论利用ANSYS的有限元分析理论可以轻松的解决一些网架的静力学问题，显示应力云图 以及位移云图。在运用ANSYS优化设计功能的基础上，为设计师提供了主动设计的依据， 在降低成本方面获得了巨大的经济相依，节省材料，实现可持续发展。在实际设计时，由于各种因数影响，截面不可能完全符合优化的结果，但是按照优化结 果去选择截面的话，用相同数量的材料，必然可以得到刚度较大的结构。通过上述简单的优化例子可知，利用ANSYS的优化设计功能进行结构的优化分析是可 行而且有效的，其方法简单，是有限元和优化分析的有机结合，为结构复杂的网架结构的设 计提供了新的方法。