发明名称 宏纤维压电复合材料压电片动力学的作动力识别方法
摘要 一种宏纤维压电复合材料压电片动力学的作动力识别方法,通过复合材料的经典层合板理论,构造了宏纤维复合材料的本构方程,利用材料力学的相关理论计算出作动力的大小,能够直观、简便的将宏纤维压电复合材料压电片的实际作用等效为压电片四个边节点载荷,方便宏纤维压电复合材料压电片的有限元动力学分析。本发明在不含压电材料单元的有限元软件中进行压电材料建模,应用于有限元软件的压电单元建模中,弥补了部分有限元方法不能进行动力学的作动力识别的不足。
申请公布号 CN106156441A 申请公布日期 2016.11.23
申请号 CN201610589044.1 申请日期 2016.07.25
申请人 西北工业大学 发明人 李斌;牛文超;杨阳
分类号 G06F17/50(2006.01)I 主分类号 G06F17/50(2006.01)I
代理机构 西北工业大学专利中心 61204 代理人 慕安荣
主权项 宏纤维压电复合材料压电片动力学的作动力识别方法,其特征在于,具体过程是:步骤1:截取压电片的分析单元;在悬臂梁的表面粘贴宏纤维压电复合材料压电片,并使该宏纤维压电复合材料压电片靠近悬臂梁的固支端一侧;将所述悬臂梁上粘贴宏纤维压电复合材料压电片的部位截取下来,得到分析单元步骤2:划分有限元单元对所述的粘贴宏纤维压电复合材料压电片的悬臂梁划分有限元单元,并清点宏纤维压电复合材料压电片长边的有限元节点数目与该宏纤维压电复合材料压电片短边的有限元节点数目;步骤3:确定宏纤维压电复合材料压电片主方向应力与应变关系;步骤4:确定宏纤维压电复合材料压电片的刚度矩阵;步骤5,确定单位电压下的宏纤维压电复合材料压电片四个边的力与力矩;步骤6,将单位电压下的宏纤维压电复合材料压电片四个边所受的力与力矩等效为节点力:所述的节点力由公式(8)计算的x、y、z三方向的力与x、y、z三方向的力矩组成;<maths num="0001"><math><![CDATA[<mrow><mtable><mtr><mtd><mrow><mi>N</mi><mi>p</mi><mi>w</mi><mo>=</mo><mo>&lsqb;</mo><mtable><mtr><mtd><mrow><mi>N</mi><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow><mo>*</mo><mi>w</mi><mi>p</mi><mo>/</mo><mi>a</mi><mo>/</mo><mn>2</mn><mi>N</mi><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow><mo>*</mo><mi>w</mi><mi>p</mi><mo>/</mo><mi>a</mi><mo>/</mo><mn>2</mn></mrow></mtd><mtd><mn>0</mn></mtd></mtr></mtable><mo>&rsqb;</mo></mrow></mtd></mtr><mtr><mtd><mrow><mi>M</mi><mi>p</mi><mi>w</mi><mo>=</mo><mo>&lsqb;</mo><mtable><mtr><mtd><mn>0</mn></mtd><mtd><mrow><mi>M</mi><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow><mo>*</mo><mi>w</mi><mi>p</mi><mo>/</mo><mi>a</mi><mo>/</mo><mn>2</mn><mo>-</mo><mi>M</mi><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow><mo>*</mo><mi>w</mi><mi>p</mi><mo>/</mo><mi>a</mi><mo>/</mo><mn>2</mn></mrow></mtd></mtr></mtable><mo>&rsqb;</mo></mrow></mtd></mtr><mtr><mtd><mrow><mi>N</mi><mi>p</mi><mi>l</mi><mo>=</mo><mo>&lsqb;</mo><mtable><mtr><mtd><mrow><mi>N</mi><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow><mo>*</mo><mi>l</mi><mi>p</mi><mo>/</mo><mi>b</mi><mo>/</mo><mn>2</mn><mi>N</mi><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow><mo>*</mo><mi>l</mi><mi>p</mi><mo>/</mo><mi>b</mi><mo>/</mo><mn>2</mn></mrow></mtd><mtd><mn>0</mn></mtd></mtr></mtable><mo>&rsqb;</mo></mrow></mtd></mtr><mtr><mtd><mrow><mi>M</mi><mi>p</mi><mi>l</mi><mo>=</mo><mo>&lsqb;</mo><mtable><mtr><mtd><mrow><mo>-</mo><mi>M</mi><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow><mo>*</mo><mi>l</mi><mi>p</mi><mo>/</mo><mi>b</mi><mo>/</mo><mn>2</mn><mi>M</mi><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow><mo>*</mo><mi>l</mi><mi>p</mi><mo>/</mo><mi>b</mi><mo>/</mo><mn>2</mn></mrow></mtd><mtd><mn>0</mn></mtd></mtr></mtable><mo>&rsqb;</mo></mrow></mtd></mtr></mtable><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>8</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001058130840000011.GIF" wi="982" he="287" /></maths>式中:Npw为宏纤维压电复合材料压电片宽边所受到的x、y、z三方向的力矩阵;Npl为宏纤维压电复合材料压电片长边所受到的x、y、z三方向的力矩阵;Mpw为宏纤维压电复合材料压电片宽边所受到的x、y、z三方向的力矩矩阵;Mpl为宏纤维压电复合材料压电片长边所受到的x、y、z三方向的力矩矩阵;wp为宏纤维压电复合材料压电片的宽度;a为宽度方向的节点数;lp为宏纤维压电复合材料压电片的长度;b为长度方向的节点数;步骤7:将对节点力进行坐标转换;对得到的等效的节点力进行坐标转换;坐标转化后的节点力为:<maths num="0002"><math><![CDATA[<mrow><mi>F</mi><mo>=</mo><mfenced open = "[" close = "]"><mtable><mtr><mtd><mo>&lsqb;</mo><mtable><mtr><mtd><mrow><mi>N</mi><mi>p</mi><mi>w</mi></mrow></mtd><mtd><mrow><mi>M</mi><mi>p</mi><mi>w</mi></mrow></mtd></mtr></mtable><mo>&rsqb;</mo><mo>*</mo><mi>T</mi></mtd></mtr><mtr><mtd><mo>&lsqb;</mo><mtable><mtr><mtd><mrow><mi>N</mi><mi>p</mi><mi>l</mi></mrow></mtd><mtd><mrow><mi>M</mi><mi>p</mi><mi>l</mi></mrow></mtd></mtr></mtable><mo>&rsqb;</mo><mo>*</mo><mi>T</mi></mtd></mtr><mtr><mtd><mo>-</mo><mo>&lsqb;</mo><mtable><mtr><mtd><mrow><mi>N</mi><mi>p</mi><mi>l</mi></mrow></mtd><mtd><mrow><mi>M</mi><mi>p</mi><mi>l</mi></mrow></mtd></mtr></mtable><mo>&rsqb;</mo><mo>*</mo><mi>T</mi></mtd></mtr><mtr><mtd><mo>-</mo><mo>&lsqb;</mo><mtable><mtr><mtd><mrow><mi>N</mi><mi>p</mi><mi>w</mi></mrow></mtd><mtd><mrow><mi>M</mi><mi>p</mi><mi>w</mi></mrow></mtd></mtr></mtable><mo>&rsqb;</mo><mo>*</mo><mi>T</mi></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>11</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001058130840000021.GIF" wi="685" he="303" /></maths>步骤8,悬臂梁动力学试验;步骤9,修正粘贴宏纤维压电复合材料压电片的悬臂梁有限元模型的模态阻尼系数:使通过悬臂梁动力学试验得到的设置的各幅值电压的中间值的试验值与该设置的各幅值电压的仿真值的误差在1%内;步骤10,任一幅值电压宏纤维压电复合材料压电片动力学的压电效应常数识别:通过公式(12)对所述的各幅值电压作用下的宏纤维压电复合材料压电片的压电效应常数d<sub>a</sub>进行修正,压电效应常数d的下标a表示信号发射器输出信号为交流信号;<maths num="0003"><math><![CDATA[<mrow><msub><mi>d</mi><mi>a</mi></msub><mo>=</mo><mfrac><msub><mi>w</mi><mn>3</mn></msub><msub><mi>w</mi><mn>4</mn></msub></mfrac><mo>&times;</mo><msub><mi>d</mi><mn>33</mn></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>12</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001058130840000022.GIF" wi="422" he="118" /></maths>w<sub>3</sub>为试验测得的所述幅值电压中间值的悬臂梁自由端位移幅值,w<sub>4</sub>为仿真计算的该幅值电压中间值的悬臂梁自由端位移幅值;对得到步骤九所述的已设置的各幅值电压作用下的宏纤维压电复合材料压电片的压电效应常数d<sub>a</sub>进行曲线拟合,得到任一幅值电压下宏纤维压电复合材料压电片动力学的压电效应常数d<sub>an</sub>;步骤11,任一幅值电压下的宏纤维压电复合材料压电片动力学的作动力识别:将得到的任一幅值电压下的宏纤维压电复合材料压电片的压电效应常数d<sub>an</sub>代入式(6)、式(7)中,并重复步骤5~7,获得任一幅值电压下宏纤维压电复合材料压电片动态作动力。
地址 710072 陕西省西安市友谊西路127号