//平面三角單元有限元//
#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個文件信息