基于有限元法的拱结构优化设计..doc

基于有限元法的拱结构优化设计 黎杰 张娅丽 （武桥重工集团股份有限公司桥梁设备研究所，430056） 摘要：本文利用有限元法对承受竖直均布载荷的静定拱（三铰拱）和超静定拱进行了高跨比优化设计，得到了如下结论：在保证强度的基础上，拱的高跨比取为1/3能使拱的自重最小。 关键词： 静定拱 超静定拱 高跨比 有限元法 优化设计 Abstract: The finite element method was used to execute the optimal design of statically determinate arch and statically indeterminate arch loaded by a vertical uniform pressure. The conclusion is as follow: on the basis of enough strength, the height-span ratio of arch is taken to be 1/3 can minimize the weight. Key words: statically determinate arch, statically indeterminate arch, height-span ratio, finite element method, optimal design. 0 引言 拱为常见建筑结构之一，型态定义为中央上半成圆弧或抛物线曲线。拱结构比桁架结构具有更大的力学优点在外荷作用下，拱主要产生压力，使构件摆脱了弯曲变形。 式中：——拱高 ——跨度 图1 竖向均布载荷作用下的静定和超静定拱计算模型 2.2 设计变量 由文献[2]可知，计算模型中主要的设计标量为拱高f和拱截面面积A。由于截面面积A由截面的轮廓尺寸（长a,宽b）以及板厚t决定，因此本文给定截面的轮廓尺寸，将板厚t和高跨比f/L作为设计变量。其中约束条件为： 1mm≤t≤10mm 0.1≤f/L≤1 2.3 约束函数 本文将满足强度要求作为约束条件，即计算所得的最大Von-mises应力应小于材料的许用应力180MPa。 2.4 优化目标 优化设计的目标函数为拱结构的自重，即拱结构的体积。当设计变量在约束条件区域内变化时满足约束函数的最小拱结构的体积为最优解。 3 命令流方法 *create,arch ! 创建优化分析文件 q=100000 ! 压力 s=180e6 ! 许用应力 L=1 ! 跨度为1m f=0.2 ! 拱高，等于高跨比 h=0.045 ! 截面宽 b=0.045 ! 截面高 t=0.002 ! 板厚 !******建模******** /PREP7 *do,I,0,50,1 x1=I/50*1 y1=4*F/(L*L)*x1*(L-x1) k,I+1,x1,y1 *enddo bsplin,all ! 创建抛物线 *do,I,1,49,1 x1=I/50*1 y1=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,HREC SECDATA,h,b,t,t,t,t ! 定义截面 lmesh,all ! 划分网格 EMODIf,all,SECNUM,1 ! 定义单元截面属性 !******加载******** D,all,uz d,1,ux d,1,uy d,2,ux, d,2,uy ! 施加位移约束 f,all,fy,-q/51 nsel,s,node,,28 ESLN,S ! （仅静定拱需要）释放中间铰点旋转自由度 ENDRELEASE,,1,ROTZ nsel,all esel,all !******求解******** /solu solve !******后处理******** ! -------------- /post1 *get,smax,SECR,,s,eqv,max !令smax为最大应力 etable,evol,volu ssum *get,vtot,ssum,,item,evol !令vtot为体积 finish *end *use,arch !