首页/文章/ 详情

10.阻尼介绍

5小时前浏览2

LS-DYNA是一款高级的有限元分析软件,主要用于模拟各种复杂的动力学行为。在LS-DYNA中阻尼的设置对于确保模拟的准确性和稳定性至关重要。


能量耗散可以发生在*damping以外的其他方式。例如沙漏力引起的能量、刚墙力引起的能量、接触摩擦力产生的能量、阻尼器产生的内部能量等。


   

1.体积粘性阻尼 (Bulk Viscosity Damping)      

   
  • 作用:主要用于控制单元在高速变形(如冲击波传播)中的高频数值振荡,避免显式算法中可能出现的数值不稳定。

  • 类型:

① 线性体积粘性阻尼:与体积应变率成正比,通过参数 QB 控制。

② 二次体积粘性阻尼:与体积应变率的平方成正比,通过参数QC控制。

  • 设置方法

*CONTROL_BULK_VISCOSITY 关键字中调整QB 和QC 的默认值(通常无需修改,除非特殊需求)。

  • 注意事项过大的体积粘性阻尼会引入非物理能量耗散,影响结果精度。


     

2. 材料阻尼 (Material Damping)      

     
  • 作用模拟材料本身的粘弹性行为或能量耗散特性。

  • 常见形式

① 粘性阻尼通过阻尼系数(如 DAMP)定义,与应变率相关。

② 瑞利阻尼 (Rayleigh Damping)与质量和刚度矩阵成比例。    


  • 设置方法:在材料定义中直接指定阻尼参数。


image.png    


       

3. 结构阻尼 (Structural Damping)        

       
  • 作用通过全局或局部阻尼系数模拟结构的能量耗散。

  • 类型

【1】质量阻尼 :与质量矩阵成比例。

质量阻尼包括*damping_global*damping_part_mass,用于低频的结构阻尼模式,但它具有阻尼刚体模式的附加效果。因此,应将经受显著刚体运动的零件排除在质量阻尼之外,否则应在该零件经受刚体运动时关闭质量阻尼或使用*damping_relative。


通过使用*damping_relative,仅对特定刚体的运动/振动起阻尼作用。*damping_relative和常规质量阻尼一样,除了相对于RB(旋转)运动的运动被阻尼。在确定相对运动时考虑刚体的平移和旋转运动。


当没有阻尼*damping_relative 时,在一个可变形单元中绘制有效应力的时间历程曲线,会出现一个不随时间衰减的振荡,如下图。



image.png

单元有效应力的时间历程图(无阻尼)


当加入阻尼*damping_relative 时,振荡会迅速消失,如下图。


image.png    

单元有效应力的时间历程图(有阻尼)


【2】刚度阻尼:与刚度矩阵成比例,在显式分析中,刚度阻尼可以减小响应的高频振荡来提高模型的稳定性。


*damping_part_stiffness用于高频模式,COEF变量(作为正值输入)近似为部分临界阻尼。COEF的典型值为0.1,约为高频模式的临界阻尼的10%。


但是增加刚度阻尼会发生不稳定,则可以通过:

(1)减小TSSFAC。如果使用*damping_part_stiffness,则应相应减小*control_timestep中的TSSFAC,因为阻尼会影响临界时间步长。

(2)降低COEF到恢复稳定性。


【3】瑞利阻尼:同时包含质量和刚度阻尼项:

image.png    
  • 设置方法

① 使用 *DAMPING 关键字定义全局瑞利阻尼参数:


image.png    


 ② 通过过 *DAMPING_PART_MASS 和 *DAMPING_PART_STIFFNESS 对特定部件定义阻尼。


         

4. 局部阻尼 (Local Damping)

         
         
  • 作用通过人工阻尼力稳定特定区域的运动(例如接触界面。

  • 设置方法接触力会在反应中引入噪声,在这种情况下通过*contact卡2上的VDC参数添加黏性阻尼有助于降低噪声。VDC 输入为临界阻尼的百分比,一个典型的值在 10~20范围。


image.png      




       
           

5. 临界阻尼比 (Critical Damping Ratio)

           
           
  • 作用:在模态分析或频域响应中,定义系统阻尼相对于临界阻尼的比例。

  • 公式

image.png      

其中           ζ 为临界阻尼比,          ω 为固有频率。

  • 设置方法过 *DAMPING_FREQUENCY_RANGE 定义频率相关的阻尼比。



       
         
参数选择建议            
         
Choose Suggestions        
1、体积粘性阻尼:保持默认值(QB=1.5, QC=0.06)除非出现数值振荡。    

2、瑞利阻尼:选择主导频率范围,计算  和 

3、材料阻尼:参考实验数据或材料手册,临界阻尼比范围:1%-10%。





       
         
注意事项          
         
precautions        

1、阻尼设置需平衡能量耗散与计算精度,过度阻尼会导致非物理结果。

2、显式分析中,材料阻尼和结构阻尼通常对低频问题影响较小,但对高频噪声抑制有效。

3、对于瞬态冲击问题,优先使用体积粘性阻尼;对于振动问题,使用瑞利阻尼更合适。


通过合理设置这些阻尼模型,LS-DYNA能够有效地模拟各种动力学行为,提高模拟的准确性和稳定性。



         
END          
       
       


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

08.MPP 和 SMP 介绍

LS-DYNA 支持两种并行计算的模式,分别为SMP和MPP;不同的计算模式适用不同的对象,各有不同的优劣,跟着我了解它们的不同吧。 一、SMP 介绍SMP: Shared Memory Parallelism,共享式多核计算。1、原理:指所有的处理器(CPU共享内存),而每一个CPU内可能含有不同的Core数。因此当接到计算工作时,系统会根据各个CPU的负载平均分配工作量给每一个 CPU。SMP串行计算所有CPU共同计算一个模型的单元应力、应变、接触、边界条件、节点加速度等。SMP2、特点:SMP 只能在一个计算机上计算,且计算内存共享,可以调用的计算内核最多为16个CPU,超过16个CPU计算效能很难提升。3、 适用类型:隐式计算只能用SMP求解器。4、 计算效率:显式计算CAE模型越大SMP求解器的求解效率越低。5、 缺点:若是计算工作量太大,内存容易不足,造成系统可能无法执行或是执行速度缓慢。二、MPP 介绍MPP: Massively Parallel Processing,分布式多核计算。1、计算原理:分布式多核计算的每一个CPU 都有自己的Core数,而在接受计算工作以前,计算工作及资料被提前分割好,并可以指定要将那一步部分的计算交由指定的CPU去执行,最后再根据网路的设定,将每一部分的资料串联在一起,最终完成计算任务。MPP并行计算依据占用的CPU数量将一个模型划分成若干个区域,每个CPU"包产到户",在每一个计算周期内单独计算各自区域内的单元应力/应变、接触、边界条件、节点加速度等信息,然后与周围其他相关区域进行同步数据交换,从而使区域边界的信息统一。MPP2、特点:MPP可以在单台计算机或者集群计算机上计算,且计算内存相互独立,可以调用的CPU核数远远超过16核。MPP需要在系统中安装一个MPI或信息传递接口(例如,Platform MPI、Open MPI、Intel MPI等)。MPI软件和动态库需要安装在本地机器上,且LS-DYNA程序可以引用它们。MPP 计算时任务的目录是在本地计算节点还是管理节点都可以,如果集群有 NFS 网络和IB,那么系统就可以在NFS磁盘上运行。但是为了获得最佳性能,我们建议使用本地磁盘。3、适用类型:显式计算既可以使用SMP求解器,也可以使用MPP求解器。4、计算效率:显式计算CAE模型越大,使用的计算机CPU数量越多,则MPP求解器相对于SMP求解器的优势就越明显。因为MPP求解器是各CPU 对模型"分片包干"后并行计算的,CPU数量翻倍,计算效率翻倍。分块模式可以是求解器自动分块,也可以由用户自定义分块,用户自定义分块可以根据需要将CPU 集中分配给计算量大的区域(如高度非线性区域),这可以说是MPP求解器相对于SMP 求解器的一个独特优势。5、缺点:MPP 求解随着 CPU 数量的增多,不同的 CPU 数量、不同的分块方式,会导致计算结果产生差异。三、SMP 和 MPP 对比 四、MPP 控制的关键字(1) *CONTROL_MPP_IO_NODUMP:抑制所有转储文件和完整卡组重启文件的输出。没有输入参数。此关键字的存在导致跳过 d3dump 和 runrsf 文输出例程。它还抑制完整卡组重启文件 d3full 的输出。(2) *CONTROL_MPP_DECOMPOSITIN :定义 MPP 的分布方法。(3) *CONTROL_MPP_DECOMPOSITION_TRANSFORMATION:通过指定的方法修改分块方式。(4) *CONTROL_MPP_DECOMPOSITION_AUTOMATIC:使求解器应用一个简单的启发式方法来对模拟过程进行适当分布。(5) *CONTROL_MPP_DECOMPOSITION_CHECK_Speed:根据从处理器的相对计算速度修改模拟计算分布。(6) *CONTROL_MPP_DECOMPOSITION_CONTACT_DISTRIBUTE:保障接触分布在所有计算处理器上,还可以为大型接触带来更好的负载平衡。(7) *CONTROL_MPP_DECOMPOSITION_CONTACT_ISOLATE:使得接触在单个处理器上隔离,从而减少相互通信。(8) *CONTROL_MPP_DECOMPOSITION_FILE:允许定义预分布、或者后续运行,以及无定义分布运行。(9) *CONTROL_MPP_DECOMPOSITION_NUMPROC:定义分布的处理器数量做整车碰撞分析时,一般采用MPP并行计算,此时LS-DYNA默认会根据网格数量均匀切分模型,保证每个计算节点计算负载均衡。但是此种切分方式不能考虑材料复杂度,因此当被分配到采用高复杂度材料单元的节点计算效率将明细低于其他节点,拖慢计算效率。所以,可以将应用高复杂度材料(如Crachfem失效材料)的零件定义成一个Part_set,通过关键字*CONTROL_MPP_DECOMPOSITION_PARTSET_DISTRIBUTE在DYNA并行计算中平均分配至所有计算节点上,保证计算效率。 END 来源:CAE碰撞仿真指导

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