空间桁架结构程序设计(Fortran).doc

实用文案

标准文档

空间桁架静力分析程序及算例

1、变量及数组说明

输入数据

控制数据

NF

单个节点的自由度数

NP

结构离散节点的总数

NE

结构离散单元的总数

NM

结构中单元不同的特征数类的总数

NR

结构受约束节点的总数

NCF

结构受外荷载作用的节点总数

ND

一个单元的节点总数

几何数据

X(NP)

节点X坐标数组

Y(NP)

节点Y坐标数组

Z(NP)

节点Z坐标数组

ME(ND,NE)

单元节点信息存储矩阵

ME(1,NE)存储杆件始端节点号

ME(2,NE)储存杆件末端节点号

RR(2,NR)

结构约束信息矩阵

RR(1,NR)存放受有约束的节点号

RR(2,NR)存放节点位移约束情况

单元特征数据

AE(2,IN)

单元特征数类数组

AE(1,IN)单元的弹性模量

AE(2,IN)单元的横截面面积

NAE(NE)

单元特征类信息存储数组

荷载数据

PF(4,NCF)

外荷载信息数组

PF(1,NCF)存放外荷载作用的节点号

PF(2,NCF)存放X方向的外荷载

PF(3,NCF)存放Y方向的外荷载

PF(4,NCF)存放Z方向的外荷载

输出数据

位移

DIST(NPF)

节点位移数组

DIST(NF*I-2)存放I节点X方向的位移

DIST(NF*I-1)存放I节点Y方向的位移

DIST(NF*I)存放I节点Z方向的位移

力

SG(NE)

单元内力数组

SM(NE)

单元截面应力数组

FL(NF*NR)

支座反力数组

FL(NF*I-2)存放受约束的I节点X方向的反力

FL(NF*I-1)存放受约束的I节点Y方向的反力

FL(NF*I)存放受约束的I节点Z方向的反力

中间变量

NPF=NF*NP

二维总刚度矩阵的最大行数

NDF=ND*NF

一个单元的自由度总数（2*3=6）

IN

单元特征类总数

AKE(2,2)

单元在局部坐标系中的刚度局矩阵

BL

杆件单元长度

T(2,6)

坐标转换矩阵

TAK(6,6)

单元在总体坐标系中的刚度矩阵

IT(NF,NP)

节点联系数组

LMT(NDF,NE)

单元联系数组

MAXA(NPF)

结构二维总刚度矩阵主对角元地址数组

NWK

结构一维总刚度矩阵的总容量

CKK(NWK)

结构一维总刚度矩阵

NN

结构矩阵方程的方程总数（去掉约束）

NNM

NNM=NN+1

V(NN)

已知节点荷载列阵数组，回代完成后为存放结构位移

PP(NPF)

所有节点荷载列阵数组

2、空间桁架结构有限元分析程序源代码

!主程序（读入文件，调用总计算程序，输出结果）

CHARACTERIDFUT*20,OUTFUT*20

WRITE(*,*)InputDataFilename:

READ(*,*)IDFUT

OPEN(11,FILE=IDFUT,STATUS=OLD)

WRITE(*,*)OutputFilename:

READ(*,*)OUTFUT

OPEN(12,FILE=OUTFUT,STATUS=UNKNOWN)

WRITE(12,*)*****************************************

WRITE(12,*)*ProgramforAnalysisofSpaceTrusses*

WRITE(12,*)*SchoolofCivilEngineeringCSU*

WRITE(12,*)*2012.6.25DesignedByMuZhaoxiang*

WRITE(12,*)*****************************************

WRITE(12,*)

WRITE(12,*)*****************************************

WRITE(12,*)*************TheInputData****************

WRITE(12,*)*****************************************

WRITE(12,100)

READ(11,*)NF,NP,NE,NM,NR,NCF,ND

WRITE(12,110)NF,NP,NE,NM,NR,NCF,ND

100FORMAT(6X,TheGeneralInformation/2X,NF,5X,NP,5X,NE,5X,NM,5X,NR,

5X,NCF,5X,ND)

110FORMAT(2X,I2,6I7)

NPF=NF*NP

NDF=ND*NF

CALLANALYSE