内容发布更新时间 : 2024/12/23 23:22:44星期一 下面是文章的全部内容请认真阅读。
文档来源为:从网络收集整理.word版本可编辑.欢迎下载支持.
GUI:main menu>preprocessor>Element type>Add/Edit/Delete 实常数和材料特性
在定义了单元类型之后,需要选择正确的实常数集。一个接触对中的接触面和目标面必须有相同的实常数号。每个接触对必须有不同的实常数号。 ANSYS 使用下伏单元的材料特性来计算一个合适的接触(或罚)刚度。在下层单元有用 TB 命令定义的塑性材料特性(不论激活与否)的情况下,接触的法向刚度可能按照系数100降低。ANSYS 自动为切向(滑动)刚度定义一个与 MU 和法向刚度成正比的缺省值。如果下伏单元是一个超单元,接触单元的材料必须与超单元形成时的原始结构单元相同。 生成接触单元
既可以通过直接生成法生成接触单元,也可以在下层单元的外表面上自动生成接触单元。推荐采用自动生成法,这种方法更为简单和可靠。可以通过下面三个步骤,来自动生成接触单元。
1、选择节点
选择已分网的柔体表面的节点。对每一个面,检查节点排列。如果用户确定某一部分节点永远不会接触到目标面,用户可以忽略它以便减少计算时间。然而用户必须保证设有漏掉可能会接触到目标面的节点。
命令: NSEL
GUI:Utility Menu>Select>Entities 2、生成接触单元 命令: ESURF
GUI:Main Menu>Preprocessor>Create>Elements>Surf to Surf
如果接触单元是附在已用实体单元划分网格的面或体上,程序会自动决定接触计算所需的外法向。如果下层单元是梁或壳单元,则必须指明哪个表面(上表面或下表面)是接触面。
命令: ESURF ,,TOP (或 BOTIOM)
GUI:Main Menu>Preprocessor>Create>Elements>Surf to Surf
使用 TOP(缺省)生成接触单元,它们的外法向与梁或壳单元的法向相同;使用 BOTIOM 生成接触单元,则它们的外法向与梁或壳单元的法向相反。必须确保
16
文档来源为:从网络收集整理.word版本可编辑.欢迎下载支持.
梁上的单元或壳单元有一致的法向。如果下伏单元是实体单元,则 TOP 或 BOTTOM 选项不起作用
3、检查接触单元外法向。当程序进行是否接触的检查时,接触面的外法线方向至关重要。对于3D单元,按节点顺序号以右手法则来决定单元的外法向。接触面的外法向应该指向目标面。否则,在开始分析计算时,程序可能会认为是有过度穿透的面,而很难找到初始解。在这些情况下,程序一般会立即停止执行。 图5-7 说明正确和不正确的外法向。
命令:/ PSYMB ,ESYS
GUI:Utility menu>plotctrls>symbols
图5-7 定义接触单元的外法向
当发现单元的外法线方向不正确时,必须通过反转所选择的不正确单元的节点号来改变它们:
命令: ESURF ,,REVE
GUI:Main Menu>Preprocessor>Create>Elements>Surf to Surf 或重新定向单元法向: 命令: ENORM
GUI:Main Menu>Preprocessor>Create>Move/Modify>Shell Normals
5.4.9 设置实常数和单元关键选项
程序使用20个实常数和数个单元关键选项,来控制面─面接触单元的接触。参见《ANSYS Elements Reference》中对接触单元的描述。 实常数
在20个实常数中,两个(R1和R2)用来定义目标面单元的几何形状。剩下的用来控制接触面单元。
R1和R2 定义目标单元几何形状。 FKN 定义法向接触刚度因子。
FTOLN 是基于单元厚度的一个系数,用于计算允许的穿透。 ICONT 定义初始闭合因子。
17
文档来源为:从网络收集整理.word版本可编辑.欢迎下载支持.
PINB 定义“Pinball\ PMIN和PMAX 定义初始穿透的容许范围。 TAUMAR 指定最大的接触摩擦。
CNOF 指定施加于接触面的正或负的偏移值。 FKOP 指定在接触分开时施加的刚度系数。 FKT 指定切向接触刚度。 COHE 制定滑动抗力粘聚力。
TCC 指定热接触传导系数。 FHTG 指定摩擦耗散能量的热转换率。 SBCT 指定 Stefan-Boltzman 常数。 RDVF 指定辐射观察系数。
FWGT 指定在接触面和目标面之间热分布的权重系数。 FACT 静摩擦系数和动摩擦系数的比率。 DC 静、动摩擦衰减系数。 命令: R
GUI:main menu> preprocessor>real constant
对实常数 FKN, FTOLN, ICONT, PINB, PMAX, PMIN, FKOP 和 FKT,用户既可以定义一个正值,也可以定义一个负值。程序将正值作为比例因子,将负值作为绝对值。程序将下伏单元的厚度作为ICON,FTOLN,PINB,PMAX 和 PMIN 的参考值。例如 ICON = 0.1 表明初始闭合因子是“0.1*下层单元的厚度”。然而,ICON = -0.1 则表示真实调整带是 0.1 单位。如果下伏单元是超单元,则将接触单元的最小长度作为厚度。参见 图5-8 。
图5-8 下层单元的厚度
在模型中,如果单元尺寸变化很大,而且在实常数如 ICONT, FTOLN, PINB, PMAX, PMIN 中应用比例系数,则可能会出现问题。因为从比例系数得到的实际结果,取决于下层单元的厚度,这就可能引起大、小单元之间的重大变化。如果出现这一问题,请用绝对值代替比例系数。
18
文档来源为:从网络收集整理.word版本可编辑.欢迎下载支持.
TCC, FHTG, SBCT, RDVF 和 FWGT 仅用于热接触分析[KEYOPT(1)=1]。 单元关键选项
每种接触单元都包括数个关键选项。对大多的接触问题,缺省的关键选项是合适的。而在某些情况下,可能需要改变缺省值。下面是可以控制接触行为的一些关键选项:
自由
度 KEYOPT(1)
接触算法(罚函数+拉格朗日乘子或罚函数) KEYOPT(2) 存在超单元时的应力状态(仅2D) KEYOPT(3) 接触检测点的位置(仅低阶接触单元) KEYOPT(4) CNOF自动调
整 KEYOPT(5)
时间步控制 KEYOPT(7)
伪接触预防 KEYOPT(8) 初始穿透或间隙的影响 KEYOPT(9)
法向和切向接触刚度修正方法控制 KEYOPT(10) 壳的厚度影响 KEYOPT(11) 接触面行为(粗糙、绑定等) KEYOPT(12) 命令: KEYOPT ET
GUI:main menu>preprocessor>Elemant Type>Add/Edit/Delete 选择接触算法
Contact wizard basic>contact algorithm
对面─面接触单元,程序可以使用增进的拉格朗日方法或罚函数方法。通过单元关键字 KEYOPT(2)来指定。
19
文档来源为:从网络收集整理.word版本可编辑.欢迎下载支持.
增进的拉格朗日方法是为了找到精确的拉格朗日乘子(即接触力),而对罚函数进行一系列修正迭代。与罚函数的方法相比,拉格朗日方法容易得到良态条件,对接触刚度的敏感性较小。然而,在有些分析中,增进的拉格朗日方法可能需要更多的迭代,特别是在变形后网格变得太扭曲时。
使用拉格朗日方法的同时应使用实常数 FTOLN。FTOLN 为拉格朗日方法指定容许的最大穿透。如果程序发现穿透大于此值时,即使不平衡力和位移增量已经满足了收敛准则,总的求解仍被当作不收敛处理。FTLON 的缺省值为0.1。用户可以改变这个值,但要注意,如果此值太小,可能会造成太多的迭代次数或者不收敛。 确定接触刚度
Contact Wizard>basic>normal penaltycontact stiffness
所有的接触问题都需要定义接触刚度,两个表面之间穿透量的大小取决于接触刚度。过大的接触刚度可能会引起总刚矩阵的病态,从而造成收敛困难。一般来说,应该选取足够大的接触刚度以保证接触穿透小到可以接受,但同时又应该让接触刚度足够小以不致引起总刚矩阵的病态而保证收敛性。
ANSYS 程序根据下伏柔体单元的材料特性,来估计一个缺省的接触刚度值。用户可用实常数 FKN 来为接触刚度指定一个比例因子或指定一个绝对值。比例因子一般在0.01和10之间;对于大变形问题,选1是比较好的;而对于弯曲为主的问题,通常为0.01~0.1。用户应当总是检验以使穿透到达极小值,而又避免过多的迭代次数。
注意 --FTOLN 和 FKN 从一个荷载步到另一个荷载步中,都可以修改。也可以在重启动中修改。这时,必须定义KEYOPT(10)=1,2。
为了确定一个较好的接触刚度值,可能需要一些经验。用户可以按下面的步骤来进行尝试:
1、 开始时取一个较低的值。低估值要比高估值好,因为由一个较低的接触
刚度导致的穿透问题,比过高的接触刚度导致的收敛性困难,要容易解决。 2、 对前几个子步进行计算分析,直到最终荷载的一个比例(刚好完全建立
接触)。 3、 检查每一子步中的穿透量和平衡迭代次数。如果总体收敛困难是由过大
的穿透引起的(而不是由不平衡力和位移增量引起的),那么可能低估了FKN 的值,或者是将 FTOLN 的值取得大小。如果总体的收敛困难是由于不平衡力和位移增量达到收敛值时需要过多的迭代次数,而不是由于过大的穿透量引起的,那么 FKN 的值可能被高估。
20