当前位置：首页 > 有限元程序的设计大作业
disx(1)=disx1(1); disy(1)=disy1(1); spt1(1)=spt11(1); spt2(1)=spt21(1); spt3(1)=spt31(1);
xpt(2)=xpt1(2); ypt(2)=ypt1(2); disx(2)=disx1(2); disy(2)=disy1(2); spt1(2)=spt11(2); spt2(2)=spt21(2); spt3(2)=spt31(2);
xpt(3)=xpt1(3); ypt(3)=ypt1(3); disx(3)=disx1(3); disy(3)=disy1(3); spt1(3)=spt11(3); spt2(3)=spt21(3); spt3(3)=spt31(3);
xpt(4)=xpt1(4); ypt(4)=ypt1(4); disx(4)=disx1(4); disy(4)=disy1(4); spt1(4)=spt11(4); spt2(4)=spt21(4); spt3(4)=spt31(4); xpt(5)=xpt1(1); ypt(5)=ypt1(1); disx(5)=disx1(1); disy(5)=disy1(1); spt1(5)=spt11(1); spt2(5)=spt21(1); spt3(5)=spt31(1);
scale=1e1; % 画变形图 subplot(2,3,1) hold on
plot(xpt,ypt,'b-')
plot(xpt+scale*disx,ypt+scale*disy,'r--') title('Deformed shape')
axis equal % axis off % 画应力图 subplot(2,3,2) hold on
fill(xpt+scale*disx,ypt+scale*disy,spt1) colorbar axis equal axis ([0 b 0 b])
title('Stress-xx distribution') % axis off subplot(2,3,3) hold on
fill(xpt+scale*disx,ypt+scale*disy,spt2) colorbar axis equal axis ([0 b 0 b])
title('Stress-yy distribution') % axis off subplot(2,3,4) hold on
fill(xpt+scale*disx,ypt+scale*disy,spt3) colorbar axis equal axis ([0 b 0 b])
title('Stress-xy distribution') % axis off end
%画x=0轴上的x向应力变化曲线 subplot(2,3,5) hold on
A=zeros(1,m+1); B=zeros(1,m+1); C=zeros(1,m+1); D=zeros(1,m+1); for i=1:1:m+1
A(i) = gNode(i*(n+1),2); B(i) = Stress(i*(n+1),1);
C(i) = p*(1+r^2/2/gNode(i*(n+1),2)^2+3*r^4/2/gNode(i*(n+1),2)^4); D(i) = p; end
plot(B,A,'r:*'); plot(C,A,'k:x');
plot(D,A,'b--');
title('x=0轴上的x向应力变化曲线') ylabel('点距小孔圆心的距离(mm)') xlabel('应力大小/ Pa') axis fill
%画小孔边缘的剪切力力变化曲线（逆时针） subplot(2,3,6) hold on
A=zeros(1,n+1); B=zeros(1,n+1); C=zeros(1,n+1); D=zeros(1,n+1); for i=1:1:n+1
A(1,i)=(i-1)*(pi/2/n); B(1,i)=Stress(i,3);
C(1,i)=p*(sin(4*A(1,i))-sin(2*A(1,i)))/2; D(1,i)=(xy_max+xy_min)/2; end
plot(A,B,'r:*'); plot(A,C,'k:x'); plot(A,D,'b--');
title('孔边缘剪切力变化曲线（逆时针)') xlabel('角度（单位为弧度）') ylabel('应力大小/ Pa') axis fill return
附录二 10x20网格计算结果
=================================================================================
=
=
=
=
结点号 X坐标 Y坐标 x-位移 y-位移 Stress-xx Stress-yy Stress-xy
---------------------------------------------------------------------------------------------------------------------------------------
1 1.0000 0.0000 3.0494e-02 -6.7903e-19 -6.7755e+07 -1.0550e+09 -4
.
2
6
9
2
e
+
0
7
2 0.9969 0.0785 3.0399e-02 -8.1180e-04 -8.3977e+07 -1.0451e+09 1
.
0
1
4
0
e
+
0
8
3 0.9877 0.1564 3.0115e-02 -1.6189e-03 -1.0225e+08 -9.7927e+08 1
.
4
9
9
6
e
+
0
8
4 0.9724 0.2334 2.9645e-02 -2.4168e-03 -1.1795e+08 -8.6325e+08 1
.
7
8
1
4
e
+
0
8
5 0.9511 0.3090 2.8990e-02 -3.2006e-03 -1.2395e+08 -7.0667e+08 1
.
7
7
9
5
e
+
0
8
6 0.9239 0.3827 2.8156e-02 -3.9659e-03 -1.1140e+08 -5.2208e+08 1
.
4
7 0.8910 0.45407
.
6
8 0.8526 0.5225-2
.
2
9 0.8090 0.5878 -1
.
4
10 0.7604 0.6494-2
.
8
11 0.7071 0.7071-4
.
2
12 0.6494 0.7604-5
.
5
13 0.5878 0.8090-6
.
6
14 0.5225 0.8526-7
.
4
15 0.4540 0.8910-7
.
7
16 0.3827 0.9239-7
.
4
17 0.3090 0.9511-6
.
7
18 0.2334 0.9724-5
.
5
19 0.1564 0.9877-3
.
9
20 0.0785 0.9969-2
.
0
21 0.0000 1.00004
.
5
22 1.0662 0.0000-5
.
0
23 1.0629 0.08374
.
0
4
5
4
2.7147e-02 -4.7082e-037
9
4
2.5971e-02 -5.4228e-034
3
6
2.4636e-02 -6.1055e-03 6
5
4
2.3148e-02 -6.7519e-03 5
6
9
2.1519e-02 -7.3580e-03 7
8
0
1.9758e-02 -7.9197e-03 9
6
0
1.7876e-02 -8.4333e-03 7
9
8
1.5886e-02 -8.8955e-03 1
2
3
1.3800e-02 -9.3028e-03 0
2
6
1.1629e-02 -9.6526e-03 9
5
1
9.3887e-03 -9.9424e-03 7
5
4
7.0915e-03 -1.0170e-02 7
3
1
4.7515e-03 -1.0334e-02 5
9
5
2.3829e-03 -1.0432e-02 4
1
9
2.0552e-18 -1.0465e-02 7
7
9
3.0658e-02 -6.3820e-191
6
3
3.0545e-02 -6.1710e-048
0
2
e
+
-7.0862e+07 e
+
6.5381e+06 e
+
1.2792e+08 e
+
2.9760e+08 e
+
5.1625e+08 e
+
7.8041e+08 e
+
1.0824e+09 e
+
1.4106e+09 e
+
1.7502e+09 e
+
2.0842e+09 e
+
2.3945e+09 e
+
2.6637e+09 e
+
2.8761e+09 e
+
3.0192e+09 e
+
3.0845e+09 e
+
-8.8630e+07 e
+
-9.2802e+07 e
+
0
8
-3.2377e+08 0
7
-1.2650e+08 0
7
5.5826e+07 0
8
2.1143e+08 0
8
3.3162e+08 0
8
4.1152e+08 0
8
4.5032e+08 0
8
4.5125e+08 0
8
4.2104e+08 0
8
3.6911e+08 0
8
3.0650e+08 0
8
2.4463e+08 0
8
1.9404e+08 0
8
1.6329e+08 0
8
1.5798e+08 0
6
-7.6429e+08 0
7
-7.5624e+08 0
7
本文作者：...共分享92篇相关文档
