织物复合材料的建模
理论与物理
什么是织物复合材料
老师，“织物复合材料”和UD材料（单向增强材料）有什么区别？
UD材料的纤维全部朝一个方向排列，但织物复合材料的纤维是编织（交织）在一起的。具有平纹、斜纹、缎纹等织造结构。
织造有什么优点？
织物的种类
| 种类 | 结构 | 特点 |
|---|---|---|
| 平纹（Plain Weave） | 1根纤维交替交织 | 稳定。铺覆性低 |
| 斜纹（Twill Weave） | 2/1, 2/2等图案 | 铺覆性高 |
| 缎纹（Satin Weave） | 5HS, 8HS（长浮线） | 铺覆性最高。纤维几乎呈直线 |
| NCF（Non-Crimp Fabric） | 纤维不编织，用缝线固定 | 无屈曲。机械性能最佳 |
“屈曲”是什么？
指在织造结构中纤维上下起伏的现象。存在屈曲时，由于纤维是弯曲的，拉伸刚度和强度会比UD材料低10~20%。NCF没有屈曲，因此具有接近UD材料的性能。
织物的FEM建模
织物复合材料的FEM建模有三个层次：
| 层次 | 方法 | 精度 |
|---|---|---|
| 宏观 | 等效均质壳（基于CLT） | 低（整体行为） |
| 细观 | RVE（纤维束+基体）建模 | 高（局部应力） |
| 微观 | 单根纤维建模 | 最高（研究用） |
细观尺度的RVE模型实用吗？
细观尺度下，用实体单元对织造图案的一个单元（Unit Cell）进行建模，并通过周期边界条件计算等效特性。可以使用TexGen或WiseTex等专用工具自动生成Unit Cell形状。
总结
我来梳理一下织物复合材料的理论。
要点：
- 纤维交织结构 — 平纹、斜纹、缎纹、NCF
- 屈曲导致刚度/强度降低 — 比UD材料低10~20%
- 铺覆性和耐损伤性优异 — 适合曲面成形
- 三个层次的建模 — 宏观（CLT）、细观（RVE）、微观（单根纤维）
- 细观尺度的RVE分析实用 — 使用TexGen/WiseTex生成Unit Cell
织物复合材料的弹性特性与单位胞理论
织物复合材料由经纱（Warp）和纬纱（Weft）交叉的单位胞重复构成。其弹性特性通过单位胞的FEM均质化分析计算，比ISO 527的简单测试值更能准确反映纤维屈曲效应。仅5%的屈曲就可使强度降低10~20%，且其面内各向同性优于单向预浸料。
各项的物理意义
- 惯性项（质量项）：$\rho \ddot{u}$，即“质量×加速度”。您有没有过急刹车时身体被向前甩出的经历？那种“被带走的感觉”正是惯性力。物体越重越难启动，一旦启动也越难停止。地震时建筑物摇晃，也是因为地面突然移动，而建筑物的质量“被落下”。静力分析中此项设为零，这是假设“缓慢施力，加速度可忽略”。但在冲击载荷或振动问题中绝对不能省略。
- 刚度项（弹性恢复力）：$Ku$ 或 $\nabla \cdot \sigma$。拉弹簧时能感觉到“想要恢复的力”吧？那就是胡克定律 $F=kx$，也是刚度项的本质。那么提问——铁棒和橡皮筋，用相同的力拉，哪个伸得更长？当然是橡皮筋。这种“难拉伸的程度”就是杨氏模量 $E$，它决定了刚度。常见的误解：“刚度高=强度高”是不对的。刚度是“不易变形的程度”，强度是“不易破坏的程度”，是不同的概念。
- 外力项（载荷项）：体积力 $f_b$（重力等）和表面力 $f_s$（压力、接触力等）。可以这样想——桥上卡车的重量是“作用在整个内部上的力”（体积力），轮胎压路面的力是“只作用在表面上的力”（表面力）。风压、水压、螺栓紧固力…这些都是外力。这里容易犯的错误：弄错载荷方向。本想“拉伸”却成了“压缩”——听起来像笑话，但在3D空间坐标系旋转时确实会发生。
- 阻尼项：瑞利阻尼 $C\dot{u} = (\alpha M + \beta K)\dot{u}$。弹一下吉他弦试试。声音会一直响吗？不，会逐渐变小。这是因为振动能量通过空气阻力和弦的内部摩擦转化成了热。汽车的减震器也是同样原理——特意吸收振动能量来提高乘坐舒适性。如果阻尼为零会怎样？建筑物在地震后会一直摇晃不停。实际上不会这样，所以设置合适的阻尼很重要。
假设条件与适用范围
量纲分析与单位制
| 变量 | SI单位 | 注意事项·换算备忘 |
|---|---|---|
| 位移 $u$ | m（米） | 输入mm时，载荷·弹性模量也需统一为MPa/N系 |
| 应力 $\sigma$ | Pa（帕斯卡）= N/m² | MPa = 10⁶ Pa。与屈服应力比较时注意单位制不一致 |
| 应变 $\varepsilon$ | 无量纲（m/m） | 注意工程应变与对数应变的区别（大变形时） |
| 弹性模量 $E$ | Pa | 钢：约210 GPa，铝：约70 GPa。注意温度依赖性 |
| 密度 $\rho$ | kg/m³ | mm系中为tonne/mm³（钢约为 10⁻⁹ tonne/mm³） |
| 力 $F$ | N（牛顿） | mm系用N，m系也用N统一 |
数值解法与实现
细观尺度RVE分析
请讲解一下RVE（代表体积单元）分析的步骤。
1. 生成Unit Cell的几何形状 — 使用TexGen, WiseTex等生成织造图案的3D形状
2. 网格生成 — 用TET10单元对Unit Cell进行网格划分
3. 应用周期边界条件 — 相对面的位移呈线性关系
4. 6个载荷工况 — 依次施加 $\varepsilon_{11}, \varepsilon_{22}, \varepsilon_{33}, \gamma_{12}, \gamma_{23}, \gamma_{13}$
5. 均质化 — 根据各载荷工况的平均应力计算等效弹性常数
从6个载荷工况确定9个弹性常数（正交各向异性）对吧。
是的。可以得到 $E_1, E_2, E_3, G_{12}, G_{23}, G_{13}, \nu_{12}, \nu_{23}, \nu_{13}$。将其用作宏观尺度壳/实体单元的材料特性。
专用工具
| 工具 | 特点 |
|---|---|
| TexGen | Unit Cell几何生成。开源（诺丁汉大学） |
| WiseTex | 织物几何+力学。KU Leuven开发 |
| DIGIMAT | 多尺度材料建模。eXstream/Hexagon |
| MicroMechanics | 纤维-基体RVE分析。MCT（多连续介质理论） |
TexGen免费真不错。
TexGen可自动生成织造结构的3D几何，并直接输出为Abaqus输入文件。在研究用的织物RVE分析中是事实上的标准工具。
宏观尺度下织物的处理
在宏观尺度（常规FEM分析）中如何处理织物？
作为等效均质材料处理。将RVE分析得到的等效弹性常数作为CLT中的一层，进行常规的层合板分析。
注意事项：
- 屈曲效应 — 确认等效特性是否包含屈曲影响
- 破坏判据 — 直接应用UD材料的Tsai-Wu/Hashin准则不准确。需要针对织物进行修正
- 铺覆 — 在FEM中反映成形时的纤维角变化
总结
我来梳理一下织物复合材料的数值方法。
要点：
- 通过RVE分析计算等效特性 — 从6个载荷工况得到9个弹性常数
- 使用TexGen（免费）生成Unit Cell — 研究的标准工具
- 使用DIGIMAT（商用）进行多尺度耦合 — 面向产业
- 宏观尺度作为等效均质材料处理 — CLT中的一层
- 不直接使用UD材料的破坏准则 — 需要针对织物进行修正
织物复合材料的均质化FEM分析
单位胞的均质化分析按以下顺序进行：①构建胞元的三维FEM模型②应用周期边界条件③根据6个分量的单位载荷对应的反力计算等效弹性矩阵。计算成本约为一个胞元FEM（约10万个单元）求解6次，在数小时以内。Software TexComp（KU Leuven）、WiseTex、TexGen等专用的unit cell meshing工具已公开，只需输入织物参数（纱线宽度·厚度·屈曲率）即可自动生成网格。
线性单元（1次单元）
节点间线性插值。计算成本低，但应力精度低。注意剪切锁定（可通过减缩积分或B-bar法缓解）。
2次单元（带中间节点）
可表现曲线变形。应力精度大幅提高，但自由度约增加2~3倍。推荐：应力评估重要时使用。
完全积分 vs 减缩积分
完全积分：有过约束（锁定）风险。减缩积分：有沙漏模式（零能模式）风险。根据情况选择。
自适应网格
基于误差指标（ZZ估计量等）的自动细化。有效提高应力集中部位的精度。有h法（单元细分）和p法（增加阶次）。
牛顿·拉弗森法
非线性分析的标准方法。每次迭代更新切线刚度矩阵。在收敛半径内具有二次收敛性，但计算成本高。
修正牛顿·拉弗森法
切线刚度矩阵使用初始值或每隔数次迭代更新。每次迭代成本低，但收敛速度为线性。
收敛判定准则
力残差范数: $||R|| / ||F_{ext}|| < \epsilon$（通常 $\epsilon = 10^{-3}$〜$10^{-6}$）。位移增量范数: $||\Delta u|| / ||u|| < \epsilon$。能量范数: $\Delta u \cdot R < \epsilon$
载荷增量法
不一次性施加全部载荷，而是分小步增加。弧长法（Riks法）可以越过载荷-位移关系的极值点进行追踪。
直接法 vs 迭代法的比喻
直接法是“用笔算精确解联立方程”的方法——可靠但大规模问题耗时过长。迭代法是“反复猜测逼近正确答案”的方法——最初答案粗略，但每次迭代精度提高。就像查字典时，从第一页开始顺序查找（直接法）不如先估计位置翻开，再前后调整（迭代法）来得高效。
网格阶次与精度的关系
1次单元是“用直尺近似曲线”——用直线折线表现，精度有限。2次单元是“柔性曲线”——能表现曲线变化，在相同网格密度下……
相关主题
なった
詳しく
報告