《ANSYS材料非线性分析》由会员分享，可在线阅读，更多相关《ANSYS材料非线性分析（31页珍藏版）》请在金锄头文库上搜索。
1、【分享】ANSYS7.0 超弹材料的定义新的曲线拟合功能-摘自 ansys 用户专区几何非线性几何非线性不受敛主要原因1.网格质量，特别是 warpage2.约束方程，少用刚性连接3.收敛准则，可适当加大容差4.荷载步设置，可适当加大步数最近碰到一个对我来说很意外的问题：如果确实如此希望大家以后小心大家知道定义接触后会自动生成一组实常数，前几天我碰到一个问题，需定义超过 10 组实常数，接触对很多，好像有 20 多处，按照常规步骤划分完所有网格，当时因为有一个实常数参数没确定，便预留了最后一组（第 10 组）实常数里面的参数为空，接下来就定义了所有的接触对，由于所有接触对里的设置一样，ANSY
2、S 在我保存 db 完重新打开后便把我所有的接触对综合成一个了！接下来我就把第十组实常数里面的参数补上了，但在求解时却提示我该实常数同时被两种单元（包括 CNTACT 单元）同时占用，出现错误！检查了半天才发现自动生成的接触对实常数把第 10 组实常数也占用了！我实在没找到什么好的解决办法，只得把接触对删除了重新定义，那可是上百多个面的选取过程，痛苦不堪简直！ANSYS 里接触对面的选取时还不能针对 Component 操作！ANSYS7.0 超弹材料的定义新的曲线拟合功能 ANSYS7.0 中的超弹材料模拟能力得到了很大的加强，在 ANSYS6.1 的超弹材料模型的基础上又增加了 Gent,
3、 Yeoh, Blatz-Ko, and Ogden (Foam)四种超弹性材料模型，使得其超弹模拟能力得到了进一步扩展。 ANSYS7.0 中对超弹能力最吸引人的增强还不在于此，而是在于其曲线拟合能力的大幅度扩展，不再像ANSYS6.1 以前的版本一样曲线拟合仅仅局限于 Mooney-Rivlin 模型，而是将其扩展到所有的超弹模型，这样，用户可以利用实验得到的应力应变数据直接让程序自己拟合出任意一种超弹材料模型的参数，大大方便了用户的使用。以下就 ANSYS7.0 的超弹拟合功能做一简单介绍。 在 ANSYS6.1 中，*mooney 命令或对应的 GUI 菜单可以根据试验数据拟合 Moo
4、ney-Rivlin 超弹材料参数，但只局限于这一种模型，而且其使用也受到了限制，其数据只能用于 hyperXX 单元，若要用于 18x 单元还需要将拟合出的参数直接输入 tb hyper 对应的超弹模型中，由于其使用的不方便性，Ansys7.0 的超弹模型中剔除了 Ansys6.1 版本中的*mooney 命令对应的曲线拟合菜单而增加了新的可适用于多个超弹模型的拟合功能，但*mooney 命令仍然存在，也就是说如果您需要使用经典的拟合方法，可以定义相应的数组参数，然后用命令流的方式进行拟合。 其实，在 Ansys7.0 中不再推荐使用 *mooney 命令来进行拟合，而建议采用新的拟合技术,
5、 可以进行多个超弹模型的拟合，而且可以直接用于 18x 单元。其拟合菜单的路径为： Material-Nonlinear-Elastic-Hyperelastic-Curve Fitting 然后按照 wizard 的提示一步步输入单轴、双轴、剪切、体积试验数据文本文件名称，如果没有任何一种试验数据，只需将该名称处空置即可，最后选择需要拟合数据的超弹模型，程序就会自动计算出相应的参数，并立刻在图形窗口显示拟合曲线与试验曲线的比较图，如果不理想，可以点击 prev 回到前面的步骤重新选取模型，如果拟合结果满意，则点击 update，拟合出来的材料参数就会被输入激活的材料号中，使用起来非常方便。
6、试验数据的文件格式需要进行说明： 单轴、等双轴、剪切数据为应变应力数据，依次输入应变、应力值，第一列为应变，第二列为应力，每一行两个数之间用空格隔开（空格数目不限） ，代表一个数据点。数据符号：拉为正，压为负。体积试验数据有所不同，每一行两个数据同样用空格隔开，第二个数为静水压力，但第一个数不是体积应变 e，而是相对体积，即 1+e，例如体积应变为-0.01，则应该在第一列输入 0.99。而且程序假定体积试验为线性关系，拟合时也只拟合出一条直线，所以通常输入两个值即可。数据符号：使体积减小的静水压力为正，这需要注意，如下例： 0 0 0.99 20 在选用 mooney 模型的情况下，单轴、等
7、双轴、剪切试验数据用来计算参数 c1-c9，而体积数据用来计算最后一个数据 d（不可压缩性因子，与泊松比有关） ，如果不给出体积试验曲线，则 d=0，这种情况下可能需要手工修改 d 值，如果知道泊松比 u，则可用如下公式计算： d=(1-2u)/(c1+c2) 此式是建立在几乎不可压缩（u 接近于或等于 0.5）的前提下的。因为准确的公式应该如下： 剪切模量 G=2(c1+c2) 体积模量 k=E/(3(1-2u) G=E/(2(1+u) k=2/d 于是 d=2/k=6(1-2u)/(4(1+u)(c1+c2) 可以看到，如果 u 约等于 0.5，则上式可以简化为前面的式子。 mooney
8、还有其他一些推导基于几乎不可压缩，所以对 mooney 模型而言，输入参数时一定要注意其泊松比应该接近 0.5，一般大于 0.49。4.2.4.4 求解步骤（命令流方法）fini/clel=60l1=15w=10ri=50ro=160routn=200pp=600pp1=30pex=5000/prep7et,1,42et,2,14keyopt,1,3,2keyopt,2,3,2确定 Mooney-Rivlin 常数的个数 概略地说，数据点的个数(即上面的 N )应至少为 Mooney-Rivlin 常数个数的两部。常数个数越多，曲线的统计量越与真实值相接近(即拟合得更好)，但曲线的形状可能会比
9、常数个数少的曲线要差。鉴于此，用户可以按顺序尝试 2 项、5 项、9 项函数，并检查其所生成的应力-应变曲线以确定到底哪一个函数在综合曲线形状以及拟合质量两方面做得最好。表 4-1 建议的 Mooney-Rivlin 常数应力-应变曲线中的点数 建议的 Mooney-Rivlin 函数无拐点(即单个曲线) 2 项 一个拐点(即 2 条曲线) 5 项 2 个拐点 9 项 图 4-14 典型的超弹性应力-应变曲线输出应力数组(CALC) ：其大小为 N 3， N 与前面的一样，该数组保存计算出来的应力值，这些应力值保存的顺序与其相应的应变保存值的顺序一致(后者以升序排列) 。排序应变数组(SORT
10、SN) ：其大小为 N 3，它保存经过排序的输入应变。排序应变数组(SORTSN) ：其大小为 N 3，它保存经过排序的输入应力。例如，如果任一个类型的测试数据包含直到 20 个数据点，而想生成 5 项的 Mooney-Rivlin 常数，则可以应用下面的命令来定义所需要的数组( 记住：用户可以代入任意有效的参数名)：*DIM,STRAIN,20,3 ! Dim. array (STRAIN) for 20 input strain-data points*DIM,STRESS,20,3 ! Dim. array (STRESS) for input stress data (20 pts.)
11、*DIM,CONST,5,1 ! Dim. array (CONST) for 5-term M-R constants*DIM,CALC,20,3 ! Dim. array (CALC) for sorted calculated stresses*DIM,SORTSN,20,3 ! Dim. array (SORTSN) for sorted input strain data*DIM,SORTSS,20,3 ! Dim. array (SORTSS) for sorted input stress data参见 * DIM 命令的说明。第二步：填充输入数据数组 当数组定义完成之后，就可以
12、用 * SET 命令(GUI：Utility Menu ParametersArray Parameters)将实验数据填入 STRAIN 数组和 STRESS 数组。请再一次记住，用户可以给这些数组任意有效的参数名；在这里所采用的参数名，只是为了讨论方便。 注意-*MOONEY 命令将所有输入应力和应变都解释为工程应力和工程应变。 这些数组大小都为 N *3，其每一列各自表示一类测试数据，顺序是：第一列：单轴拉伸和/或单轴压缩；第二列：等双轴拉伸和/或等双轴压缩；第三列：剪切(平面拉伸或压缩 )。注意这并不是说，变形模式及其等同存在 1：1 的关系。第一个变形模式-单轴拉伸- 与等效双轴压缩
13、等同，但数组第一列包含从单轴拉伸和/或单轴压缩得到的数据。类似地，第二个变形模式-等双轴拉伸-与单轴压缩等同，但数组第二列包含从等双轴拉伸和/ 或等双轴压缩得到的数据。表 4-2 应力- 应变输入数组中的数据位置变形模态 等效测试类型 测试数据的在数组中的位置单轴拉伸 单轴拉伸等双轴压缩 第一列第二列等双轴拉伸 等双轴拉伸单轴压缩 第二列第一列 剪切 平面拉伸平面压缩 第三列第三列如果只作了一种或两种测试，则须将未做的测试的相应的列置空。可用 图 4-15 来说明。图 4-15 在应力和应变输入数组中的数据位置现在考虑一个得到了单轴拉伸和剪切测试数据时的情况。在数组中存储应变和应力的命令可能象
14、下面那样(当然，数组可以有任意的名字，而且在本例中用 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 PreprocessorMaterial PropsMooney-RivlinEvaluate Const 和 Utility MenuPlotArray Parameters)来以图形的方式显示输入的和计算出的应力 -应变曲线
15、，以直观地检查计算曲线与实验数据的匹配程度。在比较这些曲线时，应比较那些代表同一变形模式的数据。也就是说，所计算的单轴拉伸曲线形状(在 * EVAL 命令中的 EVPARM = 1 )，仅应当与单轴拉伸数据( 在排序的 STRAIN 和 STRESS 数组第一列)比较。类似地，所计算的单轴压缩曲线形状，仅应当与单轴压缩数据比较；而所计算的剪切曲线形状，仅应当与剪切数据比较。当用图形显示计算的应力-应变曲线时，用户可以将显示曲线扩展到那些没有实验数据的区域。这样就可以让用户对模型在超出实验数据区域以外的地方的响应有一个定性的认识。但是，应认识到当将显示曲线扩展到一个表示另一个不同的变形模式的区域
16、时，在那个区域的显示就是没有意义的。例如，用户只能在正应变区域显示单轴拉伸曲线。总的来说，要得到一个好的结果，所作的实验数据应能代表所分析模型的所有的变形模式和响应(应变)范围。* MOONEY 命令自动把 Mooney-Rivlin 常数写到 CONST 数组中。因为 * EVAL 命令从 CONST 常数读入相同的常数，所以可以在同一个 ANSYS 阶段中，在 * MOONEY 命令后，跟着用 * EVAL 命令。如已经有了 Mooney-Rivlin 常数( 这时不必进行 * MOONEY 计算)，则必须在计算曲线前，首先定义 CONST 数组 * DIM 并用 Mooney-Rivlin 常数填充该数组 * EVAL 。可以方程容易地填充这一数组，通常给 CONST 数组12、15、19。用户也可以在 Jobname.TB 文件中，添加 * DIM 和数组填充命令，以方便操作。要
