杨氏模量数据处理：拉伸法测杨氏模量数据处理代码（山东科技大学）

拉伸法测杨氏模量数据处理代码

代码解释：
先输入标尺读数 x+ ,那一行的六个数；（单位cm）
再输入标尺读数 x- ,那一行的六个数；（单位cm）
再输入测量值 di ,那一行的六个数；（单位mm）
再输入L0,D0,b0;
得出全部结果。

#include <bits/stdc++.h>using namespace std;int main(){ double E; double m=3.0; double k=3.1415926; double a[6];// double b[6 double b[6]; cout<<"please cin 6+ 6- :"<<endl; for(int i=0 ; i<6 ; i++) { cin>>a[ i ]; } for(int i=0 ; i<6 ; i++) { cin>>b[ i ]; } //----------------------------------------------------------------------------------- cout<<"+++ "<<endl; double c[3]; for(int i=0;i<3;i++) { c [ i ] = a[ i + 3 ] - a [ i ] ; cout<< c [ i ] <<" "; } cout<<endl; //----------------------------------------------------------------------------------------------- cout<<"---- "<<endl; double d[3]; for(int i=0;i<3;i++) { d [ i ] = b [ i+3 ] - b [ i ]; cout<< d[ i ]<<" "; } cout<<endl; //---------------------------------------------------------------x的平均值------cm------------- double pingjun=0; for(int i=0;i<3;i++) { pingjun=pingjun+c[ i ]+d [ i ]; } pingjun=pingjun/6.0; cout<<"pingjun "<<pingjun; //----------------------------------------------------------求Sxm------------------cm---------------- double temp; for(int i=0;i<3;i++) { temp = (c[ i ]-pingjun)*(c[ i ]-pingjun)+ (d[ i ]-pingjun)*(d[ i ]-pingjun); } double S; S=temp/5; S=sqrt(S); cout<<"Sxm = "<<S<<endl; //---------------------------------------------------------求Ux----------------cm---------- double Ubx=0.01;//单位换算，换成了Cm double Ux=S*S+Ubx*Ubx; Ux=sqrt(Ux); cout<<"Ux = "<<Ux<<endl;//-------------------------------------------------------------钢丝-----dping------mm------Sd--------Ud----------//钢丝 double di[6]; cout<<"please cin di "<<endl; for(int i=0;i<6;i++) { cin>>di[i]; } double dping=0.547; temp=0; for(int i=0;i<6;i++) { temp=temp+(di[i]-dping)*(di[i]-dping); } double Sd=temp/5; Sd=sqrt (Sd); cout<<"Sd = "<<Sd<<endl; double UBd=0.004; //换算单位，mm double Ud=UBd*UBd+Sd*Sd; Ud=sqrt (Ud); cout<<"Ud = "<<Ud<<endl; //------------------------------------------------------------------ L0 ,D0,b0的测量 ----------------- // double UL=0.002; double UD=0.002; double Ub=0.00002; double L0 ,D0 , b0; dping=0.547; cout<<"cin: L0 D0 b0"<<endl; cin>>L0>>D0>>b0; double g=9.794; E=8*D0*m*g*L0/k/dping/dping/b0/pingjun; E=E*10000*10000; cout<<"E= "<<E<<endl; double Ue; Ue=(UL/L0)*(UL/L0) +(UD/D0)*(UD/D0) +4*(Ud/dping)*(Ud/dping) +(Ub/b0)*(Ub/b0) +(Ux/pingjun)*(Ux/pingjun); Ue=sqrt(Ue); Ue=E*Ue; cout<<"Ue= "<<Ue<<endl; cout<<"Ue/E= " <<Ue/E<<endl; return 0;}