keyopt的命令解释!.docVIP

终于明白keyopt是什么意思了！ 这几天我做ansys时，一直被一个小问题困扰着，就是设置接触时书中讲到keyopt（1－12），我不知道到底到哪里区设置，今天看到命令流 时，终于恍然大悟了，原来在GUI：main menupreprocessorElemant TypeAdd/Edit/Delete时候，options选择要更改的单元，即可看到每一行后有简写K1，k2等等，这就是对应的 keyopt1，2 （（单元关键字每种接触单元都饭知好几个关键字，对 大多的接触问题缺省的关键字是合适的，而在某些情况下，可能需要改变缺省值，来控制接触行为。接触算法（罚函数+拉格郎日或罚函数）（KEYOPT（2））出现超单元时的应力状态（DEYOPT（3））接触方位点的位置 （KEYOPI（4））刚度矩阵的选择 （KEYOPT（6））时间步长控制 （KEYOPT（7））初始渗透影响 （KEYOPT（9））接触表面情况 （KEYOPT（12））命令：KEYOPT）） 同理对于实常数的更改是GUI：main menu preprocessorreal constant时候，点击edit，选择需要更改的单元，然后会看到选项 实常数9 个实常数中，两个（R1 和 R2）用采定义目标面单元的几何形状，乘下的7 个用来控制接触行为。ANSYS 非线形分析指南 接触分析第8页R1 和R2 定义目标单元几何形状FKN 定义法向接触刚度因子FTOLN 定义最大的渗透范围ICONT 定义初始靠近因子PINB 定义“Pinball区域PMIN 和PMAX 定义初始渗透的容许范围TAUMAR 指定最大的接触摩擦命令：RGUI：main menu preprocessorreal constant对实常数 FKN，FTOLN，ICONT，PINB，PMAX，和 PMIN，你既可以定义一个正值也可以定义一个负值，程 序将正值作为比例因子，将 负值作为真实值，程序将下面覆盖原单元的厚度作为 ICON，FTOLN，PINB，PMAX 和 PMIN 的参考值，例如对 ICON，0.1 表明初始间隙因子是 0.1*下面覆盖层单元的厚度。然而，0.1表明真实缝隙是 0.1，如果下面覆盖层单元是超单元，则将接触单元的最小长度作为厚度。 真的是小问题也搞坏人 呵呵 选择接触算法：对面─面的接触单元，程序可以使用扩增的拉格朗日算法或罚函数方法，通过使用单元关键字KETOPT（2）来指定。扩张的拉格朗日算法是为了找到精确的拉格朗日乘子而对罚函数修正项进行反复迭代，与罚函数的方法相比，拉格朗日方法不易引起病态条件，对接触刚度的灵敏度较小，然而，在有些分析中，扩增的拉格朗日方法可能需要更多的迭代，特别是在变形后网格变得太扭曲时。使用拉格朗日算法的同时应使用实常数FTOLNFTOLN 为搠格朗日算法指定容许的最大渗艉，如果程序发现渗透大于此值时，即使不平衡力和位移增量已经满足了收敛准则，总的求解仍被当作不收敛处理，FTLON 的缺省值为0.1，你可以改变这个值，但要注意如果此值太小可能会造成太多的迭代次数或者不收敛 决定接触刚度所有的接触问题都需要定义接触刚度，两个表面之间渗 量的大小取决了接触刚度，过大的接触刚度可能会引起总刚矩阵的病态，而造成收敛困难，一般来说，应该选取足够大的接触刚度以保证接触渗透小到可以接受，但 同时又应该让接触刚度足够小以使不会引起总刚矩阵的病态问题而保证收敛性。程序会根据变形体单元的材料特性来估计一个缺省的接触刚度值，你能够用实常数 FKN 来为接触刚度指定一个比例因子或指定一个真正的值，比例因子一般在0.01 和10 之间，当避免过多的迭代次数时，应该尽量使渗透到达极小值。为了取得一个较好的接触刚度值，又可需要一些经验，你可以按下面的步骤过行。1、开始时取一个较低的值，低估些值要比高估些值好因为由一个较低的接触刚度导致的渗透问题要比 过高的接触刚度导致的收敛性困难，要容易解决。2、对前几个子步进行计算3、检 查渗透量和每一子步中的平衡迭代次数，如果总体收敛困难是由过大的渗透引起的（而不是由不平衡力和位移增量引起的），那么可能低估了FKN 的值或者是将FTOLN 的值取得大小，如果总体的收敛困难是由于不平衡力和位移增量达到收敛值需要过多的迭代次数，而不是由于过大的渗透量，那么FKN 的值可能被高估。 选择摩擦类型。在基本的库仑摩擦模型中，两个接触面在开始相互滑动之前，在它们的界面上会有达到某一大小的剪应力产生，这种状态则作粘合状态（stick)库仑摩擦模型定义了一个等效剪应力。）一旦剪应力超过此值后，两个表面之间将开始相互滑动，这种状态，叫作滑动状态（Sliding）粘合\滑动计算决定