基于压力隧洞模型的复合材料横向热残余应力分析

考虑碳纤维横向刚度情况下,将复合材料横向热残余应力问题简化为压力隧洞模型.利用该模型推导出复合材料固化成型过程中形成的横向热残余应力,其中分析得出了单纤维与树脂的接触压力以及该压力传递到纤维和树脂后的分布情况.结果表明接触压力传递到纤维内部后成一固定值,传递到树脂后以正比于r.函数衰减.在其基础上提出场叠加方法,得出纤维之间相互耦合的接触压力与残余应力场.通过有限元模拟,理论模型和数值模拟基本一致.     

光纤模斑谱传感器复合材料固化监测研究

提出利用一种新型光纤传感器，通过测量光纤末端近场模斑谱反映光纤埋置周围树脂折射率的变化．计算了随着树脂固化过程中折射率变化的功率谱密度分布，给出截止频率与折射率的近似关系．给出了利用这种传感器进行复合材料树脂固化监测的实验结果．发现这种光纤传感器的信号可以用来反映固化凝胶点和固化结束点．由于固化后光纤仍保持光波导特性，该光纤还可以用来监测结构受环境的扰动如温度变化、振动等．     

固化温度对碳纤维/环氧基体界面行为影响的分析

对界面粘结性能及热残余应力影响下的单纤维复合材料的界面行为进行了分析。采用界面的弹性-软化内聚力模型,用解析法对单纤维复合材料由固化引起的热残余应力、以及单纤维碎断过程纤维的轴向应力分布进行了模拟,得到了碳纤维/环氧树脂在常温和高温固化两种情况的界面粘结性能。结果表明：与常温固化相比,高温固化后,界面的剪切强度增幅不大,界面的断裂韧性显著增加;高温固化后形成的界面,使界面的软化提前、界面的脱粘延迟;高温固化产生的纤维轴向和界面径向热残余应力对界面的软化均有延迟作用;界面径向热残余应力还对界面的脱粘有延迟作用。     

三维粘弹性复合材料层合板固化残余应力研究

应用改进的有限元算法,综合考虑了碳纤维复合材料固化过程中树脂的化学收缩、应力松弛、材料的热收缩,建立了一个基于三维线性粘弹性的积分模型对复合材料层合板热压罐成型工艺过程中产生的残余应力进行了研究。结果表明:本文提出的方法能够完整地模拟树脂的固化过程并计算层合板的固化残余应力;对于文中的正交对称层合板而言,热收缩是影响残余应力大小的主要因素;固化完成后,0°铺层内的中面中心区域层间残余应力值接近于0,边缘区域会迅速增加,容易发生分层破坏;任意点处,轴向应力比其他两个方向的应力高出1~2个数量级;残余应力值沿厚度方向有较大变化;树脂的应力松弛行为会在一定程度上降低残余应力。     

界面特性对短纤维金属基复合材料蠕变行为的影响

基于短纤维增强金属基复合材料(MMC)的单纤维三维模型(三相)，利用粘弹性有限元分析方法对影响金属基复合材料的蠕变行为的因素进行了较为系统的分析。研究中主要讨论了界面特性和纤维取向角对金属基复合材料的蠕变性能的影响。研究结果发现，界面特性诸如厚度、模量和应力指数都对纤维最大轴应力和稳定蠕变率产生影响：稳态蠕变率随界面模量的增大而逐渐减小，当高于基体模量时基本保持不变；纤维轴应力的变化与蠕变率正好相反。稳态蠕变率随界面厚度、应力指数的增加而增大；而轴应力则随之减小。同时不同的纤维取向也影响金属基复合材料蠕变时的轴应力分布和稳态蠕变率。     

含界面相复合材料热残余应力分析

本文分析了含有界面相纤维增强复合材料热残余应力的空间分布。针对材料实际微结构几何特点,建立含有界面相的三维三相单丝模型,用均匀和梯度函数描述界面相模量随空间变化规律,由轴对称体弹性力学理论得到单丝热残余应力分布,结果表明梯度界面降低了残余应力。通过碳纤维电阻法测出T300/环氧树脂单丝体系固化后的纤维轴向应变,与梯度界面的分析结果基本一致。用叠加方法得到密排六方结构代表性体积元(RVE)中纤维间相互偶合的应力场,同时应用有限元法分析RVE中纤维间的残余应力分布,两者结果相互验证。     

缝纫复合材料层合板面内弹性模量分析

基于对缝纫孔附近局部细观结构的分析,提出了一种预测缝纫复合材料层合板面内力学性能的理论模型．从分析缝孔单胞的纤维弯曲几何特征入手,最终得到单向板及层合板的弹性常数．通过有限元分析研究了缝纫参数对复合材料层合板面内等效模量的影响．研究结果表明,缝纫造成单向板及层合板面内材料性质的不均匀,随着缝纫密度和缝纫线直径的增加,层合板的等效模量逐渐降低．     

纺织复合材料和结构多尺度耦合的数值分析

研究了纺织复合材料和结构多尺度耦合的数值分析模型。建立了微、细观单胞,给出了纺织复合材料平均弹性常数的逐级分析方法,着重研究了由宏观结构、到细观纤维束、再到微观纤维三个尺度耦合的应力分析方案。对于常用的板壳状纺织复合材料结构,在面内载荷下,假设每层细观单胞的平均面内应变是一致的,在弯曲、横向剪切及扭曲等非面内载荷下,在内力等效条件下将沿厚度方向连续分布的宏观应力简化为阶梯状分布,忽略了每层细观单胞范围内宏观应力沿厚度方向的梯度变化,由此利用细观单胞模型实现宏观应力与细观应力之间的传递,再利用微观单胞可得到纤维尺度的微观应力。最后以一种三维机织复合材料为例,用上述多尺度耦合的模型逐级分析了材料的平均弹性常数,并沿相反方向,由宏观结构分析逐级计算出纤维束尺度和纤维尺度的细、微观应力的局部波动。     

NiMnGa颗粒/环氧树脂复合材料的力学性能研究

将通过球磨法获得的Ni2MnGa颗粒与双酚A二缩水甘油醚环氧树脂混合均匀并常温固化,制备出相应的复合材料试件;通过拉伸-卸载试验研究了NiMnGa颗粒对树脂基体的增强效果;以试验为基础,从细观力学的角度分析了NiMnGa颗粒对树脂基复合材料弹性特性的影响;采用商用有限元软件Abaqus,对NiMnGa/树脂复合材料单胞模型进行了仿真计算。分析结果表明:随着NiMnGa颗粒体积含量的增加,试件的拉伸模量和剪切模量增大、泊松比减小;随着椭圆NiMnGa颗粒长短半轴比的减小,复合材料纵向泊松比增加、纵向拉伸模量减小,横向拉伸模量呈先减后增的趋势,剪切模量则呈先增后减的趋势。通过对比可知试验值、理论值、有限元模拟值数值相近且变化趋势一致。     

高温环境下纤维复合材料蠕变损伤的细观机理研究

首先利用复合材料纤维断裂单胞模型，编制蠕变损伤子程序，对单胞模型进行蠕变损伤分析。分析了纤维／基体弹性模量比对蠕变变形、蠕变损伤以及应力场的影响。从计算结果发现，蠕变损伤首先在纤维断裂尖端起始，然后沿着一定的角度向基体外围延伸，直至完全损伤，而且纤维／基体模量比对高温环境下的复合材料蠕变损伤产生很大的影响；纤维与基体的模量相差越大，复合材料越容易变形，抵抗蠕变变形的能力就越小，蠕变损伤越严重。经过对不同韧性的基体材料进行研究，发现基体韧性低的复合材料蠕变损伤明显高于高韧性基体复合材料，表明低韧性基体复合材料抵抗蠕变破坏的能力较低。     

Analysis of Composite Molding Using the Homogenization Method

In this paper, we present an application of the homogenization method to the analysis of Resin Transfer Molding (RTM) and Structural Reaction Injection Molding (SRIM). RTM and SRIM are relatively new molding processes for manufacturing continuous fiber reinforced polymer composites. First, the mold flow is analyzed. In the molding process, the resin experiences significant temperature changes as it fills the mold and forms a free boundary with air as it pushes out the air. In addition, the flow domain is the mold cavity packed with fiber perform, which is a porous medium. Here, the homogenization method is used to model the non-isothermal flow through porous media with free boundaries. A computer program is developed which is capable of simulating a three-dimensional mold flow using the finite element approximation. An example is provided for a three-dimensional part. Then, an analysis of the residual stress developed in the curing stage is given. The curing stage starts when the mold is completely filled and it involves chemical reaction and large temperature variation. In curing, the resin part undergoes larger volume shrinkage than the fiber part, and the residual stresses are developed due to this volume mismatch. In some cases, these stresses are large enough to cause micro-cracking and to exhaust the strength of the material. Here, a brief discussion of the application of the homogenization method to a residual stress analysis is given and one example is provided.     

A Dynamic Method for the Residual Stress Measurement During Polymer Crystallization

A new dynamic method based on bilayer system is proposed to characterize the residual stress formation during the crystallization of semi-crystalline polyethylene glycol 10000 (PEG10000). The resin is coated on a solid polymeric film to create a filmsubstrate compound. Its temperature field and dynamic deflection are monitored by synchronized optical and thermography cameras. The crystallization kinetics is first characterized from the former information. Then a simple dynamic model is proposed to relate the dynamic deflection with crystallization process. Residual stresses are established and in the range of 0–2.1 MPa. The generation of residual stresses is due to the edge constraints of the cantilever beam and to the increase of viscosity during solidification that allows the polymer to carry tensile loadings. The spherulite impingement is found to be important for this period from a microscopic view. Boundary condition should be well controlled to steer residual stresses. Such method is promising to measure residual stresses at the micro-scale for polymers to be spread on a flexible substrate and can mimic different mechanical situations of interest.     

X状Z向碳pin增强泡沫夹层结构剪切刚度预报

综合考虑了面板对横向增强Z-pin的不同约束情况,结合空间网架结构和等效夹杂方法,提出了X状Z-pins增强泡沫夹层结构剪切刚度的预报模型,经实际计算后与已有的实验值和有限元模拟值比较,证明该方法具有足够的工程精度.计算表明,X状Z-pins能大幅度地提高泡沫夹层结构的剪切刚度,并具有良好的可设计性,可以通过改变Z-pin材料,Z-pin的体积分数、角度、直径等参数改变其力学性能.其中,Z-pin体积分数越大、拉伸模量越高、直径越大,Z-pins增强泡沫夹层结构的剪切刚度越高.     

Viscoelastic analysis of residual stress in quenched thermosetting resin beams

The residual stress generated by the molding process of thermosetting resins exerts serious influences upon their mechanical properties. This residual stress is generally classified by two groups: one produced by shrinkage in the curing reaction of monomers, the other produced by the nonuniformity of the temperature distribution in the cooling process. This paper is concerned with the theoretical and experimental analysis of the generation of residual stress of the latter type, using examples of rectangular beams of thermosetting resins quenched on both the upper and lower surfaces. First, a viscoelastic model is applied to make a qualitative prediction of the residual stress in quenched beams. Second, using linear-viscoelastic theory, fundamental equations are derived for the residual stress in a viscoelastic rectangular beam, where an unsteady and nonuniform temperature distribution is assumed in the direction of depth. The theoretical values of the residual stress in rectangular beams are calculated under various quenching conditions for two resins having different viscoelastic characteristics, i.e., epoxy and unsaturated polyester. The theoretical residual-stress distributions agree fairly well with the residual stress measured experimentally at every quenching condition for both resins. The qualitative prediction that the residual stress in quenched beams is compressive in the vicinity of the upper and lower surfaces and is tensile in the inner parts is confirmed. The relaxation modulus of epoxy resin changes more greatly with time and temperature than that of unsaturated polyester resin. The theoretical and experimental analysis shows that the residual stress for the former resin is larger than that for the latter. Therefore, it is concluded that the generation of residual stress is more significant where the relaxation modulus of resin changes greatly with time and temperature.     

A micromechanical study of residual stress and its effect on transverse failure in polymer–matrix composites

Cure residual stress and its effect on damage in unidirectional fibre-reinforced polymer–matrix composites under transverse loading were studied using a micromechanical unit cell model and the finite element method. The overall residual stress introduced from curing was determined by considering two contributions: volume shrinkage of matrix resin from the crosslink polymerization during isothermal curing and thermal contraction of both resin and fibre as a result of cooling from the curing temperature to room temperature. To examine the effect of residual stress on failure, a model based on the Maximum Principal Stress criterion and stiffness degradation technique was used for damage analysis of the unit cell subjected to mechanical loading after curing. Predicted damage initiation and evolution are clearly influenced by the inclusion of residual stress. Residual stress is always detrimental for transverse compressive loading and pure shear loading. For transverse tensile loading, residual stress is detrimental for relatively low resin strength and beneficial for relatively high resin strength. Failure envelopes were obtained for both biaxial normal loading and combined shear and normal loading and the results show that residual stress results in a shifting and contraction of the failure envelopes.     

磁电复合材料中拓扑磁结构的力学调控

磁性斯格明子是在一些铁磁材料中存在的一种重要拓扑磁结构,由于其具有独特的磁-电-力-热多场耦合特性,在未来新型自旋电子器件中有着广泛的应用前景。然而,磁性斯格明子一般需要在外加磁场下才能稳定存在,极大地限制了其在自旋电子器件中的实际应用。本文基于实空间下磁电材料的相场模拟,发现铁电和铁磁复合薄膜中铁电斯格明子可以通过界面变形来稳定铁磁斯格明子。由于力电耦合效应,铁电层中铁电斯格明子的非均匀分布极化在界面产生周期性的非均匀界面变形。界面变形通过力磁耦合效应,使铁磁层中的磁性斯格明子在没有外加磁场的条件下能够稳定存在。本文的研究结果表明,基于磁电复合材料中的力-电-磁耦合效应,通过优化设计复合材料中不同组元的结构,可以实现拓扑磁结构的力学调控,从而为设计基于拓扑磁结构的新型自旋电子器件提供了新思路。     

空心玻璃微珠及玻璃纤维混合填充环氧树脂力学性能实验研究

制备了一种空心玻璃微珠质量比、四种玻璃纤维质量比的空心玻璃微珠及玻璃纤维混合填充环氧树脂复合泡沫材料。通过一系列准静态压缩、拉伸实验,研究了玻璃纤维填充量对复合材料的密度、强度、弹性模量等性能的影响,并分析了其影响原因。通过三点弯曲实验,研究了复合材料弯曲弹性模量、弯曲强度等力学性能与玻璃纤维填充量的关系。通过拉伸断口SEM照片分析了玻璃纤维增强复合材料力学性能的原因。研究表明:随着玻璃纤维填充量增加,复合材料的拉伸和弯曲强度、压缩屈服极限及弹性模量等力学性能都有较大的提高,且仍然保持了较好的弹性性能和塑性特征。     

Design and fabrication of submerged cylindrical laminates—II. Effect of fiber pre-stress

This is the sequel to the first part of this paper (Dvorak et al., 1999, Int. J. Solids Structures 36, 3917–3943) , concerned with modeling and analysis of laminated composite cylinder fabrication procedures, such as filament winding or fiber placement, which involve fiber pre-stress for waviness reduction as well as overall or local heating to and cooling from matrix curing temperatures. The fiber pre-stress applied in individual plies is shown to cause a self-stress in the respective plies, and relaxation stresses in the already completed plies and in the supporting mandrel. The final residual stress state is reached after mandrel removal. Influence functions that relate the ply stresses to the applied pre-stress forces are derived. Direct problems are solved for ply stresses caused by prescribed constant or linearly or parabolically changing pre-stress magnitudes in the layers. A superposition of the constant and parabolic distributions is shown to lead to nearly uniform stresses through the cylinder wall. The magnitudes depend on the radial stiffness of the mandrel that supports the structure during fabrication. Inverse problems are formulated as nonlinear optimizations and solved by quadratic programming. The goal is to determine fiber pre-stress distributions through the wall thickness such that the total stresses due to external hydrostatic pressure and fiber pre-stress are as uniform as possible through the wall thickness and confined by the ply strength magnitudes.