ANSYS结构非线性分析指南 - 第四章(6)

ANSYS结构非线性分析指南

图4-15 在应力和应变输入数组中的数据位置

现在考虑一个得到了单轴拉伸和剪切测试数据时的情况。在数组中存储应变和应力的命令可能象下面那样(当然，数组可以有任意的名字，而且在本例中用N1和N2表示的数据点数，可以是任意整数)：

! Uniaxial Tension Data

*SET,STRAIN(1,1), ... ! First 10 strain data points

*SET,STRAIN(11,1), ...! Strain data points 11 through N1 (if N1<21) *SET,STRESS(1,1), ... ! First 10 stress data points

*SET,STRESS(11,1), ...! Stress data points 11 through N1 ! Shear Data

*SET,STRAIN(1,3), ... ! Strain data points 1 through N2 (if N2<11) *SET,STRESS(1,3), ... ! Stress data points 1 through N2

参见 *SET 命令的说明。

第三步：计算Mooney-Rivlin常数

要自动生成Mooney-Rivlin常数，首先应执行 TB 命令，并使该命令的Lab=MOONEY，TBOPT=1。然后，执行 *MOONEY 命令，将已生成好的数组名填入其中(用户可以给这些数组任意有效的参数名，在这里所采用的参数名，只是为了讨论方便)：

TB,MOONEY,MAT,NTEMP,,1

*MOONEY,STRAIN(1,1),STRESS(1,1),,CONST(1),CALC(1),SORTSN(1), SORTSS(1),Fname,Ext

程序自动计算出Mooney-Rivlin常数，将它们写入数据库和 CONST 数组(可以是任何有效的数组名)中，并还以 TB 和 TBDATA 命令的格式写入一个ASCII文件Fname.Ext中(缺省是Jobname.TB)。

单轴公式将用于第一列的数据，而等双轴公式将用于第二列，平面(纯剪)公式将用于第三列。

注意--在 STARIN 和 STRESS 中输入的所有试验数据，将用于确定Mooney-Rivlin 超弹性材料常数。

第四步：估计Mooney-Rivlin常数的质量

4-26

ANSYS结构非线性分析指南

在输出文件中(Jobname.OUT)，检查“ROOT—MEAN—SQUARE ERROR(PERCENTAGE)”(均方根误差)和“COEFFICIENT OF DETERMINATION”(确定系数)两个输出信息，这两个值给出所计算出的应力-应变曲线与测试数据点相拟合的好坏程度的统计度量。以百分数表示的均方根误差(即2.5就表示2.5%应接近于零，确定系数接近于1.0(通常要大于0.99)。

另外，用户还应使用 *VEAL 和 *VPLOT 命令(GUI:Main Menu> Preprocessor>Material Props>Mooney-Rivlin>Evaluate Const 和 Utility Menu>Plot>Array Parameters)来以图形的方式显示输入的和计算出的应力-应变曲线，以直观地检查计算曲线与实验数据的匹配程度。在比较这些曲线时，应比较那些代表同一变形模式的数据。也就是说，所计算的单轴拉伸曲线形状(在 *EVAL 命令中的 EVPARM = 1 )，仅应当与单轴拉伸数据(在排序的STRAIN 和 STRESS数组第一列)比较。类似地，所计算的单轴压缩曲线形状，仅应当与单轴压缩数据比较；而所计算的剪切曲线形状，仅应当与剪切数据比较。

当用图形显示计算的应力-应变曲线时，用户可以将显示曲线扩展到那些没有实验数据的区域。这样就可以让用户对模型在超出实验数据区域以外的地方的响应有一个定性的认识。但是，应认识到当将显示曲线扩展到一个表示另一个不同的变形模式的区域时，在那个区域的显示就是没有意义的。例如，用户只能在正应变区域显示单轴拉伸曲线。总的来说，要得到一个好的结果，所作的实验数据应能代表所分析模型的所有的变形模式和响应(应变)范围。 *MOONEY 命令自动把Mooney-Rivlin常数写到 CONST 数组中。因为 *EVAL 命令从 CONST 常数读入相同的常数，所以可以在同一个ANSYS 阶段中，在 *MOONEY 命令后，跟着用 *EVAL 命令。如已经有了Mooney-Rivlin常数(这时不必进行 *MOONEY 计算)，则必须在计算曲线前，首先定义CONST 数组 [*DIM]并用Mooney-Rivlin常数填充该数组[*EVAL]。可以方程容易地填充这一数组，通常给 CONST 数组1×2、1×5、1×9。用户也可以在Jobname.TB文件中，添加 *DIM 和数组填充命令，以方便操作。

要检查曲线的形状，首先还必须定义[*DIM]两个表数组向量(*EVAL命令中将这两个数组名定义为 XVAL 和 ECALC，但可以应用任何有效的参数名)。这两个表数组的大小都是P维，此时P为用户想要在曲线中绘制的点的数目(通常应使用一个相对较大的P值，以使得所绘制的曲线尽可能光滑)。其次，定义变形模式、定义应变范围，用 *EVAL 命令将工程应变和计算出的工程应力数据填入数组中。最后，用 *VPLOT 命令来绘制计算出的应力-应变曲线。下面例子说明对单轴压缩变形模式的计算曲线绘图：

! Dimension strain and stress arrays for the calculated curve: ! (Any valid parameter names can be used) *DIM,XVAL,TABLE,1000 *DIM,ECALC,TABLE,1000

! Specify the mode of deformation (EVPARM), define the strain range

! (XMIN,XMAX), and use the M-R constants (CONST) to fill the strain (XVAL) ! and stress (ECALC) arrays with calculated data: *EVAL,1,2,CONST(1),XMIN,XMAX,XVAL(1),ECALC(1) ! Label the graph axes: /AXLAB,X,Engineering Strain /AXLAB,Y,Engineering Stress

! Plot the calculated uniaxial compression curve: *VPLOT,XVAL(1),ECALC(1)

参见 *DIM, *EVAL, /AXLAB, 和 *VPLOT 等命令的说明。

4-27

ANSYS结构非线性分析指南

图4-16 典型的计算出的超弹性应力-应变曲线

第五步：使用Mooney-Rivlin常数

如果用户对曲线拟合的统计数值和整个曲线的形状都非常满意，则可使用生成的Mooney-Rivlin材料性质作后续的分析(*MOONEY 命令将在数据库中保存这些常数)。在将来的另外一些分析中，如果使用同样的材料模式，则只需用/INPUT 命令读入文件“Jobname.TB”，即可将常数加载到新的数据库中。但别忘了定义材料泊松比的值[MP,NUXY,...]。用户始终要记住的一点是，应对分析结果作仔细检查，以确定原始测试数据是否涵盖了模型的变形模式和最大应变。

含有超弹性单元的分析，有时对材料性质定义和载荷施加方式非常敏感。通常是某些Mooney-Rivlin常数会导致非常稳定的刚度矩阵而其它的常数则相反。因此，应当根据经验，并特别仔细地选取常数。

ANSYS内部提供超弹性材料(基于用户输入的Mooney-Rivlin 常数)稳定检查的功能。这些检查分两个级别：

? 在分析前进行第一次稳定检查。对6个典型的应力路径(单轴拉伸和压缩，等双轴拉伸和

压缩，平面拉伸和压缩)，延伸率范围0.1～10来进行检查。如果材料在这一范围不稳定，则出现一个提示，统计在材料变得不稳定时，名义应变的临界值，列出用户输入的Mooney-Rivlin 常数。如果材料在这一范围稳定，则不出现提示。下面的警告信息，列出材料1变得不稳定时的名义应变，然后列出用户输入的Mooney-Rivlin 常数：

*** WARNING *** CP= 1.110 TIME= 16:59:52 Material 1 can become unstable under certain loading.

The strain (nominal) limits where the material becomes unstable are: UNIAXIAL TENSION 0.645E+00 UNIAXIAL COMPRESSION -0.565E+00 EQ …… 此处隐藏：4587字，全部文档内容请下载后查看。喜欢就下载吧 ……