- 积分
- 0
- 注册时间
- 2006-6-5
- 仿真币
-
- 最后登录
- 1970-1-1
|
最近,我研究了一下《有限元方法编程》第二版(Smith与Griffiths合编)中的例程5.2。对于书中的问题,我采用了ANSYS重新进行了分析,以下是分析使用的命令流:
finish
/clear
/prep7
et,1,plane42,,,2 !平面应变问题
mp,ex,1,1.E6
mp,prxy,1,0.3
n, 1,0.0000E+00,0.0000E+00
n, 2,0.0000E+00,-0.5000E+01
n, 3,0.0000E+00,-0.1000E+02
n, 4,0.1000E+02,0.0000E+00
n, 5,0.1000E+02,-0.5000E+01
n, 6,0.1000E+02,-0.1000E+02
n, 7,20.0,0.0
n, 8,0.2000E+02,-0.5000E+01
n, 9,0.2000E+02,-0.1000E+02
n, 10,0.3000E+02,0.0000E+00
n, 11,0.3000E+02,-0.5000E+01
n, 12,0.3000E+02,-0.1000E+02
e,2,5,4,1
e,3,6,5,2
e,5,8,7,4
e,6,9,8,5
e,8,11,10,7
e,9,12,11,8
/solu
antype,static
f,1,fy,-5.0
f,4,fy,-5.0
d,1,ux,0
d,2,ux,0
d,3,all,0
d,6,all,0
d,9,all,0
d,12,all,0
d,10,ux,0
d,11,ux,0
iswrite,on
solve
finish
/post1
presol,s,x
presol,s,y
presol,s,xy
其计算的结点位移结果为:
NODE UX
1 0.0000
2 0.0000
3 0.0000
4 -0.27011E-06
5 0.15050E-05
6 0.0000
7 0.31518E-06
8 0.56193E-06
9 0.0000
10 0.0000
11 0.0000
12 0.0000
NODE UY
1 -0.90221E-05
2 -0.42485E-05
3 0.0000
4 -0.37089E-05
5 -0.18870E-05
6 0.0000
7 0.79676E-06
8 0.26713E-06
9 0.0000
10 -0.10849E-07
11 0.59749E-07
12 0.0000
单元积分点应力为:
!INITIAL STRESS RECORD FOR ELEMENT 2
! SX SY SZ SXY SYZ SXZ
eis, 2,1
-0.300421 -0.932356 -0.369833 0.103297 0.00000 0.00000
-0.300421 -0.632702 -0.279937 0.103297 0.00000 0.00000
-0.204939 -0.632702 -0.251292 0.103297 0.00000 0.00000
-0.204939 -0.932356 -0.341189 0.103297 0.00000 0.00000
!
!INITIAL STRESS RECORD FOR ELEMENT 1
! SX SY SZ SXY SYZ SXZ
eis, 1,1
-0.241079 -1.03949 -0.384171 0.793198E-01 0.00000 0.00000
-0.241079 -0.664949 -0.271808 0.793198E-01 0.00000 0.00000
-0.353698 -0.664949 -0.305594 0.793198E-01 0.00000 0.00000
-0.353698 -1.03949 -0.417956 0.793198E-01 0.00000 0.00000
!
!INITIAL STRESS RECORD FOR ELEMENT 4
! SX SY SZ SXY SYZ SXZ
eis, 4,1
-0.127013 -0.381934 -0.152684 0.120921 0.00000 0.00000
-0.127013 -0.108594 -0.706821E-01 0.120921 0.00000 0.00000
-0.186843 -0.108594 -0.886311E-01 0.120921 0.00000 0.00000
-0.186843 -0.381934 -0.170633 0.120921 0.00000 0.00000
!
!INITIAL STRESS RECORD FOR ELEMENT 3
! SX SY SZ SXY SYZ SXZ
eis, 3,1
-0.147111 -0.333461 -0.144171 0.503106E-01 0.00000 0.00000
-0.147111 -0.350822E-01 -0.546579E-01 0.503106E-01 0.00000 0.00000
-0.501469E-01 -0.350822E-01 -0.255687E-01 0.503106E-01 0.00000 0.00000
-0.501469E-01 -0.333461 -0.115082 0.503106E-01 0.00000 0.00000
!
!INITIAL STRESS RECORD FOR ELEMENT 6
! SX SY SZ SXY SYZ SXZ
eis, 6,1
-0.113808E-02 0.409506E-01 0.119437E-01 0.176247E-01 0.00000 0.00000
-0.113808E-02 0.146360E-01 0.404937E-02 0.176247E-01 0.00000 0.00000
-0.367899E-01 0.146360E-01 -0.664618E-02 0.176247E-01 0.00000 0.00000
-0.367899E-01 0.409506E-01 0.124819E-02 0.176247E-01 0.00000 0.00000
!
!INITIAL STRESS RECORD FOR ELEMENT 5
! SX SY SZ SXY SYZ SXZ
eis, 5,1
-0.403818E-01 0.745721E-01 0.102571E-01 -0.290092E-01 0.00000 0.00000
-0.403818E-01 -0.159042E-02 -0.125917E-01 -0.290092E-01 0.00000 0.00000
-0.247269E-01 -0.159042E-02 -0.789518E-02 -0.290092E-01 0.00000 0.00000
-0.247269E-01 0.745721E-01 0.149536E-01 -0.290092E-01 0.00000 0.00000
而Smith的程序计算得到的结果为:
Node Displacement
1 0.0000E+00 -0.8601E-05
2 0.0000E+00 -0.4056E-05
3 0.0000E+00 0.0000E+00
4 -0.5472E-07 -0.3771E-05
5 0.1225E-05 -0.1910E-05
6 0.0000E+00 0.0000E+00
7 0.1792E-06 0.5860E-06
8 0.6280E-06 0.1707E-06
9 0.0000E+00 0.0000E+00
10 0.0000E+00 0.1129E-06
11 0.0000E+00 0.1066E-06
12 0.0000E+00 0.0000E+00
The Gauss Point stresses for element 1 are :
Point 1
-0.4299E+00 -0.1058E+01 0.1431E+00
Point 2
-0.2511E+00 -0.6411E+00 0.8629E-01
Point 3
-0.3304E+00 -0.1016E+01 0.8352E-01
Point 4
-0.1516E+00 -0.5985E+00 0.2667E-01
The Gauss Point stresses for element 2 are :
Point 1
-0.2856E+00 -0.9141E+00 0.8499E-01
Point 2
-0.1427E+00 -0.5806E+00 0.1394E+00
Point 3
-0.3808E+00 -0.9549E+00 0.3735E-01
Point 4
-0.2379E+00 -0.6215E+00 0.9177E-01
The Gauss Point stresses for element 3 are :
Point 1
-0.1513E+00 -0.3680E+00 0.6412E-01
Point 2
0.2766E-03 -0.1429E-01 0.1010E+00
Point 3
-0.2159E+00 -0.3957E+00 0.1359E-01
Point 4
-0.6432E-01 -0.4197E-01 0.5050E-01
The Gauss Point stresses for element 4 are :
Point 1
-0.2331E+00 -0.4231E+00 0.1477E+00
Point 2
-0.9447E-01 -0.9964E-01 0.1211E+00
Point 3
-0.1867E+00 -0.4032E+00 0.1015E+00
Point 4
-0.4805E-01 -0.7974E-01 0.7493E-01
The Gauss Point stresses for element 5 are :
Point 1
0.1061E-02 0.7274E-01 -0.4210E-01
Point 2
-0.2619E-01 0.9144E-02 -0.2217E-01
Point 3
-0.3382E-01 0.5779E-01 -0.3302E-01
Point 4
-0.6108E-01 -0.5806E-02 -0.1308E-01
The Gauss Point stresses for element 6 are :
Point 1
-0.4855E-01 0.1373E-01 0.3616E-01
Point 2
-0.5281E-01 0.3775E-02 0.8267E-02
Point 3
0.2640E-03 0.3465E-01 0.3758E-01
Point 4
-0.4001E-02 0.2469E-01 0.9688E-02
显然,上述两程序结果是不一样的,请问这是为什么? |
评分
-
1
查看全部评分
-
|