功能梯度轴对称压电结构的非轴对称有限单元.pdf

?Ⅸ『Ⅺ
, .. . .,. 、也 , , : 功能梯度轴对称压电结构的非轴对称有限单元 邵斌磊徐荣桥 浙江大学土木工程系 在工程上有些轴对称的物体,承受的载荷并不是轴对称的,这类问题虽不属于空间轴对称问题,但可把载荷在 圆周方向展成傅里叶级数,利用物体在几何上的对称性,把一个三维问题转化为一组二维问题求解。本文给出 了功能梯度横观各向同性的轴对称压电结构的非轴对称有限元列式,并用编制了程序,计算了圆板在 非轴对称集中载荷作用下的静力问题和压电圆板的非轴对称的自由振动,并与的结果进行了比较。 关键词:压电材料,轴对称结构,非轴对称问题,有限元 果,并与的三维单元结果进行了比较,结 .引言 果表明本文方法具有较好计算效率和精度。 压电材料由于其机电耦合特性,被广泛地应用于各 .载荷和位移的展开 个领域,例如在噪声控制、结构形状控制、振动控 我们先把作用在轴对称物体上的任意非轴对称载荷 制及振动冲击测量等领域【】。但也由于压电材料 控制方程及实际结构的复杂性,进行解析的分析变 沿柱面坐标轴,,只分解为三个分量,.,,, 得十分困难,所以在大多数情况下都有赖于数值方 ,,和,,,然后把这三个载荷分量在 法,而有限元法是其中应用最广的。 坐标方向展成傅里叶级数如下: 早在年, 】就给出了压 电材料的有限元列式。 【】利用有 ,,磁,∑磁, , 。 限元研究了功能梯度板的问题。 ∑露,., 等】基于梁理论和一阶剪 切变形理论,提出了包含压电材料质量和刚度影响 的梁、板单元,并采用等控制方法实现了结构 ,,碟¨∑露,, ? 的振动控制。丁皓江等】应用有限元方法研究了横 ∑碟,., 观各向同性压电材料的轴对称问题。 】 , 把它们推广到了非均匀的横观各向同性压电材料的 轴对称问题。但是对于轴对称结构的非轴对称问题 ,只或懈∑露, , . 、, 的单元还未见报道。 本文采用结点等参数单元,推导了非均匀轴 ∑巧,., 对称压电结构的非轴对称问题的有限元列式,并用 ,,, 编制了相应的程序。最后给出了数值结 . ??轮.//. 上式中等式右端第一项磷,/和礞都与 无关,但磁与磷为轴对称载荷,而碟为扭转载 ∽ 荷,它使物体产生扭转变形,其余各项都与有 关,为非轴对称载荷。而右端第一个级数使物体产 ,,萎笔爱耄三奄 生对称于平面的变形,我们简称这部分载荷为 对称载荷,而第二个级数使物体产生反对称于平面 的变形,称为反对称载荷。 与载荷相对应,物体上任一点的位移分量“, ,,和电势分量矽也可展成傅里叶级数: ,,瑶’∑矿,., 甘睡卜, 扛 。 ∑∥,, , ,,订懈∑,., 闰 二. 芝:∥,., 百 . ,.,,,,∑,, 扛 。 ∑订,., , 畎,.,幺菇,.,∑群,., 忙 二 ∑刃,., 式中,碚,培,菇是轴对称位移分量,订是 环向扭转位移分量,它们都与目无关;其余八个总 和项的系数矿,∥,谚,,,订,群, 刃都是,,的函数,通常式右端的第一个 级数是与对称载荷相应的与平面对称的位移分 量,式右端的第二个级数是与反对称载荷相 应的与平面口反对称的位移分量。 旧可以证明,式右端的各项一一对应于式 右端的各项。于是我们可以按照空间轴对称 。 【】 。坩三。坩。。 问题的方法对式卜中的每一项载荷求解,然后 【. 应用叠加原理得到非轴对称载荷情况下的解,使问 题得到极大的简化。 从式和式可以看出,对称载荷和 非对称载荷分别作用下的载荷函数以及相应的位移 函数,实际上只是把函数中的余弦项换为正弦项, 或把正弦项换为余弦项,所以,关于对称载荷或非 对称载荷作用下问题的求解方式是相同的。这里将 以对称载荷作用为例进行有限元公式的推导和 程序的实现。 设对称载荷的分量最,弓,巴,用单 项三角函数表示为: 式中, .单元质量矩阵 单元质量矩阵的推导过程与单元刚度的类似,但是 陋】 要注意对应于电势自由度的惯性力为零,因此其表 达式为: 塑。 。“。矿, 盟 堕螂坩 式中,Ⅳ挚ⅣⅣ 塑 一盟 掣 掣Ⅳ Ⅳ一堕警一乌Ⅳ 言咖旧 警坩 同样也是坐标的函数,我们采用与相同的 ??? 处理方法,即 一盟..? 夕∑房 警矧钊 言伽加 式中,珐为第个节点上的材料密度。 由于式中的Ⅳ存在变量,同样需要先对 .单元刚度矩阵 进行积分以求得显式的表达式,令 根据虚功原理,可以得到线性压电体的单元刚度矩 ‘矿 阵为: 同样可以使用的符号运算功能求得 的显式表达式。 置‘
。~ 式中,.,为用于坐标转换的雅克比矩阵,为线性 最后,使用高斯积分法求得单元质量矩阵为: 压电体的材料特性矩阵。对于非均匀材料,它是坐 、 ‘∑∑形%∥川 』』‘、‘, , 标的函数。在这里,我们也把它用结点的形函数进 行插值,即: .数值算例 ∑乎,叩 .算例 式中,口为第个节