層合板屈曲問(wèn)題分析.doc

《層合板屈曲問(wèn)題分析.doc》由會(huì)員上傳分享，免費(fèi)在線閱讀，更多相關(guān)內(nèi)容在工程資料-天天文庫(kù)。

1、層合板屈曲問(wèn)題分析復(fù)合材料力學(xué)一一屈曲問(wèn)題分析一、問(wèn)題的描述四邊簡(jiǎn)支的正交對(duì)稱(chēng)矩形層合板，單層厚度為0.2mm,a=800mm,b二100mm。已知各單層特性：受單向壓縮El?181GPa/E2?10.3GPa/G12?7.17GPa/?21?0.28求：臨界載荷[0/90/90/0]二、解析解1、理論分析正交對(duì)稱(chēng)層合板單向受壓的屈曲方程：?4w?4w?4w?2wD114?2?D12?2D66?22?D224+N2?0?x?x?y?y?x由Navier法設(shè)屈曲形狀為雙正弦函數(shù)：w???amnsinm?ln?l??m?xn?ysinab將w代入屈曲方程求得臨

2、界屈曲荷載為：??Ncr?2?Dll?mb/a??2?D12?2D66??D222?b?mb/a????2?2?l由上式可知N取最小值的x方向半波數(shù)m與邊長(zhǎng)比b/a及剛度有關(guān)。2、matlab編程求解El=181;v21=0.28;vl2=E2*v21/El;G12=7.17;%材料常數(shù)Qll=El/(l-vl2*v21);Q22=E2/(l-vl2*v21);Q12=E2*v21/(l-vl2*v21);Q66=G12;復(fù)合材料力學(xué)一一屈曲問(wèn)題分析%正軸剛度U1_Q=(1/8)*(3*Q11+3*Q22+2*Q12+4*Q66);U2_Q=(1/2)*(

3、Q11-Q22);U3__Q=(1/8)*(Q11+Q22-2*Q12-4*Q66);U4_Q=(1/8)*(Q11+Q22+6*Q12-4*Q66);U5_Q=(1/8)*(Q11+Q22-2*Q12+4*Q66);%單向板正軸剛度的線性組合z0=-0.4;zl=-0.2;z2=0;z3=0.2;%層合板厚度方向的坐標(biāo)thetal=O;theta2=pi/2;theta3=pi/2;theta4=0;%每層的鋪設(shè)角h=0.8;%層合板的總厚度Vl_D=(l/3)*(((zl)A3-(z0)A3)*cos(2*thetal)+((z2)A3-(zl)A3)

4、*cos(2*theta2)+((z3)A3-(z2)A3)*cos(2*theta3)+((z4)A3-(z3)A3)*cos(2*theta4));V2_D=(l/3)*(((zl)A3-(z0)A3)*cos(4*thetal)+((z2)A3-(zl)A3)*cos(4*theta2)+((z3)A3-(z2)A3)*cos(4*theta3)+((z4)A3-(z3)A3)*cos(4*theta4));V3_D=0;V4_D=0;%層合板的幾何因子Dll=Ul_Q*hA3/12+Vl_D*U2_Q+V2_D*U3_Q;D22=Ul_Q*hA3/

5、12-Vl_D*U2__Q+V2_D*U3__Q;?2?復(fù)合材料力學(xué)一一屈曲問(wèn)題分析D66=U5_Q*hA3/12-V2_D*U3_Q;%彎曲剛度a=0.8;b=0.1;%層合板的邊長(zhǎng)m=l;Nl=piA2/bA2*(Dll*(m*b/a)A2+2*(D12+2*D66)+D22/(m*b/a)A2);m=2;N2=piA2/bA2*(Dll*(m*b/a)A2+2*(D12+2*D66)+D22/(m*b/a)A2);m=3;N3=piA2/bA2*(Dll*(m*b/a)A2+2*(D12+2*D66)+D22/(m*b/a)A2);m=4;N4=pi

6、A2/bA2*(Dll*(m*b/a)A2+2*(D12+2*D66)+D22/(m*b/a)A2);m=5;N5=piA2/bA2*(Dll*(m*b/a)A2+2*(D12+2*D66)+D22/(m*b/a)A2);m=6;N6=piA2/bA2*(Dll*(m*b/a)A2+2*(D12+2*D66)+D22/(m*b/a)A2);m=7;N7=piA2/bA2*(Dll*(m*b/a)A2+2*(D12+2*D66)+D22/(m*b/a)A2);M=[l:7]N=[N1N2N3N4N5N6N7]為半波數(shù)，N為臨界荷載結(jié)果：M=1234567N=

7、1.0e+004*8.72042.32861.19180.84930.75160.76300.8370由N可知,當(dāng)半波數(shù)m=5時(shí),最小臨界載荷為7516No三、用ANSYS求解數(shù)值解網(wǎng)格劃分為：160*20邊界條件為：四邊簡(jiǎn)支，約束分別為：x、z方向位移，y、z方向位移,z?3?復(fù)合材料力學(xué)一一屈曲問(wèn)題分析方向位移，y、z方向位移數(shù)值解最小臨界載荷：7181N?4?復(fù)合材料力學(xué)一一屈rtti問(wèn)題分析四、誤差分析7516?7181?100%=4.46%,解析解與數(shù)值解的誤差在5%以內(nèi)，正確！7516