哈工大材料的力学上机实验课作业.docxVIP

Harbin Institute of Technology 材料力学电算大作业 题目名称：各种条件下的任意截面的弯矩，挠度 作者班号：1108108 作者学号：1110810812 作者姓名：刘涛 指导教师：王春香 完成时间：2013/5/23 一、 问题的提出 该程序可以解决的问题：一个力作用下（1）圆柱面杆（2）矩形截面杆（3）工字梁截面杆（4）空心圆柱截面杆 的弯矩和挠度 二、 程序流程图 输入杆长（mm），力F（kN），力作用点到左端面A的距离a（mm） 输入梁的弹性模量（Gpa） 输入所求端面到左端面的距离X 输入梁的形状 计算惯性矩 计算并输出弯矩 根据梁的形状不同输入不同的数值 计算并输出挠度 三、不同力作用下的弯矩，挠度 1.一个力作用下的任意截面的弯矩，挠度 程序： #include stdio.h #include stdlib.h #include math.h #define PI 3.1415926 int get_I(int kind); int main() { float f=100, l=100, a=50, b, x=60; float M, v; float E = 200, I; int kind=1; printf(请输入杆长l(mm)，力F（kN），力作用点到左端面A的距离a(mm)\n); scanf(%f%f%f,l,f,a); l /= 1000; f *= 1000; a /= 1000; printf(请输入梁的弹性模量(GPa)\n); scanf(%f,E); E *= 1.e9; printf(请输入所求端面到左端面A的距离X\n); scanf(%f,x); x /= 1000; printf(请输入梁的截面形状（圆柱面=1,矩形梁=2,工字梁=3,空心圆柱梁=4）\n); scanf(%d,kind); I = get_I(kind);//计算惯性矩 b = l - a; M = f * x;//计算弯矩 if (0 = x x = a) { v = (-1) * f * b * x *(l*l - x*x - b*b) / (6 * E * I * l); } else if (a x x = l) { v = (-1) * f * b * x * (l/b * pow((x-a), 3) + (l*l - b*b)*x -pow(x, 3)) / (6 * E * I * l); } else { v = 0; printf(error!); exit (0); } printf(转矩M = %gN?m\n, M); printf(挠度V = %gm\n, v); system(pause); return 0; } int get_I(int kind)//计算截面的I惯性矩 { float d,b, h;//矩形梁,工字梁，圆柱梁 float d2, alpha;//空心圆柱梁 float I; //圆柱面=1,矩形梁=2,工字梁=3,空心圆柱梁=4 if (1 == kind)//圆柱面=1 { printf (请输入圆柱梁的直径d\n); scanf(%f, d); I = PI * pow(d, 4) / 64; } else if (2 == kind)//矩形梁=2 { printf (请输入矩形截面的b和h\n); scanf (%f%f,b, h); I = b * pow(h, 3) / 12; } ///* else if (3 == kind)//工字梁=3 { printf(|——————-—————|\n); printf(|————| |——————|t\n); printf( |d|\n); printf( | |\n); printf(|————| |——————|\n);