ANSYS静力分析简支梁跨中挠度？

daphne_wj

好久不用ANSYS，本来就不熟，现在更生疏了，就是一个很简单的例子，我用LS-DYNA隐式算的结果与手算的估计值相差太大，想通过ANSYS隐式计算比较下。可是不会用，找资料要花不少时间，而我手边又没有ANSYS的书，之好求助论坛了。先谢了。
附有命令流文件
模型是简支梁，跨中施加集中力，想观察梁跨中的静态挠度。手算公式是PL^3/48EI.
截面矩形10mm*120mm
杨氏模量：2.754e5
密度：8.333e-8

/BATCH                                           !calcutate deflection of simply surported beam at the cencoid
/FILNAME,SimpleBeamMidDisplacement
/TITLE,SimpleBeamMidDisp
JPGPRF,500,100,1                                 ! MACRO TO SET PREFS FOR JPEG PLOTS
/SHOW,JPEG

! UNIT: mm, newton, tonne, second, g: 9810
/PREP7                                            !Enter preprocessor
FC=-200                                           !concentrated force at the mid span
ET,1,BEAM4                                        !choose element BEAM4 for type 1
R,1                                               !define the parameter for the Beam type 1, BEAM4
RMODIF,1,1,1200,1.0e4,1.44e6,10,120               !area:10*120;IYY,IZZ:1.0e4,1.44e6;thickness of Y,Zaxis:10,120
RMODIF,1,6,0                                      !Orientation about X axis:0 degree
MP,EX,1,2.754e5                                   !DEFINE MATERIAL PROPERTIES, Yong's Modulus
MP,DENS,1,8.333e-8                                !Define material property, density
MP,NUXY,1,0.3                                     !Define material property, possions' ratio
N,1,0,0,0
N,27,6000,0,0
FILL
TYPE,1                                             !BEAM4, element type 1
REAL,1                                             !real constant of element type 1
MAT,1                                              !Sets the element material attribute pointer,default 1
E,1,2
EGEN,26,1,1
D,1,UX,0                                           !Constrain deflections at the two ends
D,1,UY,0
D,1,UZ,0
D,27,UX,0  
D,27,UY,0
D,27,UZ,0  
SAVE, ModelData
FINISH

/SOLU
ANTYPE,0                            !perform a STATIC analysis
LUMPM,0                             !use the element-dependent mass matrix formulation maybe off
OUTRES,ALL,ALL                      !write solution of all the solution results for all substeps
AUTOTS,OFF                          !not to use automatic time stepping or load stepping.

F,14,FY,FC                          ! specify X-direction initial force at node 14 solve
FINISH

/POST1                             !Enter the general post-processor.
NSOL,2,14,U,Y,CENP_UY              !Define variable number 2 as the node 4 deflection in the Y-direction.
PLVAR,2                            !display variable 2 in the form of graph
EXTREM,2                           !list the extream value for variable 2

!!!!!!!!write all results to the file"xxx.dat"
*DIM,CENP,ARRAY,DIMT0,1,1                           !ÎòμÄàí½a￡oóÃDIM¶¨òåêy×éÃû×ÖèÎòa￡¬2»ò»¶¨òaoíÇ°Ãæìáè¡μıäá¿í¬Ãû
VGET,CENP(1),2                                     !óÃVGET°Ñ½á1û±äá¿2ìáè¡μ½Ö®Ç°óÃDIM¶¨òåμıäá¿êy×é,±ØDëêÇArray parameter vector
                                                   !¼′′æ′￠μ½Ç°Ã涨òåμÄCENPμÄμú1λÏòá¿ÖD

*GET,RESTIME,ACTIVE,,TIME,CPU                      !GETìá衽á1ûμ½±êá¿￡¬Ãû3ÆÎaRESTIME
*GET,RESMIN,VARI,2,EXTREM,VMIN, ,             !retrieve No.2 variable's EXTREM's VMIN as name RESMIN
*GET,RESMAX,VARI,2,EXTREM,VMAX, ,             !retrieve No.2 variable's EXTREM's VMAX as name RESMAX

*DIM,RES,,1,3
*VFILL,RES(1,1),DATA,RESTIME
*VFILL,RES(1,2),DATA,RESMIN
*VFILL,RES(1,3),DATA,RESMAX

*CFOPEN,midspandisp,dat   
*VWRITE,RES(1,1),RES(1,2),RES(1,3)
(F10.3,'       ',F20.15,'     ',F20.15)                  
*VWRITE,CENP(1)                    
(F14.6)                              
*CFCLOS
FINISH

daphne_wj

我知道后处理有问题，我的后处理是针对post26的，可是删除后处理部分运行仍然有错。请帮忙看一下，很简单的。

woomii

可以算啊，有错指的是还是和手算结果不付吗？手算结果是多少？

daphne_wj

是啊，差两个数量级

daphne_wj

哪位高手，帮忙看看