分销赏收藏举报申诉 / 55

立即下载开通VIP

当前位置：首页 > 学术论文 > 其他 > 有限元分析课程设计.docx

有限元分析课程设计.docx

上传人：可****

文档编号：2686439

上传时间：2024-06-04

格式：DOCX

页数：55

大小：420.22KB

《有限元分析课程设计.docx》由会员分享，可在线阅读，更多相关《有限元分析课程设计.docx（55页珍藏版）》请在咨信网上搜索。

1、有限元分析课程设计学院：土木建筑工程学院专业：工程力学班级：力学131 学号： 3121631023 姓名：崔晨露一、单元划分与节点编号二、 C语言程序 #include#include#define NE 32 /单元数#define NJ 45 /节点数#define NZ 10 /支承数#define NPJ 7 /节点荷载数#define NJ2 90 /节点位移数#define DD 14 /半带宽int LXM=0; /类型判别码double EO=170e9; /杨氏模量double MU=0.3; /泊松比double LOU=0; /容重doub

2、le TE=0.1; /厚度doubleAJZNJ+13= 0,0,0,0,0,0,0,0,0.75,0,0,1.5,0,0,2.25,0,0,3,0,1.875,0,0,1.875,0.75,0,1.875,1.5,0,1.875,2.25,0,1.875,3,0,3.75,0,0,3.75,0.75,0,3.75,1.5,0,3.75,2.25,0,3.75,3,0,5.625,0,0,5.625,0.75,0,5.625,1.5,0,5.625,2.25,0,5.625,3,0,7.5,0,0,7.5,0.75,0,7.5,1.5,0,7.5,2.25,0,7.5,3,0,9.375,0

3、,0,9.375,0.75,0,9.375,1.5,0,9.375,2.25,0,9.375,3,0,11.25,0,0,11.25,0.75,0,11.25,1.5,0,11.25,2.25,0,11.25,3,0,13.125,0,0,13.125,0.75,0,13.125,1.5,0,13.125,2.25,0,13.125,3,0,15,0,0,15,0.75,0,15,1.5,0,15,2.25,0,15,3;/共36个节点intJMNE+15= 0,0,0,0,0,1,6,7,2,0,2,7,8,3,0,3,8,9,4,0,4,9,10,5,0,6,11,12,7,0,7,12,

4、13,8,0,8,13,14,9,0,9,14,15,10,0,11,16,17,12,0,12,17,18,13,0,13,18,19,14,0,14,19,20,15,0,16,21,22,17,0,17,22,23,18,0,18,23,24,19,0,19,24,25,20,0,21,26,27,22,0,22,27,28,23,0,23,28,29,24,0,24,29,30,25,0,26,31,32,27,0,27,32,33,28,0,28,33,34,29,0,29,34,35,30,0,31,36,37,32,0,32,37,38,33,0,33,38,39,34,0,34

5、,39,40,35,0,36,41,42,37,0,37,42,43,38,0,38,43,44,39,0,39,44,45,40;/共32个单元intNZCNZ+1=0,1,2,3,4,5,6,7,8,9,10;/1-5号节点的x,y被约束doublePJNPJ+12+1=0,0,0,0,-3.75e4,50,0,-7.5e4,60,0,-7.5e4,70,0,-7.5e4,80,0,8e4,83,0,8e4,87,0,-3.75e4,90;doubleAE,JSNJ*44,jsNJ+14,KZNJ2+1DD+1,PNJ2+1,S3+18+1,KE8+18+1,SZ3+132+1;intIE

6、,JE,ME,LE;void DUGD(int,int);/生成S矩阵，KE矩阵void main()intNJ1,k,IN,IM,jn,m,i,j,z,JO,ii,jj,h,dh,E,l,zl,dl,n;doublePE,c,SIG1,SIG2,SIG3,PYL,RYL,MAYL,MIYL,CETA;doubleWY8+1,YL3+1;if(LXM!=0)/平面应力问题与平面应变问题的判别EO=EO/(1.0-MU*MU);MU=MU/(1.0-MU);for(i=0;i=NJ2;i+)for(j=0;j=DD;j+)KZij=0.0;for(E=1;E=NE;E+)DUGD(E,3);fo

7、r(i=1;i=4;i+)for(ii=1;ii=2;ii+)h=2*(i-1)+ii;dh=2*(JMEi-1)+ii;for(j=1;j=4;j+)for(jj=1;jj0)KZdhdl=KZdhdl+KEhl;/*形成P矩阵*for(i=1;i0)for(i=1;i0)for(E=1;E=NE;E+)DUGD(E,1);PE=-LOU*(AE)*TE*9.8/4;P2*IE=P2*IE+PE;P2*JE=P2*JE+PE;P2*ME=P2*ME+PE;P2*LE=P2*LE+PE;/*边界条件*for(i=1;i=NZ;i+)z=NZCi;KZz1=1.0;for(j=2;jDD)JO=

8、DD;elseJO=z;for(j=2;j=JO;j+)KZz-j+1j=0.0;Pz=0.0;/*求解方程*NJ1=NJ2-1;for(k=1;kk+DD-1)IM=k+DD-1;elseIM=NJ2;IN=k+1;for(i=IN;i=IM;i+)l=i-k+1;c=KZkl/KZk1;jn=DD-l+1;for(j=1;j=1;i-)if(DDNJ2-i+1)JO=NJ2-i+1;elseJO=DD;for(j=2;j=JO;j+)h=j+i-1;Pi=Pi-KZij*Ph;Pi=Pi/KZi1;printf(n);printf(JDUVn);for(i=1;i=NJ;i+)printf

9、(%d%-9.9f%-9.9fn,i,P2*i-1,P2*i);/*求单元应力和主应力*for(i=0;i=NJ*4;i+)for(j=0;j=3;j+)JSij=0.0;k=0;for(E=1;E=NE;E+)DUGD(E,2);for(i=1;i=4;i+)for(j=1;j=2;j+)h=2*(i-1)+j;dh=2*(JMEi-1)+j;WYh=Pdh;for(n=1;n=4;n+) for(i=1;i=3;i+) YLi=0;for(j=1;j=8;j+)YLi=YLi+SZi8*(n-1)+j*WYj;SIG1=YL1;SIG2=YL2;SIG3=YL3;PYL=(SIG1+SIG

10、2)/2;RYL=sqrt(pow(SIG1-SIG2)/2.0,2)+pow(SIG3,2);MAYL=PYL+RYL;MIYL=PYL-RYL;if(SIG2=MIYL)CETA=0;elseCETA=90-57.29578*atan2(SIG3,(SIG2-MIYL);/将节点应力值记入节点应力矩阵k+;JSk0=JMEn;JSk1=SIG1;JSk2=SIG2;JSk3=SIG3;printf(n);printf(E=%d JD=%dn,E,JMEn);printf(sx=%-9.2fsy=%-9.2ftou=%-9.2fn,SIG1,SIG2,SIG3);printf(s1=%-9.

11、2fs3=%-9.2ftheta=%-9.2fn,MAYL,MIYL,CETA);for(i=0;i=NJ;i+)for(j=0;j=3;j+)jsij=0.0;printf(n);printf(n);printf(平均节点应力:n);for(ii=1;ii=NJ;ii+)k=0;for(jj=1;jj1)for(w=1;w=4;w+)for(i=0;i=3;i+) for(j=0;j=8;j+) Bij=0.0;x=xyw1*a;y=xyw2*b; B11=(y-b)/AE; B13=(b-y)/AE; B15=(b+y)/AE; B17=-(b+y)/AE; B22=(x-a)/AE; B

12、24=-(a+x)/AE; B26=(a+x)/AE; B28=(a-x)/AE; B31=(x-a)/AE; B32=(y-b)/AE; B33=-(a+x)/AE; B34=(b-y)/AE; B35=(a+x)/AE; B36=(b+y)/AE; B37=(a-x)/AE; B38=-(b+y)/AE; for(i=1;i=3;i+) for(j=1;j=8;j+) Sij=0;for(k=1;k2)C=EO*TE/(1-MU*MU);KE11=C*(b/(3*a)+(1-MU)*a/(6*b);KE12=C*(1+MU)/8;KE13=C*(-b)/(3*a)+(1-MU)*a/(12

13、*b);KE14=C*(3*MU-1)/8);KE15=C*(-b)/(6*a)-(1-MU)*a/(12*b);KE16=(-C)*(1+MU)/8;KE17=C*(b/(6*a)-(1-MU)*a/(6*b);KE18=C*(1-3*MU)/8);KE21=KE12;KE22=C*(a/(3*b)+(1-MU)*b/(6*a);KE23=KE18;KE24=C*(a/(6*b)-(1-MU)*b/(6*a);KE25=KE16; KE26=C*(-a)/(6*b)-(1-MU)*b/(12*a); KE27=KE14; KE28=C*(-a)/(3*b)+(1-MU)*b/(12*a);K

14、E31=KE13;KE32=KE23;KE33=KE11;KE34=KE16;KE35=KE17;KE36=KE14; KE37=KE15; KE38=KE12;KE41=KE14;KE42=KE24;KE43=KE34;KE44=KE22;KE45=KE23;KE46=KE28;KE47=KE12; KE48=KE26; KE51=KE15;KE52=KE25;KE53=KE35;KE54=KE45;KE55=KE11;KE56=KE12;KE57=KE13;KE58=KE14; KE61=KE16;KE62=KE26;KE63=KE36;KE64=KE46;KE65=KE56;KE66=K

15、E22;KE67=KE18;KE68=KE24;KE71=KE17;KE72=KE27;KE73=KE37;KE74=KE47;KE75=KE57;KE76=KE67;KE77=KE11;KE78=KE16; KE81=KE18;KE82=KE28;KE83=KE38;KE84=KE48;KE85=KE58;KE86=KE68;KE87=KE78;KE88=KE22;三、C语言运行结果JD U V1 0.000000000 0.0000000002 0.000000000 0.0000000003 0.000000000 0.0000000004 0.000000000 0.000000000

16、5 0.000000000 0.0000000006 -0.000188086 -0.0001649097 -0.000085891 -0.0001456768 0.000005546 -0.0001402369 0.000097076 -0.00014716110 0.000199777 -0.00016801411 -0.000342848 -0.00053002512 -0.000162895 -0.00051856813 0.000011559 -0.00051491914 0.000186049 -0.00052001815 0.000365862 -0.00053287316 -0

17、.000462248 -0.00107802617 -0.000219051 -0.00106895718 0.000017352 -0.00106648519 0.000253834 -0.00107022920 0.000497292 -0.00108057821 -0.000545495 -0.00175892922 -0.000258019 -0.00175325823 0.000023292 -0.00175205824 0.000304770 -0.00175555825 0.000592898 -0.00176351126 -0.000595214 -0.00252754127

18、-0.000280247 -0.00252471428 0.000029780 -0.00252479229 0.000339881 -0.00252799930 0.000655116 -0.00253418331 -0.000617969 -0.00334364332 -0.000289101 -0.00334275033 0.000036320 -0.00334344334 0.000361985 -0.00334591835 0.000690729 -0.00334982336 -0.000622391 -0.00417781837 -0.000289316 -0.0041784833

19、8 0.000043378 -0.00417982839 0.000374731 -0.00418206340 0.000708887 -0.00418509041 -0.000620797 -0.00501600742 -0.000282489 -0.00501633543 0.000047901 -0.00501722444 0.000385153 -0.00501884245 0.000715651 -0.005020841E=1 JD=1sx=-18739706.59 sy=-5621911.98 tou=-5750681.88s1=-3457881.11 s3=-20903737.4

20、6 theta=110.62E=1 JD=6sx=-17302482.56 sy=-831165.23 tou=3158628.43s1=-246222.60 s3=-17887425.19 theta=79.51E=1 JD=7sx=-7120413.64 sy=2223455.45 tou=3829332.97s1=3592279.42 s3=-8489237.61 theta=70.33E=1 JD=2sx=-8557637.66 sy=-2567291.30 tou=-5079977.34s1=334757.93 s3=-11459686.89 theta=119.74E=2 JD=2

21、sx=-8557637.66 sy=-2567291.30 tou=-5079977.34s1=334757.93 s3=-11459686.89 theta=119.74E=2 JD=7sx=-8151125.39 sy=-1212250.39 tou=2891482.49s1=-165310.18 s3=-9198065.59 theta=70.10E=2 JD=8sx=959114.42 sy=1520821.56 tou=3081188.22s1=4333929.78 s3=-1853993.81 theta=47.60E=2 JD=3sx=552602.14 sy=165780.64

22、 tou=-4890271.61s1=5253286.21 s3=-4534903.43 theta=136.13E=3 JD=3sx=552602.14 sy=165780.64 tou=-4890271.61s1=5253286.21 s3=-4534903.43 theta=136.13E=3 JD=8sx=35129.71 sy=-1559127.47 tou=3089199.22s1=2428387.59 s3=-3952385.35 theta=37.77E=3 JD=9sx=9154524.94 sy=1176691.10 tou=2847712.09s1=10066722.42

23、 s3=264493.63 theta=17.76E=3 JD=4sx=9671997.38 sy=2901599.21 tou=-5131758.74s1=12434523.13 s3=139073.46 theta=151.71E=4 JD=4sx=9671997.38 sy=2901599.21 tou=-5131758.74s1=12434523.13 s3=139073.46 theta=151.71E=4 JD=9sx=8113710.99 sy=-2292688.73 tou=3821656.61s1=9366388.14 s3=-3545365.88 theta=18.15E=

24、4 JD=10sx=18346185.68 sy=777053.67 tou=3094456.29s1=18875279.52 s3=247959.84 theta=9.70E=4 JD=5sx=19904472.06 sy=5971341.62 tou=-5858959.06s1=22040677.42 s3=3835136.26 theta=159.97E=5 JD=6sx=-13982262.17 sy=164900.89 tou=-3822943.62s1=1131869.94 s3=-14949231.21 theta=104.19E=5 JD=11sx=-14563359.43 s

25、y=-1772089.99 tou=2955975.92s1=-1122021.26 s3=-15213428.16 theta=77.60E=5 JD=12sx=-6816022.82 sy=552110.99 tou=2684797.20s1=1426605.81 s3=-7690517.63 theta=71.96E=5 JD=7sx=-6234925.55 sy=2489101.88 tou=-4094122.34s1=4109479.11 s3=-7855302.78 theta=111.59E=6 JD=7sx=-7265637.30 sy=-946603.96 tou=-5031

26、972.82s1=1835537.06 s3=-10047778.32 theta=118.94E=6 JD=12sx=-7399471.66 sy=-1392718.48 tou=2205311.70s1=-670015.27 s3=-8122174.86 theta=71.86E=6 JD=13sx=871710.65 sy=1088636.22 tou=2142855.67s1=3125772.32 s3=-1165425.45 theta=46.45E=6 JD=8sx=1005545.01 sy=1534750.73 tou=-5094428.85s1=6371443.79 s3=-

27、3831148.05 theta=133.51E=7 JD=8sx=81560.30 sy=-1545198.29 tou=-5086417.85s1=4419222.89 s3=-5882860.88 theta=139.54E=7 JD=13sx=218004.27 sy=-1090385.06 tou=2146082.96s1=1807387.71 s3=-2679768.51 theta=36.52E=7 JD=14sx=8483719.49 sy=1389329.50 tou=2209756.81s1=9115713.11 s3=757335.88 theta=15.96E=7 JD

28、=9sx=8347275.52 sy=934516.27 tou=-5022744.00s1=10883107.62 s3=-1601315.83 theta=153.21E=8 JD=9sx=7306461.57 sy=-2534863.56 tou=-4048799.48s1=8758058.99 s3=-3986460.98 theta=160.28E=8 JD=14sx=7904177.59 sy=-542476.81 tou=2673789.95s1=8679414.72 s3=-1317713.94 theta=16.17E=8 JD=15sx=15587136.94 sy=176

29、2410.99 tou=2952724.10s1=16191378.35 s3=1158169.58 theta=11.57E=8 JD=10sx=14989420.92 sy=-229975.76 tou=-3769865.34s1=15872036.37 s3=-1112591.20 theta=166.82E=9 JD=11sx=-11040187.35 sy=-715138.37 tou=-3421533.32s1=315765.06 s3=-12071090.78 theta=106.77E=9 JD=16sx=-11218646.41 sy=-1310001.90 tou=2092094.30s1=-886390.77 s3=-11642257.55 theta=78.55E=9 JD=17sx=-4917357.71 sy=580384.71 tou=2008813.40s1=1236160.94 s3=-5573133.94 theta=71.92E=9 JD=12sx=-4738898.65 sy=1175248.24 tou=-3504814.21s1=2803805.18 s3=-6367455.59 theta=114.92E=10 JD=12sx=-532234

文档加载中……请稍候！
如果长时间未打开，您也可以点击刷新试试。

下载文档保存到电脑，查找使用更方便

10 金币 0人已下载

申诉本文档由用户提供并上传，收益归属内容提供方，若内容存在侵权，请申请举报、认领或删除 立即下载

配套讲稿：: 如PPT文件的首页显示word图标，表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
特殊限制：: 部分文档作品中含有的国旗、国徽等图片，仅作为作品整体效果示例展示，禁止商用。设计者仅对作品中独创性部分享有著作权。
关键词：: 有限元分析课程设计

咨信网温馨提示：
1、咨信平台为文档C2C交易模式，即用户上传的文档直接被用户下载，收益归上传人（含作者）所有；本站仅是提供信息存储空间和展示预览，仅对用户上传内容的表现方式做保护处理，对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿，我们不确定上传用户享有完全著作权，根据《信息网络传播权保护条例》，如果侵犯了您的版权、权益或隐私，请联系我们，核实后会尽快下架及时删除，并可随时和客服了解处理情况，尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确)，网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据，个别因单元格分列造成显示页码不一将协商解决，平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺，下载前须认真查看，确认无误后再购买，务必慎重购买；若有违法违纪将进行移交司法处理，若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传，付费前请自行鉴别，如您付费，意味着您已接受本站规则且自行承担风险，本站不进行额外附加服务，虚拟产品一经售出概不退款（未进行购买下载可退充值款），文档一经付费（服务费）、不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印，是因预览和防盗链等技术需要对页面进行转换压缩成图而已，我们并不对上传的文档进行任何编辑或修改，文档下载后都不会有水印标识（原文档上传前个别存留的除外），下载后原文更清晰；试题试卷类文档，如果标题没有明确说明有答案则都视为没有答案，请知晓；PPT和DOC文档可被视为“模板”，允许上传人保留章节、目录结构的情况下删减部份的内容；PDF文档不管是原文档转换或图片扫描而得，本站不作要求视为允许，下载前自行私信或留言给上传者【可****】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权，请谨慎使用；网站提供的党政主题相关内容(国旗、国徽、党徽－－等)目的在于配合国家政策宣传，仅限个人学习分享使用，禁止用于任何广告和商用目的。
6、文档遇到问题，请及时私信或留言给本站上传会员【可****】，需本站解决可联系【微信客服】、【 QQ客服】，若有其他问题请点击或扫码反馈【服务填表】；文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“【版权申诉】”（推荐），意见反馈和侵权处理邮箱：1219186828@qq.com；也可以拔打客服电话：4008-655-100；投诉/维权电话：4009-655-100。

关于本文

本文标题：有限元分析课程设计.docx
链接地址：https://www.zixin.com.cn/doc/2686439.html

可****

内容提供者

实名认证

查看上传人更多文档

部分上传会员的收益排行 01、路***（￥15400+），
02、曲****（￥15300+），
03、wei****016（￥13200+）,
04、大***流（￥12600+），
05、Fis****915（￥4200+），
06、h****i（￥4100+），
07、Q**（￥3400+），
08、自******点（￥2400+），
09、h*****x（￥1400+），
10、c****e（￥1100+）,
11、be*****ha（￥800+），
12、13********8（￥800+）。

相似文档

自信AI助手