轮胎股价大涨！高性能子午线轮胎技术改造之Abaqus轮胎建模仿真

• 2020年11月18日，赛轮轮胎年产330万套高性能智能化全钢载重子午线轮胎项目在赛轮轮胎沈阳工厂启动。 
• 2021年1月，赛轮轮胎曾公告，为进一步提升公司国际竞争力，公司全资子公司赛轮（越南）有限公司拟投资建设三期项目，具体为年产300万条半钢子午线轮胎、100万条全钢子午线轮胎及5万吨非公路轮胎项目。项目投资总额30.1亿元。 
• 2021年2月1日，赛轮轮胎又发布公告，拟投资赛轮（潍坊）高性能子午线轮胎技术升级改造项目，项目投资总额7.2亿元。山东财经报道注意到，技术升级改造项目建设完成后，赛轮(潍坊)将具备年产120万套高性能全钢载重子午线轮胎和600万套高性能半钢子午线轮胎的生产能力。

ABAQUS JOB=[JOB] USER＝ [.FOR]

Abaqus子程序的二次开发需要配置IVF和VS,网上有各种配置版本，在此不再赘述。本次课程所采用的配置为：ABAQUS2020+Visual Studio 2019+Intel Parallel Studio XE 2020。接下来，便开始讲述Abaqus 2020& Visual Studio 2019 & Intel Parallel Studio XE 2020配置

• 首先建立一空白文件命名为aba-vs-for.bat

• 打开后，输入Call+空格，在桌面右下角开始菜单栏找到右图所示，右键点击属性，对其目标处进行复制，如下所示：

• 将其复制到aba-vs-for.bat刚刚输入的后面:call %comspec% /k "D:\VS2019\VC\Auxiliary\Build\vcvarsamd64_x86.bat“
• 然后打开如下：右键点击属性，同样复制到aba-vs-for.bat文件

• 将其复制到aba-vs-for.bat刚刚输入的后面:call %comspec% /k "D:\VS2019\VC\Auxiliary\Build\vcvarsamd64_x86.bat“
• call C:\Windows\System32\cmd.exe /E:ON /V:ON /K ""D:\intel parallel studio xe 2020\compilers_and_libraries_2020.4.311\windows\bin\ipsxe-comp-vars.bat" intel64 vs2019“
• 最后输入pause

• 点击保存，以管理员身份运行aba-vs-for.bat

• 接下来关联ABAQUS，以下顺序按自己电脑存储位置进行：
• C:\SIMULIA\Commands\abq2020，将刚刚复制的两个目标复制其内即可。

• 以管理员身份运行abaus verification,查看verify.log文档：出现全PASS即成功：

• 在随意空白处cmd：abaqus information=system，可以看到Fortran配置成功

1. 程序中的变量名，不分大小写；
2. 变量名称是以字母开头再加上1到5位字母或数字构成，即变更名字串中只有前6位有效；
3. 一行只能写一个语句；
4. 程序的第一个语句固定为PROGRAM 程序名称字符串
5. 某行的第１个字符至第５个字符位为标号区，只能书写语句标号或空着或注释内容；
6. 某行的第１个字符为Ｃ或＊号时，则表示该行为注释行，其后面的内容为注释内容；
7. 某行的第６个字符位为非空格和非０字符时，则该行为上一行的续行，一个语句最多可有19个续行；
8. 某行的第７至７２字符位为语句区，语句区内可以任加空格以求美观；

9. 某行的第７３至８０字符位为注释区，８０字符位以后不能有内容。

1. GETVRN:获得节点信息,如节点温度、坐标、位移等;

2. GETVRMAVGAINODE:获得节点处平均材料积分点的信息,如应力等;

3. GETNODETOELEMMCONN:获得指定节点的单元列表;

4. GETPARTINF0 :获得节点处的部件信息。

• nStreamlines：施加磨损的参考部分（完整模型）上的节点数。

• nGenElem：模型中沿流线的扇区或元素划分的数量。

• nRevOffset：作为旋转对称模型定义的一部分而指定的节点偏移量。

• nReflOffset：作为反射对称模型的定义的一部分而指定的节点偏移量。（如果未反映模型，则将此参数设置为零）。

• jslnodes：包含所有可能在参考部分磨损的节点的必需节点信息的数组。这是一个大小为（2，nStreamlines）的数组 。对于每个流线，第一个组件是“根节点”的节点号（以下讨论中的节点a），它是参考部分上该特定流线上的节点。第二个组件是提供磨损方向的节点（节点b在下面的讨论中）。仅胎面拐角节点才需要第二部分。将其设置为等于定义磨损方向的参考部分上的节点编号。对于不在胎面拐角处的节点，将第二个数组分量设置为零。对于远离胎面拐角的这些节点，磨损将与局部3方向相反。

*HEADING*restart,read*STEP,INC=300,NLGEOM=YES,UNSYMM=YES 2:  ABLATION*STEADY STATE TRANSPORT,INERTIA=YES 3.6E5, 3.6E6, , 3.6E5*RESTART,WRITE,FREQ=1*print,contact=yes,adapt=yes*adaptive mesh,elset=TREAD,FREQ=1,MESH=4*adaptive mesh constraint,type=velocity,userNADAPT,1,,5.555555E-14*adaptive mesh constraint,constraint=lagrangianNADAPT_LAGR,*END STEP

以上为子午线轮胎磨损分析的所有内容，

注意：磨损分析前计算inp文件以及UMESHMOTION子程序FO件、注释版FO件在附件中。