有限元C++程序.doc

下载文档 降价啦

203
0
约1.68万字
约 15页
2017-11-27 发布于江西
举报
版权申诉
保障服务

有限元C++程序.doc

1、本文档共15页，可阅读全部内容。
2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击这里二次下载。
4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。

有限元C程序

#include iostream.h #include fstream.h #include iomanip.h #include math.h class Q8 { public: int NN, NF, NE, KU, KV, KRX, KRY, KQ, LK; int JEW[9], JLL[17], JE[9][501], JD[3006], KPQ[3][50], JRX[50], JRY[50], JU[50], JV[50]; double H[3+1], D[3+1], U[3007], FF[17], FS[6][6][9], FC[6][6][9], FA[6][6][9], B[5][17], C[5][17], XO[1504], YO[1504], X[50],Y[50],SK[5000], PQ[5][50], RX[50], RY[50], US[50], VS[50],EK[137]; double PO, EO, RU, XC, XA, YC, YA, YN, QQ, VY, RPM; public: /** * calculate */ void calculate() { NF = NN + NN; //NN-节点个数 NF-节点位移（载荷）个数 //弹性矩阵[D] if(D[0]==0){ D[1] = EO / (1.0 - PO * PO); D[2] = D[1] * PO; D[3] = D[1] * (1.0 - PO) / 2; } if(D[0]==2){ D[1] = EO * (1.0 - PO) / (1.0 + PO) / (1.0 - PO - PO); D[2] = D[1] * PO / (1.0 - PO); D[3] = D[1] * (1.0 - PO - PO) / (1.0 - PO) / 2; } SKDD(); NCNA(); VY = RU * (3.1415926 * RPM / 30) / 9.80; if(H[0]==2){ H[1]=1; H[2]=1; H[3]=1; } if(H[0]==3){ H[1] = 5.0 / 9.0; H[2] = 8.0 / 9.0; H[3] = H[1]; } //初始化矩阵 for(int i=0;i=NF;i++){ U[i]=0; } for (i = 0; i =LK; i++) { SK[i]=0; } for (i = 1; i = NE; i++) { int m, l1, l2; for (int j = 1; j = 8; j++) { m = JE[j][i]; l1 = j + j; l2 = l1 - 1; JLL[l2] = m + m - 1; JLL[l1] = m + m; } FEKP(i); SKKE(); } RIGHT(); for (i = 1; i = KU; i++) { FIXD(2 * JU[i] - 1, US[i]); } if (KV != 0) { for (int j = 1; j =KV; j++) { FIXD(2 * JV[j], VS[j]); } } solve(); } /** * XYCA */ void XYCA(int i, int j,int n) { double R, S, F, D; int L; XC = 0; XA = 0; YC = 0; YA = 0; YN = 0;