sph算法简介

A. autodyn并行计算有什么作用

Autodyn使用起来比较方便，与Ls-dyna相比，其主要的优缺点为：
（1）AUTODYN相对于LS-DYNA，使用相对简单，学习起来比较容易。
（2）AUTODYN自带材料库，材料参数比较齐全。
（3）AUTODYN二维不支持多cpu并行运算，无法发挥多核处理器的效能，LS-DYNA不存在该问题。
（4）LS-DYNA适合于做碰撞，但是作爆炸，效果不如AUTODYN理想。如果要做接触爆炸和近场爆炸，那么就选AUTODYN，如果想做远场，那就选ABAQUS，精度比DYNA要好。
（5）LS-DYNA以LAGRANGE算法为主，兼有和Euler和ALE算法，而AUTODYN则不同，该程序擅长计算多物质
场，其LAGRANGE算法和Euler算法的功能都很强。
（6）AUTODYN在linu系统下没有图形用户界面，但可以求解，这对随时查看计算结果有影响；LS-DYNA不
存在该问题。（ANSYS13.0开始AUTODYN可以在linux下运行）
（7）AUTODYN并行计算设置较烦，也可能是为了更有效地利用cpu的并行计算能力，而LS-DYNA适合并行计
算设置。
(8) AUTODYN多核处理器支持不理想。在windows平台，四核处理器，运算单个任务CPU利用率可达25%，而
运行3个任务时往往40%都不到。
（9）Autodyn的sph算法只能为1个part，而Ls-dyna不存在该问题。
（10）Autodyna比较贵，Ls-dyna相对便宜，Autodyn国内只有安世亚太公司在卖，Ls-dyna有不少代理公司。

B. 浅水波波动可以用什么软件模拟ansys是否可以

浅水波动现在SPH算法还是挺好的 也有用Lattice 玻尔兹曼算法也不错 有个软件RealFlow你可以参考。
Ansys主要是结构静力分析，基本不用与流动分析。Fluent到还是可以考虑，这方面的论文也很多。

C. 在欧拉法中，加速度可以分别为当地加速度和时变加速度对吗

不对，当地加速度就是时变加速度。正确的说法应该是欧拉法中加速度可以分解为时变加速度（又名当地加速度）和位变加速度（又名迁移加速度）。

（1）时变加速度（当地加速度）（localacceleration）——流动过程中流体由于速度随时间变化而引起的加速度；

（2）位变加速度（迁移加速度）（connectiveacceleration）——流动过程中流体由于速度随位置变化而引起的加速度。

(3)sph算法简介扩展阅读：

欧拉法一般用于工程分析，代表有很多，比如CFX，fluent等。

关于说的“欧拉法计算加速度的推导过程是如何保证dt内的dv是同一质点的”的疑问说明你没有明白欧拉的观点，实际上当我们已知流体中某一点的状态的时候，那么其附近质点的状态也是与这一点关联的（关联量基本的就是流体的压力分布、速度分布、密度分布等）。

所以根据这个原理，划分出许多的网格后根据流体边界条件就可以从边界网格开始逐渐计算流体内部的流动状态，直到所有网格点的数据变化小于一定的值的时候（工程上称之为解的收敛），那么我们得到的这所有网格点的状态起来就是流体的近似流动状态。

拉格朗日法现在用于电影、动画和游戏特效制作，使用的有Realflow。其中Realflow就是用粒子计算的，计算方法是SPH算法，这个算法也是近似求解NS方程的。

D. 计算机图形学中的流体模拟，国内有这个课吗

基础CFD的课倒是不少，但是流体仿真方面的公开课真没有。
现在计算机图形学的流体仿真效果比较好的是SPH算法的，入门可以读一下muller那篇经典的论文《Particle-based Fluid Simulation for Inter active Applications》。
另外Robert Bridson教授在sigraph上做的<<fluid simulation for computer graphics>>讲义也已经出版了，讲的比较系统，可以读一下。

E. LS-DYNA的分析能力

¨ 非线性动力学分析
¨ 多刚体动力学分析
¨ 准静态分析（钣金成型等）
¨ 热分析
¨ 结构-热耦合分析
¨ 流体分析：
欧拉方式
任意拉格郎日-欧拉（ALE）
流体-结构相互作用
不可压缩流体CFD分析
¨ 有限元-多刚体动力学耦合分析 （MADYMO，CAL3D）
¨ 水下冲击
¨ 失效分析
¨ 裂纹扩展分析
¨ 实时声场分析
¨ 设计优化
¨ 隐式回弹
¨ 多物理场耦合分析
¨ 自适应网格重划
¨ 并行处理（SMP和MPP）
2.材料模式库(140多种)
¨ 金属
¨ 塑料
¨ 玻璃
¨ 泡沫
¨ 编制品
¨ 橡胶（人造橡胶）
¨ 蜂窝材料
¨ 复合材料
¨ 混凝土和土壤
¨ 炸药
¨ 推进剂
¨ 粘性流体
¨ 用户自定义材料
3.单元库
¨ 体单元
¨ 薄/厚壳单元
¨ 梁单元
¨ 焊接单元
¨ 离散单元
¨ 束和索单元
¨ 安全带单元
¨ 节点质量单元
¨ SPH单元
4.接触方式(50多种)
¨ 柔体对柔体接触
¨ 柔体对刚体接触
¨ 刚体对刚体接触
¨ 边-边接触
¨ 侵蚀接触
¨ 充气模型
¨ 约束面
¨ 刚墙面
¨ 拉延筋
5.汽车行业的专门功能
¨ 安全带
¨ 滑环
¨ 预紧器
¨ 牵引器
¨ 传感器
¨ 加速计
¨ 气囊
¨ 混合III型假人模型
6.初始条件、载荷和约束功能
¨ 初始速度、初应力、初应变、初始动量（模拟脉冲载荷）；
¨ 高能炸药起爆；
¨ 节点载荷、压力载荷、体力载荷、热载荷、重力载荷；
¨ 循环约束、对称约束（带失效）、无反射边界；
¨ 给定节点运动（速度、加速度或位移）、节点约束；
¨ 铆接、焊接（点焊、对焊、角焊）；
¨ 二个刚性体之间的连接－球形连接、旋转连接、柱形连接、平面连接、万向连接、平移连接；
¨ 位移/转动之间的线性约束、壳单元边与固体单元之间的固连；
¨ 带失效的节点固连。
7.自适应网格剖分功能
自动剖分网格技术通常用于薄板冲压变形模拟、薄壁结构受压屈曲、三维锻压问题等大变形情况，使弯曲变形严重的区域皱纹更加清晰准确。
对于三维锻压问题，LS-DYNA主要有两种方法：自适应网格剖分和任意拉格朗日-欧拉网格（ALE）网格进行Rezoning），三维自适应网格剖分采用的是四面体单元。
8. ALE和Euler列式
ALE列式和Euler列式可以克服单元严重畸变引起的数值计算困难，并实现流体-固体耦合的动态分析。在LS-DYNA程序中ALE和Euler列式有以下功能：
? 多物质的Euler单元，可达20种材料；
? 若干种Smoothing算法选项；
? 一阶和二阶精度的输运算法；
? 空白材料；
? Euler边界条件：滑动或附着条件；
? 声学压力算法；
? 与Lagrange列式的薄壳单元、实体单元和梁单元的自动耦合。
9.SPH算法
SPH（Smoothed Particle Hydrodynamics）光顺质点流体动力算法是一种无网格Lagrange算法，最早用于模拟天体物理问题，后来发现解决其它物理问题也是非常有用的工具，如连续体结构的解体、碎裂、固体的层裂、脆性断裂等。SPH算法可以解决许多常用算法解决不了的问题，是一种非常简单方便的解决动力学问题的研究方法。由于它是无网格的，它可以用于研究很大的不规则结构。
SPH算法适用于超高速碰撞、靶板贯穿等过程的计算模拟，下图是泰勒杆冲击试验模拟。
10.边界元法
LS-DYNA程序采用边界元法BEM（Boundary Element Method）求解流体绕刚体或变形体的稳态或瞬态流动，该算法限于非粘性和不可压缩的附着流动。
11.隐式求解
用于非线性结构静动力分析，包括结构固有频率和振型计算。LS-DYNA中可以交替使用隐式求解和显式求解，进行薄板冲压成型的回弹计算、结构动力分析之前施加预应力等。
12.热分析
LS-DYNA程序有二维和三维热分析模块，可以独立运算，也可以与结构分析耦合，可进行稳态热分析，也可进行瞬态热分析，用于非线性热传导、静电场分析和渗流计算。
热传导单元：8节点六面体单元（3D），4节点四边形单元（2D）；
材料类型：各向同性、正交异性热传导材料，可以与温度相关，以及各向同性热传导材料的相变；
边界条件：给定热流flux边界，对流convection边界，辐射radiation边界，以及给定温度边界，它们可随时间变化；给定初始温度，可计算二个物体接触界面的热传导和热辐射，给定材料内部热生成（给定热源）；
热分析采用隐式求解方法，过程控制有：
? 稳态分析还是瞬态分析；
? 线性问题还是非线性问题；
? 时间积分法：Crank-Nicholson法（a=0.5）和向后差分法（ a=1）；
? 求解器：直接法或迭代法；
? 自动时步长控制。
13.不可压缩流场分析
LS-DYNA不可压缩流求解器是960版新增加的功能，用于模拟分析瞬态、不可压、粘性流体动力学现象。求解器中采用了超级计算机的算法结构，在确保有限元算法优点的同时计算性能得到大幅度提高，从而在广泛的流体力学领域具有很强的适用性。
14.多功能控制选项
多种控制选项和用户子程序使得用户在定义和分析问题时有很大的灵活性。
输入文件可分成多个子文件；
用户自定义子程序；
二维问题可以人工控制交互式或自动重分网格（REZONE）；
重启动；
数据库输出控制；
交互式实时图形显示；
开关控制-可监视计算过程的状态；
对32位计算机可进行双精度分析。
15.前后处理功能
LS-DYNA利用ANSYS、LS-INGRID、ETA/FEMB及LS-POST强大的前后处理模块，具有多种自动网格划分选择，并可与大多数的CAD/CAE软件集成并有接口。
后处理：结果的彩色等值线显示、梯度显示、矢量显示、等值面、粒子流迹显示、立体切片、透明及半透明显示；变形显示及各种动画显示；图形的PS、TIFF及HPGL格式输出与转换等。
16.支持的硬件平台
LS-DYNA 960版的SMP版本和MPP版本是同时发行的。MPP版本使一项任务可同时在多台分布计算机上进行计算，从而最大限度地利用已有计算设备，大幅度减少计算时间。计算效率随计算机数目增多而显着提高。
LS-DYNA 960版的SMP版本和MPP版本可以在PC机（NT、LINUX环境）、UNIX工作站、超级计算机上运行。

F. 用欧拉法计算出来的加速度怎么理解

欧拉法一般用于工程分析，代表有很多，比如CFX，fluent等。
关于你说的“欧拉法计算加速度的推导过程是如何保证dt内的dv是同一质点的”的疑问说明你没有明白欧拉的观点，实际上当我们已知流体中某一点的状态的时候，那么其附近质点的状态也是与这一点关联的（关联量基本的就是流体的压力分布、速度分布、密度分布等），所以根据这个原理，我们划分出许多的网格后根据流体边界条件就可以从边界网格开始逐渐计算流体内部的流动状态，直到所有网格点的数据变化小于一定的值的时候（工程上称之为解的收敛），那么我们得到的这所有网格点的状态起来就是流体的近似流动状态。这里并不是去求解某个质点在流体中的速度和位置，所以你的问法表现出你还没有理解欧拉法的计算思想。
拉格朗日法现在用于电影、动画和游戏特效制作，使用的有Realflow。其中Realflow就是用粒子计算的，计算方法是SPH算法，这个算法也是近似求解NS方程的。具体的可以问度娘。

G. 导入autodyn的图形应该是什么格式

刚开始接触autodyn和truegird，为了将tg的模型导入到autodyn看了两天的文档，查了好多帖子。
将收获与大家分享一下！！
《TrueGrid®User’s Guide For AUTODYN》一文中的例子如下

第一句 autodyn 要写了。不然就不知输出的格式。
很重要的就是后几句。
只有在merge环境中才能使用write命令。

在tg目录下出现一个trugrdo文件。如下

将此文件名的后缀改为.zon