γ型框架结构之最小挫曲载重-淡江大学机械与机电工程学系.docVIP

機械與機電系畢業專題報告 Γ型框架結構之最小挫曲載重 指導教授：劉昭華 淡江大學機械與機電工程學系教授 專題人員：范峻榮、林俊宏 淡江大學機械與機電工程學系大學生 摘要： 本專題是針對Γ型框架之挫曲，來探討最小挫曲的載重；其方法是將Γ型框架簡化之後，推導出其最小挫曲載重之非線性方程式，再利用數值分析（bisection method），求出挫曲載重之數值。研究動機 框架（frame）是土木建築經常使用的結構，例如室內停車場的樓房或是商業大樓。框架結構中直立的桿件（柱）通常承受軸向應力，而此軸向應力可能造成桿件挫曲（buckling）。而桿件的材質、粗細、長度和形狀都可能影響其所能承受之最大軸向應力，本專題是針對最簡單之Γ型框架結構討最小的挫曲載重。 貳.Γ型框架結構 2.1式 Γ型框架由構件AB和BC組成，並B點，且A和C以插銷。AB承受壓力P，且忽略橫桿BC產生之變形BC桿可以簡化為旋轉彈簧，如圖二所示( [1], 818頁)。 圖三顯示此結構之自由體圖，力平衡 (1) 其中H為支撐點反作用力，若旋轉彈簧之彈性係數為(R，由上式可知 (2) 以下是推導出柱位移曲線方程式[1] (3) 上式整理後可 (4) 其中 (5) (4)式之通解 (6) 由於A點為插銷端，v(0)為0，將此條件代入通解後可得。因不考慮橫桿BC之變形，故v(L)為0，代入通解後得 由於在B點柱子之旋轉角為， 將這些條件代入(6)式後， (7) 將上式中的消代入整理後即可得 (8) 2-2數值分析 βL，若是=3，我們畫出的圖形如圖四所示，由此圖可知函數在x=nπ (n=1,2,3,…) 時會趨近於無窮大，因此我們每次只尋找nπ與(n+1)π。。 kL kL kL 0.5 3.2860 3.5 3.7810 6.5 3.9991 1 3.4056 4 3.8289 7 4.0236 1.5 3.5059 4.5 3.8711 7.5 4.0459 2 3.5909 5 3.9086 8 4.0662 2.5 3.6636 5.5 3.9420 8.5 4.0848 3 3.7264 6 3.9720 9 4.1020 表一顯示在不同係數時所尋找到最小的根，其圖形如圖五所示 2.3結果與討論 由圖五(或表一)的結果可知若其他參數(E, I, L)皆為固定，係數βR愈高，則最小挫曲載重數值亦隨之增高，這是因為當旋轉彈簧之彈性係數愈大時，在B點此直立柱愈難旋轉，這使住子愈不容易產生側向位移，而此側向位移正是在挫曲狀態下之位移，因此β數值增加使得直立柱愈不容易挫曲，故需要更大的壓力P。Γ型框架結構之挫曲載重A點是其他形式支撐的狀況。參考文獻： Gere, J. M., Mechanics of Materials, 6th ed., Brooks/Cole, 2004. [2] Curtis F. Gerald、Patrick O. Wheatly：Applied Numerical Analysis7th ed. 附錄：主程式 clc,clear k=1 for i=0.1:0.1:10 j=1; U0=0.001; Umax=pi-0.001; c=i; [Uhalf,fhalf]=bisection(U0,Umax,c) ; while abs(fhalf)0.001 j=j+1; U0=(j-1)*pi+0.001; Umax=j*pi-0.001; [Uhalf,fhalf]=bisection(U0,Umax,c); end disp([c Uhalf]); x(k)=c; y(k)=Uhalf; k=k+1; end plot(x,y) xlabel(x=(Br*L)/(EI)) ylabel(y=kL) 二分法之函數程式 function [Uhalf,fhalf]=bisection(U0,Umax,c) syms u f=inline(c*(u*cot(u)-1)-u^2); for i=1:50; Uhalf=(U0+Umax)/2; fhalf=f(Uhalf); f0=f(U0); if fhalf==0 break lseif (fhal