22#材料模型的坐标定义

sdwangq

使用22#材料模型时，在AOPT=2.0时，如何使用一个矢量定义材料方向啊？
*MAT_COMPOSITE_DAMAGE
         1 0.170E+04 0.527E+11 0.158E+11 0.158E+11  0.280000    0.0       0.0
0.546E+10 0.546E+10 0.546E+10 0.02     2.0   
0.0       0.0       0.0       0.0       0.0       1.0
0.0       0.0       0.0       0.0       0.0       0.0       0.0
1.730E+07 0.812E+09 0.214E+08 0.119E+09
我使用上述参数进行计算时，提示说矢量定义错误？该如何定义该矢量？

gdyu_yu

我看了你传给我的dyna2.k，有两个问题可以探讨下：
1、自动单面接触没有定义ssid和sstyp。（这个是否需要设置，有没有设置究竟影响大不大，你可以自己验算一下）
2、沙漏控制选项IHQ=3不合适于壳单元。
试算后并没有你所说的矢量出错的提示。（我用lsdyna 971R5.1.1单精度求解器算的，可以正常计算下去）
针对以上两点，我做了相应的改动，并改变了*CONTROL_SHELL中的一个选项 lamsht的默认值，设为=1。
*CONTACT_AUTOMATIC_SINGLE_SURFACE
$#     cid                                                                 title
$#    ssid      msid     sstyp     mstyp    sboxid    mboxid       spr       mpr
            1           0           3             0            0             0          1           1
$#      fs          fd          dc           vc        vdc    penchk        bt           dt
      0.10      0.10        0.0          0.0         0.0             0       0.0   1.0E+7
$#     sfs       sfm         sst         mst         sfst        sfmt        fsf          vsf
        0.0        0.0        0.0          0.0         0.0          0.0       0.0         0.0

*HOURGLASS
$#    hgid       ihq        qm       ibq        q1        q2    qb/vdc        qw
             1         4      0.10          0     1.50     0.06          0.0        0.0

*CONTROL_SHELL
$#  wrpang     esort     irnxx    istupd    theory       bwc     miter      proj
          20.0            1          -1          1             2           2           1          0
$#  rotascl      intgrd    lamsht    cstyp6    tshell    nfail1    nfail4   psnfail
           1.0            0             1            1          0         0         0            0
$# psstupd    irquad     cntco    itsflg    irquad
              0             0            0         0         2

LAMSHT:=For composite shells with material types:
*MAT_COMPOSITE_DAMAGE
*MAT_ENHANCED_COMPOSITE_DAMAGE.
If this flag is set laminated shell theory is used. Lamination theory is applied to correct for the assumption of a uniform constant shear strain through the thickness of the shell. Unless this correction is applied, the stiffness of the shell can be grossly incorrect if there are drastic differences in the elastic constants from ply to ply, especially for sandwich type shells. Generally, without this correction the results are too stiff. For the discrete Kirchhoff shell elements, which do not consider transverse shear, this option is ignored.
EQ.0: do not update shear corrections,
EQ.1: activate laminated shell theory.

gdyu_yu

可以看关键字用户手册2#材料MAT_ORTHOTROPIC_ELASTIC的说明。

sdwangq

gdyu_yu 发表于 2012-5-3 15:02
可以看关键字用户手册2#材料MAT_ORTHOTROPIC_ELASTIC的说明。

看了，不过两者的参数设置不太一样啊，望指点

gdyu_yu

AOPT的定义方法是一样的。

sdwangq

gdyu_yu 发表于 2012-5-3 17:49
AOPT的定义方法是一样的。

是一样的，不怕你笑话，我是新手，不明白的是该如何定义矢量，如果是a矢量是由其端点坐标与原点来定义，那我上述定义应当是正确的，我定义a矢量是笛卡尔坐标系的Z轴，如果不是这么定义，那又当如何定义矢量a啊，望指点

gdyu_yu

sdwangq 发表于 2012-5-4 16:35
是一样的，不怕你笑话，我是新手，不明白的是该如何定义矢量，如果是a矢量是由其端点坐标与原点来定义， ...

AOPT=2时，需要按照*DEFINE_COORDINATE_VECTOR的方法来定义。这个需要两个矢量。
你用的是壳单元还是实体单元？
参考资料如下：
1）关键字用户手册 *DEFINE_COORDINATE_VECTOR部分；
2）*DEFINE---->最底下EXAMPLES关于*DEFINE_COORDINATE_VECTOR
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
$$$$ *DEFINE_COORDINATE_VECTOR
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
$ Define local coordinate system number 4 using two vectors.
$ Vector 1 is defined from (0.0, 0.0, 0.0) to (1.0, 1.0, 0.0)
$ Vector 2 is defined from (0.0, 0.0, 0.0) to (1.0, 1.0, 1.0)
$ See the corresponding keyword command for a description.
$
*DEFINE_COORDINATE_VECTOR
$
$...>....1....>....2....>....3....>....4....>....5....>....6....>....7....>....8
$ cid Xx Yx Zx Xv Yv Zv
4 1.0 1.0 0.0 1.0 1.0 1.0
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
3）*MAT_002/*MAT_{OPTION}TROPIC_ELASTIC 图2.1

sdwangq

gdyu_yu 发表于 2012-5-4 20:21
AOPT=2时，需要按照*DEFINE_COORDINATE_VECTOR的方法来定义。这个需要两个矢量。
你用的是壳单元还是实体 ...

谢谢你的指点，我用的是壳单元，按照*MAT_002/*MAT_{OPTION}TROPIC_ELASTIC 图2.1的说明，我只需定义一个矢量a即可，所以在单元22中，需要在关键字文件中增加以下输入：按照实例定义两个矢量，然后在在矢量a的参数A1A2A3中输入相应的矢量ID？？而不是默认一个原点后在参数中输入矢量a的终点坐标，是不是能这么理解啊？？

gdyu_yu

由于你用的是壳单元，所以当AOPT=2时，在已经明确壳单元法向（法向为矢量c）前提下，确实只需要定义一个矢量a。矢量d无需定义。
矢量a,b,c的定义方式分别为：
c=n
b=cxa
a=bxc
A1,A2,A3 Components of vector a for AOPT = 2.
D1,D2,D3 Components of vector d for AOPT = 2.
在*MAT_22中关键字卡片中对应的只需要输入矢量a的分量(A1,A2,A3)的值，矢量d的分量(D1,D2,D3)都为0。
前面所说的，AOPT=2时，需要按照*DEFINE_COORDINATE_VECTOR的方法来定义。并不是说要定义*DEFINE_COORDINATE_VECTOR这个关键字，而是按这个关键字所述的定义方法来填写*MAT_22中对应矢量的分量。
你设置的K文件中，a为指向Z正轴的矢量。如果这个矢量不与壳单元法向(即c矢量)相冲突的话，是可以的。

sdwangq

gdyu_yu 发表于 2012-5-15 01:09
我看了你传给我的dyna2.k，有两个问题可以探讨下：
1、自动单面接触没有定义ssid和sstyp。（这个是否需要设 ...

谢谢你详细的回复，跟你的交流真是受益匪浅，希望能给你长期交流下去，谢谢！！常交流

gdyu_yu

ls625rayallen 发表于 2012-11-14 16:45
】也就是说 如果我做的是一个U形结构，那底面和侧面需要设置不同的MAT，不同的AOPT,不同的a向量，对吗 ...

我熟悉的是冷轧薄板。轧制方向（R0方向）已经是定下来的，也就是料卷的方向。如果a为轧向，则设a为(1,0,0)，此时b则是R90的方向。
同一材料，只有一个轧制方向。

ls625rayallen

你设置的K文件中，a为指向Z正轴的矢量。如果这个矢量不与壳单元法向(即c矢量)相冲突的话，是可以的。

】也就是说 如果我做的是一个U形结构，那底面和侧面需要设置不同的MAT，不同的AOPT,不同的a向量，对吗

gdyu_yu

ls625rayallen 发表于 2012-11-14 12:17
高手啊，我想问下，定义的abc材料的坐标系是初始状态的时候的吗？单元变形后材料的坐标系是自动跟着单元变 ...

abc是初始状态。

ls625rayallen

高手啊，我想问下，定义的abc材料的坐标系是初始状态的时候的吗？单元变形后材料的坐标系是自动跟着单元变形而改变的吗？

roczhp

ls625rayallen 发表于 2012-11-14 16:45
】也就是说 如果我做的是一个U形结构，那底面和侧面需要设置不同的MAT，不同的AOPT,不同的a向量，对吗 ...

曲面结构有点复杂。你用shell还是solid，一般用shell。你可以在lspp中的ident下选element，然后选shell，选mat dir查看材料方向，abc显示出来，很直观，连不同积分点处的材料方向都可以查看。看是不是和你实际的材料方向一致，不一致就要想办法重新定义了。如果用solid，一层复合材料一层单元，也可以在lspp的ident下看。

ccang33

楼上的大神们，小弟刚接触符合材料，望大神指点！qq：549037089，谢谢