nCode:DesignLife 案例十一——使用Python定义一个SN方法
在本案例中,我们将使用Python脚本语言来:
- 定义包括平均应力校正方法(mean stress correction)的非标准SN曲线公式(例如在航空航天应用中可能使用的公式);
- 自定义DesignLife以使用新公式;
- 然后将其应用于飞机起落架的一部分分析中;
本案例中概述的相同过程也可用于将自定义SN方法添加到GlyphWorks中的Stress Life图形符号中。
有关Python的安装指南,请参见Glyph Reference Guide > GlyphBuilder Glyphs > 109 Scripting Glyph > 110.2 Python and NumPy。
我们不建议剪切和粘贴这些文档中的代码。这是因为PDF格式可能会更改代码示例的格式,例如缩进和换行,或者在复制的文本中引入隐藏的无效字符。
11.1 案例文件
此案例的文件位于nCode安装文件夹中的\demo\designlife\11_CustomSNPython下。需要的文件是:
- AeroSN.npd
- AeroSN.py
- CustomSNMethods.sys
- LandingGear.flo
- loadcases.lcs
- lower.op2
- mxd_ncode_Aero_SN.xdef
- spectrum.spe
- stress_history.s3t
将上述文件复制到工作文件夹中,并使用副本文件完成该实例。
11.2 使用Python
要运行此示例,应在运行DesignLife的同一台计算机上安装Python。有关安装Python的详细信息请参考《Python For nCode》。
本案例假设你不了解使用Python,但如果你有使用其他语言的编程经验,则会有所帮助。本示例旨在作为一个入门,介绍如何使用Python脚本添加新的疲劳分析方法。这不是一个关于Python的教程;相反,我们将关注Python如何与DesignLife交互并扩展它。
DesignLife的标准安装提供了七种定义SN材料数据的不同方法:
- Standard——在对数应力与对数寿命图上定义线性线段的参数;
- MultiMeanCurve——不同平均应力值下的数字化SN曲线;
- MultiRRatioCurve——不同负载比下的数字化SN曲线;
- Haigh——数字化的恒定寿命曲线;
- Bastenaire——使用Bastenaire公式定义SN曲线的参数;
- MultiTemperatureCurve——不同温度下的数字化SN曲线;
- Custom——用户定义的S-N方法。该选项允许用户使用Python脚本和材料数据的自定义定义来自定义软件。
此外,有不同的方法将平均应力效应纳入计算:
- None——不校正平均应力;
- Goodman——用Goodman方法校正平均应力;
- GoodmanTensionOnly——使用Goodman方法校正平均应力仅适用于拉应力;
- Gerber——用Gerber方法修正平均应力;
- GerberTensionOnly - 用Gerber方法修正平均应力仅适用于拉应力;
- Interpplate——在不同SN或恒定寿命曲线之间进行插值;
- FKM——FKM指南中描述的方法;
所有这些方法在DesignLife Theory Guide指南中都有详细描述。
尽管这些选项将涵盖大多数要求,但一些用户可能希望使用不同于其中任何一种的SN公式和/或平均应力校正方法。这可能是因为:
- 他们的公司已经建立或拥有使用不同于所提供的SN公式的方法;
- 他们需要使用手册数据,这些数据以不同于任何常用设计寿命方法的方式进行参数化;
- 他们使用的材料具有不同的建模方式
因此,DesignLife提供了定义自定义SN方法的能力,这样用户可以使用自己的材料数据和SN公式,同时仍然可以从DesignLife的其他功能中获益。使用自定义SN方法的能力已在低库级别的软件中实现。由于这些库是在DesignLife和GlyphWorks之间共享的,因此新的自定义SN方法(如下所述)将在两种产品中都可用。
自定义SN方法的语法和属性在《DesignLife User Guide》(第1.13节自定义方法)中有详细描述。
每个自定义SN方法有效地替换了代码中基于每个循环的范围和平均应力值(包括任何平均应力或其他影响)计算每个循环疲劳损伤的部分。一旦实现,自定义SN方法可以像通常一样使用,例如,自定义SN方法可以与随机PSD加载一起使用。注意,在使用自定义SN方法时可能会有性能损失。
11.3 自定义DesignLife
我们现在将自定义DesignLife以使用新的SN曲线公式,包括四个平均应力校正选项。
本节中描述的步骤向你展示了如何自定义DesignLife以添加新的SN方法。完成后,新方法作为CustomSNMethod属性的选项可用。在DesignLife中,当SN Analysis engine上的SNMethod = Custom时,CustomSNMethod处于激活状态。
在GlyphWorks中,通过设置属性CalculationMethod=Custom,然后从CustomSNMethod选择列表中选择所需的方法,可以在Stress Life fatigue图形符号中使用相同的方法。
这些步骤可以总结为:
-
修改材料数据库:
● 创建新的材料格式
● (可选)添加新的绘图定义(new plot definition,简称npd)
● 创建所需的数据库项 -
添加一个新的自定义方法:
通过创建CustomSNMethods.sys系统文件。
● 通过创建Python脚本(.py)文件来编写新方法
11.3.1 修改材料数据库
SN曲线方程:
其中S∞,A,B和C是材料参数。
- 要添加在MXD材料数据库中存储适当数据的功能,我们必须创建包含以下文本的xdef文件。我们把这个文件叫做mxd_ncode_Aero_SN.xdef。它必须放在DesignLife/GlyphWorks安装目录下的mats文件夹中(类似于\ncode2019\GlyphWorks\mats),因此我们将需要对此文件夹的写入权限。xdef文件的名称并不重要,只要它的前缀是“mxd_ncode_”,并且不覆盖现有的xdef文件,但是文件第二行的id很重要,这定义了方法的名称。
- 使用编辑器创建包含以下内容的XML格式文件,或者将提供的文件放在正确的位置。下次打开材料数据库时,可用的数据类型将自动更新以包含新格式。
<?xml version="1.0" encoding="UTF-8"?>
<object name="Custom Aero SN data set" id="Aero_SN" version="1">
<definition>
<child_types/>
</definition>
<common>
<property name="MaterialType" label="Material Type" type="9" help="The material
type - an integer that refers to a list of known types available in the help documentation.">
<integer default_value="99"/>
</property>
<property name="UTS" label="Ultimate Tensile Strength" type="1">
<real min="100" max="4000" default_value="500.0" default_unit="MPA"/>