TexGen导出ABAQUS的inp文件没有材料？三步补齐材料定义实战指南

1. 先说结论：材料属性为什么会在导出时"消失"

做织物流变学分析、复合材料细观力学建模的朋友，十有八九在TexGen和ABAQUS之间倒腾过inp文件。模型建得漂漂亮亮，纱线路径、横截面、层叠结构一个不落，结果导出inp之后一提交计算，要么报错，要么算完应力全是为零的荒谬结果。打开inp文件一查，发现里面连一个*Material关键字都搜不到，整个人直接愣在原地。

这个问题的根源其实不在"导出"这个动作本身，而在于我们对TexGen的定位理解有一点偏差。TexGen从骨子里来说是一个几何建模工具，它的核心数据模型里存的是纱线路径、截面形状、织物层叠关系这些几何信息。虽然TexGen的界面里也有材料属性相关的面板，但那个更多是为了一些内置的有限元分析流程服务的——包括周期性边界条件施加、织物体积分数计算这些模块。当你选择"导出ABAQUS inp"这个功能时，TexGen默认做的是把当前织物几何转化为网格节点和单元，同时根据纱线名称生成一系列Element Set（单元集合），但它并不会负责任地把材料本构、截面属性替你写进去。

我日常用到的TexGen版本是4.4.x，导出的inp文件结构大致长这样：开头是Node定义节点坐标，然后是Element定义单元连接关系，接着是Elset或Solid Section的雏形，偶尔会带Orientation。但请注意，这里的Solid Section通常是残缺的，它可能只定义了单元集合和截面类型，但Material名称是空的，或者干脆整段缺失。如果你用ABAQUS/CAE直接打开这个inp，会看到每个纱线对应的Set下面没有关联任何Section。你提交Job的时候ABAQUS就会报出类似"The following elements have no section assigned"之类的信息。就算你强行让它算，没有材料本构的单元在结果里也是一堆垃圾数据。

修复思路其实很简单：既然TexGen负责几何，那材料本构就由用户在inp里手动补齐，或者通过ABAQUS/CAE的界面来补。听起来是两步走，实际上掌握方法之后，5分钟真的够用。这篇就围绕这个核心思路展开，把修复流程、避坑经验和常见报错一次性讲清楚。

2. 三步补齐材料定义：不需要重建模

2.1 第一步：看懂TexGen导出inp的文件结构

动手修之前，先把TexGen导出的inp文件拆开看一遍。不同TexGen版本的导出格式略有差异，但核心结构是一致的。这里以我手头一个平纹织物2×2的代表性格子模型为例，导出的inp片段大致长这样：

*Node 1, 0.0, 0.0, 0.0 2, 0.25, 0.0, 0.0 …… *Element, type=C3D8R 1, 1, 2, 4, 3, 5, 6, 8, 7 …… *Elset, elset=Yarn_1 1, 2, 3, 4, 5, 6, 7, 8 …… *Elset, elset=Yarn_2 9, 10, 11, 12, 13, 14, 15, 16 ……

注意几个关键信息：节点坐标单位通常是毫米（TexGen默认mm制）；单元类型默认是C3D8R减缩积分六面体单元；每个纱线对应一个*Elset，名字就是你在TexGen里给纱线定义的名称前缀加编号。如果你是给整个织物基体单独建模的话，TexGen也会生成一个Matrix或Resin的实体集合。

整个inp文件中，和材料相关的部分——*Material、*Solid Section、*Density、*Elastic——通常一个都没有。所以你拿到的其实是一个只有几何外壳的"半成品"文件。把这个文件的逻辑捋清楚之后，修复就有明确的抓手了：你只需要在合适的位置插入两段话——材料本构定义和截面赋值定义。

有一点提醒：如果你在TexGen里勾选了"Export yarn orientations"之类的选项，inp里会带上*Orientation定义。这个对后面的材料方向校准很有用，建议导出时保留。还有一个小习惯：导出之前把每个纱线的名称改得规范一点。默认的Yarn1、Yarn2倒是也能用，但如果你在ABAQUS里要做后处理，看到"Yarn_1_weft_90deg"这样的名字，提取结果时会舒服得多。

2.2 第二步：准备一份通用材料定义模板

补材料本质上就是往inp里写标准的ABAQUS关键字。我的建议是准备一份通用的材料参数模板，放在手边，随用随抄。这里用单向玻璃纤维纱线（E-glass / Epoxy）作为示例，材料参数是复合材料手册里常见的典型值。

首先要有密度：

*Material, name=GlassEpoxy *Density 1.9e-09,

注意单位。TexGen导出节点坐标是毫米，时间单位是秒，那么对应的密度单位是吨每立方毫米，应力单位是兆帕。很多人在这里栽跟头：用kg/m³的密度直接写进去，结果模态分析和屈曲分析的频率差了三个数量级。1.9e-09这个量级对应的就是把1.5~2.1 g/cm³（也就是1500~2100 kg/m³）的树脂基复合材料换算成吨/立方毫米。默认的密度写不写其实不影响纯静力分析，但对动力学、模态、冲击分析至关重要，建议规规矩矩写上。

然后是弹性本构。单向纱线是典型的横观各向同性材料，沿着纤维方向的模量远大于横向。如果只用各向同性弹性的两个参数，算出来的应力分布基本没有参考价值。ABAQUS里定义横观各向同性弹性有两种方式：一种是Elastic, type=ENGINEERING CONSTANTS，直接填入工程常数；另一种是Elastic, type=ANISO，需要给完整刚度矩阵。对大多数织物细观模型来说，工程常数法最直观，也不容易填错。模板如下：

*Elastic, type=ENGINEERING CONSTANTS 45000., 9500., 9500., 0.3, 0.3, 0.4, 5200., 3800., 3800.

这一行九个数字依次对应E1、E2、E3、Nu12、Nu13、Nu23、G12、G13、G23。E1是纤维方向的杨氏模量，E2、E3是垂直于纤维方向的模量。注意对于横观各向同性材料，E2和E3理论上是相等的，Nu12和Nu13相等，G12和G13相等。很多新手会把九个参数抄得七零八落，最后ABAQUS提示材料刚度矩阵不是正定的，或者出现负特征值警告。我的建议是：拿不准的时候先保证E2=E3，Nu12=Nu13，G12=G13，这样即使参数不是完美准确，至少数值上是自洽的。

如果还要做纱线内部损伤，可以进一步加上Hashin失效准则：

*Damage Initiation, criterion=HASHIN 4000., 3500., 200., 180., 260., 250.

这行参数分别对应纤维拉伸强度、纤维压缩强度、基体拉伸强度、基体压缩强度、纵向剪切强度、横向剪切强度。如果只是做线弹性细观应力分析，Damage Initiation这段可以先不加，等后续扩展做渐进损伤时再补。

如果是基体材料（树脂），材料定义简单得多，各向同性弹性就够了：

*Material, name=EpoxyMatrix *Density 1.2e-09, *Elastic 3500., 0.35

基体材料的参数我一般会根据实际树脂体系调整，3500 MPa和0.35泊松比是常见环氧树脂的典型值。

2.3 第三步：把材料写回inp并关联截面

材料定义准备好了，接下来就是把它安放到inp的正确位置，并通过Solid Section把材料和单元集合关联起来。这里有一个关键点：ABAQUS中材料必须通过Section才能赋给单元。也就是说，光有Material没有*Solid Section，单元依然不知道用哪份材料算。

把材料定义放在inp的哪个位置有讲究。我的一般做法是放在Part定义之前，或者紧跟Part定义结束之后、*Assembly之前。ABAQUS支持在part级别定义section，也支持在assembly级别直接给单元集合赋section。TexGen导出的inp通常没有严格的part/instance层级，所以最省事的方式是直接在单元集合上定义截面。

需要插入的关键字块长这样：

*Solid Section, elset=Yarn_1, material=GlassEpoxy 1.,

这一行的含义是：给名为Yarn_1的单元集合赋予截面属性，材料使用前面定义的GlassEpoxy。后面的"1."是截面厚度或者积分点相关的参数，对于三维实体单元，ABAQUS会忽略这个数值，但你得保留这个占位符，否则语法报错。如果是针对基体集合：

*Solid Section, elset=Matrix, material=EpoxyMatrix 1.,

搞定这两步，整个inp就有了完整的几何加材料定义。保存后直接用ABAQUS Command提交或者导入CAE都能正常计算。注意一个细节：在inp里，材料定义关键字和截面定义关键字的顺序不能乱。Material必须出现在引用它的Solid Section之前，否则ABAQUS解析时找不到材料名称直接报错。

我习惯的做法是：导出的inp开头是Node，我先不动它，直接拉到文件最前面，在Node之前插入一大段材料定义和截面定义。这样逻辑最清晰，后续检查也方便。如下所示：

*Material, name=GlassEpoxy *Density 1.9e-09, *Elastic, type=ENGINEERING CONSTANTS 45000., 9500., 9500., 0.3, 0.3, 0.4, 5200., 3800., 3800. *Solid Section, elset=Yarn_1, material=GlassEpoxy 1., *Solid Section, elset=Yarn_2, material=GlassEpoxy 1.,

这样改完之后，把文件另存为新的inp，再打开ABAQUS/CAE，进入Property模块就可以看到每个section都有对应的材料了。计算出来的应力分布、纱线间载荷传递规律，至少从量级上是对的。

3. 不想手改inp？走ABAQUS/CAE的替代方案

3.1 导入inp后在Property模块批量赋材料

手改inp并不复杂，但不是所有人都习惯编辑文本文件。另一个常用做法是：先正常导入TexGen导出的inp，然后在ABAQUS/CAE的界面里补材料。操作路径是：打开ABAQUS/CAE，File → Import → Model，选择inp文件。导入后进入Property模块，左侧Model树里能看到一个Part或者直接看到各个Set（取决于TexGen导出的结构）。如果没有Part，模型树里的Sets会直接挂在Assembly下面，单元集合一目了然。

接下来创建材料：Module下拉框切到Property，点击Material Manager，新建材料，弹出材料编辑器后分别添加Density和Elastic。材料参数和前面一样，照着填入即可。创建完材料之后，点击Assign Section图标，在视口里选择纱线对应的单元集合（或者直接在工具栏的Set拾取器里面勾选），然后指定Section名称、材料名称。对第二个纱线集合重复操作。全部赋完之后，可以在Mesh模块查看单元属性确认：点击Assign Section检查每个集合的截面状态。

这个路线适合不喜欢动手改文本的人，尤其是刚开始接触inp的初学者。它和手改inp本质上是同一个操作，一个用文本表达，一个用图形界面表达。我个人的习惯是：如果模型复杂、纱线数量多，手改inp加上Excel批量生成材料块，反而比在CAE里一个个点更快。如果纱线就三五根，CAE界面操作也完全能接受。

3.2 材料方向校验：这步别省

在ABAQUS里给纱线赋完材料只是第一步。纱线是横观各向同性材料，它有一个"主方向"——沿着纤维方向。如果材料方向定义错误，载荷传递路径就全错了。TexGen导出inp时往往会生成*Orientation信息，这就是用来指导材料方向的关键。

如果你走的是CAE界面路线，导入inp之后，应该能在某个Set或者单元输出里看到方向相关的信息。如果CAE没有自动读入，建议打开inp确认*Orientation段落是否存在。*Orientation的定义格式通常长这样：

*Orientation, name=Ori_Yarn_1, system=RECTANGULAR 1.0, 0.0, 0.0, 0.0, 1.0, 0.0 3, 0.

这个定义的含义是：用第一个向量和第二个向量确定一个局部坐标系，第三个数字3指定绕某轴的旋转方式。在TexGen导出的文件里，这个方向通常是根据纱线路径的切向量自动生成的。如果你在inp里看到类似的段落，恭喜你，方向信息保住了。接下来只要在查看单元属性时打开Visualization里的Section剖切显示，就能直观地看到每个纱线的主方向箭头。我每次计算前都会花两分钟检查这个方向：如果纱线主方向大致沿着纱线的走向，那说明几何导出时方向没有丢。如果方向箭头一团乱麻，尤其是弯曲剧烈的区域，就要考虑回到TexGen里增大纱线分段数重新导出。

如果你的inp里压根没有*Orientation，也有补救办法：在ABAQUS/CAE里手动定义材料方向。选中对应纱线集合，点击Assign Material Orientation，指定一个柱坐标系或者通过三点定义方向。但说实话，手动定义几十根纱线的方向非常痛苦，我更推荐回到TexGen里确认导出时勾选了方向输出选项，一劳永逸。

3.3 Python脚本批量补材料的思路

如果你经常干织物流变分析这活儿，强烈建议学会用ABAQUS的Python脚本接口来做批量赋材料。这样每次TexGen导出inp后，只需要一个脚本就把所有材料、截面、方向一次性搞定，连手改inp的时间都省了。核心脚本思路如下：

from abaqus import * from abaqusConstants import * # 创建材料 mdb.models['Model-1'].Material(name='GlassEpoxy') mdb.models['Model-1'].materials['GlassEpoxy'].Density(table=((1.9e-09,),)) mdb.models['Model-1'].materials['GlassEpoxy'].Elastic( table=((45000.0, 9500.0, 9500.0, 0.3, 0.3, 0.4, 5200.0, 3800.0, 3800.0),), type=ENGINEERING_CONSTANTS) # 遍历所有集合，自动创建截面并赋材料 for setName in mdb.models['Model-1'].rootAssembly.sets.keys(): if setName.startswith('Yarn'): mdb.models['Model-1'].HomogeneousSolidSection( name='Sec_%s' % setName, material='GlassEpoxy') region = mdb.models['Model-1'].rootAssembly.sets[setName] mdb.models['Model-1'].rootAssembly.SectionAssignment( region=region, sectionName='Sec_%s' % setName)

这段脚本的思路很简单：先建材料，然后遍历所有以Yarn开头的集合，为它们批量创建同名的截面并指派。你可以把这段Python脚本保存成一个.py文件，每次TexGen导出inp后，在CAE的File → Run Script里跑一遍，材料赋值直接完成。实测在几十根纱线的模型中，整个流程几秒钟搞定，比手动操作快了不止一个量级。

4. 常见报错与排查速查表

4.1 高频报错对照表

实际操作中大家遇到的报错，其实就那么几类。我整理了一个速查表，基本覆盖了从TexGen到ABAQUS这条路上最常见的坑：

这些里面，我自己踩得最多的坑有两个：一个是"应力全为0"，一个是"单元面积为负"。

4.2 "应力全为0"的几种可能性

很多人第一次从TexGen导出inp算细观应力，结果拿到的云图是一片均匀的蓝（也就是应力为零），当场心态炸裂。这个现象的原因大概率是以下三种：材料没有成功赋上、单元类型不适合当前分析、或者载荷和边界条件本身没有正确传递到模型上。材料赋没赋上很好查，Visualization里查看单元Section状态一目了然。单元类型方面，TexGen默认的C3D8R在弯曲问题里可能出现沙漏现象，导致应力分布不合理，必要时可以换成C3D8I或C3D8H试试。如果材料已经正常赋上、单元类型也没问题，那就要回过头检查边界条件——周期性边界条件有没有完全约束住模型刚体位移，加载面有没有正确选到节点集合。很多时候模型一条边的节点没有绑定约束，结构出现机构位移，应力和应变自然算不出来。

还有一种比较多见的情况：TexGen导出小尺寸织物件时，模型只有一个单胞，体积很小，而载荷又给得不合理，导致应力在数值上小到几乎接近零。这不算bug，而是量纲和载荷量级的匹配问题。单胞模型建议通过施加均匀应变或周期性位移边界条件来获得应力场，而不是直接给力载荷。这个问题在细观力学分析中比较隐蔽，建议新手优先采用施加位移边界条件的做法。

4.3 网格畸变与负体积问题

TexGen生成的网格质量，在大多数情况下是可靠的，但有个前提：织物几何不能太极端。如果纱线截面高宽比过大，或者纱线路径曲率半径小于截面尺寸，网格划分时就会出现畸变甚至负体积单元。我遇到过一次在纱线交错位置出现单元自交叉的情况，导出的inp在ABAQUS里一运行，直接报"The area of 3 elements is zero or negative"。

解决办法分两步。第一个优先回到TexGen里优化几何：增大纱线分段数（默认的分段可能不够密）、调整截面插值方式、减少纱线间的穿透。第二个是在ABAQUS里用Mesh模块的Verify功能定位畸变单元，然后用Edit Mesh手动修复或者用remesh思路重新划分。如果畸变单元数量不多，手动修复还比较可行；如果大面积畸变，那基本上就是TexGen端几何模型的问题，不要强行在ABAQUS里硬修，回头改几何更高效。

4.4 多相模型的材料指定：不止纱线，还有基体和界面

有些织物件不止纱线和基体，还有界面相、涂层层、夹杂相等。TexGen导出时会把不同几何实体生成为不同的单元集合。在补材料的时候，一个典型的错误是只给纱线赋了材料，基体和界面全忘掉了，或者把纱线材料和基体材料赋反了。我见过有朋友做三相交织复合材料，结果把树脂基体的模量写成45000 MPa，整个应力分布完全失真。建议在赋完材料之后，进入Mesh模块，逐个集合点击查看Assign Section状态，确保每个集合对应的材料名称和实际想要的材料完全一致。这是最简单的验证方式，也是最后一道防错网。

4.5 进阶：cohesive界面单元的材料与方向

如果你要做织物的界面脱粘分析，需要在纱线表面插入cohesive单元，这时候材料缺失和方向错误的问题会更明显。cohesive单元需要定义牵引-分离本构，用Cohesive Section配合Material中的*Elastic, type=TRACTION指定刚度。方向方面，cohesive单元的初始法向必须和界面法向一致。TexGen导出的inp如果没有专门处理界面网格，那你是看不到cohesive单元的，需要先在TexGen里生成界面网格或者用ABAQUS的Embedded/cohesive技术另行构建。这部分已经不是单纯的材料缺失修复，而是网格策略问题。如果有相关需求，建议在TexGen的"Yarn Contact"或"Interface Mesh"模块里提前设置好，再导出网格。

5. 一些实操之后的个人习惯

搞TexGen和ABAQUS这套流程也有好几年了，踩过的坑不少，最后分享几个我一直在沿用的动作。第一个是每次TexGen导出inp后，在模型树里新建一个"材料接口.txt"文件，把本次模型用到的材料参数、纱线名称、单位制全部记录进去。虽然看起来是多了一步，但三个月后再翻出这个模型时，不用重新猜参数。第二个是保留TexGen的工程文件（tg3或json格式），只把inp当作临时交换格式。因为inp是单向导出的，想改几何还得回TexGen。这也意味着，如果项目中期要调整纱线间距、截面形状，直接改TexGen重新导出，比在inp里硬编辑安全得多。

最后再讲一个我常用的"土办法"：TexGen导出的inp文件，一拿到手就Ctrl+S另存为带材料版本的bak，然后在备份文件里补材料。这样万一补错了，随时能回到原始几何版本重新来。不要小看这个备份习惯。网格重导出虽然不麻烦，但重复劳动真的会消磨耐心，尤其在交期紧的时候，一个能回退的备份文件能帮你省下大把时间。希望这篇能帮你把TexGen和ABAQUS之间的"最后一公里"走顺。