|
- 2015
树脂基复合材料曲面结构件固化变形数值模拟
|
Abstract:
为了研究树脂基复合材料曲面结构件的固化变形过程, 首先分析了碳纤维增强树脂基复合材料在固化过程中密度、模量、热膨胀系数、比热容及热传导系数等材料物性的变化, 并将这些变化引入到数值模拟当中。接着, 针对复合材料复杂曲面结构件, 提出了利用定常流动的流线方程构建曲线坐标系的新方法。然后, 根据建立的曲线坐标系, 运用有限元法计算了某轻型飞机机翼上蒙皮板在固化过程中内部温度、固化度和内应力的分布情况以及材料物性随固化度的变化情况。最后, 计算了由于内部温度场和固化度场的不均匀、热膨胀系数的各向异性和固化引起的树脂体积收缩而导致的结构变形。结果表明:引入材料物性变化使固化过程的数值模拟更加合理、模拟结果更加精确, 利用定常流动的流线方程构建的曲线坐标系适用于复合材料曲面结构件的有限元分析。所得结论对研究树脂基复合材料的固化变形过程和各向异性复合材料复杂曲面构件的三维实体建模均具有指导意义。 In order to investigate the curing deformation process of resin matrix composite curved structure, first, the evolutions of physical properties such as density, modulus, thermal expansion coefficient, specific heat capacity and thermal conductivity et al during the curing process of carbon fiber reinforced resin matrix composites were analyzed, and those evolutions were introduced into the numerical simulation. Second, the new approach of constructing the curvilinear coordinate system by streamline equation of steady flow was proposed for the composite complex curved structure. Third, the distribution of internal temperature, curing degree and internal stresses, and evolution situations of material physical properties along with curing degree during the curing process of the skin plate of a light aircraft wing were calculated by finite-element method based on the curvilinear coordinate system constructed. Finally, the deformations of the structure caused by the uneven distributions of internal temperature field and cure degree field, anisotropy of thermal expansion coefficient and volumetric shrinkage of resin result from curing were calculated. The results show that the numerical simulation of curing process will be more reasonable and the simulation results will be more accurate when the evolutions of material physical properties are adopted, and the curvilinear coordinate system constructed by streamline equation of steady flow is applicable to the finite element analysis of composite curved structure.The conclusions have instruction significance to the research of curing deformation process of resin matrix composites and the establishment of three dimensional solid model of composite complex curved structure which is anisotropic. 国家"973"计划(2014CB046500)
[1] | Costa V, Sousa A. Modeling of flow and thermo-kinetics during the cure of thick laminated composites[J].International Journal of Thermal Sciences, 2003, 42(1):15-22. |
[2] | Michaud D J, Beris A N, Dhurjati P S. Thick-sectioned RTM composite manufacturing, part II. Robust cure cycle optimization and control[J]. Journal of Composite Materials, 2002, 36(10): 1201-1230. |
[3] | Chern B C, Moon T J, Howell J R, et al. New experimental data for enthalpy of reaction and temperature-and degree-of-cure-dependent specific heat and thermal conductivity of the Hercules 3501-6 epoxy system[J]. Journal of Composite Materials, 2002, 36(17): 2061-2072. |
[4] | Bogetti T A, Gillespie J W. Process-induced stress and deformation in thick-section thermoset composite laminates[J]. Journal of Composite Materials, 1992, 26(5):626-660. |
[5] | Bogetti T A, Gillespie J W. Two-dimensional cure simulation of thick thermosetting composites[J]. Journal of Composite Materials, 1991, 25(3): 239-273. |
[6] | Prasatya P, McKenna G B, Simon S L. A viscoelastic model for predicting isotropic residual stresses in thermosetting materials: Effects of processing parameters[J]. Journal of Composite Materials, 2001, 35(10): 826-849. |
[7] | Yan X. Finite element modeling of curing of epoxy matrix composites[J]. Journal of Applied Polymer Science, 2007, 103(4): 2310-2319. |
[8] | Svanberg J M, Holmberg J A. Prediction of shape distortions Part I. FE-implementation of a path dependent constitutive model[J]. Composites Part A: Applied Science and Manufacturing, 2004, 35(6): 711-721. |
[9] | Liang G, Chandrashekhara K. Cure kinetics and rheology characterization of soy-based epoxy resin system[J]. Journal of Applied Polymer Science, 2006, 102(4): 3168-3180. |
[10] | Yan X. Consolidation simulation of composite laminates: Formulation andvalidation verification[J]. Polymer Composites, 2005, 26(6): 813-822. |
[11] | Svanberg J, Holmberg J. Prediction of shape distortions. Part II. Experimental validation and analysis of boundary conditions[J]. Composites Part A: Applied Science and Manufacturing, 2004, 35(6): 723-734. |
[12] | Kim Y K, White S R. Process-induced stress relaxation analysis of AS4/3501-6 laminate[J]. Journal of Reinforced Plastics and Composites, 1997, 16(1): 2-16. |
[13] | Loos A C, Springer G S. Curing of epoxy matrix composites[J]. Journal of Composite Materials, 1983, 17(2): 135-169. |
[14] | Scott E P, Beck J V. Estimation of thermal properties in epoxy matrix/carbon fiber composite materials[J]. Journal of Composite Materials, 1992, 26(1): 132-149. |
[15] | Tavman I, Akinci H. Transverse thermal conductivity of fiber reinforced polymer composites[J]. International Communications in Heat and Mass Transfer, 2000, 27(2): 253-261. |
[16] | Tan H, Yan S L. Three-dimensional simulation of curing process for thermoset composites[J]. Acta Materiae Compositae Sinica, 2004,21(6):167-172 (in Chinese). 谭华, 晏石林. 热固性树脂基复合材料固化过程的三维数值模拟[J]. 复合材料学报, 2004,21(6): 167-172. |
[17] | Shin D D, Hahn H T. Compaction of thick composites: Simulation and experiment[J]. Polymer Composites, 2004, 25(1): 49-59. |