1)计算效率方面:发现其计算时间比单纯用SPH快了16倍(SPH算法耗时约16小时,而采用FEM-SPH耦合算法仿真耗时1小时13分钟左右),仿真结果文件D3PLOT文件平均十秒一个,仿真过程极度丝滑。 2)SPH边界难以处理问题:FEM-SPH耦合模型不用考虑用虚粒子约束法定义SPH边界,在耦合模型中直接将外层SPH粒子固连在FEM网格上即可。
不同于FE算法,本文中SPH模型以粒子代替网格,相当于有限网格的磨粒与粒子化工件的不连续加工过程。因此,FE的面面接触算法已经不再适用。本文对于有限元单元与光滑粒子接触界面的相互作用(磨粒与工件),则是通过罚函数算法来定义,耦合接触算法采用自动点面接触算法,主面设为磨粒,从面设置为SPH工件(MSTYP=3,SSTYP=4),...
(1)计算效率方面:发现其计算时间比单纯用SPH快了16倍(SPH算法耗时约16小时,而采用FEM-SPH耦合算法仿真耗时1小时13分钟左右),仿真结果文件D3PLOT文件平均十秒一个,仿真过程极度丝滑。 (2)SPH边界难以处理问题:FEM-SPH耦合模型不用考虑用虚粒子约束法定义SPH边界,在耦合模型中直接将外层SPH粒子固连在FEM网格上即可。
3、删除掉wb中接触,需要在prepost中定义。 4、划分网格如下图所示,(sph中的部分节点需要与FEM中的对齐) 4、设置边界条件和初始条件。设定FEM-SPH界面的segment,设置周围为固定边界。 5、用prepost打开K文件 6、将FEM网格变为SPH网格。 7、完成SPH网格的建立,删除掉元Part。 7、设置set_node,用于定义接触. 8、...
步骤一,当sph粒子位于sph-fem耦合界面时或者在耦合界面附近时,如图2其光滑函数w的支持域就会与sph的问题域相交,传统的sph方法无法适用于整个求解域。此时如果在求解域的边界处仍然采用传统的sph积分格式就不可避免地会出现误差,降低计算精度,因此在sph-fem固结算法中要消除耦合边界对sph积分的影响。
为解决爆炸与冲击荷载下有限元的大变形问题,FEM-SPH耦合法被广泛应用于数值模拟中。 本课程全面介绍FEM-SPH耦合算法的建模细节及工程应用。 主要内容 1. 复杂模型的SPH粒子生成方法 2. SPH-FEM耦合关键字设置 3. 弹丸侵彻靶板案例讲解 4. 爆炸与冲击...
FEM/SPH耦合算法的射弹冲击起爆仿真 张忠,陈卫东,陈浩 (哈尔滨工程大学航天与建筑工程学院,哈尔滨150001) 摘要:为了说明有限元与耦合算法的优越性,利用多种算法对高速射弹撞击屏蔽炸药起爆问题进 行了数值仿真,对各种算法的优劣性进行讨论及对比研究。利用升降法得到了不同算法的临界起爆速度,结果表明有限元与SPH耦合算...
respectively.Underwhentheparticlespeedaftercollisionprocessisbetween5~25m/s.Keywords:aluminumliquidgranulation;SPH;FEM;contactingalgorithm;numericalsimulation XefXY@ghc-ijklgm nbXe( o-pqne) restuv 1600℃ wx,yz 2000kJ/kg {|, }~{| ,Xe 68kgce/t z...
FEM-SPH耦合算法岩片采用改进的有限元-光滑粒子流体动力学耦合(FEM-SPH)的方法,建立了TBM单滚刀和双滚刀破岩的三维数值耦合模型,研究了TBM破岩机理,研究结果表明,压碎区形成阶段,在单滚刀推力的作用下,微裂纹向下扩展形成扇形破坏区;在切向力的作用下,微裂纹向前扩展形成另一个扇形破坏区,两个扇形破坏区组成了锥形...
基于FEM-SPH耦合算法的磨粒仿真研究.pdf,题目:基于FEM-SPH耦合算法的磨粒仿真研究 目录 1.仿真目的及意义 2.磨粒加工模型 3.磨粒加工模型简化 4.有限元算例实现基本流程 5.软件操作难点 6.仿真工件-工具参数设置 7.分析与讨论 8.总结与展望 1.仿真目的及意义 ⚫ 目的及意