首页/文章/ 详情

04.接触的搜索方式

4小时前浏览4

    在接触算法中,为了判定一个从节点是否发生接触,首先应当对它周围的主段进行搜索,发现最近的主段进行接触判定。


LS-DYNA有三种接触搜索方式

● 基于节点增量的搜索方式

● 基于段的bucket搜索方式

● 基于面段的搜索方式。


【1】基于节点增量的搜索方式

    增量搜索方式是最早采用的方法,对每一个从节点,搜索最近的主节点,找到最近的主点后,再搜索主节点周围最近的主段,检查若的渗透发生,施加接触力。


增量法搜索过程



所有的非自动接触方式都采用增量搜索方式。

优点:该方法简单高效;

缺点:要求主面的网格必须是连续的,对于不连续的的主面将会出现漏搜索。

    如下图,主面不连续导致段2和段3将不会被搜索,而可能实际情况就是段2或段3与从点发生接触。


主面不连续


    另外,对于网格差别较大和形状复杂的主面来说,增量搜索方式可能会产生错误的接触行为,如下图所示,搜索到的最近的主点并不包含最近的主段。



【2】基于段的bucket分类搜索方式

    基于段的bucket分类搜索方式是新的接触搜索方式,主要是为了解决增量搜索方式的缺陷。

原理是对于每个从点,搜索最近的主段,首先把接触空间划分成很多小的bucket,对于每个从节点,按照它的空间位置指定于某个bucket中,然后在该范围内搜索最近的主段,运行一定的时间后(比如100个时间循环),根据该节点空间位置的变化再重新指定新的bucket。


一维bucket分类搜索:


    上图中的黑点代表接触段的质心点(通过该质心点定位主段),对于落入某bucket的从节点而言,搜素该bucket中的段,找到最近的主段,若没有找到,再搜索附近的bucket,直到找到最近的主段。这种方法由于不断更新bucket分类,所以要耗费更多的CPU时间。


二维bucket分类搜索:

    搜索原理同一维bucket分类搜索方式;



三维bucket分类搜索:

    搜索原理同一维bucket分类搜索方式;



增量法和Bucket:

    如下图所示,若应用基于节点的增量搜索方式,找到最近的主节点,与该节点相关联的主段是1,2,3,但都不包含该从点,所以接触搜索失效;但若用基于段的bucket 分类搜索方式,可以很方便的找到包含该从点的主段,建立接触关系。



【3】基于面段的搜索方式

    SOFT=2调用了基于面段的接触算法,这种算法起源于Belytschko提出的弹球接触(Pinball contact)。

    在这种算法下,程序会计算面段之间的接触问题,而不是传统的节点-面段接触问题。当两个四节点的面段相互接触时,程序向八个节点施加接触力以避免面段之间的穿透。这种计算方法可以实现更合理的接触力分布,而且对一些特别疑难的接触问题也很有效。



常用接触的搜索方式:

    接触搜索方式由定义的接触类型(Automatic)和参数选择(SOFT)决定,在接触计算中不会发生变化。

① 若定义的接触类型为非Automatic接触,则采用基于节点增量的搜索方式;

② 若定义的接触类型为Automatic接触,则采用基于段的bucket搜索方式;

③ 若在接触定义中选择SOFT=2,则采用基于面段的搜索方式。



     

     
END
     

    

来源:CAE碰撞仿真指导
ACTLS-DYNA
著作权归作者所有,欢迎分享,未经许可,不得转载
首次发布时间:2025-06-24
最近编辑:4小时前
CAE碰撞仿真指导
本科 | CAE高级仿真工... CAE 碰撞仿真指导
获赞 6粉丝 5文章 26课程 0
点赞
收藏
作者推荐

07.应力波传递描述

显式算法在仿真分析中比较常见,经常应用在瞬态动力学、准静态问题分析中;显式求解算法基于动力学方程,当前时刻的位移只与前一时刻的速度和位移相关,求解过程中无需迭代。本节试图让读者从概念上理解应用显式动态算法时应力是如何在模型中传播的。应力波传播过程: 为了简单且明了的描述应力波在显式算法中如何传递的,本文以杆件为受力对象进行分析。 如图1所示,描述应力波沿着一个包含三个单元(三个节点)的杆件模型传播的过程,同时研究杆件随着时间增长的各个状态,以及受力过程;左侧加集中力P,右侧约束。图1 自由端作用有集中力的杆件的初始状态单元①应力求解过程: 在第一个时间增量段,施加在节点1的集中力P的作用结果是使节点1具有一个加速度。这个加速度引起节点1产生速度,接着在单元1内引起应变速率。沿第一个时间增量段内对应变速率进行积分获得了单元1的应变增量。总应变,是初始应变和应变增量的和。此问题中初始应变为零。一旦单元的应变求出之后,可以通过材料的本构关系求出单元的应力。求解过程 对于线弹性材料,应力就是弹性模量与总应变的乘积。这个过程如图2所示。在第一个时间增量段,节点2和3因为没有力作用其上,所以没有移动。图2 自由端作用有集中力的杆件在第一时间增量段结束时的状态单元①计算过程:单元②应力求解过程: 在第二个时间增量段,由单元1的应力得到内力,施加到与单元1相连的节点上,如图3所示。这些单元应力随后用于计算节点1和节点2的动力平衡方程。图3 杆件在第二个时间增量段开始时的状态单元②计算过程:单元③应力求解过程: 这个过程继续下去,到第三个时间增量段开始时,单元1和2已存在应力,以及节点1、2、和3存在了作用力,如图3.4所示。这个过程继续下去,直到总的分析时间结束。图4 杆在第三时间增量段开始时的状态 END 来源:CAE碰撞仿真指导

未登录
还没有评论
课程
培训
服务
行家
VIP会员 学习计划 福利任务
下载APP
联系我们
帮助与反馈