有限元编程的c++实现算例.doc

有限元编程的c++实现算例http://read./downloads76/doc/fileformat/290377/  1.#include<>    2.#include<>    3.     4.     5.#definene3                                                  //单元数    6.#definenj4                                                  //节点数    7.#definenz6                                                  //支撑数    8.#definenpj0                                                 //节点载荷数    9.#definenpf1                                                 //非节点载荷数   10.#definenj312                                                //节点位移总数   11.#definedd6                                                  //半带宽   12.#                                             //弹性模量   13.#                                              //截面积   14.#-4                                            //单元惯性距    15.#                                   16.    17.    [ne+1][3]={{0,0,0},{0,1,2},{0,2,3},{0,4,3}};            /*gghjghg*/                [ne+1]={,,,};                                [ne+1]={,,,};   [ne+1]={,a0,a0,a0};   [ne+1]={,i0,i0,i0};   [nz+1]={0,1,2,3,10,11,12};   [npj+1][3]={{,,}};   [npf+1][5]={{0,0,0,0,0},{0,-20,,,}};   [nj3+1][dd+1],p[nj3+1];   [7],f[7],f0[7],t[7][7];   [7][7],kd[7][7];   29.    30.    31.//**kz[][]—整体刚度矩阵   32.//**ke[][]—整体坐标下的单元刚度矩阵   33.//**kd[][]—局部坐标下的单位刚度矩阵   34.//**t[][]—坐标变换矩阵   35.    36.//**这是函数声明   (int);   (int);   (int);   (int);   41.    42.    43.//**主程序开始   ()   45.{   46. inti,j,k,e,dh,h,ii,jj,hz,al,bl,m,l,dl,zl,z,j0;   47. doublecl,wy[7];   48. intim,in,jn;   49.    50.//***********************************************   51.//<功能:形成矩阵P>   52.//***********************************************   53.    54. if(npj>0)   55. {   56.     for(i=1;i<=npj;i++