ABAQUS+ALE自适应网格关键技术.docVIP

ABAQUS ALE自适应网格技术 为了方便了解，先整体介绍一下ALE网格自适应方法基础过程，一个完整ALE过程能够分为若干个网格remesh子过程，而每一次remesh过程能够分为两步： 1生成一个新网格（create a new mesh），利用多种算法和控制策略生成一个良好网格，关键包含划分频率和算法。 2环境变量转换（advection variales），也就是将旧网格中变量信息利用remapping技术转换到新网格中，也有不一样算法，其中包含静变量（应力场，应变场等）转换和动变量（速度场，加速度场等）转换。 上面两步在软件设置上面，可认为是对网格划分区域控制（ALE Adaptive Mesh Domain）和算法控制（ALE Adaptive Mesh Controls）。 1 ALE区域控制 （1）几何区域选择（set） ※ No ALE adaptive mesh domain for this step 该分析步没有使用ALE技术。 ※ Use the ALE adaptive mesh domain below 将以下区域定义为ALE区域。 （2）ALE Adaptive Mesh Controls 自适应技术控制选项，后面介绍 （3）Frequency 频率控制，关键是对整个step time中网格remesh次数进行控制。Remesh次数n能够由n=Increment number /Frequency来表示其意义，当frequency值为i时，表示每i个增量步进行一次remesh。 一个经典ALE过程，在每5-100个增量步就需要一次remesh，对于拉格朗日问题，改参数默认值为10，若变形实在太大，可合适调高，以增加网格重画强度，对于爆炸，碰撞等变形时间极短问题求解，则在每一个增量步全部需要一次remesh，这时Frequency值需要设置得很小，比如设为1，当然，adaptive remesh过程强度也很高，也会很废时。对于其它变形不是很猛烈问题求解，该参数值能够合适调高。对于欧拉问题，默认值为1。 图1 Frequency设置 （4）Remeshing sweeps per increment 一个频率下迭代次数，当该参数值为n时，每一个remesh过程将对网格进行n次sweep，其实这个参数能够了解为对整个adaptive remesh过程每一个子过程（remesh过程）强度进行控制。 那么，我们先来了解一下sweep概念，每sweep一次，abaqus将利用我们设置好算法（体积算法，拉普拉斯算法或等位算法）生成一套新网格，但这个网格不一定是符合要求，所以，需要在生成新网格基础上用一样方法再进行sweep，就像我们求解方程时迭代概念是一样。就这么一直sweep下去直到sweep次数达成mesh sweeps参数值，这么就完成了一个remesh过程中新网格生成。一样，mesh sweeps参数值越高，adaptive remesh过程强度越高，网格优化情况良好机率也就越大。 图2 Remeshing sweeps per increment设置 （5）Initial remeshing sweeps 也就是ALE过程开始之前对网格一个优化，概念和mesh sweeps类似，因为我们有可能利用已经变形很厉害网格进行分析，这时，在分析开始之前，就需要对网格进行重画。 图3 Initial remeshing sweeps设置 2 ALE过程控制 算法控制包含两部分；一为网格算法控制；其二为变量转换算法控制。 （1）Priority 也就是指网格梯度控制（是否保持初始网格梯度，若需要保持初始网格梯度，则对网格质量将会有影响）。 Improve aspect ratio 在计算过程中将考虑到网格单元高宽比改善，不考虑对初始网格梯度保持。 Preserve initial mesh grading 在计算过程中确保初始网格梯度，但不会考虑到网格宽高比改善。 图4 Priority设置界面 （2）smoothing algorithm ※ Use enhanced algorithm based on evolving element geometry 关键是在几何学方面对我们定义网格sweep算法（前面提到三种算法）进行增强，目标是为了确保adaptive remesh过程健壮性，为推荐选项，选它就行了 ※ conventional smoothing 利用我们定义好算法进行计算，无几何增强。 图5 smoothing algorithm界面控制 （3）Meshing Predictor 也就是网格节点位置控制（理想网格节点位置控制，将会降低需要网格sweeps次数，降低资源浪费）。 Cur