c语言数据文件流程图,C语言程序设计等参元程序流程图.ppt

程序需要的主要数据 StructFE COORDS 定义结构体类型intn 结点总数double x x 2n 结点坐标 x n y n 注 前n个为结点的x坐标 后n个为结点的y坐标 structFE COORDSF X 定义结构体变量F X 来存放结点信息 设为全局变量 程序需要的主要数据 structFE TOPOLOGY 定义类型longn 单元总数long t t 4n 单元的节点编码 structFE TOPOLOGYF T 定义结构体变量F T 来存放单元信息 可设为全局变量 如第5单元的结点编号为F T t 16 F T t 19 structFE PKMATRIX 定义结构体类型intn 自由度总数 即总刚方程的阶数int i i n 总刚矩阵的索引数组 即主元素在k中的位置double k k i n 1 1 变半带宽一维压缩存储的总刚矩阵 structFE PKMATRIXF K 定义结构体变量F X 来存放总刚矩阵的信息 如 F K i i 代表下标为i行主元素在压缩存储总刚F K k 中的位置 nd F X i i 1 F K i i 代表下标为i行非0元素的个数 structFE BOUNDARY 定义结构体类型intn 约束总数或结点荷载数int node 约束对应的结点编号int attr 约束的方向值 1 x方向 2 y方向 3 转动方向double value 支座沉降值或荷载值 structFE BOUNDARYF U F P 定义结构体变量F U F P 分别来存放边界约束和结点荷载的信息 structFE STRESS F Slongn 输出的应力 应变 点总数double x x 2n 相应应力点的坐标double s s 3n 应力 应变 值 structFE STRESSF S 记录的单元应力 应变 信息 包括每个应力点的位置坐标和三个应力 应变 值 doubleMAT 5 定义材料信息doubleE 材料弹性模量doublev 泊松比doublet 厚度doublepx x方向的容重doublepy y方向的容重 平面刚架静力计算程序流程图 单元循环 单元循环 平面刚架静力计算程序流程图 单元循环 单元循环 程序模块设计 shortReaddata 4 char file 读取平面四结点等参元分析需要的所有数据 如结点信息 单元信息 截面信息 边界支座信息 荷载信息 材料信息等 在此函数中定义全局变量F X F T F U F P MAT 并从文件中读取相应的数据 函数所需参数 file为数据文件名 成功读取数据 函数的返回值为0 数据文件不存在 返回 1 内存出错 返回1 shortKindex 4 建立采用4结点等参元计算平面问题的变带宽一维压缩存储的总刚矩阵的索引数组F K i 即总刚矩阵的主元素k i i 被压缩存储后 在F K k 中的位置 定义存储总刚信息的结构体变量F K 并计算出总刚矩阵F K k 元素的个数 并为之分配存储单元 函数需要的数据 结点信息F X 单元信息F T 均为全局变量 返回总刚矩阵的最大半带宽值 如果内存出错 则返回0 voidEkm2d 4 计算变带宽一维压缩存储的总刚度矩阵F K k 所需数据 结点信息F X 单元信息F T 截面信息F AI 材料信息MAT 总刚索引信息F K 均为全局变量 需调用的函数 voidEvdt2d double d 计算弹性矩阵d 5 voidEk2002 double xy double d double ek 计算4结点等参元的单元刚度矩阵ek 8 8 voidP load 4 float p 计算结构总的荷载向量 p 所需数据 结点信息F X 单元信息F T 荷载信息F P 结构自重 均为全局变量 需调用的函数 voidPbd2002 double xy double pxy double ep 计算单元自重引起的等效结点荷载 voidZuv2d structFE BOUNDARYu float p 根据边界条件修改刚度方程 主元素乘大数法 主元素置1法 u 记录边界支座条件的结构体变量 p 结构总的的荷载向量 高斯法解方程组 intSolvc1 long n float ak longjn Solvc1 总刚矩阵三角分解 ak L U Uji Lij Lii 正确返回1 当主元素ifk i i 0 返回 i voidSolvc2 float p long n float ak intjn Solvc2 回代求解 ak 一维压缩存储的总刚矩阵F K k n jn 总刚矩阵索引数组 jn 总刚方程的阶数p jn 荷载向量 函数返回时 存放结点位移 voidStr2d 4 double w structFE STRESS s 计算所有单元高斯点的应力向量F S 所需数据 结点位移信息 w 单元信息F T 节点信息 F X 均为全局变量 s 所有单元高斯点的应力向量 调用函数 voidStr2002 double xu double d double p double sr 计算4节点等参元的单元高斯点应力向量 8 voidStn2d 4 double w structFE STRESS s 计算所有单元高斯点的应变向量F S 所需数据 结点位移信息 w 单元信息F T 节点信息 F X 均为全局变量 s 所有单元高斯点的应变向量 调用函数 voidStn2002 double xu double d double p double sr 计算4节点等参元的单元高斯点应变向量

展开阅读全文