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

有限元编程的C++实现算例
n. com/dow nl oads76/doc/fileformat/290377/ganjia

#in clude<>
#in clude<>
3.
4.
#defi ne ne
3
//单元数
#defi ne nj
4
//节点数
#defi ne nz
6
//支撑数
#defi ne npj
0
//节点载荷数
#defi ne npf
1
//非节点载荷数
#defi ne nj3
12
//节点位移总数
#defi ne dd
6
//半带宽
#defi ne e0

//弹性模量
#defi ne a0

//截面积
#defi ne i0
-4
//单元惯性距
#defi ne pi

16.
17.
i nt
/*gghjghg*/
jm[ ne+1][3]={{0,0,0},{0,1,2},{0,2,3},{0,4,3}};
double
gc[ ne+1]={,,,};
double gj[ ne+1]={,,,};
double mj[ ne+1]={,a0,a0,a0};
double gx[ ne+1]={,i0,i0,i0};
i nt zc[ nz+1]={0,1,2,3,10,11,12};
double pj[ npj+1][3]={{,,}};
double pf[n pf+1][5]={{0,0,0,0,0},{0,-20,,,}};
double kz[nj3+1][dd+1],p[nj3+1];
double pe[7],f[7],f0[7],t[7][7];
double ke[7][7],kd[7][7];
29.
30.
//**kz [][] —整体刚度矩阵
//**ke[][] —整体坐标下的单元刚度矩阵
//**kd[][] —局部坐标下的单位刚度矩阵
//**t[][] —坐标变换矩阵
//** 这是函数声明
void jdugd(i nt);
void zb(i nt);
void gdn 1( in t);
void dugd(i nt);
//** 主程序开始
void mai n()
{
int i,j,k,e,dh,h,ii,jj,hz,al,bl,m,l,dl,zl,z,jO;
double cl,wy[7];
int im,i n,jn;
//*
**********************************************
//< 功能：形成矩阵P>
//*
**********************************************
if(n pj>0)
{
for(i=1;i<=n pj;i++)
{
//
j=pj[i][2];
p[j]=pj[i][1];
}
}
if(n pf>0)
把节点载荷送入P
for(i=1;i<=n pf;i++)
35.
36.
37.
38.
39.
40.
41.
42.
43.
44.
45.
46.
47.
48.
49.
50.
51.
52.
53.
54.
55.
56.
57.
58.
59.
60.
61.
62.
63.
64.
66.
67.
68.
69.
70.
71.
72.
73.
74.
75.
76.
77.
78.
79.
80.
81.
82.
83.
84.
85.
86.
87.
88.
89.
90.
91.
//*
//<功能：生成整体刚度矩阵kz[][]>
//求固端反力F0
hz=i;
gdn l(hz);
e=(i nt)pf[hz][3];
zb(e);
//求单元号码
for(j=1;j<=6;j++)
//求坐标变换矩阵T
{
pe[j]=;
for(k=1;k<=6;k++)
//求等效节点