有限元分析课程设计.docx

有限元分析 课程设计 学 院：土木建筑工程学院 专 业： 工程力学 班 级： 力学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; /容重double 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,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,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,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,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);for(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=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(%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+SIG2)/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.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; B24=-(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*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);KE31=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=KE22;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.0000000005 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.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 -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.00417848338 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.46 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=2sx=-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 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 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=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 sy=-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=-5031972.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=-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=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=1762410.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