资源描述：
《基于有限元法拱结构优化设计》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、基于有限元法拱结构优化设计 摘要：本文利用有限元法对承受竖直均布载荷的静定拱（三铰拱）和超静定拱进行了高跨比优化设计，得到了如下结论：在保证强度的基础上，拱的高跨比取为1/3能使拱的自重最小。关键词：静定拱超静定拱高跨比有限元法优化设计Abstract:Thefiniteelementmethodwasusedtoexecutetheoptimaldesignofstaticallydeterminatearchandstaticallyindeterminatearchloadedbyaverticaluniformpressure.Theconc
2、lusionisasfollow:onthebasisofenoughstrength,theheight-spanratioofarchistakentobe1/3canminimizetheweight.Keywords:staticallydeterminatearch,staticallyindeterminatearch,height-spanratio,finiteelementmethod,optimaldesign.中图分类号：S611文献标识码：A文章编号：0引言8拱为常见建筑结构之一，型态定义为中央上半成圆弧或抛物线曲线。拱结构比桁
3、架结构具有更大的力学优点。在外荷作用下，拱主要产生压力，使构件摆脱了弯曲变形。在拱结构形状设计中，最重要的是确定合理的拱轴线方程，使拱结构的强度得到充分的利用的前提下具有最好的经济性。在拱结构设计中，一般载荷和跨度是已知量，拱轴线方程主要由拱高决定。拱的高跨比究竟取何值，才能使设计出的拱结构最优，造价最低？对于这个问题很多专家和学者都做出过研究，如文献[1]给出了拱的高跨比的一些推荐值，但对高跨比的最佳取值没有提出建议。文献[2]对静定拱的最优高跨比给出了理论上的推导公式，但对超静定拱未能进行理论推导计算。本文利用有限元软件ANSYS对静定拱（三铰拱）
4、和超静定拱进行了高跨比优化设计，对文献[2]的理论推导结果进行验证和补充。1ANSYS优化设计简介ANSYS优化设计大致有如下七步：1有限元参数化建模；2定义优化变量；3选择优化方法；4执行优化分析；5后处理。ANSYS提供了两种求最小目标函数的方法——零阶方法和一阶方法，另外还允许用户使用自己开发的算法，每种算法的特点如下：零阶方法：这种方法利用曲线拟合的方式逼近状态变量和目标函数，对硬件要求低，计算耗时少，但计算精度相对较低。8一阶方法：这种方法使用状态变量和目标函数对设计变量的一阶偏导数。在每次迭代中，梯度计算(用最大斜度法或共轭方向法)确定搜索
5、方向，并用线搜索法对无约束问题进行最小化此方法精度很高，适合独立变量变化范围很大和较大的设计空间的情况，计算耗时较多。2拱结构优化设计模型2.1计算模型本文模型与文献[2]大体一致，拱轴线方程为抛物线方程，竖向均布沿水平轴均匀变化的压力。拱结构形式为静定和超静定两种，如图一所示。拱轴线方程为：式中：——拱高——跨度图1竖向均布载荷作用下的静定和超静定拱计算模型2.2设计变量由文献[2]可知，计算模型中主要的设计标量为拱高f和拱截面面积A。由于截面面积A由截面的轮廓尺寸（长a,宽b）以及板厚t决定，因此本文给定截面的轮廓尺寸，将板厚t和高跨比f/L作为设
6、计变量。其中约束条件为：1mm≤t≤10mm0.1≤f/L≤12.3约束函数8本文将满足强度要求作为约束条件，即计算所得的最大Von-mises应力应小于材料的许用应力180MPa。2.4优化目标优化设计的目标函数为拱结构的自重，即拱结构的体积。当设计变量在约束条件区域内变化时满足约束函数的最小拱结构的体积为最优解。3命令流方法*create,arch!创建优化分析文件q=100000!压力s=180e6!许用应力L=1!跨度为1mf=0.2!拱高，等于高跨比h=0.045!截面宽b=0.045!截面高t=0.002!板厚!******建模******
7、**/PREP7*do,I,0,50,1x1=I/50*1y1=4*F/(L*L)*x1*(L-x1)k,I+1,x1,y1*enddo8bsplin,all!创建抛物线*do,I,1,49,1x1=I/50*1y1=4*F/(L*L)*x1*(L-x1)HPTCREATE,LINE,1,,COORD,x1,y1,0*enddo!创建硬点et,1,188!定义单元类型mp,ex,1,2.06e11!定义杨氏模量mp,nuxy,1,0.3!定义泊松比SECTYPE,1,beam,HRECSECDATA,h,b,t,t,t,t!定义截面lmesh,all!
8、划分网格EMODIf,all,SECNUM,1!定义单元截面属性!******加载******
