桥梁工程CAD学习教案.pptx

《桥梁工程CAD学习教案.pptx》由会员分享，可在线阅读，更多相关《桥梁工程CAD学习教案.pptx（35页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。

1、会计学1桥梁工程桥梁工程CAD第一页，共35页。教材教材教材教材(jioci)(jioci)及参考资料及参考资料及参考资料及参考资料桥梁工程桥梁工程CADCAD，李丽芬编写，李丽芬编写(binxi)(binxi)，校内讲义校内讲义Midas/civilMidas/civil操作操作(cozu)(cozu)，李丽芬编写，校内，李丽芬编写，校内讲义讲义桥梁工程桥梁工程CADCAD课程实验指导书课程实验指导书，李丽芬编写，李丽芬编写，校内使用校内使用桥梁结构电算程序设计桥梁结构电算程序设计，颜东煌，田仲初，李学文编著，湖南大学出版社，颜东煌，田仲初，李学文编著，湖南大学出版社桥梁工程桥梁工程，罗旗帜

2、主编，华南理工大学出版社，罗旗帜主编，华南理工大学出版社Fortran90Fortran90学习教程学习教程，赵翔龙编著，北京大学出版社，赵翔龙编著，北京大学出版社Fortran PowerStationFortran PowerStation基础教程基础教程，徐明编著，清华大学出版社，徐明编著，清华大学出版社第1页/共35页第二页，共35页。第第1 1章章 编程基础知识编程基础知识 Fortran Fortran语言语言(yyn)(yyn)编程基础编程基础 1.1 Fortran 1.1 Fortran 语言语言(yyn)(yyn)程序设计初步程序设计初步 源程序的书写源程序的书写(shxi

3、)(shxi)格式格式 一个简单的一个简单的Fortran90 Fortran90 程序基本结构如下：程序基本结构如下：program program 程序名称程序名称 声明语句部分声明语句部分 执行语句部分执行语句部分 end program end program程序名称程序名称 的内容是可选的，只有end才是一个Fortran程序所必须的end 意味着程序的编译到此为止，即程序的结束 第2页/共35页第三页，共35页。源程序的书写源程序的书写(shxi)(shxi)格式格式 Fortran 90Fortran 90的语句的语句(yj)(yj)行可以是行可以是0 0到到132132个字符个

4、字符 Fortran 90Fortran 90允许允许(ynx)(ynx)出现空语句行出现空语句行 在在Fortran 90Fortran 90中，除赋值语句外，每个语句都要使用关键字开头中，除赋值语句外，每个语句都要使用关键字开头 一般而言，每行有一条语句，较短的语句可以放在一行使用，语句间一定要用分号隔开一般而言，每行有一条语句，较短的语句可以放在一行使用，语句间一定要用分号隔开 Fortran 90 Fortran 90 的注释语句是以感叹号为标志的，注释语句可以单独占一行，也可在程序的其他语句后面出现。的注释语句是以感叹号为标志的，注释语句可以单独占一行，也可在程序的其他语句后面出现。

5、如果遇到一条语句的长度超过了如果遇到一条语句的长度超过了Fortran 90Fortran 90所允许的行最大长度，需要移到下一行，则在语句末尾添加续行符（所允许的行最大长度，需要移到下一行，则在语句末尾添加续行符（&），编译器就会把下一行作为续行来处理。），编译器就会把下一行作为续行来处理。第3页/共35页第四页，共35页。常量常量(chngling)(chngling)Fortran 90Fortran 90提供了提供了5 5种内部数据类型：整型、实型、复型、字符型和逻辑种内部数据类型：整型、实型、复型、字符型和逻辑(lu j)(lu j)型。其中，整型、实型、复型为数值型数据类型，字符型

6、和逻辑型。其中，整型、实型、复型为数值型数据类型，字符型和逻辑(lu j)(lu j)型为非数值型数据类型。型为非数值型数据类型。uu整型常量实际上是整型常数或简称整型常量实际上是整型常数或简称(jinchng)(jinchng)为整数。整型常量的描述形式为：为整数。整型常量的描述形式为：nn。例如：。例如：+25+25，-339-339，0 0，+9999+9999uu实型常量实际上是实型常数或简称为实数。实型常量实际上是实型常数或简称为实数。小数形式实型常量的描述形式有如下小数形式实型常量的描述形式有如下3 3种：种：n.mn.m；n.n.；.m.m。合法的实型常量如。合法的实型常量如+1

7、.34+1.34，-0.25-0.25，9.9.，.3.3，-.567-.567指数形式的实型常量描述形式有如下指数形式的实型常量描述形式有如下4 4种：种：n.mESn.mES；n.ESn.ES；nESnES；.mES.mES。下面的表示方式都是合法的：。下面的表示方式都是合法的：2.0E22.0E2，3E-53E-5，-5.23E+4-5.23E+4，1.E61.E6，4E44E4，.45E3.45E3。第4页/共35页第五页，共35页。常量常量(chngling)(chngling)uu字符字符(zf)(zf)型常量，最基本的字符型常量，最基本的字符(zf)(zf)型常量是字符型常量是字

8、符(zf)(zf)串，它是由单引号或双引号括起来的若干个字符串，它是由单引号或双引号括起来的若干个字符(zf)(zf)构成。字符构成。字符(zf)(zf)型常量中空格是有用的，大小写也是有区别的。型常量中空格是有用的，大小写也是有区别的。uu逻辑型常量是表示条件逻辑型常量是表示条件(tiojin)(tiojin)成立与否的逻辑量，逻辑型常量只有两个，即条件成立与否的逻辑量，逻辑型常量只有两个，即条件(tiojin)(tiojin)成立为真，不成立为假。在成立为真，不成立为假。在Fortran90Fortran90中表示为如下形式：中表示为如下形式：.True.True.，表示，表示“真真”，即

9、满足逻辑条件，即满足逻辑条件(tiojin)(tiojin)。.False.False.，表示，表示“假假”，即不满足逻辑条件，即不满足逻辑条件(tiojin)(tiojin)。uu复型常量由实部和虚部组成。其描述形式如下：（复型常量由实部和虚部组成。其描述形式如下：（a a，b b），其中），其中a a是复数的实部，是复数的实部，b b是复数的虚部。例如：（是复数的虚部。例如：（1.01.0，-1.0-1.0）表示）表示1.0-1.0i1.0-1.0i；（；（4 4，5.3E25.3E2）表示）表示4+530i4+530i；（；（-3-3，4 4）表示）表示-3+4i-3+4i。第5页/共3

10、5页第六页，共35页。变量变量(binling)(binling)变量名：变量名：Fortran90Fortran90规定，变量、程序、函数、数组等的名字的大小为规定，变量、程序、函数、数组等的名字的大小为1 13131个字符，可以是字母或数字，下划线也可使用个字符，可以是字母或数字，下划线也可使用(shyng)(shyng)，但开头字符必须是字母。名字中的大写和小写字母是等价的。为了避免混淆，建议不要使用，但开头字符必须是字母。名字中的大写和小写字母是等价的。为了避免混淆，建议不要使用(shyng)Fortran(shyng)Fortran中已有的具有特定意义的名字作为变量名。中已有的具有特

11、定意义的名字作为变量名。变量类型：变量类型也分为变量类型：变量类型也分为(fnwi)5(fnwi)5类：整型变量、实型变量、复型变量、字符型变量和逻辑型变量。变量类型声明符分别为：类：整型变量、实型变量、复型变量、字符型变量和逻辑型变量。变量类型声明符分别为：integerinteger，realreal，complexcomplex，charactercharacter和和logicallogical。Fortran 90Fortran 90中的中的I-NI-N规则：除非特别声明外，在程序中的变量名，凡以规则：除非特别声明外，在程序中的变量名，凡以I I，J J，KK，L L，MM，NN这这

12、6 6个字母开关的都被默认为整型变量，以其他字母开头的变量名表示实型变量。个字母开关的都被默认为整型变量，以其他字母开头的变量名表示实型变量。第6页/共35页第七页，共35页。内部内部(nib)(nib)函数函数 函数的自变量可以是常量(chngling)、变量或表达式。如sin（2.0）、sin（x）（x是实型变量）和sin（）都是合法的。注意：三角函数中的角度单位是“弧度(hd)”而不是“度”，在实际编程时要注意转化。第7页/共35页第八页，共35页。算术算术(sunsh)(sunsh)表达式表达式 内部内部(nib)(nib)算术运算符是算术运算符是Fortran90Fortran90内

13、置的，共有加（内置的，共有加（+）、减（）、减（-）、乘（）、乘（*）、除（）、除（/）和乘方（）和乘方（*）5 5种算术运算符。种算术运算符。在进行整数除法在进行整数除法(chf)(chf)运算的时候，如果是整数之间的除法运算的时候，如果是整数之间的除法(chf)(chf)，那么结果也是整数。如：，那么结果也是整数。如：12/512/5的结果为的结果为2 2；7/87/8的结果为的结果为0 0。Fortran 90Fortran 90允许不同类型的数值型数据类型之间进行算术运算，但不允许在数值型数据与非数值型数据类型之间进行算术运算。允许不同类型的数值型数据类型之间进行算术运算，但不允许在数

14、值型数据与非数值型数据类型之间进行算术运算。同类型之间的算术运算的结果仍然保持原数据类型不变。如：同类型之间的算术运算的结果仍然保持原数据类型不变。如：2*52*5的值为的值为1010，是整数；的值为，是整数；的值为10.010.0，是实数。，是实数。不同类型的数值型数据类型之间的算术运算遵循一定的规则，即低级类型转换成高级类型，这种转换是系统自动进行的。整型数据和实型数据之间进行算术运算的时候，整型数据将转换为实型数据。不同类型的数值型数据类型之间的算术运算遵循一定的规则，即低级类型转换成高级类型，这种转换是系统自动进行的。整型数据和实型数据之间进行算术运算的时候，整型数据将转换为实型数据。

15、第8页/共35页第九页，共35页。赋值语句赋值语句(yj)(yj)一般格式：变量一般格式：变量(binling)=(binling)=常量或表达式常量或表达式 uu表控输入语句的一般格式为：表控输入语句的一般格式为：Read*Read*，变量，变量(binling)(binling)列表列表“*”“*”符号表示符号表示“表控输入表控输入”，不必指定输入数据的格式，只需将数据按其合法形式依次输入即可，所以又称自由格式输入。变量，不必指定输入数据的格式，只需将数据按其合法形式依次输入即可，所以又称自由格式输入。变量(binling)(binling)列表中的变量列表中的变量(binling)(bi

16、nling)之间用逗号分隔开。之间用逗号分隔开。简单的表控输入输出语句简单的表控输入输出语句 uu 表控输出语句的一般格式为：表控输出语句的一般格式为：Print*Print*，列表，列表 列表内容可以是常量、变量、表达式或字符串，它们之间用逗号分隔开。列表内容可以是常量、变量、表达式或字符串，它们之间用逗号分隔开。第9页/共35页第十页，共35页。1.2 Fortran 901.2 Fortran 90数据结构数据结构(sh j ji u)(sh j ji u)整型数据整型数据(shj)(shj)整型数据定义的一般形式整型数据定义的一般形式(xngsh)(xngsh)为：为：IntegerI

17、nteger（kind=kind_valuekind=kind_value），属性列表，属性列表变量列表变量列表 Kind定义了整型数据类别 属性列表中可以是数组维数说明符Dimension和指针类型说明符Pointer等实实型数据型数据 实型数据有单精度和双精度之分，一般形式为：实型数据有单精度和双精度之分，一般形式为：RealReal（kind=kind_valuekind=kind_value），属性列表，属性列表 变量列表变量列表 Kind定义了实型数据类别 第10页/共35页第十一页，共35页。1.2 Fortran 901.2 Fortran 90数据结构数据结构(sh j ji

18、u)(sh j ji u)复型数据复型数据(shj)(shj)复型数据与实型数据一样也有单精度和双精度之分，一般形式如下：复型数据与实型数据一样也有单精度和双精度之分，一般形式如下：ComplexComplex（kind=kind_valuekind=kind_value），属性列表，属性列表(li bio)(li bio)变量列表变量列表(li bio)(li bio)Kind定义了复型数据类别 复数由两部分组成，实部和虚部。字符型数据字符型数据 字符型数据类型声明的一般形式为：字符型数据类型声明的一般形式为：CharacterCharacter类型参数类型参数，属性列表，属性列表 变量列表

19、变量列表 类型参数有两种形式：类型参数有两种形式：Len=Len=字符串长度字符串长度 例：例：charactercharacter（len=4len=4）*字符串长度字符串长度 例：例：character*4character*4 第11页/共35页第十二页，共35页。1.2 Fortran 901.2 Fortran 90数据结构数据结构(sh j ji u)(sh j ji u)逻辑逻辑(lu j)(lu j)型数据型数据 逻辑型数据逻辑型数据(shj)(shj)定义的一般形式为：定义的一般形式为：LogicalLogical（kind=kind_valuekind=kind_value

20、），属性列表，属性列表变量列表变量列表 Kind定义了逻辑型数据类别 逻辑型常量只有两个：.True.和.False.使用使用DataData语句初始化变量语句初始化变量 Fortran 90Fortran 90提供的提供的DataData语句进行初始化变量，语句进行初始化变量，DataData语句可以同时为多个变量赋初值，其一般形式为：语句可以同时为多个变量赋初值，其一般形式为：Data Data 变量变量1/1/数据列表数据列表/，变量变量2/2/数据列表数据列表/第12页/共35页第十三页，共35页。第第2 2章章 截面截面(jimin)(jimin)几何特性几何特性计算计算 2.1 2

21、.1 截面几何截面几何(j h)(j h)特性计算一般方法特性计算一般方法 截面几何特性包括内容很多，桥梁计算中通常用到截面几何特性包括内容很多，桥梁计算中通常用到截面面积截面面积 A A、形心位置、形心位置(wi zhi)yc(wi zhi)yc、面积矩、面积矩 SxSx、抗弯惯性、抗弯惯性IxIx、抗扭惯矩、抗扭惯矩 IT IT和截面模量和截面模量WxWx等。等。第13页/共35页第十四页，共35页。第第2 2章章 截面几何特性截面几何特性(txng)(txng)计算计算 1.2 1.2 概述概述(i sh)(i sh)第14页/共35页第十五页，共35页。第第2 2章章 截面截面(jim

22、in)(jimin)几何特性几何特性计算计算 如果区域如果区域A A可剖分成为可剖分成为(chngwi)n(chngwi)n个子区域个子区域AiAi假设每个假设每个AiAi的形状很规则，其面积矩和惯性矩很容易用公式的形状很规则，其面积矩和惯性矩很容易用公式(gngsh)(gngsh)表示，表示，则累加各个子区域和面积矩和惯性矩就可以求得整个则累加各个子区域和面积矩和惯性矩就可以求得整个A A的面积矩和惯性矩。有的面积矩和惯性矩。有的截面的截面A A，可以剖分成为梯形之和，则用梯形分块法；有的截面，可以剖分成为梯形之和，则用梯形分块法；有的截面A A，可以剖分，可以剖分为三角形之和，则用三角形分

23、块法。如果区域的边界可剖分为左、右两条曲为三角形之和，则用三角形分块法。如果区域的边界可剖分为左、右两条曲线，见图线，见图2-12-1，即，即f(y)f(y)和和g(y)g(y)，则：，则：第15页/共35页第十六页，共35页。第第2 2章章 截面截面(jimin)(jimin)几何特性几何特性计算计算 封闭薄壁截面封闭薄壁截面(jimin)(jimin)的抗扭惯矩公式为：的抗扭惯矩公式为：倘若遇到封闭薄壁截面上带有倘若遇到封闭薄壁截面上带有“翅翼翅翼”的一般情的一般情况况(qngkung)(qngkung)，如图，如图2-22-2所示，则其总抗扭惯矩所示，则其总抗扭惯矩可近似地叠加计算：可近

24、似地叠加计算：第16页/共35页第十七页，共35页。第第2 2章章 截面几何截面几何(j(j h)h)特性计算特性计算 式中，式中，cici矩形截面矩形截面(jimin)(jimin)抗扭刚度系数，根据抗扭刚度系数，根据t/bt/b比值计算比值计算第17页/共35页第十八页，共35页。第第2 2章章 截面几何特性截面几何特性(txng)(txng)计算计算 2.2 2.2 梯形梯形(txng)(txng)分块法分块法 梯形梯形(txng)(txng)上底、下底和高分别为上底、下底和高分别为a a、b b和和h h 梯形计算公式梯形计算公式 面积：面积：形心位置：形心位置：对形心轴的惯性矩：对形

25、心轴的惯性矩：第18页/共35页第十九页，共35页。图示截面，可以图示截面，可以(ky)(ky)剖分成若干个梯形从下至上顺序连续排列。用节线（分割梯形的线段）处的宽度和高度（节线的纵坐标）来描述，设：剖分成若干个梯形从下至上顺序连续排列。用节线（分割梯形的线段）处的宽度和高度（节线的纵坐标）来描述，设：顺序排列梯形顺序排列梯形(txng)(txng)分块法分块法 N N：节线总条数；：节线总条数；B B（N N）：节线宽度）：节线宽度(kund)(kund)数组；数组；H H（N N）：节线高度数组。）：节线高度数组。对于第对于第i i个梯形，个梯形，其下底宽其下底宽b b=B B(i i)，

26、上底宽上底宽a a=B B(i i+1)+1)，高高h h=HH(i+1i+1)-)-H H(i i)对对x x轴的面积矩轴的面积矩A Aii对对x x轴的惯性矩轴的惯性矩 整个截面的形心轴至的距离整个截面的形心轴至的距离 整个截面对形心轴的惯性矩整个截面对形心轴的惯性矩 第19页/共35页第二十页，共35页。图示截面，虽然图示截面，虽然(surn)(surn)也可分割成若干个梯形，但各梯形并非顺序连续排列，不能用节线来描述了，可以对各个梯形单独描述。由于只计算截面竖向平面内的几何特性，故梯形的横向位置可以任意，描述单个梯形用如下参数：也可分割成若干个梯形，但各梯形并非顺序连续排列，不能用节线

27、来描述了，可以对各个梯形单独描述。由于只计算截面竖向平面内的几何特性，故梯形的横向位置可以任意，描述单个梯形用如下参数：任意排列任意排列(pili)(pili)梯形分块法梯形分块法 N N：梯形数；：梯形数；A A（N N）梯形上底宽；梯形上底宽；B B（N N）梯形下底宽；梯形下底宽；H1H1（N N）梯形高度；梯形高度；H H（N N）梯形下底至截面梯形下底至截面(jimin)(jimin)底缘的距离（高度）。底缘的距离（高度）。求几何特性的方法与顺序排列梯形分块法完全相同求几何特性的方法与顺序排列梯形分块法完全相同 第20页/共35页第二十一页，共35页。梯形梯形(txng)(txng)

28、分块法程序分块法程序 转转Fortran PowerstationFortran Powerstation讲解运行讲解运行(ynxng)(ynxng)程序程序:截面抗弯截面抗弯.f90.f90例题：梯形例题：梯形(txng)(txng)分块法计算截面特性分块法计算截面特性 右图：全截面分右图：全截面分5 5个梯形，个梯形，6 6条节线条节线按顺序排列梯形分块法计算按顺序排列梯形分块法计算截面抗弯截面抗弯.f90.f90程序中的输入语句程序中的输入语句OPEN(3,FILE=datainputKW.dat,STATUS=OLD)OPEN(3,FILE=datainputKW.dat,STATUS

29、=OLD)READ(3,*)IYREAD(3,*)IYREAD(3,*)NREAD(3,*)N DO 100 J=1,N DO 100 J=1,NREAD(3,*)NU(J)READ(3,*)NU(J)READ(3,*)(B(I),H(I),I=1,NU(J)READ(3,*)(B(I),H(I),I=1,NU(J)100 CONTINUE100 CONTINUE对应输入语句，建立数据文件对应输入语句，建立数据文件文件名：文件名：inputKW.datinputKW.dat文件内容：文件内容：1 11 16 60.6,0.00.6,0.00.6,0.50.6,0.50.18,0.650.18,

30、0.650.18,2.110.18,2.112.2,2.332.2,2.332.2,2.452.2,2.45变量变量IYIY、N N、NU(N)NU(N)、B(NU),H(BU)B(NU),H(BU)的含义的含义参见源程序参见源程序截面抗弯截面抗弯.f90.f90中的注释语句。中的注释语句。第21页/共35页第二十二页，共35页。梯形梯形(txng)(txng)分块法程序分块法程序 例题：梯形分块法计算截面例题：梯形分块法计算截面(jimin)(jimin)特性特性 截面截面(jimin)(jimin)抗弯抗弯.f90.f90程序中的输出语句程序中的输出语句OPEN(4,FILE=dataKW

31、.dat,STATUS=UNKNOWN)OPEN(4,FILE=dataKW.dat,STATUS=UNKNOWN)DO 100 J=1,NDO 100 J=1,NWRITE(4,80)AA(J),HX,RI(J)WRITE(4,80)AA(J),HX,RI(J)8080FORMAT(F15.5,1X,F15.5,1X,F15.5,1X)FORMAT(F15.5,1X,F15.5,1X,F15.5,1X)100 CONTINUE100 CONTINUE对应输出语句，生成的计算结果数据文件对应输出语句，生成的计算结果数据文件文件名：文件名：KW.datKW.dat文件内容：文件内容：1.1471

32、0 1.47399 .934251.14710 1.47399 .93425操作步骤操作步骤（1 1）启动）启动Fortran Powerstation/Microsoft Developer StudioFortran Powerstation/Microsoft Developer Studio（2 2）新建）新建Project WorkspaceProject Workspace。例如：。例如：NameName为为kwgjkwgj，路径；，路径；E:E:桥梁桥梁CADkwgjCADkwgj（3 3）建立数据文件。文件名：）建立数据文件。文件名：inputKW.datinputKW.dat

33、，存放路径：存放路径：E:E:桥梁桥梁CADkwgjdata inputKW.datCADkwgjdata inputKW.dat（4 4）将源程序文件）将源程序文件“截面抗弯截面抗弯.90.90”复制到复制到E:E:桥梁桥梁CADkwgjCADkwgj（5 5）将源程序）将源程序“E:E:桥梁桥梁CADkwgjCADkwgj截面抗弯截面抗弯.90.90”insert files into projectinsert files into project（6 6）Compiling ,Limking Compiling ,Limking （7 7）Build/execute kwgj.exeB

34、uild/execute kwgj.exe。程序运行结束，自动生成数据文件程序运行结束，自动生成数据文件E:E:桥梁桥梁CADkwgjdataKW.datCADkwgjdataKW.dat变量变量AA(J),HX,RI(J)AA(J),HX,RI(J)的含义参见源程序的含义参见源程序截面抗弯截面抗弯.f90.f90中的注释语句。中的注释语句。第22页/共35页第二十三页，共35页。梯形梯形(txng)(txng)分块法程序分块法程序 例题例题(lt)(lt)：梯形分块法计算截面特性：梯形分块法计算截面特性 右图：全截面右图：全截面(jimin)(jimin)分分5 5个梯形，个梯形，6 6条节

35、条节线线按任意排列梯形分块法计算按任意排列梯形分块法计算截面抗弯截面抗弯.f90.f90程序中的输入语句程序中的输入语句OPEN(3,FILE=datainputKW.dat,STATUS=OLD)OPEN(3,FILE=datainputKW.dat,STATUS=OLD)READ(3,*)IYREAD(3,*)IYREAD(3,*)NREAD(3,*)N DO 100 J=1,N DO 100 J=1,NREAD(3,*)NU(J)READ(3,*)NU(J)READ(3,*)(B(I),A(I),H1(I),H(I),I=1,NU(J)READ(3,*)(B(I),A(I),H1(I),

36、H(I),I=1,NU(J)100 CONTINUE100 CONTINUE对应输入语句，建立数据文件对应输入语句，建立数据文件文件名：文件名：inputKW.datinputKW.dat文件内容：文件内容：2 21 15 50.6,0.6,0.5,0.0,0.6,0.6,0.5,0.0,0.6,0.18,0.15,0.5,0.6,0.18,0.15,0.5,0.18,0.18,1.46,0.65,0.18,0.18,1.46,0.65,0.18,2.2,0.22,2.11,0.18,2.2,0.22,2.11,2.2,2.2,0.12,2.332.2,2.2,0.12,2.33计算结果与按顺

37、序排列梯形计算结果与按顺序排列梯形分块法的结果一致。分块法的结果一致。变量变量IYIY、N N、NU(N)NU(N)、B(NU)B(NU)、A(NU)A(NU)、H1(NU)H(BU)H1(NU)H(BU)的含义参见源程序的含义参见源程序截面抗弯截面抗弯.f90.f90中的注释语句。中的注释语句。第23页/共35页第二十四页，共35页。例题：梯形例题：梯形(txng)(txng)分块法计算截面特分块法计算截面特性性 右图：全截面右图：全截面(jimin)(jimin)分分8 8个梯个梯形，形，9 9条节线条节线按顺序排列梯按顺序排列梯形形(txng)(txng)分分块法计算块法计算数据文件名：

38、数据文件名：inputKW.datinputKW.dat1 11 19 96.0,0.06.0,0.06.0,0.256.0,0.251.4,0.251.4,0.251.0,0.451.0,0.451.0,2.01.0,2.06.514,2.226.514,2.228.514,2.228.514,2.2210.0,2.3510.0,2.3510.0,2.510.0,2.5按任意排列梯形按任意排列梯形分块法计算分块法计算数据文件名：数据文件名：inputKW.datinputKW.dat2 21 19 96.0,6.0,0.25,0.06.0,6.0,0.25,0.00.7,0.5,0.2,0.

39、250.7,0.5,0.2,0.250.7,0.5,0.2,0.250.7,0.5,0.2,0.250.5,0.5,1.55,0.450.5,0.5,1.55,0.450.5,0.5,1.55,0.450.5,0.5,1.55,0.450.5,3.257,0.22,2.00.5,3.257,0.22,2.00.5,3.257,0.22,2.00.5,3.257,0.22,2.08.514,10.0,0.13,2.228.514,10.0,0.13,2.2210.0,10.0,0.15,2.3510.0,10.0,0.15,2.35生成的计算结果数据文件生成的计算结果数据文件文件名：文件名：KW

40、.datKW.dat文件内容：文件内容：6.81995 1.51387 5.962076.81995 1.51387 5.96207第24页/共35页第二十五页，共35页。截面截面(jimin)(jimin)抗扭惯矩计算程序抗扭惯矩计算程序 转转Fortran PowerstationFortran Powerstation讲解讲解(jingji)(jingji)运行程序运行程序:Jmtx-:Jmtx-knj.f90knj.f90例题例题(lt)(lt)：计算：计算T T形截面抗扭惯矩形截面抗扭惯矩 右图：全截面分右图：全截面分5 5个矩形。其中，矩形个矩形。其中，矩形和和 的的宽取上、下底平

41、均宽，分别为宽取上、下底平均宽，分别为0.390.39和和1.191.19Jmtx-knj.f90Jmtx-knj.f90程序中的输入语句程序中的输入语句OPEN(3,FILE=datainputKN.dat,STATUS=OLD)OPEN(3,FILE=datainputKN.dat,STATUS=OLD)READ(3,*)IYREAD(3,*)IYREAD(3,*)NREAD(3,*)N DO 100 J=1,N DO 100 J=1,NREAD(3,*)NU(J)READ(3,*)NU(J)READ(3,*)(B(I),H(I),I=1,NU(J)READ(3,*)(B(I),H(I),

42、I=1,NU(J)100 CONTINUE100 CONTINUE对应输入语句，建立数据文件对应输入语句，建立数据文件文件名：文件名：inputKN.datinputKN.dat文件内容：文件内容：3 31 15 50.6,0.50.6,0.50.39,0.150.39,0.151.46,0.181.46,0.181.19,0.221.19,0.222.2,0.122.2,0.12 变量变量IYIY、N N、NU(N)NU(N)、B(NU),H(BU)B(NU),H(BU)的含义参的含义参见源程序见源程序Jmtx-knj.f90Jmtx-knj.f90中的注释语句。中的注释语句。第25页/共3

43、5页第二十六页，共35页。截面抗扭惯矩计算截面抗扭惯矩计算(j sun)(j sun)程序程序 截面抗弯截面抗弯.f90.f90程序程序(chngx)(chngx)中的输出语句中的输出语句OPEN(4,FILE=dataKN.dat,STATUS=UNKNOWN)OPEN(4,FILE=dataKN.dat,STATUS=UNKNOWN)DO 100 J=1,NDO 100 J=1,NWRITE(4,80)XIT(J)WRITE(4,80)XIT(J)8080FORMAT(F12.5,1X)FORMAT(F12.5,1X)100 CONTINUE100 CONTINUE对应输出语句对应输出语句(yj)(yj)，生成的计算结果数据文件，生成的计算结果数据文件文件名：文件名：KN.datKN.dat文件内容：文件内容：.00760 .00760变量变量XIT(J)XIT(J)表示截面抗扭惯矩。表示截面抗扭惯矩。第26页/共35页第二十七页，共35页。第27页/共35页第二十八页，共35页。第28页/共35页第二十九页，共35页。第29页/共35页第三十页，共35页。第30页/共35页第三十一页，共35页。第31页/共35页第三十二页，共35页。第32页/共35页第三十三页，共35页。第33页/共35页第三十四页，共35页。第34页/共35页第三十五页，共35页。