关于压电片二维弹性系数矩阵问题

怡人

各位，在输入压电材料的“stiffness” Matrix【C】(6*6的对称矩阵)时，看ansys帮助中的例子，他所解决的是三维问题，在输入弹性矩阵时的命令流如下：
tbdata,1,12.80e10,6.8e10,6.6e10
tbdata,7,12.80e10,6.6e10
tbdata,12,11.0e10
tbdata,16,2.1e10
tbdata,19,2.1e10
tbdata,21,2.1e10
在这里请问各位，命令中的1，7，12，16，19，21这些位置是怎么定出来的？满足什么规律？如果要解决二维问题，即：我要输入一个4*4的对称矩阵，我的位置应该怎么定？我直接把矩阵的每一行写出来，不写成这种简写的形势，而是直接写成下面这种命令是否可以，我的命令如下：
tbdata,1,12.80e10,6.80e10,0,0
tbdata,5,6.80e10,12.8e10,0,0
tbdata,9,0,0,2.1e10,0
tbdata,13,0,0,0,2.1e10
这样是否可以？因为在计算结果的时候在input4*4矩阵这块会提示出错，望各位老师指点~~还有一个小问题，就是在流固耦合面或者固固耦合面上是不是不能加载位移约束？因为我的耦合点处也会出现提示出错的对话框，如果把这些点处的位移约束去掉的话就不会出现错误了，还望各位指导，不胜感激~~

criyuel

压电矩阵在材料里头都会给你，命令流里头的位置参数就对应的压电矩阵里头的各个参数，一一对应进去就好了，一目了然，比如一般都是类似于下边这种形式的6*6矩阵，1就是第一行第一列的元素，然后按行往后查就好了，7对应的就是第二行的第一个以此类推，大概就是这个意思，后三行的前几个因子都是0，就不计算在内，多以加的时候就把他们去掉就行了。
1 2 3 0 0 0
0 1 2 0 0 0
0 0 1 0 0 0
0 0 0 1 0 0
0 0 0 0 1 0
0 0 0 0 0 1

沙场勇士

回复 2 # criyuel 的帖子

关键是怎么定义4X4的矩阵，ansys中还是6X6的矩阵，只是矩阵中的参数变少了，所以运行的时候会出现错误说这个矩阵不正定