材料力学编程Word下载.docx
- 文档编号:19089033
- 上传时间:2023-01-03
- 格式:DOCX
- 页数:26
- 大小:354.60KB
材料力学编程Word下载.docx
《材料力学编程Word下载.docx》由会员分享,可在线阅读,更多相关《材料力学编程Word下载.docx(26页珍藏版)》请在冰豆网上搜索。
,&
l);
printf("
Pleaseinputa/m:
a);
b=l-a;
PleaseinputF/N:
F);
Pleaseinputx/m:
x);
PleaseinputEI:
EI);
if((x>
=a)&
&
(x<
=l))
{
v=(-F*b*((l*pow((x-a),3)/b)-(pow(l,2)-pow(b,2))*x-pow(x,3)))/(6*l);
M=(F*a*(l-x))/l;
}
else
{
v=(-F*b*x*(pow(l,2)-pow(x,2)-pow(b,2)))/(6*l);
M=(F*b*x)/l;
v(%.2lf)=%.2lfmm\nM(%.2lf)=%.2lfkN.m\n"
x,v*1000/EI,x,M/1000);
scanf("
%lf%lf%lf%lf%lf"
l,&
a,&
x,&
F,&
}
一个力偶作用下的任意截面的弯矩,挠度
编译程序
{doublea,b,l,me,x,v,M,EI;
一个力偶作用下的任意截面的弯矩,挠度\n"
Pleaseinputme/N.m:
me);
v=(me*(-pow(x,3)+3*l*pow((x-a),2)+(pow(l,2)-3*pow(b,2))*x))/(6*l);
M=(me*(l-x))/l;
v=(me*x*(pow(l,2)-3*pow(b,2)-pow(x,2)))/(6*l);
M=(-me*x)/l;
me,&
一段分布力作用下的任意截面的
doublea,b,c,l,q,x,v,M,Fa,Fb,c1,c2,c3,d3,EI;
一段分布力作用下的任意截面的弯矩,挠度\n"
Pleaseinputc/m:
c);
b=l-a-c;
Pleaseinputq/N/m:
q);
Fa=(q*c*(c+2*b))/(2*l);
Fb=(q*c*(c+2*a))/(2*l);
c3=(Fa*pow((a+c),3))/(3*l)-(q*a*pow(c,3))/(6*l)-(q*pow(c,4))/(8*l)+(Fb*pow(b,2))/2-(Fb*pow(b,3))/(3*l);
c2=c3-(Fb*pow(b,2))/2-(Fa*pow((a+c),2))/2+(q*pow(c,3))/6;
c1=c2;
d3=-c3*l;
(a+c))&
v=((Fb*pow((l-x),3))/6+c3*x+d3);
M=Fb*(l-x);
a)&
=(a+c)))
v=((Fa*pow(x,3))/6-(q*pow((x-a),4))/24+c2*x);
M=((q*c*x*(c+2*b))/(2*l))-((q*pow((x-a),2))/2);
=0)&
=a))
v=((Fa*pow(x,3))/6+c1*x);
M=Fa*x;
v(%.2lf)=%.2lfmm\n"
x,v*1000/EI);
M(%.2lf)=%.2lfkN.m\n"
x,M/1000);
n个力共同作用下的任意截面的弯矩,挠度
#defineN1000
doubleFunctionV(inti,intn,doublel,doublex,doublea[],doubleb[],doubleF[]);
doubleFunctionM(inti,intn,doublel,doublex,doublea[],doubleb[],doubleF[]);
inti,n;
doublev[N],M[N],a[N],b[N],F[N],sumv=0,sumM=0,l,x,EI;
n个力共同作用下的任意截面的弯矩,挠度\n"
Pleaseinputn:
%d"
n);
for(i=0;
i<
n;
i++)
Pleaseinputa[%d]/m:
i);
a[i]);
b[i]=l-a[i];
printf("
PleaseinputF[%d]/N:
F[i]);
v[i]=FunctionV(i,n,l,x,a,b,F);
M[i]=FunctionM(i,n,l,x,a,b,F);
sumv=sumv+v[i];
sumM=sumM+M[i];
x,sumv*1000/EI);
x,sumM/1000);
doubleFunctionV(inti,intn,doublel,doublex,doublea[],doubleb[],doubleF[])
doublev;
=a[i]))
v=(-F[i]*b[i]*x*(pow(l,2)-pow(x,2)-pow(b[i],2)))/(6*l);
v=(-F[i]*b[i]*((l*pow((x-a[i]),3)/b[i])+(pow(l,2)-pow(b[i],2))*x-pow(x,3)))/(6*l);
returnv;
doubleFunctionM(inti,intn,doublel,doublex,doublea[],doubleb[],doubleF[])
doubleM;
M=(F[i]*b[i]*x)/l;
M=(F[i]*a[i]*(l-x))/l;
returnM;
h个力偶共同作用作用下的任意截面的弯矩,挠度
doubleFunctionV(inti,intn,doublel,doublex,doublea[],doubleb[],doubleme[]);
doubleFunctionM(inti,intn,doublel,doublex,doublea[],doubleb[],doubleme[]);
inti,h;
doublel,x,EI,v[N],M[N],a[N],b[N],me[N],sumv=0,sumM=0;
h个力偶共同作用下的任意截面的弯矩,挠度\n"
Pleaseinputh:
h);
h;
Pleaseinputme[%d]/N.m:
me[i]);
v[i]=FunctionV(i,h,l,x,a,b,me);
M[i]=FunctionM(i,h,l,x,a,b,me);
doubleFunctionV(inti,inth,doublel,doublex,doublea[],doubleb[],doubleme[])
v=(me[i]*x*(pow(l,2)-3*pow(b[i],2)-pow(x,2)))/(6*l);
v=(me[i]*(-pow(x,3)+3*l*pow((x-a[i]),2)+(pow(l,2)-3*pow(b[i],2))*x))/(6*l);
doubleFunctionM(inti,inth,doublel,doublex,doublea[],doubleb[],doubleme[])
M=(-me[i]*x)/l;
M=(me[i]*(l-x))/l;
m段分布力共同作用下的任意截面的弯矩,挠度
doubleFunctionV(inti,intm,doublel,doublex,doublea[],doubleb[],doublec[],doubleq[],doubleFa[],doubleFb[],doublec1[],doublec2[],doublec3[],doubled3[]);
doubleFunctionM(inti,intm,doublel,doublex,doublea[],doubleb[],doublec[],doubleq[],doubleFa[],doubleFb[],doublec1[],doublec2[],doublec3[],doubled3[]);
inti,m;
doublel,x,EI,v[N],M[N],a[N],b[N],c[N],q[N],Fa[N],Fb[N],c1[N],c2[N],c3[N],d3[N],sumv=0,sumM=0;
m段分布力共同作用下的任意截面的弯矩,挠度\n"
Pleaseinputm:
m);
m;
Pleaseinputc[%d]/m:
c[i]);
b[i]=l-a[i]-c[i];
Pleaseinputq[%d]/N/m:
q[i]);
Fa[i]=(q[i]*c[i]*(c[i]+2*b[i]))/(2*l);
Fb[i]=(q[i]*c[i]*(c[i]+2*a[i]))/(2*l);
c3[i]=(Fa[i]*pow((a[i]+c[i]),3))/(3*l)-(q[i]*a[i]*pow(c[i],3))/(6*l)-(q[i]*pow(c[i],4))/(8*l)+(Fb[i]*pow(b[i],2))/2-(Fb[i]*pow(b[i],3))/(3*l);
c2[i]=c3[i]-(Fb[i]*pow(b[i],2))/2-(Fa[i]*pow((a[i]+c[i]),2))/2+(q[i]*pow(c[i],3))/6;
c1[i]=c2[i];
d3[i]=-c3[i]*l;
v[i]=FunctionV(i,m,l,x,a,b,c,q,Fa,Fb,c1,c2,c3,d3);
M[i]=FunctionM(i,m,l,x,a,b,c,q,Fa,Fb,c1,c2,c3,d3);
doubleFunctionV(inti,intm,doublel,doublex,doublea[],doubleb[],doublec[],doubleq[],doubleFa[],doubleFb[],doublec1[],doublec2[],doublec3[],doubled3[])
v=((Fa[i]*pow(x,3))/6+c1[i]*x);
a[i])&
=(a[i]+c[i])))
v=((Fa[i]*pow(x,3))/6-(q[i]*pow((x-a[i]),4))/24+c2[i]*x);
(a[i]+c[i]))&
v=((Fb[i]*pow((l-x),3))/6+c3[i]*x+d3[i]);
returnv;
doubleFunctionM(inti,intm,doublel,doublex,doublea[],doubleb[],doublec[],doubleq[],doubleFa[],doubleFb[],doublec1[],doublec2[],doublec3[],doubled3[])
M=Fa[i]*x;
M=((q[i]*c[i]*x*(c[i]+2*b[i]))/(2*l))-((q[i]*pow((x-a[i]),2))/2);
M=Fb[i]*(l-x);
returnM;
n个力,h个力偶,m段分布力共同作用下的任意截面的弯矩,挠度
doubleFunctionV1(inti,intn,doublel,doublex,doublea[],doubleb[],doubleF[]);
doubleFunctionM1(inti,intn,doublel,doublex,doublea[],doubleb[],doubleF[]);
doubleFunctionV2(inti,intn,doublel,doublex,doublea[],doubleb[],doubleme[]);
doubleFunctionM2(inti,intn,doublel,doublex,doublea[],doubleb[],doubleme[]);
doubleFunctionV3(inti,intm,doublel,doublex,doublea[],doubleb[],doublec[],doubleq[],doubleFa[],doubleFb[],doublec1[],doublec2[],doublec3[],doubled3[]);
doubleFunctionM3(inti,intm,doublel,doublex,doublea[],doubleb[],doublec[],doubleq[],doubleFa[],doubleFb[],doublec1[],doublec2[],doublec3[],doubled3[]);
inti,n,h,m;
doublel,x,EI,a[N],b[N],v[N],M[N],sumv=0,sumM=0;
doubleF[N],sumv1=0,sumM1=0;
doubleme[N],sumv2=0,sumM2=0;
doublec[N],q[N],Fa[N],Fb[N],c1[N],c2[N],c3[N],d3[N],sumv3=0,sumM3=0;
请输入:
n,h,m,l/m,x/m,EI:
"
%d%d%d%lf%lf%lf"
n,&
h,&
m,&
n=%d,h=%d,m=%d,l=%lfm,x=%lfm,EI%lf\n"
n,h,m,l,x,EI);
计算%d个力共同作用下的任意截面的弯矩,挠度\n"
n);
v[i]=FunctionV1(i,n,l,x,a,b,F);
M[i]=FunctionM1(i,n,l,x,a,b,F);
sumv1=sumv1+v[i];
sumM1=sumM1+M[i];
力共同作用结果为:
v(%.2lf)=%.2lfmmM(%.2lf)=%.2lfkN.m\n"
x,sumv1*1000/EI,x,sumM1/1000);
计算%d个力偶共同作用下的任意截面的弯矩,挠度\n"
h);
v[i]=FunctionV2(i,h,l,x,a,b,me);
M[i]=FunctionM2(i,h,l,x,a,b,me);
sumv2=sumv2+v[i];
sumM2=sumM2+M[i];
力偶共同作用结果为:
x,sumv2*1000/EI,x,sumM2/1000);
计算%d段分布力共同作用下的任意截面的弯矩,挠度\n"
m);
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 材料力学 编程
![提示](https://static.bdocx.com/images/bang_tan.gif)