子模型是得到模型部分区域中更加精确解的有限单元技术。在有限元分析中往往出现这种情况，即对于用户关心的区域，如应力集中区域，网格太疏不能得到满意的结果，而对于这些区域之外的部分，网格密度已经足够了。
要得到这些区域的较精确的解，可以采取两种办法：
- 用较细的网格重新划分并分析整个模型；
- 只在关心的区域细化网格并对其分析。
显而易见，方法a太耗费机时，方法b即为子模型技术。
子模型方法又称为切割边界位移法或特定边界位移法。切割边界就是子模型从整个较粗糙的模型分割开的边界。整体模型切割边界的计算位移值即为子模型的边界条件。
子模型基于圣维南原理，即如果实际分布载荷被等效载荷代替以后，应力和应变只在载荷施加的位置附近有改变。这说明只有在载荷集中位置才有应力集中效应，如果子模型的位置远离应力集中位置，则子模型内就可以得到较精确的结果。
除了能求得模型某部分的精确解以外，子模型技术还有几个优点：
- 它减少甚至取消了有限元实体模型中所需的复杂的传递区域。
- 它使得用户可以在感兴趣的区域就不同的设计(如不同的圆角半径)进行分析。
- 它帮助用户证明网格划分是否足够细。
使用子模型的限制如下：
- 只对体单元和壳单元有效。
- 子模型的原理要求切割边界应远离应力集中区域。
生成并分析较粗糙的模型。
第一个步骤是对整体建模并分析。分析类型可以是静态或瞬态的，其操作与各分析的步骤相同。单元类型——子模型技术只能使用块单元和壳单元。分析模型中可以有其他单元类型(如梁单元作为加强筋)，但切割边界只能经过块和壳单元。
生成子模型。
子模型是完全依靠粗糙模型的。使用与粗糙模型中同样的单元类型。同时应指定相同的单元实参(如壳厚)和材料特性。
- 提供切割边界插值。
用户定义切割边界的结点，ANSYS程序用粗糙模 型结果插值方法计算这些点上的自由度数值(位移等)。对于子模型切割边界上的所有结点，程序用粗糙模型网格中相应的单元确定自由度数值，然后这些数值用单元形状功能插值到切割边界上。
- 分析子模型。
用户指定分析类型和分析选项，加入插值的DOF数值(和温度数值)，施加其他的载荷和边界条件，指定载荷步选项，并对子模型求解
- 验证切割边界和应力集中区域的距离应足够远。最后一步是验证子模型切割边界是否远离应力集中部分。可以通过比较切割边界上的结果(应力，磁通密度等)与粗糙模型相应位置的结果是否一致来验证。如果结果符合得很好，证明切割边界的选取是正确的。如果不符合的话，就要重新定义离感兴趣部分更远一些的切割边界重新生成和计算子模型。
子模型切分边界的方法与注意事项
WB可以较容易的实现子模型的边界切分，其分析流程如下图。
需要注意的是B3(全模型的几何)和C3(子模型的几何)必须基于统一的坐标系。
网格无关性的评价标准
数值计算中之所以需要网格，是由所采取的算法密切相关的。当前的主流偏微分方程数值离散方法都是先计算节点上的物理量，然后通过插值在方式求得节点间的值。因此，从理论上讲，网格点布置得越密集，所得到的计算结果也越精确。
网格越密，计算量越大，计算周期也越长。而我们的计算资源总是有限的。其次，随着网格的加密，计算机浮点运算造成的舍入误差也会增大。因此在实际应用中，使用者总是在计算精度与计算开销间寻求一个比较合适的点，这个点所处的位置就是达到网格无关的阈值。实际上，大部分的情况下，网格加密到一定程度后，计算结果的变化已经相当小了。
所谓网格无关性验证，实际上就是验证计算结果对于网格密度变化的敏感性。也就是不断的改变网格的疏密，观察计算结果的变化，若其变化幅度在允许的范围之内，我们就可以说计算值已经与网格无关了。对于瞬态分析，需验证时间步长的独立性。
来都来了，点个赞再走吧~~~
扫一扫
被折叠的 条评论
为什么被折叠?
到【灌水乐园】发言
