孫云

摘要：在《結(jié)構(gòu)力學(xué)》課程教學(xué)中，三鉸拱的內(nèi)力計(jì)算是個(gè)難點(diǎn)。本文應(yīng)用MATLAB編制了通用程序，可以解決三鉸拱在復(fù)雜荷載作用下截面內(nèi)力的計(jì)算問題。

關(guān)鍵詞：結(jié)構(gòu)力學(xué)；MATLAB；三鉸拱；內(nèi)力

中圖分類號(hào)：G642.0？搖 文獻(xiàn)標(biāo)志碼：A 文章編號(hào)：1674-9324（2014）05-0096-02

三鉸拱是《結(jié)構(gòu)力學(xué)》課程其中的一個(gè)章節(jié)，由于拱的截面法線與水平軸的夾角φ隨截面位置發(fā)生改變，截面內(nèi)力又與φ有關(guān)，因而求解三鉸拱，尤其在復(fù)雜荷載作用下截面的內(nèi)力是件煩瑣的事情。為了讓師生能從煩瑣、重復(fù)的手工計(jì)算的勞動(dòng)中解放出來，把更多的時(shí)間用于基本概念、基本理論的理解及問題的解決方法上[1]，本文利用MATLAB軟件的符號(hào)運(yùn)算、數(shù)值計(jì)算、編程等功能，編制了三鉸拱在復(fù)雜荷載作用下截面內(nèi)力計(jì)算的通用程序，不僅提高了教學(xué)效率，也激發(fā)了學(xué)生學(xué)習(xí)《結(jié)構(gòu)力學(xué)》和MATLAB的極大興趣。

三鉸平拱在只有豎向荷載作用下任一截面的內(nèi)力與同跨度、同荷載的對(duì)應(yīng)簡支梁對(duì)應(yīng)截面的內(nèi)力有如下關(guān)系[2]：

M=M0-FHy，F(xiàn)s=F0scosφ-FHsinφ，F(xiàn)N=-F0Ssinφ-FHcosφ （1）

因此，先介紹簡支梁任一截面的內(nèi)力計(jì)算，再介紹三鉸拱任一截面的內(nèi)力計(jì)算。

一、簡支梁的內(nèi)力

若簡支梁[3]上作用有均布荷載qi（i=1，nq），集中荷載 Fpi（i=1，nF）和集中力偶Mi（i=1，nM）。均布荷載的左、右端點(diǎn)到支座A的距離分別為ai（i=1，nq）、bi（i=1，nq），集中荷載和集中力偶到支座A的距離分別為ci（i=1，nF）、di（i=1，nM）（見圖1）。設(shè)支座反力FAy、FBy已經(jīng)由平衡方程求出。

求任一截面x的內(nèi)力。根據(jù)截面法，取x截面以左部分為隔離體，考慮隔離體上的荷載是任意、復(fù)雜的，可能有均布荷載，可能有集中荷載，也有可能有集中力偶，且荷載數(shù)目可能是多個(gè)。小變形條件下，可用疊加法計(jì)算內(nèi)力[4]。根據(jù)截面上的剪力和彎矩分別等于隔離體上所有外力沿截面切向的投影代數(shù)和以及對(duì)所求截面力矩代數(shù)和的計(jì)算法則，分別考慮隔離體上支座反力FAy、均布荷載qi、集中荷載Fpi和集中力偶Mi等對(duì)x截面內(nèi)力的影響。

（一）反力FAy的影響

M1=FAy·x，F(xiàn)S1=FAy （2）

（二）均布荷載qi的影響

均布荷載與所求截面的位置有如下3種情況：（1）均布荷載qi全部作用在隔離體上（圖2a），此時(shí)考慮全部均布荷載的作用；（2）均布荷載qi部分作用在隔離體上，部分在隔離體外（圖2b），此時(shí)只考慮隔離體上部分荷載的作用，則 bi=x；（3）均布荷載qi全部作用在隔離體外（圖2c），則不考慮均布荷載的作用，即qi=0，ai=bi=0。

nq個(gè)均布荷載產(chǎn)生的x截面內(nèi)力為：M2=-■qi（bi-ai）xi-■，F(xiàn)S2=-■qi（bi-ai） （3）

對(duì)于不同的均布荷載情況，公式中的qi、ai、bi取不同值。

（三）集中荷載Fpi的影響

集中荷載與所求截面的位置有如下2種情況：（1）集中荷載Fpi作用在隔離體上，應(yīng)考慮Fpi的作用；（2）集中荷載 Fpi作用在隔離體外，不考慮Fpi的作用，則Fpi=0，ci=0。

nF個(gè)集中荷載產(chǎn)生的x截面內(nèi)力為：M3=-■Fpi（xi=ci）， FS3=-■Fpi （4）

對(duì)于不同的集中荷載情況，F(xiàn)pi、ci取不同值。

（四）集中力偶Mi的影響

集中力偶與所求截面的位置有兩種情況：（1）集中力偶Mi作用在隔離體上，應(yīng)考慮Mi的作用，則di

nM個(gè)集中力偶產(chǎn)生的x截面內(nèi)力為：M4=-■Mi，F(xiàn)S4=0（5）

對(duì)于不同的集中力偶情況，Mi取不同值。

則反力FAy、均布荷載qi、集中荷載FPi和集中力偶Mi共同作用下x截面內(nèi)力為：

M=M1+M2+M3+M4，F(xiàn)s=Fs1+Fs2+Fs3+Fs4 （6）

利用MATLAB編制通用的程序，只需分別將均布荷載、集中荷載和集中力偶的個(gè)數(shù)：nq、nF、nM、大小qi、Fpi、Mi，以及分別到左端支座A的距離ai、bi、ci、di，梁長L以及截面位置x等輸入，就可以利用程序快速求解出復(fù)雜荷載作用下簡支梁任一截面的內(nèi)力。數(shù)據(jù)輸入時(shí)，規(guī)定荷載向下為正；集中力偶順時(shí)針轉(zhuǎn)為正。結(jié)果輸出時(shí)，規(guī)定彎矩使桿件下側(cè)受拉為正，剪力使桿件順時(shí)針轉(zhuǎn)為正（這里MATLAB程序未給出）。

例1：求如圖3所示簡支梁C截面的內(nèi)力。

圖3中nq=2，q=[10 15]，a=[4 0]，b=[9 13]；nF=1，F(xiàn)=50，c=4；nM=1，M=-30，d=11；L=13，x=6。運(yùn)行程序后得到，M=566.54kN·m，F(xiàn)s=-0.58kN。

二、三鉸拱的內(nèi)力

考慮拱高為f的三鉸平拱，設(shè)在豎向荷載作用下的水平反力FH已經(jīng)求出，其對(duì)應(yīng)簡支梁的任一截面的彎矩和剪力用上述方法已可求出，則根據(jù)三鉸平拱的內(nèi)力計(jì)算公式（1），只要確定截面的豎向坐標(biāo)y和夾角φ，即可求出三鉸拱任一截面的內(nèi)力。截面的豎向坐標(biāo)y和夾角φ的確定過程如下：先對(duì)拱軸線方程y=f（x）求導(dǎo)，得到tanφ=y'，這可以通過MATLAB中的diff（y）實(shí)現(xiàn)，然后求出φ=arctan（y'）。最后計(jì)算指定x1截面的y1和φ1值，可通過y1=subs（y，x1）和φ1=subs（φ，x1）實(shí)現(xiàn)。

例2：求如圖4所示三鉸拱D截面的內(nèi)力。已知拱軸線方程為y=4fx（L-x）/L2。

輸入nq=1，q=10，a=6，b=12；nF=2，F(xiàn)p=[10 20]，c=[2 4]；nM=0；M=0，d=0；f=4，L=12，x=9。運(yùn)行程序后得D截面的內(nèi)力為：M=10kN·m，F(xiàn)s=0，F(xiàn)N=-42.06kN。同樣得到C截面的M=0，F(xiàn)s=6.67kN，F(xiàn)N=-35kN。

三、結(jié)語

用MATLAB編制了通用程序，可快速簡捷地求解簡支梁和三鉸拱在復(fù)雜荷載作用下任一截面的內(nèi)力。仿照上述編程思想，也可解決其他類型簡單梁的內(nèi)力計(jì)算和其他問題。

參考文獻(xiàn)：

[1]何結(jié)兵，顧愛軍，楚海建.Matlab在工程力學(xué)課程教學(xué)中的實(shí)踐[J].力學(xué)與實(shí)踐，2003，25（6）：73-75.

[2]曾又林，周建波，蔣寅軍，等.結(jié)構(gòu)力學(xué)題解[M].武漢：華中科技大學(xué)出版社，2005.

[3]夏健明.用Excel繪制三鉸拱的內(nèi)力圖[J].力學(xué)與實(shí)踐，2010，（4）：104-106.

[4]周博，李喜斌.基于Matlab的材料力學(xué)電算法[J].應(yīng)用科技， 2002，29（005）：37-39.

摘要：在《結(jié)構(gòu)力學(xué)》課程教學(xué)中，三鉸拱的內(nèi)力計(jì)算是個(gè)難點(diǎn)。本文應(yīng)用MATLAB編制了通用程序，可以解決三鉸拱在復(fù)雜荷載作用下截面內(nèi)力的計(jì)算問題。

關(guān)鍵詞：結(jié)構(gòu)力學(xué)；MATLAB；三鉸拱；內(nèi)力

中圖分類號(hào)：G642.0？搖 文獻(xiàn)標(biāo)志碼：A 文章編號(hào)：1674-9324（2014）05-0096-02

三鉸拱是《結(jié)構(gòu)力學(xué)》課程其中的一個(gè)章節(jié)，由于拱的截面法線與水平軸的夾角φ隨截面位置發(fā)生改變，截面內(nèi)力又與φ有關(guān)，因而求解三鉸拱，尤其在復(fù)雜荷載作用下截面的內(nèi)力是件煩瑣的事情。為了讓師生能從煩瑣、重復(fù)的手工計(jì)算的勞動(dòng)中解放出來，把更多的時(shí)間用于基本概念、基本理論的理解及問題的解決方法上[1]，本文利用MATLAB軟件的符號(hào)運(yùn)算、數(shù)值計(jì)算、編程等功能，編制了三鉸拱在復(fù)雜荷載作用下截面內(nèi)力計(jì)算的通用程序，不僅提高了教學(xué)效率，也激發(fā)了學(xué)生學(xué)習(xí)《結(jié)構(gòu)力學(xué)》和MATLAB的極大興趣。

三鉸平拱在只有豎向荷載作用下任一截面的內(nèi)力與同跨度、同荷載的對(duì)應(yīng)簡支梁對(duì)應(yīng)截面的內(nèi)力有如下關(guān)系[2]：

M=M0-FHy，F(xiàn)s=F0scosφ-FHsinφ，F(xiàn)N=-F0Ssinφ-FHcosφ （1）

因此，先介紹簡支梁任一截面的內(nèi)力計(jì)算，再介紹三鉸拱任一截面的內(nèi)力計(jì)算。

一、簡支梁的內(nèi)力

若簡支梁[3]上作用有均布荷載qi（i=1，nq），集中荷載 Fpi（i=1，nF）和集中力偶Mi（i=1，nM）。均布荷載的左、右端點(diǎn)到支座A的距離分別為ai（i=1，nq）、bi（i=1，nq），集中荷載和集中力偶到支座A的距離分別為ci（i=1，nF）、di（i=1，nM）（見圖1）。設(shè)支座反力FAy、FBy已經(jīng)由平衡方程求出。

求任一截面x的內(nèi)力。根據(jù)截面法，取x截面以左部分為隔離體，考慮隔離體上的荷載是任意、復(fù)雜的，可能有均布荷載，可能有集中荷載，也有可能有集中力偶，且荷載數(shù)目可能是多個(gè)。小變形條件下，可用疊加法計(jì)算內(nèi)力[4]。根據(jù)截面上的剪力和彎矩分別等于隔離體上所有外力沿截面切向的投影代數(shù)和以及對(duì)所求截面力矩代數(shù)和的計(jì)算法則，分別考慮隔離體上支座反力FAy、均布荷載qi、集中荷載Fpi和集中力偶Mi等對(duì)x截面內(nèi)力的影響。

（一）反力FAy的影響

M1=FAy·x，F(xiàn)S1=FAy （2）

（二）均布荷載qi的影響

均布荷載與所求截面的位置有如下3種情況：（1）均布荷載qi全部作用在隔離體上（圖2a），此時(shí)考慮全部均布荷載的作用；（2）均布荷載qi部分作用在隔離體上，部分在隔離體外（圖2b），此時(shí)只考慮隔離體上部分荷載的作用，則 bi=x；（3）均布荷載qi全部作用在隔離體外（圖2c），則不考慮均布荷載的作用，即qi=0，ai=bi=0。

nq個(gè)均布荷載產(chǎn)生的x截面內(nèi)力為：M2=-■qi（bi-ai）xi-■，F(xiàn)S2=-■qi（bi-ai） （3）

對(duì)于不同的均布荷載情況，公式中的qi、ai、bi取不同值。

（三）集中荷載Fpi的影響

集中荷載與所求截面的位置有如下2種情況：（1）集中荷載Fpi作用在隔離體上，應(yīng)考慮Fpi的作用；（2）集中荷載 Fpi作用在隔離體外，不考慮Fpi的作用，則Fpi=0，ci=0。

nF個(gè)集中荷載產(chǎn)生的x截面內(nèi)力為：M3=-■Fpi（xi=ci）， FS3=-■Fpi （4）

對(duì)于不同的集中荷載情況，F(xiàn)pi、ci取不同值。

（四）集中力偶Mi的影響

集中力偶與所求截面的位置有兩種情況：（1）集中力偶Mi作用在隔離體上，應(yīng)考慮Mi的作用，則di

nM個(gè)集中力偶產(chǎn)生的x截面內(nèi)力為：M4=-■Mi，F(xiàn)S4=0（5）

對(duì)于不同的集中力偶情況，Mi取不同值。

則反力FAy、均布荷載qi、集中荷載FPi和集中力偶Mi共同作用下x截面內(nèi)力為：

M=M1+M2+M3+M4，F(xiàn)s=Fs1+Fs2+Fs3+Fs4 （6）

利用MATLAB編制通用的程序，只需分別將均布荷載、集中荷載和集中力偶的個(gè)數(shù)：nq、nF、nM、大小qi、Fpi、Mi，以及分別到左端支座A的距離ai、bi、ci、di，梁長L以及截面位置x等輸入，就可以利用程序快速求解出復(fù)雜荷載作用下簡支梁任一截面的內(nèi)力。數(shù)據(jù)輸入時(shí)，規(guī)定荷載向下為正；集中力偶順時(shí)針轉(zhuǎn)為正。結(jié)果輸出時(shí)，規(guī)定彎矩使桿件下側(cè)受拉為正，剪力使桿件順時(shí)針轉(zhuǎn)為正（這里MATLAB程序未給出）。

例1：求如圖3所示簡支梁C截面的內(nèi)力。

圖3中nq=2，q=[10 15]，a=[4 0]，b=[9 13]；nF=1，F(xiàn)=50，c=4；nM=1，M=-30，d=11；L=13，x=6。運(yùn)行程序后得到，M=566.54kN·m，F(xiàn)s=-0.58kN。

二、三鉸拱的內(nèi)力

考慮拱高為f的三鉸平拱，設(shè)在豎向荷載作用下的水平反力FH已經(jīng)求出，其對(duì)應(yīng)簡支梁的任一截面的彎矩和剪力用上述方法已可求出，則根據(jù)三鉸平拱的內(nèi)力計(jì)算公式（1），只要確定截面的豎向坐標(biāo)y和夾角φ，即可求出三鉸拱任一截面的內(nèi)力。截面的豎向坐標(biāo)y和夾角φ的確定過程如下：先對(duì)拱軸線方程y=f（x）求導(dǎo)，得到tanφ=y'，這可以通過MATLAB中的diff（y）實(shí)現(xiàn)，然后求出φ=arctan（y'）。最后計(jì)算指定x1截面的y1和φ1值，可通過y1=subs（y，x1）和φ1=subs（φ，x1）實(shí)現(xiàn)。

例2：求如圖4所示三鉸拱D截面的內(nèi)力。已知拱軸線方程為y=4fx（L-x）/L2。

輸入nq=1，q=10，a=6，b=12；nF=2，F(xiàn)p=[10 20]，c=[2 4]；nM=0；M=0，d=0；f=4，L=12，x=9。運(yùn)行程序后得D截面的內(nèi)力為：M=10kN·m，F(xiàn)s=0，F(xiàn)N=-42.06kN。同樣得到C截面的M=0，F(xiàn)s=6.67kN，F(xiàn)N=-35kN。

三、結(jié)語

用MATLAB編制了通用程序，可快速簡捷地求解簡支梁和三鉸拱在復(fù)雜荷載作用下任一截面的內(nèi)力。仿照上述編程思想，也可解決其他類型簡單梁的內(nèi)力計(jì)算和其他問題。

參考文獻(xiàn)：

[1]何結(jié)兵，顧愛軍，楚海建.Matlab在工程力學(xué)課程教學(xué)中的實(shí)踐[J].力學(xué)與實(shí)踐，2003，25（6）：73-75.

[2]曾又林，周建波，蔣寅軍，等.結(jié)構(gòu)力學(xué)題解[M].武漢：華中科技大學(xué)出版社，2005.

[3]夏健明.用Excel繪制三鉸拱的內(nèi)力圖[J].力學(xué)與實(shí)踐，2010，（4）：104-106.

[4]周博，李喜斌.基于Matlab的材料力學(xué)電算法[J].應(yīng)用科技， 2002，29（005）：37-39.

摘要：在《結(jié)構(gòu)力學(xué)》課程教學(xué)中，三鉸拱的內(nèi)力計(jì)算是個(gè)難點(diǎn)。本文應(yīng)用MATLAB編制了通用程序，可以解決三鉸拱在復(fù)雜荷載作用下截面內(nèi)力的計(jì)算問題。

關(guān)鍵詞：結(jié)構(gòu)力學(xué)；MATLAB；三鉸拱；內(nèi)力

中圖分類號(hào)：G642.0？搖 文獻(xiàn)標(biāo)志碼：A 文章編號(hào)：1674-9324（2014）05-0096-02

三鉸拱是《結(jié)構(gòu)力學(xué)》課程其中的一個(gè)章節(jié)，由于拱的截面法線與水平軸的夾角φ隨截面位置發(fā)生改變，截面內(nèi)力又與φ有關(guān)，因而求解三鉸拱，尤其在復(fù)雜荷載作用下截面的內(nèi)力是件煩瑣的事情。為了讓師生能從煩瑣、重復(fù)的手工計(jì)算的勞動(dòng)中解放出來，把更多的時(shí)間用于基本概念、基本理論的理解及問題的解決方法上[1]，本文利用MATLAB軟件的符號(hào)運(yùn)算、數(shù)值計(jì)算、編程等功能，編制了三鉸拱在復(fù)雜荷載作用下截面內(nèi)力計(jì)算的通用程序，不僅提高了教學(xué)效率，也激發(fā)了學(xué)生學(xué)習(xí)《結(jié)構(gòu)力學(xué)》和MATLAB的極大興趣。

三鉸平拱在只有豎向荷載作用下任一截面的內(nèi)力與同跨度、同荷載的對(duì)應(yīng)簡支梁對(duì)應(yīng)截面的內(nèi)力有如下關(guān)系[2]：

M=M0-FHy，F(xiàn)s=F0scosφ-FHsinφ，F(xiàn)N=-F0Ssinφ-FHcosφ （1）

因此，先介紹簡支梁任一截面的內(nèi)力計(jì)算，再介紹三鉸拱任一截面的內(nèi)力計(jì)算。

一、簡支梁的內(nèi)力

若簡支梁[3]上作用有均布荷載qi（i=1，nq），集中荷載 Fpi（i=1，nF）和集中力偶Mi（i=1，nM）。均布荷載的左、右端點(diǎn)到支座A的距離分別為ai（i=1，nq）、bi（i=1，nq），集中荷載和集中力偶到支座A的距離分別為ci（i=1，nF）、di（i=1，nM）（見圖1）。設(shè)支座反力FAy、FBy已經(jīng)由平衡方程求出。

求任一截面x的內(nèi)力。根據(jù)截面法，取x截面以左部分為隔離體，考慮隔離體上的荷載是任意、復(fù)雜的，可能有均布荷載，可能有集中荷載，也有可能有集中力偶，且荷載數(shù)目可能是多個(gè)。小變形條件下，可用疊加法計(jì)算內(nèi)力[4]。根據(jù)截面上的剪力和彎矩分別等于隔離體上所有外力沿截面切向的投影代數(shù)和以及對(duì)所求截面力矩代數(shù)和的計(jì)算法則，分別考慮隔離體上支座反力FAy、均布荷載qi、集中荷載Fpi和集中力偶Mi等對(duì)x截面內(nèi)力的影響。

（一）反力FAy的影響

M1=FAy·x，F(xiàn)S1=FAy （2）

（二）均布荷載qi的影響

均布荷載與所求截面的位置有如下3種情況：（1）均布荷載qi全部作用在隔離體上（圖2a），此時(shí)考慮全部均布荷載的作用；（2）均布荷載qi部分作用在隔離體上，部分在隔離體外（圖2b），此時(shí)只考慮隔離體上部分荷載的作用，則 bi=x；（3）均布荷載qi全部作用在隔離體外（圖2c），則不考慮均布荷載的作用，即qi=0，ai=bi=0。

nq個(gè)均布荷載產(chǎn)生的x截面內(nèi)力為：M2=-■qi（bi-ai）xi-■，F(xiàn)S2=-■qi（bi-ai） （3）

對(duì)于不同的均布荷載情況，公式中的qi、ai、bi取不同值。

（三）集中荷載Fpi的影響

集中荷載與所求截面的位置有如下2種情況：（1）集中荷載Fpi作用在隔離體上，應(yīng)考慮Fpi的作用；（2）集中荷載 Fpi作用在隔離體外，不考慮Fpi的作用，則Fpi=0，ci=0。

nF個(gè)集中荷載產(chǎn)生的x截面內(nèi)力為：M3=-■Fpi（xi=ci）， FS3=-■Fpi （4）

對(duì)于不同的集中荷載情況，F(xiàn)pi、ci取不同值。

（四）集中力偶Mi的影響

集中力偶與所求截面的位置有兩種情況：（1）集中力偶Mi作用在隔離體上，應(yīng)考慮Mi的作用，則di

nM個(gè)集中力偶產(chǎn)生的x截面內(nèi)力為：M4=-■Mi，F(xiàn)S4=0（5）

對(duì)于不同的集中力偶情況，Mi取不同值。

則反力FAy、均布荷載qi、集中荷載FPi和集中力偶Mi共同作用下x截面內(nèi)力為：

M=M1+M2+M3+M4，F(xiàn)s=Fs1+Fs2+Fs3+Fs4 （6）

利用MATLAB編制通用的程序，只需分別將均布荷載、集中荷載和集中力偶的個(gè)數(shù)：nq、nF、nM、大小qi、Fpi、Mi，以及分別到左端支座A的距離ai、bi、ci、di，梁長L以及截面位置x等輸入，就可以利用程序快速求解出復(fù)雜荷載作用下簡支梁任一截面的內(nèi)力。數(shù)據(jù)輸入時(shí)，規(guī)定荷載向下為正；集中力偶順時(shí)針轉(zhuǎn)為正。結(jié)果輸出時(shí)，規(guī)定彎矩使桿件下側(cè)受拉為正，剪力使桿件順時(shí)針轉(zhuǎn)為正（這里MATLAB程序未給出）。

例1：求如圖3所示簡支梁C截面的內(nèi)力。

圖3中nq=2，q=[10 15]，a=[4 0]，b=[9 13]；nF=1，F(xiàn)=50，c=4；nM=1，M=-30，d=11；L=13，x=6。運(yùn)行程序后得到，M=566.54kN·m，F(xiàn)s=-0.58kN。

二、三鉸拱的內(nèi)力

考慮拱高為f的三鉸平拱，設(shè)在豎向荷載作用下的水平反力FH已經(jīng)求出，其對(duì)應(yīng)簡支梁的任一截面的彎矩和剪力用上述方法已可求出，則根據(jù)三鉸平拱的內(nèi)力計(jì)算公式（1），只要確定截面的豎向坐標(biāo)y和夾角φ，即可求出三鉸拱任一截面的內(nèi)力。截面的豎向坐標(biāo)y和夾角φ的確定過程如下：先對(duì)拱軸線方程y=f（x）求導(dǎo)，得到tanφ=y'，這可以通過MATLAB中的diff（y）實(shí)現(xiàn)，然后求出φ=arctan（y'）。最后計(jì)算指定x1截面的y1和φ1值，可通過y1=subs（y，x1）和φ1=subs（φ，x1）實(shí)現(xiàn)。

例2：求如圖4所示三鉸拱D截面的內(nèi)力。已知拱軸線方程為y=4fx（L-x）/L2。

輸入nq=1，q=10，a=6，b=12；nF=2，F(xiàn)p=[10 20]，c=[2 4]；nM=0；M=0，d=0；f=4，L=12，x=9。運(yùn)行程序后得D截面的內(nèi)力為：M=10kN·m，F(xiàn)s=0，F(xiàn)N=-42.06kN。同樣得到C截面的M=0，F(xiàn)s=6.67kN，F(xiàn)N=-35kN。

三、結(jié)語

用MATLAB編制了通用程序，可快速簡捷地求解簡支梁和三鉸拱在復(fù)雜荷載作用下任一截面的內(nèi)力。仿照上述編程思想，也可解決其他類型簡單梁的內(nèi)力計(jì)算和其他問題。

參考文獻(xiàn)：

[1]何結(jié)兵，顧愛軍，楚海建.Matlab在工程力學(xué)課程教學(xué)中的實(shí)踐[J].力學(xué)與實(shí)踐，2003，25（6）：73-75.

[2]曾又林，周建波，蔣寅軍，等.結(jié)構(gòu)力學(xué)題解[M].武漢：華中科技大學(xué)出版社，2005.

[3]夏健明.用Excel繪制三鉸拱的內(nèi)力圖[J].力學(xué)與實(shí)踐，2010，（4）：104-106.

[4]周博，李喜斌.基于Matlab的材料力學(xué)電算法[J].應(yīng)用科技， 2002，29（005）：37-39.