哈工大有限元大作业

机械有限元习题答案——哈工大

机械有限元习题答案——哈⼯⼤第⼆章习题2.1 解释如下的概念：应⼒、应变，⼏何⽅程、物理⽅程、虚位移原理。

解○1应⼒是某截⾯上的应⼒在该处的集度。

○2 应变是指单元体在某⼀个⽅向上有⼀个ΔU 的伸长量，其相对变化量就是应变。

X U Xx ??=ε表⽰在x 轴的⽅向上的正应变，其包括正应变和剪应变。

○3⼏何⽅程是表⽰弹性体内节点的应变分量与位移分量之间的关系，其完整表⽰如下：Txz yz xy z y x x w z u zv y w y u x v z w y vx u x w z u z v y w y u x v z w y v x u ??+++=+????+????+?=??????=γγγεεεε○4物理⽅程：表⽰应⼒和应变关系的⽅程某⼀点应⼒分量与应变分量之间的关系如下：=???????????????????=666564636261565554535251464545434241363534333231262524232221161514131211αααααααααααααααααααααααααααααααααααατττσσσσxz yz xy z y xxz yz xy zz yy xx γγγεεε○5虚位移原理：在弹性有⼀虚位移情况下，由于作⽤在每个质点上的⼒系，在相应的虚位移上虚功总和为零，即为：若弹性体在已知的⾯⼒和体⼒的作⽤下处于平衡状态，那么使弹性体产⽣虚位移，所有作⽤在弹性体上的体⼒在虚位移上所做的⼯就等于弹性体所具有的虚位能。

2.2说明弹性体⼒学中的⼏个基本假设。

○1 连续性假设：就是假定整个物体的体积都被组成该物体的介质所填满，不存在任何间隙。

○2 完全弹性假设：就是假定物体服从虎克定律。

○3 各向同性假设：就是假定整个物体是由同意材料组成的。

○4 ⼩变形和⼩位移假设：就是指物体各点的位移都远远⼩于物体原来的尺⼨，并且其应变和转⾓都⼩于1。

2.3简述线应变与剪应变的⼏何含义。

哈工大机械原理大作业1 连杆机构运动分析牛头刨床学号5号

大作业1 连杆机构运动分析运动分析题目图1-13 牛头刨床机构简图一、按比例画出机构运动简图见A3图纸二、机构的结构分析及基本杆组的划分（1）机构的结构分析该机构为平面连杆机构，活动构件数n=5；有A、B、C、E、F五个转动副和B、F两个移动副，没有高副，故低副数PL=7，高副数PH=0。

机构的自由度F=3n-2PL-PH=3*5-2*7=1。

机构有确定运动轨迹。

（2）基本杆组的划分图1-13中，原动件构件1为Ⅰ级杆组，移除之后按照拆杆组法对机构进行拆分，可得由滑块2和杆3组成的RPRⅡ级杆组，以及由杆4和滑块5组成的RRPⅡ级杆组。

拆分图如下：图1-13-1 Ⅰ级杆组图1-13-2 RPRⅡ级杆组图1-13-3 RRPⅡ级杆组三、各杆组的运动分析数学模型图1-13-1 Ⅰ级杆组 x B =x A +l AB cos ϕy B =y A +l AB sin ϕ图1-13-2 RPR Ⅱ级杆组 x E =x B +(l CE -s)cos ϕj y E =y B +(l CE -s)sin ϕjϕj =arctan(B 0/A 0)A 0=xB -x CB 0=y B -y CS =2020B A + 图1-13-3 RRP Ⅱ级杆组F 的位置 x F =x E -l EF cos ϕiy F =y E +l EF sin ϕi ϕi =arcsin(A 1/l EF )A 1=H-H 1-y EF 的速度F 的加速度四、建立坐标系以A 为原点建立坐标系，则A （0,0）五、编程使用MATLAB 编写程序如下lAB=108;lCD=620;lEF=300;H1=350;H=635;xC=0;yC=-350;syms t;f=(255*pi/30)*t;xB=lAB*cos(f);yB=lAB*sin(f);A0=xB-xC;B0=yB-yC;S=sqrt(A0^2+B0^2);f1=atan(B0/A0);xE=xB+(lAB-S)*cos(f1);yE=yB+(lAB-S)*sin(f1);A1=H-H1-yE;f2=asin(A1/lEF);xF=xE-lEF*cos(f2);yF=yE+lEF*sin(f2);a=0:0.001:0.5;xF=subs(xF,t,a);plot(a,xF)title('位移x随时间t变化图像') xlabel('t(s)'),ylabel(X')lAB=108;lCD=620;lEF=300;H1=350;H=635;xC=0;yC=-350;syms t;f=(255*pi/30)*t;xB=lAB*cos(f);yB=lAB*sin(f);A0=xB-xC;B0=yB-yC;S=sqrt(A0^2+B0^2);f1=atan(B0/A0);xE=xB+(lAB-S)*cos(f1);yE=yB+(lAB-S)*sin(f1);A1=H-H1-yE;f2=asin(A1/lEF);xF=xE-lEF*cos(f2);yF=yE+lEF*sin(f2);vF=diff(xF,t);a=0:0.001:0.5;vF=subs(vF,t,a);plot(a,vF)title('速度v随时间t变化图像') xlabel('t(s)'),ylabel(V')lAB=108;lCD=620;lEF=300;H1=350;H=635;xC=0;yC=-350;syms t;f=(255*pi/30)*t;xB=lAB*cos(f);yB=lAB*sin(f);A0=xB-xC;B0=yB-yC;S=sqrt(A0^2+B0^2);f1=atan(B0/A0);xE=xB+(lAB-S)*cos(f1);yE=yB+(lAB-S)*sin(f1);A1=H-H1-yE;f2=asin(A1/lEF);xF=xE-lEF*cos(f2);yF=yE+lEF*sin(f2);aF=diff(xF,t,2);a=0:0.001:0.5;aF=subs(aF,t,a);plot(a,aF)title('加速度a随时间t变化图像') xlabel('t(s)'),ylabel(A')六、位移、速度、加速度图像位移xE随时间t变化图像速度vE随时间t变化图像加速度aE随时间t变化图像七、运动分析由MATLAB编程后绘制所得图线分析可知，牛头刨床的刨头F行程约为500mm，在0.06s~0.17s之间速度和加速度波动较大，在0.17s~0.29s 之间速度和加速度波动较小，这就保证了刨头在空行程时有急回运动，而在工作行程是由较均匀的切削速度，运动连续，且运动性能良好稳定。

哈工大有限元大作业

作业一一．计算程序和结果展示1.程序clearsyms a b c x E lD=E*pi*(b*x+c)^4/64;B(:,:,1)=[-6/l^2+12*x/l^3 4/l-6*x/l^2];B(:,:,2)=[6/l^2-12*x/l^3 ,2/l-6*x/l^2];n=0;for i=1:2for j=1:2n=n+1;f=B(:,:,i)*D*transpose(B(:,:,j));k(:,:,n)=int(f,x,0,l);endendk11=k(:,:,1);k12=k(:,:,2);k21=k(:,:,3);k22=k(:,:,4);K=[k11 k12;k21 k22];K=simple(K);2.结果（1）bx+cK =[ (3*pi*E*(11*b^4*l^4 + 49*b^3*c*l^3 + 84*b^2*c^2*l^2 + 70*b*c^3*l +35*c^4))/(560*l^3), -(pi*E*(19*b^4*l^4 + 84*b^3*c*l^3 + 147*b^2*c^2*l^2 + 140*b*c^3*l + 105*c^4))/(1120*l^2), -(3*pi*E*(11*b^4*l^4 + 49*b^3*c*l^3 + 84*b^2*c^2*l^2 + 70*b*c^3*l + 35*c^4))/(560*l^3), -(pi*E*(47*b^4*l^4 + 210*b^3*c*l^3 + 357*b^2*c^2*l^2 + 280*b*c^3*l + 105*c^4))/(1120*l^2)][ -(pi*E*(19*b^4*l^4 + 84*b^3*c*l^3 + 147*b^2*c^2*l^2 + 140*b*c^3*l +105*c^4))/(1120*l^2), (pi*E*(3*b^4*l^4 + 14*b^3*c*l^3 + 28*b^2*c^2*l^2 +35*b*c^3*l + 35*c^4))/(560*l), (pi*E*(19*b^4*l^4 + 84*b^3*c*l^3 + 147*b^2*c^2*l^2 +140*b*c^3*l + 105*c^4))/(1120*l^2), (pi*E*(13*b^4*l^4 + 56*b^3*c*l^3 +91*b^2*c^2*l^2 + 70*b*c^3*l + 35*c^4))/(1120*l)][ -(3*pi*E*(11*b^4*l^4 + 49*b^3*c*l^3 + 84*b^2*c^2*l^2 + 70*b*c^3*l +35*c^4))/(560*l^3), (pi*E*(19*b^4*l^4 + 84*b^3*c*l^3 + 147*b^2*c^2*l^2 + 140*b*c^3*l +105*c^4))/(1120*l^2), (3*pi*E*(11*b^4*l^4 + 49*b^3*c*l^3 + 84*b^2*c^2*l^2 +70*b*c^3*l + 35*c^4))/(560*l^3), (pi*E*(47*b^4*l^4 + 210*b^3*c*l^3 + 357*b^2*c^2*l^2 + 280*b*c^3*l + 105*c^4))/(1120*l^2)][ -(pi*E*(47*b^4*l^4 + 210*b^3*c*l^3 + 357*b^2*c^2*l^2 + 280*b*c^3*l +105*c^4))/(1120*l^2), (pi*E*(13*b^4*l^4 + 56*b^3*c*l^3 + 91*b^2*c^2*l^2 +70*b*c^3*l + 35*c^4))/(1120*l), (pi*E*(47*b^4*l^4 + 210*b^3*c*l^3 + 357*b^2*c^2*l^2 +280*b*c^3*l + 105*c^4))/(1120*l^2), (pi*E*(17*b^4*l^4 + 77*b^3*c*l^3 +133*b^2*c^2*l^2 + 105*b*c^3*l + 35*c^4))/(560*l)]（2）ax^2+bx+c（将程序中D的直径换成“ax^2+bx+c”）K =[ (pi*E*(518*a^4*l^8+2233*a^3*b*l^7+2420*a^3*c*l^6+3630*a^2*b^2*l^6+7920*a^2*b*c*l^5 + 4356*a^2*c^2*l^4 +2640*a*b^3*l^5+8712*a*b^2*c*l^4+9702*a*b*c^2*l^3+3696*a*c^3*l^2 + 726*b^4*l^4 + 3234*b^3*c*l^3 + 5544*b^2*c^2*l^2 + 4620*b*c^3*l +2310*c^4))/(12320*l^3), -(pi*E*(938*a^4*l^8+4004*a^3*b*l^7+4290*a^3*c*l^6+6435*a^2*b^2*l^6+13860*a^2*b*c*l^ 5+7524*a^2*c^2*l^4+4620*a*b^3*l^5+15048*a*b^2*c*l^4+16632*a*b*c^2*l^3+6468*a*c^3*l ^2 +1254*b^4*l^4+5544*b^3*c*l^3+9702*b^2*c^2*l^2+9240*b*c^3*l+6930*c^4))/(73920*l^2), -(pi*E*(518*a^4*l^8 + 2233*a^3*b*l^7 + 2420*a^3*c*l^6 + 3630*a^2*b^2*l^6 +7920*a^2*b*c*l^5 + 4356*a^2*c^2*l^4 + 2640*a*b^3*l^5 + 8712*a*b^2*c*l^4 +9702*a*b*c^2*l^3 + 3696*a*c^3*l^2 + 726*b^4*l^4 + 3234*b^3*c*l^3 + 5544*b^2*c^2*l^2 + 4620*b*c^3*l + 2310*c^4))/(12320*l^3), -(pi*E*(2170*a^4*l^8 + 9394*a^3*b*l^7 +10230*a^3*c*l^6 + 15345*a^2*b^2*l^6 + 33660*a^2*b*c*l^5 + 18612*a^2*c^2*l^4 +11220*a*b^3*l^5 + 37224*a*b^2*c*l^4 + 41580*a*b*c^2*l^3 + 15708*a*c^3*l^2 +3102*b^4*l^4 + 13860*b^3*c*l^3 + 23562*b^2*c^2*l^2 + 18480*b*c^3*l +6930*c^4))/(73920*l^2)][ -(pi*E*(938*a^4*l^8 + 4004*a^3*b*l^7 + 4290*a^3*c*l^6 + 6435*a^2*b^2*l^6 + 13860*a^2*b*c*l^5 + 7524*a^2*c^2*l^4 + 4620*a*b^3*l^5 + 15048*a*b^2*c*l^4 +16632*a*b*c^2*l^3 + 6468*a*c^3*l^2 + 1254*b^4*l^4 + 5544*b^3*c*l^3 + 9702*b^2*c^2*l^2 + 9240*b*c^3*l + 6930*c^4))/(73920*l^2), (pi*E*(434*a^4*l^8 + 1848*a^3*b*l^7 + 1980*a^3*c*l^6 + 2970*a^2*b^2*l^6 + 6435*a^2*b*c*l^5 + 3564*a^2*c^2*l^4 +2145*a*b^3*l^5 + 7128*a*b^2*c*l^4 + 8316*a*b*c^2*l^3 + 3696*a*c^3*l^2 + 594*b^4*l^4 + 2772*b^3*c*l^3 + 5544*b^2*c^2*l^2 + 6930*b*c^3*l + 6930*c^4))/(110880*l),(pi*E*(938*a^4*l^8 + 4004*a^3*b*l^7 + 4290*a^3*c*l^6 + 6435*a^2*b^2*l^6 +13860*a^2*b*c*l^5 + 7524*a^2*c^2*l^4 + 4620*a*b^3*l^5 + 15048*a*b^2*c*l^4 +16632*a*b*c^2*l^3 + 6468*a*c^3*l^2 + 1254*b^4*l^4 + 5544*b^3*c*l^3 + 9702*b^2*c^2*l^2 + 9240*b*c^3*l + 6930*c^4))/(73920*l^2), (pi*E*(1946*a^4*l^8 + 8316*a^3*b*l^7 +8910*a^3*c*l^6 + 13365*a^2*b^2*l^6 + 28710*a^2*b*c*l^5 + 15444*a^2*c^2*l^4 +9570*a*b^3*l^5 + 30888*a*b^2*c*l^4 + 33264*a*b*c^2*l^3 + 12012*a*c^3*l^2 +2574*b^4*l^4 + 11088*b^3*c*l^3 + 18018*b^2*c^2*l^2 + 13860*b*c^3*l +6930*c^4))/(221760*l)][ -(pi*E*(518*a^4*l^8 + 2233*a^3*b*l^7 + 2420*a^3*c*l^6 +3630*a^2*b^2*l^6 + 7920*a^2*b*c*l^5 + 4356*a^2*c^2*l^4 + 2640*a*b^3*l^5 +8712*a*b^2*c*l^4 + 9702*a*b*c^2*l^3 + 3696*a*c^3*l^2 + 726*b^4*l^4 + 3234*b^3*c*l^3 +5544*b^2*c^2*l^2 + 4620*b*c^3*l + 2310*c^4))/(12320*l^3), (pi*E*(938*a^4*l^8 + 4004*a^3*b*l^7 + 4290*a^3*c*l^6 + 6435*a^2*b^2*l^6 + 13860*a^2*b*c*l^5 +7524*a^2*c^2*l^4 + 4620*a*b^3*l^5 + 15048*a*b^2*c*l^4 + 16632*a*b*c^2*l^3 +6468*a*c^3*l^2 + 1254*b^4*l^4 + 5544*b^3*c*l^3 + 9702*b^2*c^2*l^2 + 9240*b*c^3*l + 6930*c^4))/(73920*l^2), (pi*E*(518*a^4*l^8 + 2233*a^3*b*l^7 +2420*a^3*c*l^6 + 3630*a^2*b^2*l^6 + 7920*a^2*b*c*l^5 + 4356*a^2*c^2*l^4 +2640*a*b^3*l^5 + 8712*a*b^2*c*l^4 + 9702*a*b*c^2*l^3 + 3696*a*c^3*l^2 + 726*b^4*l^4 + 3234*b^3*c*l^3 + 5544*b^2*c^2*l^2 + 4620*b*c^3*l + 2310*c^4))/(12320*l^3),(pi*E*(2170*a^4*l^8 + 9394*a^3*b*l^7 + 10230*a^3*c*l^6 + 15345*a^2*b^2*l^6 +33660*a^2*b*c*l^5 + 18612*a^2*c^2*l^4 + 11220*a*b^3*l^5 + 37224*a*b^2*c*l^4 +41580*a*b*c^2*l^3 + 15708*a*c^3*l^2 + 3102*b^4*l^4 + 13860*b^3*c*l^3 +23562*b^2*c^2*l^2 + 18480*b*c^3*l + 6930*c^4))/(73920*l^2)][ -(pi*E*(2170*a^4*l^8 + 9394*a^3*b*l^7 + 10230*a^3*c*l^6 + 15345*a^2*b^2*l^6 +33660*a^2*b*c*l^5 + 18612*a^2*c^2*l^4 + 11220*a*b^3*l^5 + 37224*a*b^2*c*l^4 +41580*a*b*c^2*l^3 + 15708*a*c^3*l^2 + 3102*b^4*l^4 + 13860*b^3*c*l^3 +23562*b^2*c^2*l^2 + 18480*b*c^3*l + 6930*c^4))/(73920*l^2), (pi*E*(1946*a^4*l^8 +8316*a^3*b*l^7 + 8910*a^3*c*l^6 + 13365*a^2*b^2*l^6 + 28710*a^2*b*c*l^5 +15444*a^2*c^2*l^4 + 9570*a*b^3*l^5 + 30888*a*b^2*c*l^4 + 33264*a*b*c^2*l^3 +12012*a*c^3*l^2 + 2574*b^4*l^4 + 11088*b^3*c*l^3 + 18018*b^2*c^2*l^2 + 13860*b*c^3*l + 6930*c^4))/(221760*l), (pi*E*(2170*a^4*l^8 + 9394*a^3*b*l^7 + 10230*a^3*c*l^6 +15345*a^2*b^2*l^6 + 33660*a^2*b*c*l^5 + 18612*a^2*c^2*l^4 + 11220*a*b^3*l^5 +37224*a*b^2*c*l^4 + 41580*a*b*c^2*l^3 + 15708*a*c^3*l^2 + 3102*b^4*l^4 +13860*b^3*c*l^3 + 23562*b^2*c^2*l^2 + 18480*b*c^3*l + 6930*c^4))/(73920*l^2),(pi*E*(2282*a^4*l^8 + 9933*a^3*b*l^7 + 10890*a^3*c*l^6 + 16335*a^2*b^2*l^6 +36135*a^2*b*c*l^5 + 20196*a^2*c^2*l^4 + 12045*a*b^3*l^5 + 40392*a*b^2*c*l^4 +45738*a*b*c^2*l^3 + 17556*a*c^3*l^2 + 3366*b^4*l^4 + 15246*b^3*c*l^3 +26334*b^2*c^2*l^2 + 20790*b*c^3*l + 6930*c^4))/(110880*l)]（3）b=0（等截面梁)K =[ (3*pi*E*c^4)/(16*l^3),-(3*pi*E*c^4)/(32*l^2),-(3*pi*E*c^4)/(16*l^3),-(3*pi*E*c^4)/(32*l^2)] [ -(3*pi*E*c^4)/(32*l^2), (pi*E*c^4)/(16*l), (3*pi*E*c^4)/(32*l^2), (pi*E*c^4)/(32*l)][ -(3*pi*E*c^4)/(16*l^3), (3*pi*E*c^4)/(32*l^2), (3*pi*E*c^4)/(16*l^3),(3*pi*E*c^4)/(32*l^2)] [ -(3*pi*E*c^4)/(32*l^2), (pi*E*c^4)/(32*l), (3*pi*E*c^4)/(32*l^2), (pi*E*c^4)/(16*l)]总结：将结果（3）与教材等截面梁刚度矩阵比较，发现表达式一样，侧面证明了程序的正确性。

哈工大机械原理大作业连杆机构9

机械原理大作业大作业一：连杆机构运动分析学生姓名：学号：指导教师：丁刚完成时间：机电工程学院机械设计系制二〇一八年四月连杆机构运动分析1题目（9）图1 设计题目在图1所示的机构中，已知l AB=60mm，l BC=180mm，l DE=200mm，l CD=120mm，l EF=300mm，h=80mm，h1=85mm，h2=225mm，构件1以等角速度ω1=100rad/s 转动。

求在一个运动循环中，滑块5的位移、速度和加速度曲线。

2分析结构1、杆1为RR主动件，绕A以ω1 转动，自由度1.2、4杆和滑块5为RRP II级杆组.，自由度0.3、2，3杆组成II级杆组RRR，自由度0.总共自由度为F=5*3-2*7=1 .由上述的杆组类型，确认出所需运动分析数学模型：同一构件上的点、RRP、RRR。

3.杆组法对平面连杆机构进行运动分析3.1对主动件杆1 RR I级构件的分析主动杆1转角：φ= [0°,360°) δ=0°,则φ’=ω1=100 rad/s角加速度φ’’=0.已知h2=225mm, h=80mm, l AB=60mm 所以A（225mm,80mm）A点速度(0,0),加速度(0,0)B点位置(x A+l AB*cos(φ), Y A+l AB*sin(φ))B点速度(-l AB*sin(φ), l AB*cos(φ)),加速度(-l AB*cos(φ), -l AB*sin(φ))3.2RRRII 级杆组分析（模型参考教材P37-38）图3 如图所示两个构件组成II 级杆组。

已知了B 的位置（x B ，y B ）= (x A +l AB *cos(φ), Y A +l AB *sin(φ)),速度（x ’B ，y ’B ）和加速度（x ’’B ,y ’’B ）, 已知运动副D （0,0），还可知，x ’D =y ’D =0, x ’’D =y ’’D =0. l BC =180 mm, l CD = 120mm所以，x c =x D +l CD *cos(φi)= x B +l BC *cos(φj) y c =x D +l CD *sin(φi)= x B +l BC *sin(φi) 对于φ的求解： A 0=2*l CD （x B -x D ） B 0=2*l BC （y B -y D ） C 0=l CD 2+ l BD 2- l BC 2为了保证机构的装配正常：l BD ≤l CD + l BC AND l BD ≥Abs （l CD - l BC ）可求3杆的转角φi=2*arctan(（B 0±sqrt （A 02 + B 02- C 02））/（A 02+ C 02)），角速度w3=φi ’和角加速度α3= φi ’’3.3 同一构件上的点(模型参考书P35-36)Φiφjφi已知D(0,0),速度(0,0),加速度(0,0),3杆转角φi 角速度φi’角加速度φi’’，Φi和它的导数在3.2都有体现LDE= 200mm可求出E的坐标,速度,加速度.x E =x c+lCE*cos(φi)y E =x C+lCE*sin(φi)同样地，速度、加速度通过求导即可得出算式，可以编出程序。

哈工大研究生塑性体积成形大作业2

基于元胞自动机法的微观组织模拟方法及其在热加工中的应用学院：材料科学与工程学院姓名：孙巍学号：14S109063摘要：建立微观组织演变模型，模拟晶粒形核生长过程及组织分布和偏析情况在近年来已经成为国内外科研工作者的研究热点。

在微观组织的模拟方法中，常用的有蒙特卡罗（MC）法和元胞自动机（CA）法。

CA法由于其模型结构简单，在空间和时间尺度上都不受限制，适合于大型机的并行计算，而受到更多研究者的青睐。

研究表明CA法在金属凝固，静态再结晶及轧制变形所引起的动态再结晶等过程的模拟中，表现出了其特有的优越性。

本文概述了元胞自动机（CA）方法在微观组织模拟中的建模方法及应用，并综述了目前此方面研究工作的进展现状。

关键词：微观组织；模拟；元胞自动机；应用1.1CA的发展史元胞自动机方法这一思想最早由计算机创始人，著名数学家V.Neumann[1]在50年代提出。

它是一种用来描述复杂系统在离散空间一时间上演化规律的数学算法。

最初应用于生物体发育中细胞的自我复制。

为了避免当时电子管计算机的限制，他提出了一个简单的模式，把一个长方形平面划分为若干个网络，每一个格点表示一个细胞或系统的基元，把他们的状态赋值为0或1，在网格中用空格或实格来表示。

在预先设定的规则下，细胞或基元的演化就用网格中的空格或实格的变动描述。

该方法称为元胞自动机法。

1970年，剑桥大学的J.H.Conway 利用元胞自动机法编制了一个名为“生命”的游戏程序，并由M.Gardner[2]通过《科学美国人》介绍到全世界。

从80年代开始，就有人对金属凝固过程中结晶组织的形成用元胞自动机法来进行模拟研究，到1991年，Hesselbarth和Gobel[3]才首先提出将元胞自动机用于再结晶模拟的研究。

目前在材料科学领域中主要应用在如下几个方面：1、在凝固过程[4-5]中，元胞自动机方法己经用于枝晶晶粒结构的模拟，并用于处理晶粒形核与生长竞争等方面的问题：2、在再结晶过程[6-12]的模拟中，元胞自动机能模拟出再结晶核心的形成、生长，位错密度的演化，以及预测出应力-应变曲线：3、在维持动力学的模拟中，用于模拟位错样式的形成等：4、在热处理相变过程中[13]，模拟新相晶粒形核与生长过程。

哈工大机械原理大作业19齿轮.docx

机械原理大作业（三）作业名称：齿轮传动设计设计题目：院系：机电工程学院班级：1108108设计者：许彬彬学号：1110810816指导教师：林琳设计时间：2013年6月1日哈尔滨工业大学机械设计1、设计题目 1.1机构运动简图1.2机械传动系统原始参数电动机转速min /1450r n =，输出转速min/4001r n =，m in /3502r n =，m in /3003r n =，带传动的最大传动比5.2m ax =p i ,滑移齿轮传动的最大传动比4m ax =v i ，定轴齿轮传动的最大传动比4m ax =d i 。

根据传动系统的原始参数可知，传动系统的总传动比为：25.36401450011===n n i43.41351450022===n n i33.48301450033===n n i传动系统的总传动比由带传动、滑移齿轮传动和定轴齿轮传动三部分实现。

设带传动的传动比为5.2m ax =p i ，滑移齿轮的传动比为321v v v i i i 、、，定轴齿轮传动的传动比为f i，则总传动比f v p i i i i 1m ax 1= f v p i i i i 2max 2= fv p i i i i 3max 3= 令4max 3==v v i i则可得定轴齿轮传动部分的传动比为833.44*5.233.48max max 3===v p f i i i i滑移齿轮传动的传动比为000.3833.4*5.225.36max 11===fp v i i i i429.3833.4*5.243.41max 22===fp v i i i i设定轴齿轮传动由3对齿轮传动组成，则每对齿轮的传动比为4691.1833.4max 33=≤===d f d i i i 3、齿轮齿数的确定根据滑移齿轮变速传动系统中对齿轮齿数的要求，可大致选择齿轮5、6、7、8、9和10为角度变位齿轮，其齿数：37,9,35,10,33,111098765======z z z z z z ；它们的齿顶高系数1=*a h ，径向间隙系数25.0=*c ，分度圆压力角020=α，实际中心距mm a 46'=。

(完整word版)哈工大机械原理大作业凸轮DOC

H a r b i n I n s t i t u t e o f T e c h n o l o g y机械原理大作业二课程名称：机械原理设计题目: 凸轮机构设计院系：班级：设计者：学号：指导教师：哈尔滨工业大学一、设计题目如右图所示直动从动件盘形凸轮机构,选择一组凸轮机构的原始参数,据此设计该凸轮机构。

凸轮机构原始参数序号升程（mm）升程运动角升程运动规律升程许用压力角27130150正弦加速度30°回程运动角回程运动规律回程许用压力角远休止角近休止角100°余弦加速度60°30°80°二. 凸轮推杆升程、回程运动方程及推杆位移、速度、加速度线图凸轮推杆升程运动方程:)]512sin(2156[130s ϕππϕ-= )512sin(4.374)]512cos(1[156v 211ϕπϕπωω=-=a％ t 表示转角,s 表示位移t=0：0.01:5*pi/6;%升程阶段s= [(6＊t)/(5*pi ）- 1/（2＊pi ）*sin(12＊t/5）］＊130; hold on plot(t ，s ）; t= 5＊pi/6:0。

01:pi; ％远休止阶段s=130； hold on plot(t,s ）;t=pi ：0.01:14*pi/9；%回程阶段s=65*［1+cos(9*（t-pi ）/5)]； hold on plot(t ，s ）;t=14*pi/9：0.01:2＊pi ；s=0；hold onplot（t,s）；grid onhold off％t表示转角,令ω1=1t=0：0。

01:5*pi/6;%升程阶段v=156*1*［1-cos（12＊t/5）]/pi hold onplot(t，v)；t= 5＊pi/6：0。

01:pi;v=0hold onplot(t，v)；t=pi:0.01:14＊pi/9;％回程阶段v=—117*1*sin（9＊(t—pi)/5) hold onplot（t，v)；t=14＊pi/9:0。

哈尔滨工业大学 04 结构力学——平面有限元1

1 Nj 2A(aj bjxcjy)
单元形函数
1
Nk2A(akbkxcky)

哈工大土木工程学院
26
第五章平面问题有限元分析
uN iu iN jujN ku k
vN iviN jvjN kvk
u v

Ni 0
0 Ni
Nj 0
0 Nj
Nk 0
dNe
三角形单元共有6个结点位移，所以有6个自由度。故选位移模式：
u12x3y v46x6y
d

u
v

1 0
x 0
y 0
结点坐标矩阵
N0
哈工大土木工程学院
待定系数
广义坐标
1

2

0 1
0 x
0 y

3 4
5-1-2 平面问题的总势能表达式
E P1 2A T d A A F b T d d A
S F S T d d S F e T e
哈工大土木工程学院

16
第五章平面问题有限元分析 §5-2 常应变三角形单元
5-2-1 单元结点位移和结点力
三角形单元
12N2u
yk3
ve
vvij
1 1
xi xj
vk 1 xk
yi yj
54N2v
yk6
1 xi D 1 xj
1 xk
yi yj 2A yk
i j k 保证三角形面积为正
zx 0
zy 0
哈工大土木工程学院
三个应变分量只发生在 xy面上，故称平面应变问题。

5

哈工大机械原理大作业1

11.440890
420.636249
198
-25.780027
-7.408008
448.711622
19
17.496096
11.511568
389.280068
199
-27.066134
-7.329197
454.399338
20
19.511173
11.576658
356.588380
200
-28.338399
-6.837081
484.309157
26
31.792111
11.837106
133.586249
206
-35.672384
-6.752179
488.602299
27
33.860111
11.856813
92.240035
207
-36.843419
-6.666544
492.701606
28
35.930920
（4）利用导数的定义与其物理意义
利用上述公式，选取适当的步长h，利用F点位移就可以得出速度与加速度
5、用VC编程
#include <stdio.h>
#include <math.h>
#define pi 3.14159265358979323846
//定义全局变量
double Lab,Lbc,Lcd,Xf,Xd,Yd;//定义已知位置量
11.079883
-670.906136
223
-53.471235
-5.223053
533.839546
44
68.312705
10.958430

哈工大有限元大作业1-16S108346殷琪

有限元大作业1
一、题目
两对边简支、另两对边固支的矩形薄板（如图1），其上承受均布压力q的作用，已知：a = 1000mm， b = 500mm，板厚 h = 20mm， q=-1N/mm2，E=2.1×105N，μ = 0 .3 求对称截面上的挠度w及弯矩 Mx、My的分布曲线。

图1.矩形薄板
二、求解过程
1．进入ANSYS 2．设置计算类型 3．选择单元类型 4．定义材料参数5．确定薄板的厚度 6．生成几何模型 7．网格划分
8．模型施加约：
a 分别给两长边施加固支约束、给两短边施加简支约束
b 给薄板施加均布载荷
9．分析计算
11．弯矩Mx云图显示（如下图）
13．定义与YOZ平面和长边平行的对称面Mx上的路径14．显示对称截面Mx上的挠度（Mx）曲线（如下图）
16．定义与XOZ平面和长边平行的对称截面Mx上的路径17．显示对称截面My上的挠度（My）曲线（如下图）
三、计算结果及感受
在本次分析报告中,平面Mx为与YOZ平面平行的对称截面，平面My为与XOZ 平面平行的对称截面。

我根据《ANYSYS基础进阶》一书中所提供的步骤来做的分析，得到的结果如上面的图中所示，根据薄板Z方向的挠度云图可知，薄板的挠度为-1.08559mm。

通过这次作业感受到了陈老师所说的ANSYS等有限元分析软件只是一个工具，而真正的奥妙在有限元理论中，以前只是单纯的用智能划分网格（smart），现在知道了去想如何划分才能得到较合理分析。

哈工大机械原理大作业(连杆机构)

三、建立数学模型
建立坐标系：以C为原点，水平方向为X轴，CA所在直线为Y轴建立直角坐标系（如图4）。
取曲柄1水平且位于A点右侧为初始时刻，设曲柄1角速度为w，由题意知w= =8.5π rad/s………………（1）
设曲柄1转角为θ，则B点坐标：
xB=ιABcosθ=ιABcoswt
yB=H1+ιABsinθ=H1+ιABsinwt………………（2）
form=1:length(t)-1
ddxF(m)=(dxF(m+1)-dxF(m))/0.0001;
end
ddxF(length(t))=ddxF(length(t)-1);
figure
plot(t,ddxF)
title('¼ÓËÙ¶ÈÍ¼Ïñ');
xlabel('t /s'),ylabel('v /(m/s^2)');
输出图像：
xE(m)=yE(m)/k(m);
xF(m)=xE(m)-(-H^2+lEF^2-yE(m)^2+2*yE(m)*H)^(1/2)+0.1142;
end
form=1:length(t)-1
dxF(m)=(xF(m+1)-xF(m))/0.0001;
end
dxF(length(t))=dxF(length(t)-1);
∵ιEF+ιCE>H且ιCE<H
∴E点始终在F点的右下方
∴xF<xE,所以x2舍去，只取xF=x1……………（8）
∴点F坐标为（xF，H）
当t=0时，可得F点初始位置坐标，不妨设为（xo,H）。
则F点位移（通过计算，t=0时，得xo=-0.1142）

哈工大有限元大作业

机械有限元习题答案——哈工大

哈工大机械原理大作业1 连杆机构运动分析 牛头刨床 学号5号

哈工大有限元大作业

哈工大机械原理大作业连杆机构9

哈工大研究生塑性体积成形大作业2

哈工大机械原理大作业19齿轮.docx

(完整word版)哈工大机械原理大作业凸轮DOC

哈尔滨工业大学 04 结构力学——平面有限元1

哈工大机械原理大作业1

哈工大有限元大作业1-16S108346殷琪

哈工大机械原理大作业(连杆机构)

哈工大机械原理大作业1 连杆机构运动分析牛头刨床学号5号