OptFuture | 力学benchmark验证_静力学_拓扑优化_航空_航天_汽车_增材_材料

OptFuture | 力学benchmark验证

前言

本篇将以NAFEMS中静力学benchmark为准，验证OptFuture线性静力模块求解精度，涉及到的模型名称分别为Background to Benchmarks和Scordelis–Lo roof。

算例一

椭圆环在面外均布压力下的静力学分析

此benchmark中几何结构的俯视图为一个椭圆环，其上表面有均布压强。由于载荷和几何结构的对称性，仅分析椭圆环的四分之一，如图1。材料为各向同性材料，杨氏模量为210GPa，泊松比0.3。椭圆环上表面负Z方向受到大小为 0.5MPa 的均布压强载荷，在X、Y、Z的对称面处分别设置对称约束，椭圆环外表面施加位移约束，约束X和Y自由度，如图2所示。此算例将计算出模型位置D点，坐标（2, 0, 0.6）处的Y向正应力值并与参考资料1中给出的值进行比较。

图1 Background to Benchmarks简化模型（1/8）示意

图2 1/8模型约束示意（几何线框图）

总位移和Y向正应力计算结果如图3所示。

图3 位移与y向正应力

表1给出椭圆形孔内部坐标为 (2, 0, 0.6) 的点 D 处上表面的Y向正应力。采用粗化网格的情况下，该值与 NAFEMS 基准（参考资料1）基本一致，相对偏差小于 0.2%。本算例验证了OptFuture求解静力学问题的精度。

表1 OptFuture计算值与对标值比较（单位：MPa）

算例二

薄壁结构静力学分析

Scordelis–Lo roof

此算例是一个被广泛使用的基准模型，用于验证CAE软件在薄壁结构问题中的求解精度，屋面长度2L为50m，屋面半径R为25m，厚度为0.25m。外直面为自由面，外曲面受到约束，约束沿Y和Z方向平移，屋面顶部的直面和内曲面施加对称约束。在模型上表面的Z方向施加压强-90Pa。将计算出模型的最大Z方向变形并与参考资料2中给出的值进行比较。

图4 Scordelis–Lo roof模型示意

Scordelis–Lo roof模型Z向位移计算结果如图5所示，最大Z向位移解为0.311m。

图5 Z向位移云图

表2为OptFuture计算值与对标值的比较。参考资料2中引用的中面Z向位移的参考解为 -0.3086 m。实际上，在其他已发表的基准结果中该值为-0.302m，我们采用这个值作为该模型的基准解。此算例表明OptFuture在薄壁结构的求解中同样具有很高的计算精度。

表2 OptFuture计算值与对标值比较（单位：m）

结论

本次推送通过NAFEAMS中的两个经典静力算例对OptFuture静力模块计算准确性进行了评估。OptFuture线性静力计算结果显示，算例一中计算值与对标值的差异为-0.053MPa，相对偏差小于0.2%，算例二中计算值与对标值的差异为-0.009m，相对偏差小于3%，表明OptFuture计算值与NAFEAMS的对标值基本一致。

算例一验证了OptFuture在经典静力学问题中的应力求解精度，算例二验证了OptFuture求解薄壁结构问题的精度和效率。值得强调的是，OptFuture采用独创的空间网格技术，无需划分贴体网格；在薄壁结构中，也无需手动抽取中面。这大大简化了工程实际问题的求解流程，在汽车、航空航天等各工业部门中有广泛的应用前景。

参考资料

G.A.O. Davies, R.T. Fenner, and R.W. Lewis, Background to Benchmarks, NAFEMS, Glasgow, 1993.
R.H. MacNeal and R.L. Harder, Proposed Standard Set of Problems to Test Finite Element Accuracy, Finite Elements in Analysis and Design, 1, 1985.

往期推荐

OptFuture | 增材制造中的晶格结构

OptFuture | 拓扑优化制造约束（二）

OptFuture | 拓扑优化制造约束（一）

来源：OptFuture优解未来

优解案例 | 身轻如燕，基于性能驱动的机械臂优化设计

一、背景简介全方位移动机器人车（见图1，Mountain工作室），选用进口麦克纳姆轮、伺服电机、控制系统，具备无线遥控功能，广泛用于各行业的科研和生产，一台机器人可配合多台机床工作，可实现上下料、长距离搬运、喷涂等多种功能。图1 全方位移动机器人车优解未来是国内自主研发新一代云CAX智能设计工业软件的高科技公司，主营产品OptFuture软件基于云原生架构，可跨平台/地区实现团队协同设计，具有空间网格前处理器、高效CAE求解器、优化算法引擎和制造工艺等功能模块，支持多目标/多尺度/多材料/多工况等优化类型。拓扑优化/衍生式设计/创成式设计代表了事物设计方式的一种范式转变，并正在改变着事物的制造方式，OptFuture将助力设计师把任务需求更快更好的变成最终产品方案。二、机械臂优化设计上图1中机械臂的质量较大，会带来机械动作不灵活、伺服电机工作功率高等问题，故优解未来团队尝试利用拓扑优化技术对机械臂进行结构再设计，以探索满足性能指标条件下的最优轻质化构型。在实际工作环境中，机械臂具有多种典型工作状态（见图2），最大可搬运重量为500kg。通过多体动力学仿真，提取在每一种工作状态下机械臂关节处的载荷数据，以作为多工况分析的载荷边界条件。图2 机械臂典型工况利用机械臂的原有模型包络边界，可获得初始优化区域（见图3）。机械臂材质为铸铁，体分比0.10，底部关节处给定固定约束，端部施加载荷边界条件，进行基于多工况的机械臂拓扑优化设计（见图4）。图3 机械臂优化区域图4 机械臂优化设计三、分析与对比在最大搬运重量500kg下，机械臂的位移见图5，原有方案的最大位移为0.2mm，优化方案的最大位移为3.35mm（相比于臂长735mm，仍属于弹性小变形范畴内）。图5 位移云图机械臂的应力见图6，铸铁的屈服强度为370MPa，原有方案的最大应力为20.46MPa，剩余安全系数为18；优化方案的最大应力为200.17MPa，剩余安全系数为1.8。可见，优化方案既保留了一定的安全裕度，又较为合理的实现材料最佳分布。图6 Von-Mises应力云图两种机械臂方案的数据对比见图7，在满足一定强度安全裕度下，优化方案减重约75.6%，即115.5kg。产品制造成本较为接近，仍采用铸造+焊接的工艺方案。图7 数据对比四、最终产品效果经过拓扑优化设计和产品配色渲染，全方位移动机器人车的最终产品效果见图8，身轻如燕，如臂使指。图8 最终产品效果来源：OptFuture优解未来