基于Fortran90的满应力设计计算.pdfVIP

Science TechnologyVision 科 技 视 界 科技 探·索·争鸣 基于Fortran90的满应力设计计算 谢雨涵 邵雪杰 周 政 (河海大学 大禹学院，江苏 南京210098) 【摘 要】利用Fortran90的接近数学公式的 自然描述的特性，可以直接对矩阵和复数进行运算，来完成建筑结构设计中的满应力设计计 算。 关【键词】满应力；应力比法；Fortran90 在结构力学计算设计中，满应力计算是比较基本的计算设计，也 直接对矩阵和复数进行运算 ．这一点类似matlab：自诞生以来广泛地 是确定结构建筑工程最节省方案的重要步骤。由于现实生活中结构多 应用于数值计算领域，积累了大量高效而可靠的源程序：Fortran90的 为超静定结构 ，满应力设计时采用人工计算工作量较大，历时较长，精 推出使Fortran语言具备了现代高级编程语言的一些特性 度较差。因此，本文引用Fortrna90，采用应力比法设计迭代进程，从而 2．2 建立程序的基础资料数据库 可极大提高计算速度和精度。 见于满应力设计所涉及到的基础资料数据较多．本程序数据库选 用MicrosoftExcel2007作为其矩阵计算过程的同一数据库 将控制条 1 满应力计算原理和方法 件、结构资料 (包括初始截面积 、桁架杆件数和连接方式)以及工况资 1．1 满应力简介 料 (包括广义作用力、作用点和作用方向)输人矩阵中。 满应力设计 (FullyStressedDesign．FSD)是结构优化设计中准则 2．3 Fortran90程序流程 法的一种。FSD直接从结构力学的基本原理出发 ，以满应力为准则，保 对满应力设计 (应力比法)计算原理和计算过程进行程序语言化 证杆件的材料能够得到充分利用。 处理．具体如下程序流程(具体如图1)。 所谓满应力是指结构的各个杆件至少在一组确定的荷载组合下 承受容许应力或临界力。满应力设计的思路是在结构几何形状和结构 材料已经确定的情况下，通过调整杆件的截面，使其满足满应力准则。 这样得到的设计认为是处于满应力状态。由于所确定的杆件截面为其 下限值 。因此．满应力设计所得到的结果通常同时也是重量最轻的。 1．2 应力比法计算原理和方法 求解满应力采用应力比法 ．应力比法是满应力准则法中基本的一 种迭代方法．该法的主要原理是 ： 取杆件截面面积A、(各工况)轴力N、各工况最大拉力 N 、各工 况最大压力N… (o)、(各工况)应力o-=NA／、各工况最大拉应力 。=N ／ A、轴心压杆的稳定系数 (强度控制时， =i)、各工况最大等效压应 力 =N ／(A)(0)、杆件的拉应力比：o-I40-]、杆件的等效压应力比： ／【】、杆件的应力比r=max{o。，州【， ／[州)。 满应力法的主要公式： 应力比：r~)：max{0[)／[]，口[】}。 面积迭代式：A ： A 应力比ru(1表示应力大于容许应力．说明杆件超载，截面小了， ． 故A： A 大截面面积：应力比r1，表示应力小于容许应力，说明 材料还未被充分利用，采用的截面大了，故A=ru{~AtuJ减小截面面积。 具体计算步骤如下： (1)估算出结构的初始设计方案。若设计变量是杆件的截面积，则 初始设计为： A ：(A A A …A (2)对于某个杆件 i，计算各工况下杆件轴力N 计算杆件的应力 图 1 比 2．4 程序具体编制过程 (3)按 A。： 计算杆件新的面积。