结构仿真要点总结 第1篇
谈到固化变形，就很难回避固化残余应力。本质上来讲，固化过程导致的结果是复合材料构件内部产生固化残余应力，而工程上反映出来的，就是构件脱模后固化残余应力释放，产生固化变形，在此我没有把二者进行严格的区分。但需要指出的是，固化残余应力在宏细观不同尺度上，是有所不同的。宏观层面上，固化残余应力引起固化变形，带来的最主要问题是构件成型精度降低，引起装配问题；而细观尺度上的固化残余应力，也就是纤维、树脂以及界面相上的残余应力，还会进一步对构件的力学性能，例如横向拉压性能产生影响。我现在的想法是把二者放在一个统一的框架内进行计算求解，目前宏观到细观已经搞定，但是细观尺度上的残余应力如何影响宏观，还没有完全解决。
言归正传，固化变形可以说是连接复合材料结构设计、制造以及装配过程的一个核心问题之一。需要指出的是，目前在结构设计过程中，可能对于固化变形的关注并不明显，但相信随着复合材料结构设计制造一体化的进一步推动，在结构设计过程中考虑固化残余应力和固化变形，是不可避免的。
结构仿真要点总结 第2篇
复合材料固化过程内部温度场分布，一方面影响构件成型后内部固化残余应力的分布情况，另一方面过高的温度峰值以及固化不均匀对构件成型后的力学性能也要很大影响。针对固化温度场的研究，现阶段主要使用的子程序有HETVAL、FILM、DISP以及USDFLD等，在庄茁先生的著作中有部分源代码，忽忽沛幺同志也进行过介绍，在此不再赘述。
我想要强调的是，HETVAL程序最大的问题，就是无法定义随温度和固化度变化的导热系数和比热容。实际上，复合材料在固化过程中，其材料性能参数并不是固定不变的，大多数情况下，与内部温度、固化度的变化密切相关。此时，UMATHT子程序便可以解决导热系数以及比热容随温度、固化度的变化问题，这是HETVAL程序所完成不了的。
如果大家对UMATHT子程序感兴趣，回头我可以再专门写一篇HETVAL和UMATHT子程序的比较分析，二者在处理树脂固化反应放热上是一致的，但UMATHT可以定义导热以及比热随温度、固化度的变化。
仿真模型
结构模型来源于文献[1]，如图1-1所示，材料为AS4/8552，固化后单层厚度，构件铺层为[02/902]s，文献[1]中并没有进行相应温度场和固化度场的计算，在此我查了资料，找到了相应材料的热力学参数，进行温度场的仿真分析。这里我是利用HETVAL子程序进行树脂固化反应放热计算的。复合材料传热的问题，本质上可以看作是一个具有非线性内热源的热传导问题。
考虑到树脂基体的化学反应放热，并忽略树脂流动过程产生的热量传递，则复合材料制件固化过程中三维热传导方程可以用式（1-1）和（1-2）表示。在HETVAL程序中，FLUX(1)用于定义当前节点上的热流密度，即树脂的内部反应放热Q，见式（1-1）。STATEV一方面用于定义固化率（固化度的变化率，见式（1-2），以表征材料的固化动力学方程；另一方面用作子程序间传递的媒介。譬如常见的STATEV(1)表示固化度，即树脂固化程度；STATEV(2)表示固化速率。对于材料AS4/8552，其固化速率可以用公式（1-3）进行表示[2]。在温度场计算过程中，可以进一步计算出树脂玻璃化温度变化曲线，如公式（1-4）所示[2]，该曲线对于后续进行固化残余应力分析尤其重要。
（1-1）
（1-2）
（1-3）
（1-4）
式中：λx、λy、λz分别为全局坐标系下沿x，y，z三个方向复合材料制件的热传导系数；T为温度；Q为热生成率；ρc为复合材料密度；Cc为复合材料比热容；t为时间。详细参数的含义在此不细说，上述提到的论文中大家都可以找到。具体的取值见表1-1所示。温度场计算的边界条件为构件外表面均采用70W/m2·K的对流换热系数[2]，以此来表征热压罐内高温气体与复合材料构件之间的换热。固化过程中，初始温度和固化最终温度均为20℃，第一个保温平台温度为110℃，保温时长60min，第二个保温平台温度为180℃，保温时长120min，升降温速率均为2℃/min。固化工艺曲线和对流换热系数可以通过FILM子程序进行定义，也可以利用Abaqus的“Interaction”模块直接定义，如图1-2所示。
图1-1 结构模型[1]
图1-2 定义对流换热系数
表1 -1AS4/8552的热力学参数[2]
在这里需要说明两个问题。第一个问题是，我们在用实体单元做复合材料分析的时候，通常是采用Abaqus里面的复材单元模块进行铺层以及区域的选择（“Property”—“Creat Composite Layup”），而帮助文档里面明确说明“Continuum shell and solid composite layups are expected to have a single element through the entire thickness across a combination of all the regions specified in the composite layup. Each single element through the thickness contains the multiple plies that you defined in the ply table. If the region to which you assign your continuum shell or solid composite layup contains multiple elements through the thickness, each element will contain all of the plies that you defined in the ply table and the analysis results will not be as expected. If your model contains multiple continuum shell or solid elements through the thickness of a region, you can obtain correct results by defining a separate composite layup for each layer of elements. You can define a composite layup for each layer by selecting the elements of a native Abaqus/CAE mesh or the orphan elements when you specify the region of the layup. You must create a layup for each layer of elements.”这里大家都能看懂，我就不详细解释了。简单来说，就是你在厚度方向上定义一层网格和多层网格时，在复合材料建模的过程中需要注意厚度尺寸上网格数量的问题。用一层网格可能完全足够做力学分析，但是在传热过程中，更希望厚度上有更多层的网格，这样才能更好地观察到厚度方向上是否存在温度梯度或者固化不均匀。
第二个问题是，在使用Abaqus的复材模块时，是不能用来进行瞬时温度场分析的，这也是帮助文档里面提过的（The use of composite solids is limited to three-dimensional brick elements that have only displacement degrees of freedom.）。（不知道最新版本的Abaqus是否可以，如果可以了，请大家给我留言，指正我的错误，避免我误导他人。）在这里进行瞬时温度场分析时，我是采用厚度上四层网格，考虑到铺层为[02/902]s，因此我是两个连续相同的铺层作为一层网格，然后在“Property”模块人为的赋予每层网格材料方向。对于这个模型，四层网格应该足够了，但是如果构件厚度增加，厚度方向上可能需要更多的网格。（本来我想做一个4、8、12层网格的比较，但是到8层的时候，我的电脑表示它负担不了这些网格的计算了。。。）。然后在力学分析的时候，可以使用传热分析时的网格，然后采用复材模块进行铺层，将传热计算结果导入应力分析模型，进行固化变形计算。
固化过程温度和固化度计算结果如图1-3所示，可以看到，在第一个保温平台结束前，树脂实际上是几乎没有固化的，这一过程树脂黏度逐渐降低，以树脂流动为主，树脂固化主要发生在第二个升温阶段。最终固化完成后，树脂的固化度为，可以认为固化完全。以图1-3(c)中的A点，进行整个固化阶段温度、固化度以及Tg曲线的演化规律研究，如图1-4所示，可以看到由于树脂固化反应放热的原因，在第二个保温平台的初始时刻是存在温度过热现象的，但是并不明显，并且该温度峰值很快消退。因此可以认为构件在厚度方向上几乎是没有温度梯度发生的，在实际的固化残余应力计算过程中，对于该模型，可以将热-化学-力问题，转换为化学-力问题，从而减少建模和计算工作量。
（a）温度和固化度计算结果 6300s（第一个保温阶段结束）
（b）温度和固化度计算结果 8500s（第二个升温阶段结束）
（c）温度和固化度计算结果 固化完成
图1-3 复合材料结构固化过程温度和固化度计算结果
图1-4温度、固化度以及Tg曲线的演化规律
结构仿真要点总结 第3篇
可以看到，上述的计算模型与结构其实都相对简单，并且所采用的本构模型也很容易编入UMAT。但实际工程应用中的结构要比文中的结构复杂的多，图3-1是最近几年做过的一些实际结构。在此有几个问题，是我在一直思考的，希望可以集思广益。
（1）树脂的流动主要发生在液态（黏流态）阶段，这个时候树脂实际上并没有开始固化，或者固化度非常小，那这个阶段究竟对固化残余应力有多少影响？真的可以忽略吗？
（2）对于结构复杂的大尺寸构件，例如加筋壁板，采用实体单元进行计算固然准确，但计算效率如何保证？而采用壳单元，是否可以将固化变形的诱导因素都考虑在内？
（3）现有的模型多是分别针对不同类型结构的，例如L型或者U型主要考虑圆角处的回弹，平板则是考虑模具作用或者铺层角度的变化，但如果对于厚度较厚的复杂形状构件，则可能发生纤维体积分数分布不均、固化不均等其他因素，如何能在复杂模型中考虑尽可能多的因素，也是目前建模需要考虑的问题。
（4）虽然现有一些软件，例如ESI-RTM，Compro等软件，也可以做复合材料固化仿真，但如何突破国外软件的封锁，在复材领域也任重道远。
图3-1 复杂结构零件固化变形
参考文献
[1] Wucher B, Lani F, Pardoen T, et al. Tooling geometry optimization for compensation of cure-induced distortions of a curved carbon/epoxy C-spar[J]. Composites Part A, 2014, 56:27-35.
[2] Bellini C, Sorrentino L. Analysis of cure induced deformation of CFRP U-shaped laminates[J]. Composite Structures, 2018, 197:1-9.
[3] Ersoy N, Garstka T, Potter K, et al. Modelling of the spring-in phenomenon in curved parts made of a thermosetting composite[J]. Composites Part A, 2010, 41(3):410-418.
[4] Cinar K, Ersoy N . Effect of fibre wrinkling to the spring-in behaviour of L-shaped composite materials[J]. Composites Part, 2015, 69:105-114.
[5] Cinar K, Oztrk UE, Ersoy N, et al. Modelling manufacturing deformations in corner sections made of composite materials. J Compos Mater, 2014, 48(7):799–813