转自www.simwe.com badapple100 wrote: 怎么在ansys里得到任意点(不在节点上)上的应力结果。 上有个mac是专门干这个用的 http:///ansys/macros/sxyz.mac 还有,我用solid92单元,怎么结果的element list里只有4个节点上结果,边节点上的节点信息没有呢? 要显示中结点应力必须/efacet,2 要用/efacet,2必须/GRAPHICS,POWER 而用了/GRAPHICS,POWER之后只对露在表面的结点进行平均。 文件里有写: 把它拷到ansys工作目录下,然后运行ansys,在命令输入框输入 sxyz,X,Y,Z,avg where: X= x coordinate Y= y coordinate Z= z coordinate avg=0 for averaged result, 1 for unaveraged result 结果好像有点问题。我输入一个已知节点的坐标,它的计算结果和list里列出的结果有相当大的出入。如下: sxyz结果: Averaged Stress Results SX = -1.397905908E-05. SY = -8.651374752E-06. SZ = -5.144134442E-05. SXY = -3.945676892E-08. SYZ = -1.143312579E-06. SXZ = 2.10849745E-05. S1 = -4.44601699E-06. S2 = -8.692004552E-06. S3 = -6.09337567E-05. SINT = 5.648777906E-05. SEQV = 5.448900002E-05. liat结果: NODE SX SY SZ SXY SYZ SXZ 4 -0.17083E-04-0.33587E-05-0.91514E-04-0.21577E-05 0.52513E-06 0.38929E-04 我试了一下,没问题啊。 输入文件: fini /clear /prep7 et,1,82 mp,ex,1,2e11 mp,nuxy,1,0.3 rect,,1,,1 esize,0.1 amesh,1 /solu nsel,s,loc,x,0 d,all,all,0 nsel,s,loc,x,1 nsel,r,loc,y,0 f,all,fy,-1000 allsel,all solve /post1 /output,sxyz1,txt sxyz,0.5,0.5,0,0 /output 使用 sxyz,0.5,0.5,0,0 得到的结果: USE COMMAND MACRO SXYZ ARGS= 0.50000 0.50000 Averaged Stress Results at X= 0.5, Y= 0.5, Z= 0. SX = 511.178887. SY = -94.6331167. SZ = 0. SXY = -1665.15374. SYZ = 0. SXZ = 0. S1 = 1900.7828. S2 = 0. S3 = -1484.23703. SINT = 3385.01983. SEQV = 2938.90236. PRINTOUT RESUMED BY /GOP PRNSOL(list结点应力)得到的结果: NODE SX SY SZ SXY SYZ SXZ 211 511.18 -94.633 0.0000 -1665.2 0.0000 0.0000 NODE S1 S2 S3 SINT SEQV 211 1900.8 0.0000 -1484.2 3385.0 2938.9 结果是一样的。 我试了一下,这个mac对平面问题(四边形和三角单元)没问题。对三维问题只有对六面体单元是对的,对四面体单元计算有出入。我用的是solid92。 怎么办?请大家献计献策!! 与/GRAPHICS模式有关。当/GRAPHICS,power时为提高显示速度,只对露在表面的节点应力进行平均。 输入文件: fini /clear /prep7 et,1,92 mp,ex,1,2e11 mp,nuxy,1,0.3 block,,1,,1,,1 esize,0.1 MSHAPE,1,3D MSHKEY,0 vmesh,1 /solu nsel,s,loc,x,0 d,all,all,0 nsel,s,loc,x,1 nsel,r,loc,y,0 f,all,fy,-1000 allsel,all solve /post1 /output,sxyz2,txt sxyz,1,1,1,0 /output 结果: sxyz是通过定义path得出的,不受/GRAPHICS模式的影响: USE COMMAND MACRO SXYZ ARGS= 1.0000 1.0000 1.0000 Averaged Stress Results at X= 1, Y= 1, Z= 1. SX = -392.640175. SY = -435.647891. SZ = 336.343785. SXY = 554.031306. SYZ = -168.488453. SXZ = -153.169632. S1 = 485.79912. S2 = -9.20216826. S3 = -968.56425. SINT = 1454.87185. SEQV = 1281.20515. 当/GRAPHICS,full时list节点应力的结果(与sxyz的结果基本相同): 507 -392.95 -435.93 336.34 554.03 -168.50 -153.18 507 485.80 -9.2631 -969.07 1454.9 1281.2 当/GRAPHICS,power时list节点应力的结果(与sxyz的结果相差较大): 507 -452.84 -481.39 368.61 604.46 -197.24 -167.05 507 541.36 -27.676 -1079.3 1620.7 1424.8 再次感谢pjwseu !的确像你说的。 如果换一种说法的话,这是不是意味着在显示coutour plot或list结果的时候应该用/GRAPHICS,full得到的才是正确的结果,否则只能是变相的 |
|