资源描述：
《有限元法大作业》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。
1、有限元法大作业一平面刚架的程序用VisualC++编制的平面刚架的源程序如下：///////////////////////////////////////////////////////程序开始//////////////////////////////////////////////////////////////////#include"iostream.h"#include"math.h"#include"stdlib.h"#include"conio.h"//*****************//声明必要变量//*****************#definePI3.14159
2、2654intNE;//单元数intNJ;//节点数intNZ;//支承数intNPJ;//有节点载荷作用的节点数intNPF;//非节点载荷数intHZ;//载荷码intE;//单元码intfangchengshu;doubleF[303];//各节点等效总载荷数值intdym_jdm[100][2];//单元码对应的节点码:dym_jdm[][0],dym_jdm[][1]分为前后节点总码intzhichengweizhi[300];//记录支持节点作用点的数组intfjzhzuoyongdanyuan[100];//非节点载荷作用单元intfjzhleixing[100];//非节
3、点载荷类型:1-均布,2-垂直集中,3-平行集中,4-力偶,5-角度集中doublefjzhzhi[100];//非节点载荷的值doublefjzhzuoyongdian[100];//非节点载荷在各竿的作用点doublefjzhjiaodu[100];//非节点载荷作用角度intjdzhzuoyongdian[100];//节点载荷作用的节点数组doublejiedianzaihe[101][3];//节点载荷值,其jiedianzaihe[][0]--jiedianzaihe[][2]分别为U,V,Mdoublezhengtigangdu[303][303];//整体刚度数组doub
4、lechangdu[100];//各单元竿长数组doublejiaodu[100];//各单元角度数组doubletanxingmoliang[100];//各单元弹性模量数组doubleJ_moliang[100];//各单元J模量数组doublemianji[100];//各单元面积数组doubleweiyi[303];//记录各个节点位移的数组doubledy_weiyi[100][6];//各个单元在局部坐标系中的位移数组dy_weiyi[i][0]-dyweiyi[i][6]分别为第i+1单元的u1,v1,@1,u2,v2,@2doubledy_neili[100][6];//
5、各个单元在局部坐标系中的固端内力dy_weiyi[i][0]-dyweiyi[i][6]分别为第i+1单元的U1,V1,M1,U2,V2,M2doublegan_neili[100][6];//各个单元的竿端内力数组,gan_neili[i][6]表示第i+1单元的6内力.//*******************//一些处理过程的函数31//*******************voidinput();//数据的输入voidzonggang();//计算总刚度,存放于zhengtigangdu[][]数组中voidzongzaihe();//计算等效总节点载荷voidzhichengy
6、inru();//引入支承条件voidjsweiyi();//求各个节点位移voidjs_dy_weiyi();//求局部坐标系中的位移voidganduanneili();//求竿端内力voiddy_gangdu(inti,doubledg[6][6]);//求单元在局部坐标系中的单刚voidjs_T_T1(inti,doubleT[6][6],doubleT1[6][6]);//求单元的转换矩阵及其逆阵//************//主函数//************voidmain(){input();cout<<"输出结果"<7、＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝”<
