|
楼主 |
发表于 2007-1-7 12:38
|
显示全部楼层
1.如果处理成orthotropic,很明显我的矩阵是三维下的orthotropic退化到平面的orthotropic是
c11 c12 c16
c12 c22 c26
c16 c26 c66
(c16= c26=0)
2.这样的话,与我的矩阵是不符合的.(我的c16= c26不等于零,尽管趋于零,但我在另外一个算例中,这两个数为负数,并不趋于零)
3.如何处理?我这样做对不对?
i.写出三维下弹性本构方程(6个,并且注意到是正交异性,三个方向工程弹性常数不同)
ii.然后,既然我这里定的是平面问题,就有平面问题属性的定义,我定义为平面应变状态.所以令其中一个方向(3方向)的正应变为零,得出相关的关系式,重新带入本构方程,并且整理得到
iii.得到平面应变下的正交异向弹性柔度矩阵Sij:
(1-u31*u13)/E1 -(u21+u31*u23)/E2 0
-(u12+u13*u32)/E1 (1-u32*u23)/E2 0
0 0 1/G12
iv.当然里面有一定关系:
u21/E2=u12/E1;
u31/E3=u13/E1;
u32/E3=u23/E2;
v.然后建立等价关系
以上处理过程对吗? |
|