非匀质材料弹性模量预测的数值方法.pdfVIP

第 52卷第 3期 甘 肃 水 利 水 电 技 术 Vo1．52．No．3 2016年 3月 GANsuWATERRES0URCESANDHYDROPOWERTECHNOL0GY Mar．，2016 · 设计与研究 · 非匀质材料弹性模量预测的数值方法 顾 政 (中国美术学院，浙江 杭州 310000) 摘要：工程材料大多属于非匀质材料，因此在求解弹性模量的过程 中首先要解决材料非匀质性的实现问题。快速傅 里叶变换法利用多个像素点反应材料 内部结构与性质，再与迭代法相结合可以预测了非匀质材料的弹性模量。以非 匀质杆件和 多孔材料为例 ，通过数值解与解析解比较 ，初步验证 了该方法的有效性。 关键词：非匀质材料；弹性模量；快速傅里叶变换；迭代法 中图分类号：TB383．4 文献标志码：A 文章编号：2095一o144(2016)03—0015—03 1 前 言 的像素点上获得结构单元的代表性信息l8]。以非匀 工程材料大多属于非匀质材料 ，这些材料广泛 质杆件为例进行模拟 ，首先为了获得各段弹性模量 地应用于各个工程领域 ，它们不仅具有多种优异 的 E1，E2，E3，…，EN的杆件 的图像，在图像区域等距 性能 ．而且制造工艺简单。通过研究其力学性能，有 离选取像素点，根据快速傅里叶变换原理 ，每个一维 助于我们更好地 了解材料非匀质特性 ，为新型材料 图像都包含 2个像素点，这些像素点的力学性质相 设计开发提供理论依据。而在求解弹性模量的过程 互独立 ，计算步骤可以归纳如下 ： 中，首先要解决材料非均质性 的实现 问题 。快速傅 (1)将杆件划分为N=2段 ，如图 1所示 ，这里 里叶变换法通过利用图像选取多个像素点 ，反映了 的K为像素点个数 的指数 。取每段的中点为像素 材料内部结构与性质 。这个方法与有限元法 j相 比， 点 用坐标可表示为P／=(i一0．5)L／2，i的取值从 l 避免 了网格划分 的困难 ，同时也节省 了总刚度矩 阵 至0J7V 求解需要耗费的大量时间。 快速傅立 叶变换不是一种新型的变换算法 ，它 是为了能够快速实现离散傅里叶变换而进行 的改 图1 像素点的选取 进 ，被广泛应用于各个领域 。傅立叶变换最早 由法 国数学家和物理学家傅里叶于 1807年提出，之后 ， (2)对每个像素点进行判断，如果落在第 l段 ， 为了扩大傅里叶变换的使用范围，傅里叶变换的表 则弹性模量赋为 E1；如果落在第 2段 ，则弹性模量 达也派生 出了多种形式。离散博里叶变换就是最常 赋为 ￡2：同理可得像素点落在其他段时的弹性模 见的傅里叶变换形式之一 3¨]，为 了加快运算 ，Cooley 量 。 和Tukey在论文 《机器计算傅里叶级数的一种算法》 由于各段弹性模量均不同，取各点弹性模量 E 中提 出了快速傅里叶变换思路 ，后经过各 国学者 的 (z)是坐标 的函数 ，其应力～应变关系可表示成 ： 研究完善 ，形成了高效 的快速傅里叶变换 ，简称为 o-=E~e+r() (1) FFT…。根据抽取方式主要分为两大类 ：按时间抽取 式 (1)中， 为弹性模量常数 ，丁()定义为： 和按频率抽取法 。前者是将时域信号序列按偶奇分 丁()：(E()一Eo)s (2) 排 ．后者是将频域信号序列按前后对半分排 。根据 式 (1)的解可表示为： 分