作业有限元(必威体育精装版整理版).pdfVIP

作业1：求截面面积F、形心位置y0、x0，惯性矩Ix0、Iy0

一、流程图：

开始

输入各节点

的X、Y坐标

将坐标原点移至第

一节点

i从1到n-1循环

计算单个三角形的

面积、面积矩、形

心位置，惯性矩，

并叠加面积、面积

矩和惯性矩

计算截面面积、形

心坐标和惯性矩

将坐标原点移回原

来位置

数据输出

结束

二、变量说明：

n—截面上的节点个数；

X[100]、Y[100]—用来存放各节点的X、Y坐标；

f—单个三角形的面积；

y0—截面形心到x轴距离；

x0—截面形心到y轴距离；

Ix0、Iy0—截面对形心的惯性矩；

F—截面的面积；

yf—单个三角形的形心y坐标；

xf—单个三角形的形心x坐标；

jx—单个三角形对x轴的惯性矩；

jy—单个三角形对y轴的惯性矩；

Sx—对x轴的静面矩；

Sy—对y轴的静面矩。

三、源程序代码：

#includeiostream.h

#includemath.h

voidmain()

{

intn,i;

cout请输入节点个数:;

cinn;

doubleX[100],Y[100];

for(i=0;in;i++)

{

cout请输入第i+1个节点的X坐标:;

cinX[i];

cout请输入第i+1个节点的Y坐标:;

cinY[i];

}

doublef,y0,x0,Ix0,Iy0,F=0.0,yf,xf,jx,jy,Sx=0.0,Sy=0.0;

//将坐标原点移至第一点；

for(i=1;in;i++)

{

X[i]=X[i]-X[0];

Y[i]=Y[i]-Y[0];

}

for(i=1;in-1;i++)

{

f=(X[i]*Y[i+1]-Y[i]*X[i+1])/2.0;

yf=(Y[i]+Y[i+1])/3.0;

xf=(X[i]+X[i+1])/3.0;

jx=f*(Y[i]*Y[i]+Y[i]*Y[i+1]+Y[i+1]*Y[i+1])/6.0;

jy=f*(Y[i]*X[i]+X[i]*X[i+1]+X[i+1]*X[i+1])/6.0;

F=F+f;

Sx=Sx+f*yf