原创力文档- 1、本文档共13页，可阅读全部内容。
- 2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
- 3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
查看更多
作业有限元
作业1：求截面面积F、形心位置y0、x0，惯性矩Ix0、Iy0
一、流程图：
二、变量说明：
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
void main()
{
int n,i;
cout请输入节点个数:;
cinn;
double X[100],Y[100];
for(i=0;in;i++)
{
cout请输入第i+1个节点的X坐标:;
cinX[i];
cout请输入第i+1个节点的Y坐标:;
cinY[i];
}
double f,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;
Sy=Sy+f*xf;
}
y0=Sx/F;
x0=Sy/F;
Ix0=jx-F*y0*y0;
Iy0=jy-F*x0*x0;
y0=y0+Y[0];
x0=x0+X[0];
cout**********计算结果**********/n;
coutF=F;
couty0=y0;
coutx0=x0;
coutIx0=Ix0;
coutIy0=Iy0endl;
cout****************************endl;
}
四、算例：
如上图所示：节点数N=22，节点编号如图所示。
各节点的坐标如下：
计算结果如下：
作业2：将BASIC程序改为FORTRAN程序（或C程序）
流程图
变量说明：
L — 跨度；
LB — 单主梁宽；
RB — 人行道宽度；
XB — 行车道宽度；
DX — 插值函数宽度；
NCD — 车道数；
YE[350] — 用来存放插值点坐标；
N — 桥面板分块数；
MM — 主梁数；
NN — 插值点数；
WIJ[15][350] — 用来存放挠度值；
ETAIJ[15][350] — 用来存放影响线竖向坐标；
MCQ[15] — 用来存放汽车荷载横向分布系数；
MCL[15] — 用来存放履带车荷载横向分布系数；
MCG[15] — 用来存放挂车荷载横向分布系数；
MCR[15] — 用来存放人群荷载横向分布系数；
Y[350] — 用来存放主梁Y坐标；
DE — 弹性模量；
G — 剪切模量；
MC[4] — 记录四条车道的初始分布系数；
JX — 纵向单宽抗弯惯性矩；
JY — 横向单宽抗弯惯性矩；
JTXY — 纵横向抗扭惯性矩之和。
源程序代码：
#includeiostream.h
#includemath.h
#includeiomanip.h
void main()
{
int i,N,MM=7,NN,NCD,k,m,j,IEND,IK1,IK2,IK,IL,J1,J;
double YE[350],Y[15],WIJ[15][350],ETAIJ[15][350],WB[350],MC[4],MCQ[15],MCG[15],MCL[15],MCR[15];
double DX,L,LB,RB,XB,JX,JY,JTXY,DE=1.0,DB,ALPHA,A,B,THETA,OMEGA,PHI,ETA,GAMMA,DELTA,DM,DN,C1,C2,C3,C4;
double PSI,D1,D2,CC,SS,SC
文档评论（0）