C語言實訓(xùn)報告

C語言實訓(xùn)報告機械工程學(xué)院機設(shè)xxxx班Xxx一、設(shè)計題目：蝸桿傳動的優(yōu)化設(shè)計設(shè)計一普通圓柱蝸桿，已知參數(shù)：輸入功率＝8.8；蝸桿轉(zhuǎn)速；傳動比；由電動機驅(qū)動，載荷平衡。蝸桿材料，表面硬度，蝸輪材料為，離心鑄造，蝸桿減速器每日8h，全年按300個工作日計，要示工作壽命不低于10年。根據(jù)以上條件進行優(yōu)化設(shè)計，通常在滿足使用要求的前提下，以結(jié)構(gòu)尺寸是否緊湊、傳動效率是否較高作為評判設(shè)計優(yōu)劣的指標(biāo)，在此以傳動中心別具匠心為目標(biāo)函數(shù)式中：——傳動中心距；——蝸桿軸向和蝸輪端面模數(shù)；——蝸桿直徑系數(shù)；——蝸桿頭數(shù)；——蝸輪齒數(shù)；——工程傳動比；題目分析：由上式可知傳動中心別具匠心與模數(shù)、蝸桿直徑系數(shù)和蝸桿頭數(shù)有關(guān)。此三個參數(shù)可作為獨立設(shè)計變量，即，與此相對應(yīng)，目標(biāo)函數(shù)可寫為:在進行蝸桿傳動設(shè)計時，各參數(shù)應(yīng)滿足強度和剛度方面的要求，應(yīng)用網(wǎng)格法可求解這一問題。網(wǎng)格法是約束直接優(yōu)化方法中較為簡單的一種方法，它的基本思想是將可行域分為許多網(wǎng)格，求出滿足設(shè)計約束的網(wǎng)格點上的目標(biāo)函數(shù)值，比較它們的大小，從中選擇函數(shù)值最小的網(wǎng)格點。依次循環(huán)，直到網(wǎng)格之間的距離達到控制精度，即可得滿足精度要求的近似最優(yōu)解。算法步驟：網(wǎng)格法的算法步驟為：（1）給定目標(biāo)函數(shù)初值（一個足夠大的正數(shù)）、對應(yīng)各設(shè)計變量的等分?jǐn)?shù)和計算精度。（2）將區(qū)間進行等分，間距為，各分點坐標(biāo)為,式中:，共有T個分點：（3）對T個分點按順序逐一進行可行性檢查，放棄那些不滿足設(shè)計約束的網(wǎng)格點，計算滿足設(shè)計約束的網(wǎng)格點所對應(yīng)的目標(biāo)函數(shù)值，并與目標(biāo)函數(shù)初值比較，若，則。否則，判別下一個網(wǎng)格點。（4）如果則停止計算，即為所求的最優(yōu)點和最優(yōu)值。否則取轉(zhuǎn)向步驟（2）繼續(xù)計算。四、流程圖源程序#include<stdio.h>/*預(yù)處理命令*/#include<stdlib.h>#include<math.h>main()/*主函數(shù)*/{doublefun(doublex[3]);intyesorno(intng,doublex[3]);doublex[3];/*定義各變量*/doublea[3]={2.0,7.0,2.0},aa[3];/*定義網(wǎng)格區(qū)間上限*/doubleb[3]={8.0,25.0,3.0},bb[3];/*定義網(wǎng)格區(qū)間下限*/intn=3;intm[3]={3,6,5};/*定義網(wǎng)格區(qū)間的等分值*/doubleeps=0.1;doublex0[3]={0.2,0.3};doublexmin[20],fmin;doublefo=10.0e5,f;/*給f0初值用于比較*/doubleh[20];doublehmax;inti;inti0,i1,i2;intng=10;/*參數(shù)約束條件個數(shù)*/do{for(i=0;i<n;i++)/*精度使用網(wǎng)格法*/{h[i]=(b[i]-a[i])/m[i];/*將區(qū)間[ai,bi]進行mi等分*/aa[i]=a[i];bb[i]=b[i];}for(i2=0;i2<=m[2];i2++)/*網(wǎng)格法劃分區(qū)間*/{x[2]=a[2]+i2*h[2];/*確定三維網(wǎng)格坐標(biāo)*/for(i1=0;i1<=m[1];i1++){x[1]=a[1]+i1*h[1];for(i0=0;i0<=m[0];i0++){x[0]=a[0]+i0*h[0];if(yesorno(ng,x)==0)continue;/*判斷剛度與強度約束判斷*/f=fun(x);if(f>=fo)continue;for(i=0;i<n;i++){x0[i]=x[i];xmin[i]=x[i];}/*取最小值*/fo=f,fmin=f;}}}hmax=0.0;/*網(wǎng)格精度細化*/for(i=0;i<n;i++)if(h[i]>hmax)hmax=h[i];/*求最大間距hmax*/{for(i=0;i<n;i++){a[i]=xmin[i]-h[i];/*進一步使用網(wǎng)格法*/b[i]=xmin[i]+h[i];if(a[i]<aa[i])a[i]=aa[i];if(b[i]>bb[i])b[i]=bb[i];}fo=10.0e4;}}while(hmax>eps);printf("theresults:\n");for(i=0;i<n;i++)printf("x**=%10.5e\n",xmin[i]);/*輸出x[0]x[1]x[2]的最小值*/printf("fo=%10.5e\n",fmin);/*輸出f0的最小值*/printf("hmax=%10.5e\n",hmax);/*輸出hmax的最小值*/getchar();return(0);}intyesorno(intng,doublex[3])/*定義yesorno函數(shù)，用于判斷是否滿足約束條件*/{inti;doubleZp,Z1,Z2,Z3,Z4,Z5;intyesno=1;doublegx[20];gx[0]=-x[0]+2.0;/*蝸桿軸向和蝸桿端面模數(shù)m的范圍：1.0<=m<=8.0*/gx[1]=-8.0+x[0];gx[2]=-25.0+x[1];/*蝸桿直徑系數(shù)q的范圍：7.0<=q<=25.0*/gx[3]=-x[1]+7.0;gx[4]=-x[2]+2.0;/*蝸桿頭數(shù)z1的范圍：1.0<=z1<=6.0*/gx[5]=-3.0+x[2];gx[6]=-atan(x[2]/x[1])+10.0/180.0;/*蝸桿導(dǎo)程角計算約束*/gx[7]=-20.0/180.0+atan(x[2]/x[1]);Z1=-0.8881*2.0*x[1]/(x[1]+18.0*x[2]);/*z1到z5的計算均是用于校核接觸和彎曲強度*/Zp=4.02*exp(Z1);/*計算e的z1次方的值，再乘以4.02*/Z2=x[0]*(x[1]+18.0*x[2]);Z3=pow(Z2,3.0);/*計算z2的3次方的值*/Z4=0.125*Z3;gx[8]=252.0-170100.0*Zp/sqrt(Z4);/*蝸桿的接觸強度校核計算*/Z5=x[1]+1.0;gx[9]=146.0-745000.0/(pow(x[0],3.0)*x[2]*(0.5+sqrt(Z5)));/*蝸桿彎曲強度校核*/for(i=0;i<ng;i++)if(gx[i]>0.0)/*如果滿足約束條件，yesorno=0*/{yesno=0;}printf("yesorno=%d\n",yesno);returnyesno;}doublefun(doublex[20])/*定義fun函數(shù),用于求中心距*/{doublef;f=0.5*x[0]*(x[1]+18*x[2]);/*中心距計算公式*/return(f);}運行結(jié)果：六、強度校核程序運行結(jié)果是：模數(shù)m=2蝸桿直徑系數(shù)q=17.926蝸桿頭數(shù)=2中心距a=53.926mm齒頂高=3.7037×mm查機械設(shè)計手冊，選擇標(biāo)準(zhǔn)尺寸：a/min:50i:19.5①m/mm:2d1/mm:22.4z1:2z2:39x2:-0.1γ:10°07′29"普通圓柱蝸桿傳動的基本參數(shù)及其匹配(摘自GB10085-1988neqDIN3976-1980)標(biāo)準(zhǔn)傳動比蝸輪齒面接觸強度校核：根據(jù)表面材料查表：根據(jù)接觸系數(shù)：載荷平穩(wěn)，使用系數(shù)：=1.00傳動比：i=19.5轉(zhuǎn)矩查圖表，傳動效率轉(zhuǎn)速不變，轉(zhuǎn)速系數(shù)壽命系數(shù)接觸疲勞極限接觸疲勞強度的最小安全系數(shù)，取校核合格實訓(xùn)心得通過計算機c語言強化實習(xí)，感受成功喜悅的同時，也有許多的感慨，我在這方面的知識也得到了很大的提升，在這之前我們已經(jīng)對這門課程學(xué)了一個學(xué)期，對其有一定的了解，但僅僅是停留在了解的程度而已，但真的是不知道該怎樣去運用，感到很棘手。真正對c語言開始運用，進行實際操作應(yīng)該是在這次實訓(xùn)，它是一門操作性比較強的課程，只有自己親自編寫，調(diào)試，改正，到最后運行出正確的結(jié)果，才能在這個過程中對c語言有更深的了解。在老師的細心指導(dǎo)和我的不斷努力下，我順利的完