有限元法大作业

有限元法大作业一平面刚架的程序用Visual C++编制的平面刚架的源程序如下：///////////////////////////////////////////////////////程序开始//////////////////////////////////////////////////////////////////#include"iostream.h"#include"math.h"#include"stdlib.h"#include"conio.h"//*****************//声明必要变量//*****************#define PI 3.141592654int NE; //单元数int NJ; //节点数int NZ; //支承数int NPJ; //有节点载荷作用的节点数int NPF; //非节点载荷数int HZ; //载荷码int E; //单元码int fangchengshu;double F[303]; //各节点等效总载荷数值int dym_jdm[100][2]; //单元码对应的节点码:dym_jdm[][0], dym_jdm[][1]分为前后节点总码int zhichengweizhi[300]; //记录支持节点作用点的数组int fjzhzuoyongdanyuan[100]; //非节点载荷作用单元int fjzhleixing[100]; //非节点载荷类型:1-均布,2-垂直集中,3-平行集中,4-力偶,5-角度集中double fjzhzhi[100]; //非节点载荷的值double fjzhzuoyongdian[100]; //非节点载荷在各竿的作用点double fjzhjiaodu[100]; //非节点载荷作用角度int jdzhzuoyongdian[100]; //节点载荷作用的节点数组double jiedianzaihe[101][3];//节点载荷值,其jiedianzaihe[][0]-- jiedianzaihe[][2]分别为U, V, Mdouble zhengtigangdu[303][303]; //整体刚度数组double changdu[100]; //各单元竿长数组double jiaodu[100]; //各单元角度数组double tanxingmoliang[100]; //各单元弹性模量数组double J_moliang[100]; //各单元J模量数组double mianji[100]; //各单元面积数组double weiyi[303]; //记录各个节点位移的数组double dy_weiyi[100][6]; //各个单元在局部坐标系中的位移数组dy_weiyi[i][0]-dyweiyi[i][6]分别为第i+1单元的u1,v1,@1,u2,v2,@2double dy_neili[100][6];//各个单元在局部坐标系中的固端内力dy_weiyi[i][0]-dyweiyi[i][6]分别为第i+1单元的U1,V1,M1,U2,V2,M2double gan_neili[100][6];//各个单元的竿端内力数组,gan_neili[i][6]表示第i+1单元的6内力.//*******************//*******************void input(); //数据的输入void zonggang(); //计算总刚度,存放于zhengtigangdu[][]数组中void zongzaihe(); //计算等效总节点载荷void zhichengyinru(); //引入支承条件void jsweiyi(); //求各个节点位移void js_dy_weiyi(); //求局部坐标系中的位移void ganduanneili(); //求竿端内力void dy_gangdu(int i,double dg[6][6]); //求单元在局部坐标系中的单刚void js_T_T1(int i,double T[6][6],double T1[6][6]); //求单元的转换矩阵及其逆阵//************//主函数//************void main(){input();cout<<" 输出结果"<<endl;cout<<"＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝”<<endl;zonggang();zongzaihe();zhichengyinru();jsweiyi();cout<<"输出位移:"<<endl;for( int i=0;i<NJ;i++){cout<<"u["<<i+1<<"]:"<<weiyi[i*3]<<endl;cout<<"v["<<i+1<<"]:"<<weiyi[i*3+1]<<endl;cout<<"@["<<i+1<<"]:"<<weiyi[i*3+2]<<endl;}js_dy_weiyi();ganduanneili();//输出竿端内力cout<<"＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝”<<endl;cout<<"输出竿内力:"<<endl;for( i=0;i<NE;i++){cout<<"第"<<i+1<<"单元:"<<endl;cout<<" U1:"<<gan_neili[i][0]<<" V1:"<<gan_neili[i][01]<<" M1:"<<gan_neili[i][2]<<endl;cout<<" U2:"<<gan_neili[i][3]<<" V2:"<<gan_neili[i][4]<<" M2:"<<gan_neili[i][5]<<endl;}cout<<" ==========================================="<<endl;cout<<"输出完毕,按任意键退出！"<<endl;getchi();}//***************************//***************************void dy_gangdu(int i,double dg[6][6]){dg[0][0]=tanxingmoliang[i]*mianji[i]/changdu[i]; dg[0][1]=0; dg[0][2]=0;dg[0][3]=-tanxingmoliang[i]*mianji[i]/changdu[i]; dg[0][4]=0; dg[0][5]=0; dg[1][0]=0;dg[1][1]=12*tanxingmoliang[i]*J_moliang[i]/(changdu[i]*changdu[i]*changdu[i]);dg[1][2]=-6*tanxingmoliang[i]*J_moliang[i]/(changdu[i]*changdu[i]); dg[1][3]=0;dg[1][4]=-12*tanxingmoliang[i]*J_moliang[i]/(changdu[i]*changdu[i]*changdu[i]);dg[1][5]=-6*tanxingmoliang[i]*J_moliang[i]/(changdu[i]*changdu[i]); dg[2][0]=0;dg[2][1]=-6*tanxingmoliang[i]*J_moliang[i]/(changdu[i]*changdu[i]);dg[2][2]=4*tanxingmoliang[i]*J_moliang[i]/changdu[i]; dg[2][3]=0;dg[2][4]=6*tanxingmoliang[i]*J_moliang[i]/(changdu[i]*changdu[i]);dg[2][5]=2*tanxingmoliang[i]*J_moliang[i]/changdu[i];dg[3][0]=-tanxingmoliang[i]*mianji[i]/changdu[i]; dg[3][1]=0; dg[3][2]=0;dg[3][3]=tanxingmoliang[i]*mianji[i]/changdu[i]; dg[3][4]=0; dg[3][5]=0; dg[4][0]=0;dg[4][1]=-12*tanxingmoliang[i]*J_moliang[i]/(changdu[i]*changdu[i]*changdu[i]);dg[4][2]=6*tanxingmoliang[i]*J_moliang[i]/(changdu[i]*changdu[i]); dg[4][3]=0;dg[4][4]=12*tanxingmoliang[i]*J_moliang[i]/(changdu[i]*changdu[i]*changdu[i]);dg[4][5]=6*tanxingmoliang[i]*J_moliang[i]/(changdu[i]*changdu[i]); dg[5][0]=0;dg[5][1]=-6*tanxingmoliang[i]*J_moliang[i]/(changdu[i]*changdu[i]);dg[5][2]=2*tanxingmoliang[i]*J_moliang[i]/changdu[i]; dg[5][3]=0;dg[5][4]=6*tanxingmoliang[i]*J_moliang[i]/(changdu[i]*changdu[i]);dg[5][5]=4*tanxingmoliang[i]*J_moliang[i]/changdu[i];}//*************************//求单元的转换矩阵及其逆阵//*************************void js_T_T1(int i,double T[6][6],double T1[6][6]){int Tti,Tti2;for( Tti=0;Tti<6;Tti++){//计算单元的转换矩阵Tfor( Tti2=0;Tti2<6;Tti2++){if(Tti==0&&Tti2==0)T[Tti][Tti2]=cos(jiaodu[i]);else if(Tti==0&&Tti2==1)T[Tti][Tti2]=sin(jiaodu[i]);else if(Tti==1&&Tti2==0)T[Tti][Tti2]=-sin(jiaodu[i]);else if(Tti==1&&Tti2==1)T[Tti][Tti2]=cos(jiaodu[i]);else if(Tti==2&&Tti2==2)T[Tti][Tti2]=1;T[Tti][Tti2]=cos(jiaodu[i]);else if(Tti==3&&Tti2==4)T[Tti][Tti2]=sin(jiaodu[i]);else if(Tti==4&&Tti2==3)T[Tti][Tti2]=-sin(jiaodu[i]);else if(Tti==4&&Tti2==4)T[Tti][Tti2]=cos(jiaodu[i]);else if(Tti==5&&Tti2==5)T[Tti][Tti2]=1;elseT[Tti][Tti2]=0;if(fabs(T[Tti][Tti2])<0.0000001)T[Tti][Tti2]=0;if(T[Tti][Tti2]>0.99999999)T[Tti][Tti2]=1;if(T[Tti][Tti2]<-0.99999999)T[Tti][Tti2]=-1;}}//计算T的转置矩阵T1for(Tti=0;Tti<6;Tti++){for(Tti2=0;Tti2<6;Tti2++)T1[Tti2][Tti]=T[Tti][Tti2];}}//*******************//计算竿端内力的函数//*******************void ganduanneili(){int Ti,Ti2;double chengji=0;double danyuan_gangdu[6][6]; //单元刚度数组for(int i=0;i<NE;i++){dy_gangdu(i,danyuan_gangdu);for(Ti=0;Ti<6;Ti++){for(Ti2=0;Ti2<6;Ti2++)chengji=chengji+danyuan_gangdu[Ti][Ti2]*dy_weiyi[i][Ti2];gan_neili[i][Ti]=chengji+dy_neili[i][Ti];chengji=0;}}//************************************************************//计算各个单元在局部坐标系中的位移的函数,存放在dy_weiyi[][]中//************************************************************void js_dy_weiyi(){int Ti,Ti3;double T[6][6],T1[6][6];double chengji=0;for(int i=0;i<NE;i++){js_T_T1(i,T,T1);int qianjiedian,houjiedian;qianjiedian=dym_jdm[i][0];houjiedian=dym_jdm[i][1];for(Ti=0;Ti<6;Ti++){dy_weiyi[i][Ti]=T[Ti][0]*weiyi[(qianjiedian-1)*3]+T[Ti][1]*weiyi[(qianjiedian-1)*3+1]+T[Ti][2]* weiyi[(qianjiedian-1)*3+2]+T[Ti][3]*weiyi[(houjiedian-1)*3]+T[Ti][4]*weiyi[(houjiedian-1)*3+1]+T[Ti][5]* weiyi[(houjiedian-1)*3+2];}}}//***********************//计算得位移的函数//***********************void jsweiyi(){int Ti,Ti2;int i,j,k;int i_max;double max,zhongjian;double gd_kz[303][303];for(Ti=0;Ti<NJ*3;Ti++){for(Ti2=0;Ti2<NJ*3;Ti2++)gd_kz[Ti][Ti2]=zhengtigangdu[Ti][Ti2];gd_kz[Ti][NJ*3]=F[Ti];}//用高斯主列消去法解方程,先计算增广矩阵的变形矩阵for(Ti=0;Ti<NJ*3-1;Ti++){max=fabs(gd_kz[Ti][Ti]);i_max=Ti;for(Ti2=Ti;Ti2<NJ*3;Ti2++){{i_max=Ti2-1;max=gd_kz[Ti2][Ti];}}if(i_max!=Ti){for(i=Ti;i<NJ*3+1;i++){zhongjian=gd_kz[Ti][i];gd_kz[Ti][i]=gd_kz[i_max][i];gd_kz[i_max][i]=zhongjian;}}for(j=Ti+1;j<NJ*3;j++){double shang=gd_kz[j][Ti]/gd_kz[Ti][Ti];for(k=0;k<NJ*3+1;k++)gd_kz[j][k]=-shang*gd_kz[Ti][k]+gd_kz[j][k];}}//反代求值double chengji=0;weiyi[NJ*3-1]=gd_kz[NJ*3-1][NJ*3]/gd_kz[NJ*3-1][NJ*3-1];for(i=NJ*3-2;i>=0;i--){for(j=i+1;j<NJ*3;j++)chengji=chengji+gd_kz[i][j]*weiyi[j];weiyi[i]=(gd_kz[i][NJ*3]-chengji)/gd_kz[i][i];chengji=0;}}//********************//引入支承条件的函数//********************void zhichengyinru(){int i,Ti;int weizhi;fangchengshu=NJ*3-NZ;for(i=0;i<NZ;i++){weizhi=zhichengweizhi[i]-1;for(Ti=0;Ti<NJ*3;Ti++)zhengtigangdu[weizhi][Ti]=0;zhengtigangdu[Ti][weizhi]=0;zhengtigangdu[weizhi][weizhi]=1;F[weizhi]=0;}}//*************************//计算等效总节点载荷的函数//*************************void zongzaihe(){int i,Ti;double U1=0,V1=0,M1=0,U2=0,V2=0,M2=0;double c,G,l,d;double Ff[6];double T[6][6],T1[6][6]; //T为坐标转换矩阵,T1为T的转置矩阵for(i=0;i<NJ*3;i++)F[i]=0;for(i=0;i<6;i++)Ff[i]=0;for( i=0;i<NE;i++){for(int jl=0;jl<6;jl++)dy_weiyi[i][jl]=0;}for(i=0;i<NPF;i++){int dym=fjzhzuoyongdanyuan[i];U1=V1=M1=U2=V2=M2=0;for(int t=0;t<6;t++)Ff[t]=0;c=fjzhzuoyongdian[i];l=changdu[(fjzhzuoyongdanyuan[i]-1)];d=l-c;G=fjzhzhi[i];if(fjzhleixing[i]==1){V1=-0.5*G*c*(2-2*c*c/(l*l)+c*c*c/(l*l*l));V2=-G*c-V1;M1=G*c*c*(6-8*c/l+3*c*c/(l*l))/12;M2=-G*c*c*c*(4-3*c/l)/(12*l);U1=U2=0;}else if(fjzhleixing[i]==2){V1=-G*(l+2*c)*d*d/(l*l*l);M1=G*c*d*d/(l*l);M2=-G*c*c*d/(l*l);U1=U2=0;}else if(fjzhleixing[i]==3){U1=-G*d/l;U2=-G*c/l;M1=M2=V1=V2=0;}else if(fjzhleixing[i]==4){U1=U2=0;V1=-6*G*c*d/(l*l*l);V2=-V1;M1=G*d*(2-3*d/l)/l;M2=G*c*(2-3*c/l)/l;}else if(fjzhleixing[i]==5){V1=-G*sin(fjzhjiaodu[i])*(l+2*c)*d*d/(l*l*l);V2=-G*sin(fjzhjiaodu[i])*(l+2*d)*c*c/(l*l*l);M1=G*sin(fjzhjiaodu[i])*c*d*d/(l*l);M2=-G*sin(fjzhjiaodu[i])*c*c*d/(l*l);U1=U2=0;U1=-G*cos(fjzhjiaodu[i])*d/l+U1;U2=-G*cos(fjzhjiaodu[i])*c/l+U2;}//记录竿的固端内力dy_neili[dym-1][0]=U1;dy_neili[dym-1][1]=V1;dy_neili[dym-1][2]=M1;dy_neili[dym-1][3]=U2;dy_neili[dym-1][4]=V2;dy_neili[dym-1][5]=M2;js_T_T1( fjzhzuoyongdanyuan[i]-1, T,T1);for(Ti=0;Ti<6;Ti++)Ff[Ti]=-(T1[Ti][0]*U1+T1[Ti][1]*V1+T1[Ti][2]*M1+T1[Ti][3]*U2+T1[Ti][4]*V2+T1[Ti][5]*M2);//将载荷转换到整体坐标系中F[(dym_jdm[dym-1][0]-1)*3]=F[(dym_jdm[dym-1][0]-1)*3]+Ff[0];F[(dym_jdm[dym-1][0]-1)*3+1]=F[(dym_jdm[dym-1][0]-1)*3+1]+Ff[1];F[(dym_jdm[dym-1][0]-1)*3+2]=F[(dym_jdm[dym-1][0]-1)*3+2]+Ff[2];F[(dym_jdm[dym-1][1]-1)*3]=F[(dym_jdm[dym-1][1]-1)*3]+Ff[3];F[(dym_jdm[dym-1][1]-1)*3+1]=F[(dym_jdm[dym-1][1]-1)*3+1]+Ff[4];F[(dym_jdm[dym-1][1]-1)*3+2]=F[(dym_jdm[dym-1][1]-1)*3+2]+Ff[5];for(i=0;i<NPJ;i++){F[(jdzhzuoyongdian[i]-1)*3]=F[(jdzhzuoyongdian[i]-1)*3]+jiedianzaihe[i][0];F[(jdzhzuoyongdian[i]-1)*3+1]=F[(jdzhzuoyongdian[i]-1)*3+1]+jiedianzaihe[i][1];F[(jdzhzuoyongdian[i]-1)*3+2]=F[(jdzhzuoyongdian[i]-1)*3+2]+jiedianzaihe[i][2];}}//***********************************************//计算总刚度的函数,存放于zhengtigangdu[][]数组中//***********************************************void zonggang(){double gd[6][6],gd2[6][6]; // gd1为中间变量.double T[6][6],T1[6][6]; //T为坐标转换矩阵,T1为T的转置矩阵double danyuan_gangdu[6][6];int Ti,Ti2,Ti3;double chengji=0; //中间变量for(Ti=0;Ti<NJ*3;Ti++) //先将整体刚度付0{for(Ti2=0;Ti2<NJ*3;Ti2++)zhengtigangdu[Ti][Ti2]=0;}for(int i=0;i<NE;i++){dy_gangdu(i,danyuan_gangdu); //局部坐标系中的单元刚度js_T_T1( i, T,T1); //计算单元的转换矩阵T//以下为计算单元在整体坐标系中的单刚存在gd[]数组中.for(Ti=0;Ti<6;Ti++){for(Ti2=0;Ti2<6;Ti2++){gd[Ti][Ti2]=0;gd2[Ti][Ti2]=0;}}chengji=0;for(Ti=0;Ti<6;Ti++){for(Ti2=0;Ti2<6;Ti2++){for(Ti3=0;Ti3<6;Ti3++)chengji=chengji+T1[Ti][Ti3]*danyuan_gangdu[Ti3][Ti2]; //改过gd2[Ti][Ti2]=chengji;chengji=0;}chengji=0;for(Ti=0;Ti<6;Ti++){for(Ti2=0;Ti2<6;Ti2++){for(Ti3=0;Ti3<6;Ti3++)chengji=chengji+gd2[Ti][Ti3]*T[Ti3][Ti2]; //改过gd[Ti][Ti2]=chengji;chengji=0;}}chengji=0;//以下为将每个单元刚度集成到整体刚度中计算整体刚度int qianjiedian,houjiedian;qianjiedian=dym_jdm[i][0];houjiedian=dym_jdm[i][1];for(Ti=0;Ti<3;Ti++){for(Ti2=0;Ti2<3;Ti2++)zhengtigangdu[(qianjiedian-1)*3+Ti][(qianjiedian-1)*3+Ti2]=zhengtigangdu[(qianjiedian-1)*3+Ti][(qian jiedian-1)*3+Ti2]+gd[Ti][Ti2];}for(Ti=0;Ti<3;Ti++){for(Ti2=0;Ti2<3;Ti2++)zhengtigangdu[(houjiedian-1)*3+Ti][(houjiedian-1)*3+Ti2]=zhengtigangdu[(houjiedian-1)*3+Ti][(houjie dian-1)*3+Ti2]+gd[Ti+3][Ti2+3];}for(Ti=0;Ti<3;Ti++){for(Ti2=0;Ti2<3;Ti2++)zhengtigangdu[(qianjiedian-1)*3+Ti][(houjiedian-1)*3+Ti2]=zhengtigangdu[(qianjiedian-1)*3+Ti][(houji edian-1)*3+Ti2]+gd[Ti][Ti2+3];}for(Ti=0;Ti<3;Ti++){for(Ti2=0;Ti2<3;Ti2++)zhengtigangdu[(houjiedian-1)*3+Ti][(qianjiedian-1)*3+Ti2]=zhengtigangdu[(houjiedian-1)*3+Ti][(qianji edian-1)*3+Ti2]+gd[Ti+3][Ti2];}}}//**************************//输入数据的函数//**************************{char queding;//基本参数输入cout<<"请输入基本参数:单元数,节点数,支持数,有节点载荷的节点数,非节点载荷数;"<<endl<<endl;cin>>NE>>NJ>>NZ>>NPJ>>NPF;cout<<endl<<endl;//输入单元几何参数cout<<"请输入单元的几何参数:"<<endl;for(int i=0;i<NE;i++){cout<<"第"<<i+1<<"单元"<<endl;cout<<"长度:"; cin>>changdu[i];cout<<"角度:"; cin>>jiaodu[i]; jiaodu[i]=jiaodu[i]*PI/180;cout<<"弹性模量E:"; cin>>tanxingmoliang[i];cout<<"J模量:"; cin>>J_moliang[i];cout<<"面积:"; cin>>mianji[i];cout<<endl;//cout<<"第"<<i+1<<"单元长度:"<<changdu[i]<<" 角度:"<<jiaodu[i]<<" 弹性模量E:"<<tanxingmoliang[i]<<" J量:"<<J_moliang[i]<<" 面积:"<<mianji[i]<<endl<<endl;}//输入单元对应的节点总码cout<<"请输入个单元对应的节点总码"<<endl;for(i=0;i<NE;i++){cout<<"第"<<i+1<<"单元前面节点的总码:"; cin>>dym_jdm[i][0];cout<<"第"<<i+1<<"单元后面节点的总码:"; cin>>dym_jdm[i][1];}cout<<endl;//输入支承点位置cout<<"请输入支承点位置:"<<endl;for(i=0;i<NZ;i++){cout<<"第"<<i+1<<"支承位:"; cin>>zhichengweizhi[i];}cout<<endl;//输入节点载荷值cout<<"请输入节点载荷作用节点及其载荷值:"<<endl;for(i=0;i<NPJ;i++){cout<<"载荷所在节点位置:"; cin>>jdzhzuoyongdian[i];}for(i=0;i<NPJ;i++){cout<<"第"<<jdzhzuoyongdian[i]<<"节点的载荷值:"<<endl;cout<<"U:"; cin>>jiedianzaihe[i][0];cout<<"V:"; cin>>jiedianzaihe[i][1];cout<<"M:"; cin>>jiedianzaihe[i][2];}cout<<endl;cout<<"请输入非节点载荷作单元,载荷类型,载荷值,作用点(指离第一个节点的距离),作用角度:"<<endl;cout<<"其中1-均布,2-垂直集中,3-平行集中,4-力偶,5-角度集中"<<endl;for(i=0;i<NPF;i++){cout<<"第"<<i+1<<"个非节点载荷作用单元:"; cin>>fjzhzuoyongdanyuan[i];cout<<"载荷类型:"; cin>>fjzhleixing[i];cout<<"载荷值(如果是类型5的力,只输入力的绝对值) :"; cin>>fjzhzhi[i];cout<<"离第一个节点的距离:"; cin>>fjzhzuoyongdian[i];cout<<"力在该单元坐标中的作用角度:"; cin>>fjzhjiaodu[i];fjzhjiaodu[i]=fjzhjiaodu[i]*PI/180; cout<<endl;}cout<< endl<<endl;cout<<"按任意键显示输入情况:"<<endl;getch();//显示输入的情况cout<<"输入的数据如下:"<<endl;cout<<"===========================================:"<<endl;cout<<"单元数:"<<NE<<" 节点数:"<<NJ<<" 支持数:"<<NZ<<"有节点载荷的节点数:"<<NPJ<<" 非节点载荷数:"<<NPF<<endl;cout<<"============================================"<<endl;cout<<"各单元几何参数如下:"<<endl;cout<<"单元"<<"\t"<<"单元长度:"<<"\t"<<"角度: "<<"\t"<<"弹性模量E:"<<"\t"<<"J量: "<<"\t"<<"面积:"<<endl;for(i=0;i<NE;i++)cout<<i+1<<"\t"<<changdu[i]<<"\t\t"<<jiaodu[i]<<"\t\t"<<tanxingmoliang[i]<<"\t\t"<<J_moliang[i]<<"\ t\t"<<mianji[i]<<endl;cout<<"==========================================="<<endl;cout<<"各单元对应节点总码如下:"<<endl;cout<<"单元号"<<"\t"<<"前节点总码"<<"\t"<<"后节点总码"<<endl;for(i=0;i<NE;i++)cout<<i+1<<"\t\t"<<dym_jdm[i][0]<<"\t\t"<<dym_jdm[i][1]<<endl;cout<<" ==========================================="<<endl;cout<<"支承位置如下:"<<endl;cout<<"支承位置"<<"\t";for(i=0;i<NZ;i++)cout<<zhichengweizhi[i]<<" ";cout<<endl;cout<<" ==========================================="<<endl;cout<<"节点载荷输入如下:"<<endl;cout<<"节点位置"<<"\t"<<"U "<<"\t\t"<<"V "<<"\t\t"<<"M "<<endl;for(i=0;i<NPJ;i++){cout<<jdzhzuoyongdian[i]<<"\t\t"<<jiedianzaihe[i][0]<<"\t\t"<<jiedianzaihe[i][1]<<"\t\t"<<jiedianzaihe[ i][2];}cout<<" ==========================================="<<endl;cout<<"非节点载荷的作用情况如下:"<<endl;cout<<"作用单元"<<"\t"<<"载荷类型"<<"\t"<<"载荷值"<<"\t\t"<<"作用点"<<"\t\t"<<"作用角度"<<endl;for(i=0;i<NPF;i++)cout<<fjzhzuoyongdanyuan[i]<<"\t\t"<<fjzhleixing[i]<<"\t\t"<<fjzhzhi[i]<<"\t\t"<<fjzhzuoyongdian[i]<< "\t\t"<<fjzhjiaodu[i]<<endl;cout<<" ==========================================="<<endl;cout<<"按任意键继续！"<<endl;getch();cout<<endl;}///////////////////////////////////////////////////////程序结束////////////////////////////////////////////////////////////////// 用以上程序计算教材例题3-7以及3-8结果都和例题的结果一样。