有限元三角形与悬臂梁

有限元c语言程序+报告（三角形与悬臂梁）可任意改变单元个数

//平面三角单元有限元//
#include“stdio.h“
#include“math.h“
double** dmake（int mint n）
{
	double **a;
	a=new double *[m];
	for（int i=0;i		a[i]=new double[n];
	return a;
}
void ddelete（double **aint mint n）
{
	for（int i=0;i		delete a[i];
	delete a;
}

int** imake（int mint n）
{
	int **a;
	a=new int *[m];
	for（int i=0;i		a[i]=new int[n];
	return a;
}
void idelete（int **aint mint n）
{
	for（int i=0;i		delete a[i];
	delete a;
}

void main（）
{
	//初始化；
	int linerow;
	//	double k[20][20]r[20];
	//	float x[10]y[10]eee1vvhqdip[20];
	float changgaoeee1vvhqe2;
	//	float e2str[3][8]strn[7][10];
	int iee;//ie[8][3];
	int npoin;
	float b[3]c[3]area;
	//	int in1ir[20]indid[20]iffix[2];
	int in1indiffix[2];
	int ijmnnnadd;
	float att;
	int nmaxnmax1;
	float smaxsmax1;
	int ii;
	
	FILE *fp*fp1*fp2*fp3*fp4*fp5;
	//输入数据；
	printf（“dd“）;
	fp=fopen（“data.txt““w“）;
	fp1=fopen（“init.dat““r“）;
	fp2=fopen（“coord.txt““w“）;
	fp3=fopen（“ie.txt““w“）;
	fp4=fopen（“ir.txt““w“）;
	fp5=fopen（“iffix.txt““w“）;
	printf（“please input the data.\n“）;
	printf（“npoin:\n“）;
	fscanf（fp1“%i“&line）;
	fscanf（fp1“%f“&chang）;
	gao=chang;
	printf（“COORDINATION:\n“）;
	int *jd=new int[line+1];
	int *dy=new int[line];
	dy[0]=0;jd[0]=0;jd[1]=1;
	for（i=1;i	{
		//fscanf（fp2“%f%f“&x[ii]&y[ii]）;
		jd[i+1]=jd[i]+i+1;
		dy[i]=dy[i-1]+2*i-1;
	}
	npoin=jd[line];
	iee=dy[line-1];
	printf（“the number of elements:%d\n“npoin）;
	printf（“the number of elements:%d\n“iee）;

	double **k**str**strn;
	k=dmake（2*npoin2*npoin）;
	str=dmake（3iee）;
	strn=dmake（7npoin）;
	int **ie;
	ie=imake（iee3）;
	double *r=new double[2*npoin];
	float *x=new float[npoin];
	float *y=new float[npoin];
	float *dip=new float[2*npoin];
	int *ir=new int[line];
	int *id=new int[2*npoin];

	printf（“the number of elements:\n“）;
	//fscanf（fp1“%i“&iee）;
	fscanf（fp1“%f“&ee）;
	printf（“E:%f\n“ee）;
	printf（“Poisson Rate:\n“）;
	fscanf（fp1“%f“&vv）;
	printf（“Poisson Rate:%f\n“vv）;
	
	//输入结点编号
	printf（“reading ie[i][j]...\n“）;
	for（ii=1;ii<=iee;ii++）
	{
		for（i=0;ii-dy[i]>0;i++）;
		i=i-1;
		ie[ii-1][0]=（ii-dy[i]+1）/2+jd[i];
		if（（ii-dy[i]）%2==1）
		{
			ie[ii-1][1]=（ii-dy[i]+1）/2+jd[i+1];
			ie[ii-1][2]=ie[ii-1][1]+1;
		}
		else
		{
			ie[ii-1][1]=（ii-dy[i]+1）/2+jd[i+1]+1;
			ie[ii-1][2]=ie[ii-1][0]+1;
		}
		fprintf（fp3“%i %i %i\n“ie[ii-1][0]ie[ii-1][1]ie[ii-1][2]）;
	}
	for（ii=1;ii<=npoin;ii++）
	{
		for（i=0;ii-jd[i]>0;i++）;
		x[ii-1]=（ii-jd[i-1]-1）*chang/（line-1）;
		y[ii-1]=（line-i）*gao/（line-1）;
		fprintf（fp2“%f %f\n“x[ii-1]y[ii-1]）;
	}
	fclose（fp2）;

	//总刚赋零；
	for（i=0;i<2*npoin;i++）
	{
		for（j=0;j<2*npoin;j++）
		{
			k[i][j]=0;
		}
	}
	//计算iee单元面积；
	for（i=0;i	{
		float x1x2x0y1y2y0;
		x0=x[ie[i][0]-1];
		x1=x[ie[i][1]-1];
		x2=x[ie[i][2]-1];
		
		y0=y[ie[i][0]-1];
		y1=y[ie[i][1]-1];
		y2=y

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件      92007  2010-06-27 00:39  有限元\截图\3.jpg

     文件      14919  2010-06-27 01:31  有限元\截图\4.jpg

     文件      54286  2010-06-27 01:48  有限元\截图\5.jpg

     文件      61014  2010-06-27 01:32  有限元\截图\截图01.jpg

     文件     196608  2010-08-15 20:10  有限元\报告\三角形.doc

     文件     312320  2010-08-15 20:10  有限元\报告\悬臂梁.doc

     文件        114  2010-08-15 20:01  有限元\程序\C语言程序\三角形\coord.txt

     文件       1841  2010-08-15 20:01  有限元\程序\C语言程序\三角形\data.txt

     文件         28  2010-08-15 20:01  有限元\程序\C语言程序\三角形\ie.txt

     文件         18  2010-08-15 20:01  有限元\程序\C语言程序\三角形\iffix.txt

     文件         33  2010-05-30 21:40  有限元\程序\C语言程序\三角形\init.dat

     文件          3  2010-08-15 20:01  有限元\程序\C语言程序\三角形\ir.txt

     文件       8845  2010-05-30 22:24  有限元\程序\C语言程序\三角形\main.cpp

     文件       3377  2010-06-24 21:16  有限元\程序\C语言程序\三角形\main.dsp

     文件        516  2010-06-24 21:16  有限元\程序\C语言程序\三角形\main.dsw

     文件      41984  2010-08-15 20:01  有限元\程序\C语言程序\三角形\main.ncb

     文件      48640  2010-08-15 20:01  有限元\程序\C语言程序\三角形\main.opt

     文件        242  2010-08-15 20:01  有限元\程序\C语言程序\三角形\main.plg

     文件        855  2010-08-15 20:01  有限元\程序\C语言程序\三角形\Release\coord.txt

     文件      29903  2010-08-15 20:01  有限元\程序\C语言程序\三角形\Release\data.txt

     文件        606  2010-08-15 20:01  有限元\程序\C语言程序\三角形\Release\ie.txt

     文件         68  2010-08-15 20:01  有限元\程序\C语言程序\三角形\Release\iffix.txt

     文件         33  2010-06-24 21:19  有限元\程序\C语言程序\三角形\Release\init.dat

     文件          3  2010-08-15 20:01  有限元\程序\C语言程序\三角形\Release\ir.txt

     文件      69632  2010-06-24 21:17  有限元\程序\C语言程序\三角形\Release\main.exe

     文件      17045  2010-06-24 21:17  有限元\程序\C语言程序\三角形\Release\main.obj

     文件     221704  2010-06-24 21:17  有限元\程序\C语言程序\三角形\Release\main.pch

     文件      33792  2010-06-24 21:17  有限元\程序\C语言程序\三角形\Release\vc60.idb

     文件       2052  2010-05-30 22:29  有限元\程序\C语言程序\悬臂梁\coord.txt

     文件     130240  2010-05-30 22:29  有限元\程序\C语言程序\悬臂梁\data.txt

............此处省略107个文件信息