连续纤维复合材料微观结构的量化分析.docVIP

连续纤维复合材料微观结构的量化分析 １６ 连 续 纤 维 复 合 材 料 微 观 结 构 的 量 化 分 析 维普资讯 doc. ２００６年 ７ 月 连 续 纤 维 复合 材 料 微 观 结构 的量 化 分 析 孙 江， 肖 琪 （１．嘉兴学 院机 电工程学 院 ，浙江 ３１４００１；２．嘉 ***工程训 练中心 ，浙江 ３１４００１） 摘要 ：本文通过对连续 纤维 复合材料微观 图像 的量化分析 ，获取其二阶 密度 函数 Ｋ（ｒ），并以 Ｋ（ｒ）为主要 函数建 立 目标 函 数 ，利用 ｍａｆ ｌａｂ的遗传优化 工具箱得 出代表体单元 。得 出的代表 体单元在力 学性 能上基 本反映 了原微观 结构 的力学性能。 关键词 ：连某某维复合材料 ；微观 结构；二 阶密度 函数 ；代表体单元 中图类 号 ：ＴＱ３２７．１ 文献标识码 ：Ａ 文章编 号 ：１００３一Ｏ９９９ （２００６）０４—００１６—０４ １ 前 言 连 续纤 维 与树脂 可 用不 同 的工艺 制备 各种 复合 材 料 制 品。这 些制 品 由于单 位质 量 的强度 及 弹性 模 量 较传 统 金属 材料 高 和一些 其 它特殊 性 能 而被广 泛 应用于航空航天、建筑、化工 、交通运输及体育等行 业 中。连某某维复合材料的研究包括制备工艺和性 能研 究 。对其 性 能的研 究必 然 涉及 到微 观结 构形 态 的分 析 。连续 纤维 复合 材 料微 观形 态一般 都 是 随机 分布的。以往对其性能的分析为了简化计算 ，将微 观结 构形 态视 为 纤维 按 一 定规 则 排 列 ，因 而选 择 的 代表体单元 （ＲＶＥ）为纤维按正方形 排列或按六角 形排列。这样 的代表体单 元可使力学模 型简化 ，但 却与 实 际 的微 观结 构形 态 有较 大 的差距 。 早在 ８Ｏ年代 ，Ｊａｍｅｓ，Ｔｏｒｑｕａｔｏ，Ｓｔｅｌｌ等 人 采 用数 值 图像 处理 技术 获取 复合 材 料 中异质 相 的空 间相 关 函数 ¨Ｊ。之后 ，Ｐｏｖｉｒｋ也采用数值 图像 处理技术 ， 通过指示 函数及对指示 函数的频谱 密度 的分析 ，找 出具有周期特性 的复合材料代 表体单元 Ｈ］。近年 来 ，Ｚｅｍａｎ采用 Ｒｅｌｙ等人提 出二 阶密度 函数 ，通过 遗传优化算法对两相复合材料进行分析处理 ，以得 到其 代 表体 单元 ，进 而应 用 代 表 体 单 元 分 析 纤 维 复 合材料 的性能Ｊ。本文介 绍 了微观 形态 的分 析方 法 ，利用二阶密度 函数对模拟二相复合材料作分析 ， 选 出其 代表 体单 元 ，并 对 其性 能作 分 析 。 ２ 复合材料微观结构 的数值 描述 对任一 由两相介质随机组成的复合材料微观图 像 ，文献 ［５］对 复合材 料微观结构采用指示 函数， （ｘ，Ｏｔ）和 ｎ点 概率 函数作 数值描述 。指示 函数， （ｘ，仅）定义 为 点 ｘ位 于某某 Ｏｔ的 ｒ相 区域 时为 １，否 则 为零 。用 数学 式 表示 即为 ： Ｘｒ（ｘ，仅）：ｆ ∈Ｄ，（ 。 （１） ｔ０ ∈ Ｄ ， （仅） 式 中 Ｄｒ（Ｏｔ）表示 在样 本 Ｏｔ中 ｒ相 占据 的 区域 。假 设 ｒ＝ｍ，其 中 ｍ为基体相 ，为增强相。例如 ，对于 具有圆形截面的纤维浸蘸树脂所形成 的复合材料 ， 以厂Ｆ标表示纤维 ，ｍ下标表示树脂 ，则指示 函数 （ｘ，Ｏｔ）和 （ｘ，Ｏｔ）有 下述 关 系 ： ｒ（ｘ，仅）＋ （ｘ，仅）＝ １。 （２） 以指 示 函数乘 积 的样本 平 均值 表示 ｎ点概率 为 Ｓ ．．，ｒ（ｘ１．．．，ｘ） ＝Ｘｒ，（ｘｌ，Ｏｔ）… （ｘ，Ｏｔ）。 （３） 由于 ｎ点概率 在 描述 复合 材 料微 观结 构时 非常 难操作 ，实际上仅使用它 的一 、二 阶函数 ，也 即在 戈 点 发现 ｒ相 的一 点 概率 函数 Ｓ，（ｘ） ＝ （ｘ，仅） （４） 在两 点 ｘ、ｘ点 同时发 现 ｒ相 、ｓ相 的两 点概 率 函数 Ｓ（ｘ，ｘ）＝ （，Ｏｔ）置 （ｘ，Ｏｔ）。 （５） 一 般 情 况 下 要 找 出这 些 指 示 函数 还 是 相 当 困 难 ，因此 常用 的处 理 办 法 是 将 材 料 看 成是 统 计 均 质 的。这样 就有 Ｓ，（ｘ）＝ Ｓ （６） Ｓ（ｘ，ｘ）＝ Ｓ（ｘ—ｘ） （７） 进一步再将介质看 成是 统计各 向同性的。这样 Ｓ （ｘ—ｘ）可 以简化 成 Ｓ（ｘ，ｘ） ＝Ｓ（１ｌｘ—ｘｌ１）。 （８） 对材 料作 出各 态 历 经 的假 设 时 ，可 用 ｒ相 的体 平 均 ｃ，值来代替一点概率函数 ｓ，，即 Ｓ，＝ｃ， （９） 下述 用 ｓ，和 Ｓ描述 一碳 纤维 ／环 氧树 脂复 合材 料 的 收稿 日期 ：２００５—１１－０