应力分布不对称

hitboy

我加的是对称载荷，可是在后处理的应力图上，每层的应力都不对称，而且应力图上面都有网格，哪位大师帮忙看看。谢谢

hitboy

附件上是我的程序，改成TXT就可以运行了，盼望哪位牛人给解决一下，非常感谢。

ldhltw

结构是否完全满足堆成要求
最好直接把命令流贴出来

tonnyw

1# hitboy
PDF file cannot be opened.

To get symmetric results, you have to have: 1. symmetric load; 2. symmetric mesh; 3. symmetric boundary condition.

In plotting results, it is better to plot the invariant, such as effective stress, displacement magnitude. If you want to plot tensor component or vector, you gotta make sure what kind of coordinate system you are using.

hitboy

/PREP7
ET,1,SOLID46
KEYOPT,1,2,0
KEYOPT,1,8,1
R,1,20,1
RMORE
RMORE,1,0,0.1308,1,45,0.1308
RMORE,1,-45,0.1308,1,+45,0.1308
RMORE,1,-45,0.1308,1,+45,0.1308
RMORE,1,-45,0.1308,1,90,0.1308
RMORE,1,90,0.1308,1,90,0.1308

ET,2,200
KEYOPT,2,1,6


EX=146900      
EY=11380
EZ=12380
VXY=0.3
VYZ=0.4
VXZ=0.3
GXY=6185
GXZ=6185
GYZ=5780


enum=24000  
L=88.9
W=6.35
D=1.5875

K,1,-L,-W
K,2,L,-W
K,3,L,W
K,4,-L,W
K,5,10,W
K,6,-10,W
K,7,-10,-W
K,8,10,-W
L,1,7
L,7,8
L,8,2
L,2,3
L,3,5
L,5,6
L,6,4
L,4,1
L,6,7
L,8,5
L,6,8
L,5,7
K,9,0,0
CIRCLE,9,D,,,,,  
LPTN,11,14,16
LPTN,12,13,15
LCOMB,21,25
LCOMB,26,22
LCOMB,16,19
LCOMB,20,24
LDELE,23
LDELE,27
NUMCMP,ALL
AL,1,9,7,8
AL,3,4,5,10
AL,12,17,14,9
AL,2,15,18,12
AL,15,10,11,13
AL,16,11,6,14
LPLOT


TYPE,2  
LSEL,S,,,1,7,6
LESIZE,ALL,,,40
LSEL,S,,,8,9
LESIZE,ALL,,,10

LSEL,S,,,4,10,6
LESIZE,ALL,,,10
LSEL,S,,,3,5,2
LESIZE,ALL,,,40

LSEL,S,,,9,17,8
LESIZE,ALL,,,10
LSEL,S,,,12,14,2
LESIZE,ALL,,,10

LSEL,S,,,2,18,16
LESIZE,ALL,,,10
LSEL,S,,,12,15,3
LESIZE,ALL,,,10
ALLSEL

LSEL,S,,,10,13,3
LESIZE,ALL,,,10
LSEL,S,,,11,15,4
LESIZE,ALL,,,10

LSEL,S,,,6,16,10
LESIZE,ALL,,,10
LSEL,S,,,11,14,3
LESIZE,ALL,,,10
MSHKEY,1     
AMESH,ALL


TYPE,1
EXTOPT,ESIZE,20,0   
VEXT,1,6,1,0,0,2.616  
aclear,1,6
NUMMRG,NODE
NUMMRG,ELEM
numcmp,all  

*do,i,1,enum
MP,EX,i,EX1(i)
MP,EY,i,EY1(i)
MP,EZ,i,EZ1(i)
MP,GXY,i,GXY1(i)
MP,GXZ,i,GXZ1(i)
MP,GYZ,i,GYZ1(i)
MP,prXY,i,PRXY1(i)
MP,prXZ,i,PRXZ1(i)
MP,prYZ,i,PRYZ1(i)
*enddo
*do,i,1,enum
mpchg,i,i
*enddo

finish



/solu
antype,static
nsel,s,loc,x,0
d,all,UX,0
nsel,s,loc,Y,0
D,ALL,UY,0
nsel,s,loc,x,-88.9
sf,all,pres,-10
nsel,s,loc,x,88.9
sf,all,pres,-10
allsel
solve

hitboy

从新发了一下源程序，希望哪位给看看，谢谢

hitboy

自己顶一下，希望哪位牛人帮忙解决一下，非常感谢

tonnyw

7# hitboy

I cannot open your file. I am using Ansys 8.1.

hitboy

o ,呵呵，谢谢楼上的关注。希望哪位大师能给运行看看，谢谢

hitboy

怎么没消息呢？哪位楼主有时间也给看一下吧，非常感谢。实在是不知道问题出在哪里呢

hitboy

还是没人吗？哪位帮帮忙啊

tonnyw

Are you sure this part is correct?

Do you want to assign the value of EX to EX1(i)? Or the other way?

*do,i,1,enum
MP,EX,i,EX1(i)
MP,EY,i,EY1(i)
MP,EZ,i,EZ1(i)
MP,GXY,i,GXY1(i)
MP,GXZ,i,GXZ1(i)
MP,GYZ,i,GYZ1(i)
MP,prXY,i,PRXY1(i)
MP,prXZ,i,PRXZ1(i)
MP,prYZ,i,PRYZ1(i)
*enddo

XJTU.FX

应该是约束加的不对称吧？

hitboy

这部分我觉得是没有问题的，我后面单元破坏后，要修改这些特性参数，所以定义的是EX(i)等

hitboy

非常感谢tonnyw

*do,i,1,enum
MP,EX,i,EX1(i)
MP,EY,i,EY1(i)
MP,EZ,i,EZ1(i)
MP,GXY,i,GXY1(i)
MP,GXZ,i,GXZ1(i)
MP,GYZ,i,GYZ1(i)
MP,prXY,i,PRXY1(i)
MP,prXZ,i,PRXZ1(i)
MP,prYZ,i,PRYZ1(i)
*enddo
这部分我觉得是没有问题的，我后面单元破坏后，要修改这些特性参数，所以定义的是EX(i)等

hitboy

哦，谢谢tonnyw
我发现我发的程序是错误的，少些了一些东西。
/PREP7
ET,1,SOLID46
KEYOPT,1,2,0
KEYOPT,1,8,1
R,1,20,1
RMORE
RMORE,1,0,0.1308,1,45,0.1308
RMORE,1,-45,0.1308,1,+45,0.1308
RMORE,1,-45,0.1308,1,+45,0.1308
RMORE,1,-45,0.1308,1,90,0.1308
RMORE,1,90,0.1308,1,90,0.1308

ET,2,200
KEYOPT,2,1,6


EX=146900      
EY=11380
EZ=12380
VXY=0.3
VYZ=0.4
VXZ=0.3
GXY=6185
GXZ=6185
GYZ=5780

*dim,m_1,array,enum

*do,i,1,enum
m_1(i)=1
*enddo

*dim,EX1,array,enum            !!!定义单元突降刚度系数数组
*dim,EY1,array,enum                                                                                
*dim,EZ1,array,enum
*dim,GXY1,array,enum
*dim,GYZ1,array,enum
*dim,GXZ1,array,enum
*dim,PRXY1,array,enum
*dim,PRYZ1,array,enum
*dim,PRXZ1,array,enum       !!!定义单元突降刚度系数数组
*vfact,EX
*vfun,EX1(1),copy,m_1(1)     !!ex1中存储了单元的无损伤的初始刚度系数，ex1存放突降刚度系数
*vfact,EY
*vfun,EY1(1),copy,m_1(1)
*vfact,EZ
*vfun,EZ1(1),copy,m_1(1)
*vfact,GXY                                       
*vfun,GXY1(1),copy,m_1(1)
*vfact,GYZ
*vfun,GYZ1(1),copy,m_1(1)
*vfact,GXZ
*vfun,GXZ1(1),copy,m_1(1)
*vfact,VXY
*vfun,PRXY1(1),copy,m_1(1)  
*vfact,VYZ
*vfun,PRYZ1(1),copy,m_1(1)  
*vfact,VXZ
*vfun,PRXZ1(1),copy,m_1(1)

enum=24000  
L=88.9
W=6.35
D=1.5875

K,1,-L,-W
K,2,L,-W
K,3,L,W
K,4,-L,W
K,5,10,W
K,6,-10,W
K,7,-10,-W
K,8,10,-W
L,1,7
L,7,8
L,8,2
L,2,3
L,3,5
L,5,6
L,6,4
L,4,1
L,6,7
L,8,5
L,6,8
L,5,7
K,9,0,0
CIRCLE,9,D,,,,,  
LPTN,11,14,16
LPTN,12,13,15
LCOMB,21,25
LCOMB,26,22
LCOMB,16,19
LCOMB,20,24
LDELE,23
LDELE,27
NUMCMP,ALL
AL,1,9,7,8
AL,3,4,5,10
AL,12,17,14,9
AL,2,15,18,12
AL,15,10,11,13
AL,16,11,6,14
LPLOT


TYPE,2  
LSEL,S,,,1,7,6
LESIZE,ALL,,,40
LSEL,S,,,8,9
LESIZE,ALL,,,10

LSEL,S,,,4,10,6
LESIZE,ALL,,,10
LSEL,S,,,3,5,2
LESIZE,ALL,,,40

LSEL,S,,,9,17,8
LESIZE,ALL,,,10
LSEL,S,,,12,14,2
LESIZE,ALL,,,10

LSEL,S,,,2,18,16
LESIZE,ALL,,,10
LSEL,S,,,12,15,3
LESIZE,ALL,,,10
ALLSEL

LSEL,S,,,10,13,3
LESIZE,ALL,,,10
LSEL,S,,,11,15,4
LESIZE,ALL,,,10

LSEL,S,,,6,16,10
LESIZE,ALL,,,10
LSEL,S,,,11,14,3
LESIZE,ALL,,,10
MSHKEY,1     
AMESH,ALL


TYPE,1
EXTOPT,ESIZE,20,0   
VEXT,1,6,1,0,0,2.616  
aclear,1,6
NUMMRG,NODE
NUMMRG,ELEM
numcmp,all  
*do,i,1,enum
MP,EX,i,EX1(i)
MP,EY,i,EY1(i)
MP,EZ,i,EZ1(i)
MP,GXY,i,GXY1(i)
MP,GXZ,i,GXZ1(i)
MP,GYZ,i,GYZ1(i)
MP,prXY,i,PRXY1(i)
MP,prXZ,i,PRXZ1(i)
MP,prYZ,i,PRYZ1(i)
*enddo
*do,i,1,enum
mpchg,i,i
*enddo

/solu
antype,static
nsel,s,loc,x,0
d,all,UX,0
nsel,s,loc,Y,0
D,ALL,UY,0
nsel,s,loc,x,-88.9
sf,all,pres,-10
nsel,s,loc,x,88.9
sf,all,pres,-10
allsel
solve



我从新发了一下程序，写的很罗嗦啊，才开始研究，希望哪位大师帮忙看看，弄了好几天了，实在是不找不到问题。很着急，谢谢

yzzyinyin

你新发的程序里，载荷施加不对称吗？

hitboy

从新编辑了一下我的程序，可是图出来很是不对，结果不是应该0度铺层在Y方向应力集中，90度铺层X方向应力集中比较严重吗？而且应该是比较对称的吧，哪位大师帮忙看看。非常感谢！
/PREP7
ET,1,SOLID46
KEYOPT,1,2,0
KEYOPT,1,8,1
R,1,20,1   
RMORE
RMORE,1,0,0.1308,1,45,0.1308
RMORE,1,-45,0.1308,1,+45,0.1308
RMORE,1,-45,0.1308,1,+45,0.1308
RMORE,1,-45,0.1308,1,90,0.1308
RMORE,1,90,0.1308,1,90,0.1308
MP,EX,1,146900
MP,EY,1,11380
MP,EZ,1,12380
MP,GXY,1,6185
MP,GXZ,1,6185
MP,GYZ,1,5780
MP,prXY,1,0.3
MP,prXZ,1,0.3
MP,prYZ,1,0.4
ET,2,200
KEYOPT,2,1,6

enum=24000   
L=88.9
W=6.35
D=1.5875
K,1,-L,-W
K,2,L,-W
K,3,L,W
K,4,-L,W
K,5,10,W
K,6,-10,W
K,7,-10,-W
K,8,10,-W
L,1,7
L,7,8
L,8,2
L,2,3
L,3,5
L,5,6
L,6,4
L,4,1
L,6,7
L,8,5
L,6,8
L,5,7
K,9,0,0
CIRCLE,9,D,,,,,  
LPTN,11,14,16
LPTN,12,13,15
LCOMB,21,25
LCOMB,26,22
LCOMB,16,19
LCOMB,20,24
LDELE,23
LDELE,27
NUMCMP,ALL
AL,1,9,7,8
AL,3,4,5,10
AL,12,17,14,9
AL,2,15,18,12
AL,15,10,11,13
AL,16,11,6,14
LPLOT

TYPE,2  
LSEL,S,,,1,7,6
LESIZE,ALL,,,40
LSEL,S,,,8,9
LESIZE,ALL,,,10
LSEL,S,,,4,10,6
LESIZE,ALL,,,10
LSEL,S,,,3,5,2
LESIZE,ALL,,,40
LSEL,S,,,9,17,8
LESIZE,ALL,,,10
LSEL,S,,,12,14,2
LESIZE,ALL,,,10
LSEL,S,,,2,18,16
LESIZE,ALL,,,10
LSEL,S,,,12,15,3
LESIZE,ALL,,,10
ALLSEL
LSEL,S,,,10,13,3
LESIZE,ALL,,,10
LSEL,S,,,11,15,4
LESIZE,ALL,,,10
LSEL,S,,,6,16,10
LESIZE,ALL,,,10
LSEL,S,,,11,14,3
LESIZE,ALL,,,10
MSHKEY,1     
AMESH,ALL

TYPE,1
EXTOPT,ESIZE,20,0   
VEXT,1,6,1,0,0,2.616  
NUMMRG,NODE
NUMMRG,ELEM
numcmp,all  
eplot
fini
/solu
nsel,s,loc,x,-88.9
d,all,,,,,,,all,,,,,
fini
/solu
nsel,s,loc,x,88.9
sfa,all,1,pres,-100
allsel
solve
fini
/post1
avprin,0,0
shell,mid
avres,2
/efacet,1
layer,9
force,total
plnsol,s,x

layer,1
plnsol,s,x
finish

FoRd

不错啊，我试试看

hitboy

不知道楼上试的结果怎么样，我的90度铺层最明显了，应力应该集中在X轴左右，而且应该最先发生破坏的是圆孔X轴方向的单元，可是图上应力最大的在两侧的边上，很是不明白。