平面梁结构的内力计算--最终修改版..doc

平面梁结构的内力计算 有一外伸梁结构如图所示，其中 ，集中力，对该梁进行分析，画出弯矩图和剪力图。求出最大剪力和最大弯矩。（材料力学 刘鸿文 P133 -4.10） 图2-1 外伸梁简化图 2、解析法求解 由梁的平衡方程，求出支反力为： 。 梁的A C、C D、Ｄ E、E B等四段内，剪力和弯矩都不能有同一个方程来表示，所以应分为四段考虑。对每一段都可以用同一个方法计算，列出剪力方程和弯矩方程，方程中ｘ以ｍ为单位，以ｋＮ为单位，以ｋＮ为单位。 在AC段内: =-30x (); =-(); 在CD段内： =10 （）=10x-25 (); 在DE段内： =-10 （）=-10x+15 (); 在EB段内： =-30x+120 (); = (); 依照剪力方程和弯矩方程，分段做剪力图和弯矩图。 剪力图： 弯矩图： 从剪力图和弯矩图中可以看出： 最大剪力，最大弯矩为。还可以看出，在集中力截面作用的两侧，剪力有一突然变化，变化的数值就等于集中力，在集中力偶作用的截面的两侧，弯矩有一突然变化，变化的数值就等于集中力偶矩。 3、模型的简化和建立 1.梁的参数设定: 长度 l=4m； 宽度 b= 0.2m； 厚度 h=0.3m 2.材料参数 材料特性应理想弹性体的线性问题的四个假定条件，即：满足完全弹性假定，连续性假定，均匀性假定，各向同性假定的理想弹性体。所以，选择球墨铸铁，它的弹性模量为1.4e11。它的弹性模量EI=140Gpa,泊松比选择 u=0.3。 3.单元选择： 由于题中梁只受均布载荷和剪力，符合梁的2维的单元， BEAM3单元，运用于2维问题，具有拉伸，压缩，弯曲特性，对于在每个节点上有3个自由度x,y方向的平动自由度以及绕z轴的旋转自由度，符合BEAM单元家族中的2D elastic3类型。即为二维梁单元。根据梁的几何参数，所以参数定义为：截面积AREA=，面积惯量Izz=,单元高度HEIGHT=。 4、梁的边界条件 在节点C处梁受X,Y两个方向的约束；节点E受只受Y方向的约束。 5、梁所受的载荷 AC之间作用着均布载荷q=30kN/m，在节点D处作用着集中力=3kN·m，方向为-y方向，所以为负值。 6.节点和单元的设定 由以上分析可知，在ｘ＝1ｍ处，有个固定端和作用着均布载荷q=30kN/m，而x=2m处还作用着一个力. x=3m处有一个滚动支座。 支承点应该取为节点。 所以将X=0,X=1,X=3设置为节点，节点均布，采用直接方法，将梁划分为20个单元，21个节点。 具体建模以及分析过程： 创建节点 2.定义单元类型 选择BEAM单元家族中2D elastic 3类型。 3.定义材料特性 因为选材料是弹性模量prex为1.4e11,泊松比EX=0.3 4.定义几何参数 横截面积AREA=,面积参量IZZ=，单元高度HEIGHT=0.3m. 5.创建单元 i.创建一个单元再Copy使其所有节点连接成20个单元 ii.显示单元资料。 6.施加约束和载荷 i.加载节点自由度约束 ii加载载荷步 7.改变分析类型 8.后处理 i．显示梁变形结果 ii建立元素结果表 创建单元表： iii列出所有表格资料 iV画剪力图 V:画弯矩图 从ANSYS软件分析的剪力图和弯矩图可以看出： 最大剪力： 最小剪力： 最大弯矩： 最小弯矩： 与之前解析法做出的解对比，解正确。 三．机械优化设计 1.问题的提出：用鲍威尔法计算f(X)=的极小值。 2.用解析法求解： 3.程序流程图 4.用C语言的程序代码如下： #include math.h #include conio.h #include stdio.h #include stdlib.h #define MAXDEM 10 //下面以函数12*XI*X1+10*X2*X2-96*X1-100*X2+200为例进行计算 double function(double X[]) //X[]用于存放各个变量的初值，变量个数与自由度数有关 { return ( 12*X[0]*X[0]+10*X[1]*X[1]