材料力學(xué)課程設(shè)計--曲柄軸的強(qiáng)度設(shè)計、疲勞強(qiáng)度校核及剛度計算

1、<p><b>　　材料力學(xué)課程設(shè)計</b></p><p><b>　　計算說明書</b></p><p>　　設(shè)計題目：曲柄軸的強(qiáng)度設(shè)計、疲勞強(qiáng)度校核及剛度計算</p><p><b>　　目錄</b></p><p>　　一、設(shè)計目的··&#

2、183;····································

3、··········3</p><p>　　二、設(shè)計任務(wù)和要求·····················&

4、#183;··················3</p><p>　　2.1、設(shè)計計算說明書的要求···········

5、83;······························3</p><p>　　2.2、分析討論及說明書部

6、分的要求····································

7、·4</p><p>　　2.3、程序計算部分的要求·····························

8、83;···············4</p><p>　　三、設(shè)計題目················&

9、#183;·································4</p><

10、p>　　1）畫出曲柄軸的內(nèi)力圖·································

11、3;············5</p><p>　　2）設(shè)計主軸頸D和曲柄頸直徑d··················

12、;·····················8</p><p>　　3）校核曲柄臂的強(qiáng)度··········

13、;····································

14、83;·8</p><p>　　4）校核主軸頸飛輪處的疲勞強(qiáng)度····························

15、3;········ 14</p><p>　　5）用能量法計算A端截面的轉(zhuǎn)角 ，·····················&

16、#183;········· 15</p><p>　　四、分析討論及必要說明····················

17、3;·············18五、設(shè)計的改進(jìn)措施及方法·····················

18、··········19六、設(shè)計體會·························&#

19、183;······················19七、參考文獻(xiàn)············

20、3;···································20</p

21、><p><b>　　附錄</b></p><p>　　一.流程圖·························

22、3;································21</p><p>　　二

23、.C語言程序···································

24、3;··················22</p><p>　　計算輸出結(jié)果·············&

25、#183;····································

26、;25</p><p><b>　　設(shè)計目的</b></p><p>　　本課程設(shè)計是在系統(tǒng)學(xué)完材料力學(xué)課程之后，結(jié)合工程實際中的問題，運(yùn)用材料力學(xué)的基本理論和計算方法，獨(dú)立地計算工程中的典型零部件，以達(dá)到綜合利用材料力學(xué)知識解決工程實際問題的目的。同時，可以使學(xué)生將材料力學(xué)的理論和現(xiàn)代計算方法及手段融為一體，既從整體上掌握了基本理論和現(xiàn)代計算方法，又提高了分析問題、

27、解決問題的能力；既是對以前所學(xué)知識（高等數(shù)學(xué)、工程圖學(xué)、理論力學(xué)、算法語言、計算機(jī)和材料力學(xué)等）的綜合運(yùn)用，又為后續(xù)課程的學(xué)習(xí)打下基礎(chǔ)，并初步掌握工程設(shè)計思路和設(shè)計方法，使實際工作能力有所提高。具體有一下六項：</p><p>　　1、使所學(xué)的材料力學(xué)知識系統(tǒng)化、完整化。</p><p>　　2、在系統(tǒng)全面復(fù)習(xí)的基礎(chǔ)上，運(yùn)用材料力學(xué)知識解決工程實際中的問題。</p><

28、p>　　3、由于選題力求結(jié)合專業(yè)實際，因而課程設(shè)計可把材料力學(xué)與專業(yè)需要結(jié)合起來。</p><p>　　4、綜合運(yùn)用以前所學(xué)的各門課程的知識（高等數(shù)學(xué)、工程圖學(xué)、理論力學(xué)、算法語言、計算機(jī)等），使相關(guān)學(xué)科的知識有機(jī)地聯(lián)系起來。</p><p>　　5、初步了解和掌握工程實際中的設(shè)計思路和設(shè)計方法。</p><p>　　6、為后續(xù)課程的教學(xué)打下基礎(chǔ)。</p

29、><p><b>　　二、設(shè)計任務(wù)和要求</b></p><p>　　參加設(shè)計者要系統(tǒng)復(fù)習(xí)材料力學(xué)課程的全部基本理論和方法，獨(dú)立分析、判斷設(shè)計題目的已知條件和所求問題，畫出受力分析計算簡圖和內(nèi)力圖，列出理論依據(jù)并到處計算公式，獨(dú)立編制計算機(jī)程序，通過計算機(jī)給出計算結(jié)果，并完成設(shè)計計算說明書。</p><p>　　2.1 設(shè)計計算說明書的要求<

30、/p><p>　　設(shè)計計算說明書是該題目設(shè)計思路、設(shè)計方法和設(shè)計結(jié)果的說明，要求書寫工整，語言簡練，條理清晰、明確，表達(dá)完整。具體內(nèi)容應(yīng)包括：</p><p>　　設(shè)計題目的已知條件、所求及零件圖。</p><p>　　畫出結(jié)構(gòu)的受力分析計算簡圖，按比例標(biāo)明尺寸、載荷及支座等。</p><p>　　靜不定結(jié)構(gòu)要畫出所選擇的基本靜定系統(tǒng)及與之相應(yīng)的

31、全部求和過程。</p><p>　　畫出全部內(nèi)力圖，并標(biāo)明可能的各危險截面。</p><p>　　危險截面上各種應(yīng)力的分布規(guī)律圖及由此判定各危險點(diǎn)處的應(yīng)力狀態(tài)圖。</p><p>　　各危險點(diǎn)的主應(yīng)力大小及主平面位置。</p><p>　　選擇強(qiáng)度理論并建立強(qiáng)度條件。</p><p>　　列出全部計算過程的理論依據(jù)、公

32、式推導(dǎo)過程以及必要的說明。</p><p>　　對變形及剛度分析要寫明所用的能量法計算過程及必要的內(nèi)力圖和單位力圖。</p><p>　　疲勞強(qiáng)度計算部分要說明循環(huán)特性，，，r，，的計算，所查，，各系數(shù)的依據(jù)，疲勞強(qiáng)度校核過程及結(jié)果，并繪出構(gòu)件的持久極限曲線。</p><p>　　2.2、分析討論及說明部分的要求</p><p>　　分析計算

33、結(jié)果是否合理，并討論其原因、改進(jìn)措施。</p><p>　　提出改進(jìn)設(shè)計的初步方案及設(shè)想。</p><p>　　提高強(qiáng)度、剛度及穩(wěn)定性的措施及建議。</p><p>　　2.3、程序計算部分的要求</p><p><b>　　程序框圖。</b></p><p><b>　　計算機(jī)程序。&

34、lt;/b></p><p>　　打印結(jié)果（數(shù)據(jù)結(jié)果要填寫到設(shè)計計算說明書上）。</p><p><b>　　三、設(shè)計題目</b></p><p>　　某曲柄軸材料為球墨鑄鐵（QT400—10），[]=120MPa，曲柄臂抽象為矩形（圖7-12），h=1.2D,b/h=2/3(左、右臂尺寸相同)，,有關(guān)數(shù)據(jù)查下表：</p>

35、<p><b>　　要求：</b></p><p>　　1、畫出曲柄軸的內(nèi)力圖。</p><p>　　2、按強(qiáng)度條件設(shè)計主軸頸D和曲軸頸的直徑d。</p><p>　　3、校核曲柄臂的強(qiáng)度。 </p><p>　　4、安裝飛輪處為鍵槽，校核主軸頸的疲勞強(qiáng)度，取疲勞安全系數(shù)n=2。</p><

36、p>　　鍵槽為端銑加工，軸頸表面為車削加工，。</p><p>　　5、用能量法計算A端截面的轉(zhuǎn)角 。</p><p>　　圖7-2 曲柄軸力學(xué)模型</p><p>　?。ㄒ唬?、畫出曲軸的內(nèi)力圖</p><p><b>　?。?）外力分析</b></p><p>　　畫出曲軸的計算簡圖（如

37、上圖），計算F分力和扭矩。</p><p><b>　　由平衡條件計算反力</b></p><p><b>　?。?）內(nèi)力分析</b></p><p>　　內(nèi)力圖如下，不計彎曲切應(yīng)力，彎矩圖畫在纖維受壓側(cè)，根據(jù)內(nèi)力圖確定危截面。(單位：力/N 力矩/N·m)</p><p>　　1）主軸

38、頸c端面為危險截面，受到兩項彎矩和扭轉(zhuǎn)</p><p>　　2）右曲柄臂下端面為危險截面，受到兩項彎矩和扭轉(zhuǎn)和軸力</p><p>　　3）左曲柄臂上端面為危險截面，受到兩項彎矩，扭轉(zhuǎn)和軸力</p><p>　　4）曲軸頸中間截面D為危險截面，受扭轉(zhuǎn)和兩向彎曲</p><p>　　、設(shè)計曲軸頸直徑d和主軸頸D</p><p

39、>　　主軸頸的危險截面為C截面，受扭轉(zhuǎn)和兩向彎曲,可用第三強(qiáng)度理論計算：</p><p><b>　　∴</b></p><p><b>　　故D取。</b></p><p>　?。?）曲軸頸屬于彎扭組合變形，由第三強(qiáng)度理論可得：</p><p><b>　　∴ </

40、b></p><p><b>　　故d取。</b></p><p><b>　　、校核曲柄臂的強(qiáng)度</b></p><p>　　曲柄臂的危險截面為矩形截面，受扭轉(zhuǎn)、兩向彎曲及軸力的作用。為確定危險點(diǎn)的位置，畫出曲柄臂上截面應(yīng)力分布圖。</p><p>　　①左曲柄臂的強(qiáng)度計算</p>

41、;<p>　　根據(jù)應(yīng)力分布圖可判定出可能的危險點(diǎn)為，，。</p><p><b>　　對點(diǎn)進(jìn)行應(yīng)力分析。</b></p><p>　　點(diǎn)處于單向拉伸和兩項彎曲的拉應(yīng)力，此時扭轉(zhuǎn)切應(yīng)力為零</p><p>　　2）對點(diǎn)進(jìn)行應(yīng)力分析。</p><p>　　由于有扭轉(zhuǎn)切應(yīng)力，查《材料力學(xué)》表3-1，得到，。<

42、;/p><p>　　正應(yīng)力由軸力、繞Z軸的彎矩共同引起。</p><p>　　由于點(diǎn)處于二向應(yīng)力狀態(tài)，故選用第三強(qiáng)度理論：</p><p>　　3）對點(diǎn)進(jìn)行應(yīng)力分析。</p><p><b>　　同理可得：</b></p><p><b>　　應(yīng)用第三強(qiáng)度理論：</b></

43、p><p>　　故，綜上所述3個危險點(diǎn)都是安全的，所以左曲軸臂安全 </p><p>　?、谟仪S臂的強(qiáng)度計算</p><p>　　根據(jù)應(yīng)力分布圖可判定出可能的危險點(diǎn)為，，。</p><p>　　1)對點(diǎn)進(jìn)行應(yīng)力分析。</p><p>　　點(diǎn)處于單向拉伸，所以正應(yīng)力</p><p>　　2）對

44、點(diǎn)進(jìn)行應(yīng)力分析。</p><p>　　由于有扭轉(zhuǎn)切應(yīng)力，查表3-1，得到，。</p><p>　　正應(yīng)力由軸力、繞Z軸的彎矩共同引起。</p><p>　　由于點(diǎn)處于二向應(yīng)力狀態(tài)，故選用第三強(qiáng)度理論：</p><p><b>　　對點(diǎn)進(jìn)行應(yīng)力分析。</b></p><p><b>　　同

45、理可得</b></p><p><b>　　應(yīng)用第三強(qiáng)度理論：</b></p><p>　　故，綜上所述3個危險點(diǎn)都是安全的，所以右曲軸臂安全 </p><p>　　（四）、 校核主軸頸安裝飛輪的鍵槽處截面的疲勞強(qiáng)度</p><p>　　查《材料力學(xué)》得有效應(yīng)力集中系數(shù)，表面質(zhì)量系數(shù)。</p>

46、<p><b>　　已知，，，n=2。</b></p><p>　　安裝飛輪的鍵槽處只受扭轉(zhuǎn)作用，并且由于曲軸的啟停，可以看做在啟動之后正常工作時是看做為所受扭轉(zhuǎn)帶來的切應(yīng)力是最大值，而停止時是最小的，可以看作為一種應(yīng)力循環(huán)。</p><p>　　忽略鍵槽對抗扭截面系數(shù)的影響，飛輪處截面抗扭截面系數(shù)</p><p>　　曲軸工作時

47、，在不變扭矩作用下，</p><p>　　曲軸不工作時， </p><p>　　∴ 故該循環(huán)為脈動循環(huán)。</p><p><b>　　安全系數(shù)：</b></p><p>　　所以，飛輪處截面的疲勞強(qiáng)度是足夠的。</p><p>　?。ㄎ澹?、 用能量法計算A截面的轉(zhuǎn)角。</p&g

48、t;<p>　　采用圖乘法分別求解截面的轉(zhuǎn)角。</p><p><b>　?、?、求：</b></p><p>　　在截面A加一單位力偶矩，并求得支座反力如下圖所示。</p><p><b>　　由平衡方程得：</b></p><p><b>　　當(dāng)時，查表可得：</b&

49、gt;</p><p>　　方向與單位力偶矩相同。</p><p><b>　?、颉⑶螅?lt;/b></p><p>　　在截面A加一單位力偶矩，支座反力，彎矩圖和軸力圖如下圖所示。</p><p>　　方向與單位力偶矩相同。</p><p><b>　　分析討論及必要說明</b>

50、;</p><p>　　在本次設(shè)計中，做以下幾點(diǎn)說明：</p><p>　　在外力分析時，在設(shè)定未知力的時候，由于已知沒有X方向的外力，故未設(shè)。</p><p>　　在畫內(nèi)力圖時，不計彎曲切應(yīng)力，故未畫剪力圖。</p><p>　　在強(qiáng)度計算方面，由于材料是球墨鑄鐵（QT400-10），其物理性質(zhì)與剛相近，用第三強(qiáng)度理論而不用第一或第二強(qiáng)度理

51、論公式。</p><p>　　在校核曲柄臂時，畫內(nèi)力分布時，把曲柄臂的危險截面看成矩形，</p><p>　　在疲勞強(qiáng)度校核飛輪截面時，忽略鍵槽對的影響。</p><p>　　五、設(shè)計的改進(jìn)措施及方法</p><p>　　1、提高曲軸的彎曲強(qiáng)度</p><p>　　合理安排曲軸的受力情況及設(shè)計合理的截面，但對于該曲軸只

52、能采用合理安排曲軸的受力情況。在結(jié)構(gòu)允許的情況下，可采取合理設(shè)計和布置支座或?qū)⒓休d荷適當(dāng)分散。</p><p><b>　　提高曲軸的彎曲剛度</b></p><p>　　提高彎曲剛度的主要措施有：改善結(jié)構(gòu)形式，減少彎矩的數(shù)值、選擇合理的截面及合理選材等。對于該曲軸可以改善結(jié)構(gòu)形式，減少彎矩的數(shù)值并且合理選材，選擇合適的材料。</p><p>

53、;<b>　　提高曲軸的疲勞強(qiáng)度</b></p><p>　　1）減緩應(yīng)力集中。為了消除和緩解應(yīng)力集中，再設(shè)計曲軸時，應(yīng)盡量避免出現(xiàn)方形直角或帶有尖角的孔和槽，即在主軸頸和曲柄臂相連處應(yīng)采用半徑較大的過度圓角</p><p>　　2）提高構(gòu)件表面強(qiáng)度等。提高曲軸表面的強(qiáng)度可通過兩方面實現(xiàn),一是從加工入手提高表面加工質(zhì)量?？刹捎镁?xì)加工，降低表面粗糙度，尤其對高強(qiáng)度鋼更

54、重要；二是增加表層強(qiáng)度，對曲軸中應(yīng)力集中的部位如鍵槽處應(yīng)采取某些工藝措施，即表面熱處理或化學(xué)處理，如表面高頻淬火、滲碳、滾壓、噴丸等。</p><p><b>　　設(shè)計體會</b></p><p>　　本次課程設(shè)計讓我受益頗多，從剛開始摸不到頭腦，到后來自己進(jìn)行計算，用CAD畫圖，用C語言設(shè)計程序，編輯word文檔，收獲了很多，也明白了材料力學(xué)并不是想象的那樣簡單。&

55、lt;/p><p>　　在計算題目時，我感受到，理論和實際還是有一定的不同的，當(dāng)我看到題目后我不懂得如何的去計算它，對于課本上所學(xué)習(xí)到的公式，我不會運(yùn)用到實際，無法從中抽象出我們所需要的模型，隨著與同學(xué)們的溝通和查閱資料，我漸漸地有了頭緒，找到了套用公式的方法；還有在進(jìn)行計算時，有些對于計算結(jié)果并沒有太大影響的條件是可以忽略的，不是所有地方都能算的，不過還好我們能夠粗略的計算一下了。而對于c語言編程更是讓人頭疼的事情

56、，不過，當(dāng)編完之后就會覺得其實也沒什么，當(dāng)然我編的程序并不是那么的完美，或者說根本就不能夠應(yīng)用于實踐，但是這次實踐中的編程應(yīng)用還是讓我覺得c語言沒有白學(xué)，每個學(xué)科之間都是有聯(lián)系的。還有CAD，許久沒有用CAD畫圖，好多操作都不記得了，這次畫圖一邊回憶，一邊摸索，又一次加深了印象，鞏固了所學(xué)。而word這個我曾經(jīng)認(rèn)為我已經(jīng)學(xué)會了軟件，帶給我了很多的麻煩，尤其是公式編輯，排版等等，真的是不應(yīng)用永遠(yuǎn)不知道哪里是不會的。不過，現(xiàn)在這些我都學(xué)會了

57、，很高興！</p><p><b>　　七、參考文獻(xiàn)</b></p><p>　　1、材料力學(xué)/聶毓琴，孟光偉主編，北京：機(jī)械工業(yè)出版社。</p><p>　　2、材料力學(xué)實驗與課程設(shè)計/聶毓琴、吳宏主編，北京：機(jī)械工業(yè)出版社。</p><p>　　3、C程序設(shè)計/第三版，譚浩強(qiáng)主編，北京：清華大學(xué)出版社。</p&

58、gt;<p><b>　　附錄</b></p><p><b>　　流程圖：</b></p><p><b>　　二、c語言程序</b></p><p>　　#include<stdio.h></p><p>　　#include<math.h&

59、gt;</p><p>　　#define Pi 3.1415926</p><p>　　#define n 2</p><p>　　#define i 120e6</p><p>　　#define E 150e9</p><p>　　int main() </p><p>

60、;　　{int angle,d,D;</p><p>　　double temps,tempc,Fy,Fz,Mx,e,l,l1,l2,l3,l4,F,FAy,FAz,FBy,FBz,W,Mx1,My1,Mz1,Mx2,My2,Mz2,DD,g,f,dd,h,b,A,s1,s2,s3,t2,t3,a,sr2,sr3,r,s11,s22,s33,t22,t33,sr22,sr33,tmax,tmin,R,tm,ta,

61、nt,T,Kt,Et,B,pt,Q1,Q11,Q12,Q13,Q14,Q15,Q16,EI1,EI2,EI3,Q2,Q21,Q22,Q23,Q24,Q25,GI1,B1,u,G;</p><p>　　printf("enter the data:F,W,l1,l2,l3,e,angle\n");</p><p>　　scanf("%lf,%lf,%lf,%lf

62、,%lf,%lf,%d",&F,&W,&l1,&l2,&l3,&e,&angle);</p><p>　　temps=sin(angle*Pi/180);</p><p>　　tempc=cos(angle*Pi/180);</p><p>　　Fy=F*temps;</p><p

63、>　　Fz=F*tempc;</p><p><b>　　Mx=Fz*e;</b></p><p>　　printf("Fy=%lf,Fz=%lf,Mx=%lf\n",Fy,Fz,Mx);</p><p>　　FBz=Fz*l1/(l1+l2);</p><p>　　FBy=(Fy*l1-W*(

64、l1+l2+l3))/(l1+l2);</p><p>　　FAz=l2*Fz/(l1+l2);</p><p>　　FAy=(Fy*l2+W*l3)/(l1+l2);</p><p>　　printf("FAy=%lf,FAz=%lf,FBy=%lf,FBz=%lf\n",FAy,FAz,FBy,FBz);</p><p&g

65、t;<b>　　l=1.5*e;</b></p><p><b>　　l4=0.5*l;</b></p><p>　　Mx1=FAz*e+(Fz-FAz)*e;</p><p>　　My1=FAz*(l1+l4)-Fz*l4;</p><p>　　Mz1=FAy*(l1+l4)-Fy*l4;<

66、/p><p>　　f=sqrt(Mx1*Mx1+My1*My1+Mz1*Mz1);</p><p>　　g=32*f/(Pi*i);</p><p>　　DD=pow(g,1.0/3.0); </p><p>　　My2=l1*FAz;</p><p>　　Mx2=e*FAz;</p><p>　　

67、Mz2=FAy*l1;</p><p>　　f=sqrt(Mx2*Mx2+My2*My2+Mz2*Mz2);</p><p>　　g=32*f/(Pi*i);</p><p>　　dd=pow(g,1.0/3.0);</p><p>　　printf("D=%.4lfm,d=%.4lfm\n",DD,dd);</p&

68、gt;<p>　　D=(int)(1000*DD);</p><p>　　d=(int)(1000*dd);</p><p>　　if(d%2==0) </p><p><b>　　d=d+2;</b></p><p>　　else d=d+1;</p><p>　　if(D%2==

69、0) </p><p><b>　　D=D+2;</b></p><p>　　else D=D+1;</p><p>　　printf("故D取%dmm,d取%dmm \n",D,d);</p><p><b>　　a=0.231;</b></p><p>

70、;<b>　　r=0.858;</b></p><p>　　h=1.2*(double)D/1000;</p><p><b>　　b=2*h/3;</b></p><p><b>　　A=h*b;</b></p><p>　　s1=FAy/A+FAz*e*6/(h*h*b)+

71、FAy*6*(l1-l4)/(b*b*h);</p><p>　　printf("s1=%.2lfPa,",s1);</p><p>　　s2=FAy/A+FAy*(l1-l4)*6/(b*b*h);</p><p>　　t2=FAz*(l1-l4)/(a*h*b*b);</p><p>　　sr2=sqrt(s2*s2+

72、4*t2*t2);</p><p>　　printf("sr2=%lfPa,",sr2);</p><p><b>　　t3=r*t2;</b></p><p>　　s3=FAy/A+FAz*e*6/(h*h*b);</p><p>　　sr3=sqrt(s3*s3+4*t3*t3);</p&g

73、t;<p>　　printf("sr3=%lfPa\n",sr3);</p><p>　　if(s1<=i&&sr2<=i&&sr3<=i)</p><p>　　printf("左曲軸臂安全\n");</p><p>　　s11=(W+FBy)/A+Mx1*6/(

74、h*h*b)+(FAy*(l1+l4)-Fy*l4)*6/(h*b*b);</p><p>　　printf("s11=%lfPa,",s11);</p><p>　　s22=(W+FBy)/A+(FAy*(l1+l4)-Fy*l4)*6/(h*b*b);</p><p>　　t22=(FAz*(l1+l4)-Fz*l4)/(a*h*b*b);&

75、lt;/p><p>　　sr22=sqrt(s22*s22+4*t22*t22);</p><p>　　printf("sr22=%lfPa,",sr22);</p><p>　　s33=(W+FBy)/A+Mx1*6/(h*h*b);</p><p>　　t33=r*t22;</p><p>　　sr

76、33=sqrt(s33*s33+4*t33*t33);</p><p>　　printf("sr33=%lfPa\n",sr33);</p><p>　　if(s11<=i&&sr22<=i&&sr33<=i)</p><p>　　printf("右曲軸臂安全\n");<

77、/p><p><b>　　Kt=1.52;</b></p><p><b>　　Et=0.76;</b></p><p><b>　　B=0.95;</b></p><p><b>　　pt=0.05;</b></p><p><

78、b>　　T=160e6;</b></p><p>　　tmax=Mx*16/(Pi*((double)D/1000)*((double)D/1000)*((double)D/1000));</p><p><b>　　tmin=0;</b></p><p>　　R=tmin/tmax;</p><p>

79、　　tm=(tmax+tmin)/2;</p><p>　　ta=(tmax-tmin)/2;</p><p>　　nt=T/(Kt/(Et*B)*ta+pt*tm);</p><p>　　printf("nt=%lf\n",nt);</p><p><b>　　if(nt>=n)</b><

80、;/p><p>　　printf("安全\n");</p><p>　　EI1=E*Pi*((double)D/1000)*((double)D/1000)*((double)D/1000)*((double)D/1000)/64;</p><p>　　EI2=E*h*b*b*b/12;</p><p>　　EI3=E*Pi*

81、((double)d/1000)*((double)d/1000)*((double)d/1000)*((double)d/1000)/64;</p><p>　　Q11=(l1-l4)*FAy*(l1-l4)/2*(l1+2*l4+3*l2)/(3*(l1+l2))/EI1;</p><p>　　Q12=FAy*(l1-l4)*e*(l2+l4)/(l1+l2)/EI2;</p&g

82、t;<p>　　Q13=(FAy*l4*l4/2*(3*l2+l4)/(3*(l1+l2))+FAy*(l1-l4)*l4*(2*l2+l4)/(2*(l1+l2))+(Fy-FAy)*l4*l4/2*(3*l2-l4)/(3*(l1+l2))+(FAy*(l1+l4)-Fy*l4)*l4*(2*l2-l4)/(2*(l1+l2)))/EI3;</p><p>　　Q14=(l2-l4)/(l1+l

83、2)*e*(FAy*(l1+l4)-Fy*l4)/EI2;</p><p>　　Q15=(l2-l4)/(l1+l2)*(l2-l4)/2*(FAy*(l1+l4+(l2-l4)/3)-Fy*(l4+(l2-l4)/3))/EI1;</p><p>　　Q16=e/((l1+l2)*E*h*b)*(W+FBy-FAy);</p><p>　　Q1=Q11+Q12+Q

84、13+Q14+Q15+Q16;</p><p>　　printf("A點(diǎn)處轉(zhuǎn)角θz=%lf\n",Q1);</p><p><b>　　u=0.27;</b></p><p>　　G=E/(2*(1+u));</p><p><b>　　B1=0.196;</b></p&g

85、t;<p>　　GI1=G*B1*h*b*b*b;</p><p>　　Q21=FAz*(l1-l4)*(l1-l4)/2*(l1+3*l2+2*l4)/(3*(l1+l2))/EI1;</p><p>　　Q22=FAz*(l1-l4)*e*(l2+l4)/(l1+l2)/GI1;</p><p>　　Q23=(FAz*l4*l4/2*(3*l2+l

86、4)/(3*(l1+l2))+FAz*(l1-l4)*l4*(2*l2+l4)/(2*(l1+l2))+(-FAz+Fz)*l4*l4/2*(3*l2-l4)/(3*(l1+l2))+(FAz*(l1+l4)-l4*Fz)*l4*(2*l2-l4)/(2*(l1+l2)))/EI3;</p><p>　　Q24=(l2-l4)/(l1+l2)*e*(FAz*(l1+l4)-l4*Fz)/GI1;</p>

87、;<p>　　Q25=(l2-l4)/(l1+l2)*(l2-l4)*(FAz*(l1+l4)-l4*Fz)/(3*EI1);</p><p>　　Q2=Q21+Q22+Q23+Q24+Q25;</p><p>　　printf("A點(diǎn)處轉(zhuǎn)角θy=%lf\n",Q2);</p><p><b>　　return 0;<