复合悬臂曲梁的有限元分析.docxVIP

复合悬臂曲梁的有限元分析 1 曲梁和复合曲梁的应力和切应力的计算公式不同 梁和复合梁是工程中常见的结构。例如，在道桥技术中的拱桥中，重量从地面上支撑的钩子可以被视为梁。主轴承的轴承和汽车板的弹簧可以被视为由多层金属板重叠形成的复合梁。与直梁和复合直梁相比较,由于曲梁曲率的影响,使得曲梁和复合曲梁的内力和变形计算比直梁和复合直梁要复杂得多。在材料力学教材中对平面曲杆纯弯曲时的正应力计算公式已经有了介绍,另外,在《高等材料力学》中对具有对称截面的曲梁应力计算也有介绍,但这些计算公式都是针对单一材料的曲梁。而文献对于一般截面形状的曲梁和复合曲梁在复杂受力情况下的正应力和切应力的计算公式进行了推导,但到目前为止还没有相应的实验来印证这些公式。为此本文建立了有限元模型,运用ANSYS软件对复合悬臂曲梁的正应力和切应力进行了简化计算,并设计了复合悬臂曲梁的实验模型,将三者结果相互对比,对复合曲梁理论公式的正确性作了进一步验证。 2 复合曲线理论的计算 2.1 复合悬臂曲梁的材料参数 文献中给出了一般截面形状的复合曲梁(见图1)在复杂受力情况下的正应力计算公式,即 σ1=E1Η(1-y/R){[(eg-f2)Νs-(bg-df)Μz+(bf-de)Μy]-[(bg-df)Νs-(ag-d2)Μz+(af-bd)Μy]y+[(bf+de)Νs-(af-bd)Μz+(ae-b2)Μy]z}σ2=E2Η(1-y/R){[(eg-f2)Νs-(bg-df)Μz+(bf-de)Μy]-[(bg-df)Νs-(ag-d2)Μz+(af-bd)Μy]y+[(bf+de)Νs-(af-bd)Μz+(ae-b2)Μy]z}σ1=E1H(1?y/R){[(eg?f2)Ns?(bg?df)Mz+(bf?de)My]?[(bg?df)Ns?(ag?d2)Mz+(af?bd)My]y+[(bf+de)Ns?(af?bd)Mz+(ae?b2)My]z}σ2=E2H(1?y/R){[(eg?f2)Ns?(bg?df)Mz+(bf?de)My]?[(bg?df)Ns?(ag?d2)Mz+(af?bd)My]y+[(bf+de)Ns?(af?bd)Mz+(ae?b2)My]z} 式中: a=E1(A1+1R∫A1ydA1+1R2Jz1)+E2(A2+1R∫A2ydA2+1R2Jz2)b=E1(∫A1ydA1+1RJz1)+E2(∫A2ydA2+1RJz2)d=E1(∫A1zdA1+1RJyz1)+E2(∫A2zdA2+1RJyz2)e=E1Jz1+E2Jz2f=E1Jyz1+E2Jyz2g=E1Jy1+E2Jy2Η=(aeg-af2-b2g+2bdf-d2e)Jy1=∫A1z21-y/RdA1?Jyz1=∫A1yz1-y/RdA1Jz1=∫A1y21-y/RdA1?Jy2=∫A1z21-y/RdA1Jyz2=∫A1yz1-y/RdA1?Jz2=∫A1y21-y/RdA1 E1、E2为材料Ⅰ和Ⅱ的弹性模量;A1、A2为材料Ⅰ和Ⅱ的截面面积;σ1、σ2为材料Ⅰ和Ⅱ上各点处的正应力;R为复合曲梁的曲率半径。 图中s为沿曲梁轴线量取的弧长;y为指向该梁轴线曲率中心的径向坐标;z为形心轴,正交于杆件平面。py(s)pz(s)为沿y和z方向单位弧长上的外力分量;Ns、Vy、Vz、My、Mz为任意截面上的轴向力、剪力和弯矩;G为组合截面的形心;H1、H2为材料Ⅰ和Ⅱ的截面高度;h1、h2为形心轴z到上、下边缘的距离;?1(y)、?2(y)为横截面上的积分区域。 对于具有直角梯形截面和1/4圆弧轴线的复合悬臂曲梁(见图2),当自由端受集中力P=500 N作用时在45°截面处的内力为Vz=0;Vy=P×cos 45°=353.55 N;NS=P×sin 45°=353.55 N;My=0;Mz=P×R×sin 45°=55.628 MN·m。 2.2 材料的切应力 根据文献,一般截面形状的复合曲梁在复杂受力情况下的切应力计算公式如下: (h2-H2≤y≤h2) 式中: α1=1R(eg-f2)-(bg-df)+μ(bf-de)α2=1R(bg-df)-(ag-d2)+μ(af-bd)α3=1R(bf-de)-(af-bd)+μ(ae-b2)μ=VzVy?V=Vy b′为距z轴为y处的截面宽度;τ1、τ2为材料Ⅰ和Ⅱ上各点处的切应力。 由于复合悬臂曲梁正应力和切应力的计算公式比较复杂,这里应用matlab软件进行编程计算。 3 局部细化单元格 为了用有限元法分析直角梯形截面复合悬臂曲梁,利用三维实体单元SOLID92将梁的上下层分别划分成多个实体单元,再在上下层结合部位将单元格进行局部细化。图3表示利用有限元程序计算三种悬臂曲梁(纯钢、铝钢组合和铜钢组合曲梁)在自由端受集中力P=500 N时