有限元上机实验报告材料(1).docxVIP

PAGE \* MERGEFORMAT1 工程有限元与数值计算上机实验报告 姓名：熊健 学号：2150100196 班级：机械56 上机实验一：ANSYS经典训练 2.1三角桁架的ANSYS仿真分析与MATLAB编程 MATLAB仿真结果截图： MATLAB主程序： E=2.1e11; A=0.01; x1=0; y1=0; x2=sqrt(3); y2=0; x3=0; y3=1; alpha1=0; alpha2=90; alpha3=150; k1=Bar2D2Node_Stiffness(E,A,x1,y1,x2,y2,alpha1) k2=Bar2D2Node_Stiffness(E,A,x1,y1,x3,y3,alpha2) k3=Bar2D2Node_Stiffness(E,A,x2,y2,x3,y3,alpha3) KK=zeros(6,6); KK=Bar2D2Node_Assembly(KK,k1,1,2); KK=Bar2D2Node_Assembly(KK,k2,1,3); KK=Bar2D2Node_Assembly(KK,k3,2,3) % k=KK([3,4],[3,4]) % p=[10000,10000]; % u=k\p q=[0 0 0.000022538 0.000077135 0 0]; P=KK*q u1=[q(1);q(2);q(3);q(4)] stress1=Bar2D2Node_Stress(E,x1,y1,x2,y2,alpha1,u1) u2=[q(1);q(2);q(5);q(6)] stress2=Bar2D2Node_Stress(E,x1,y1,x3,y3,alpha1,u1) u3=[q(3);q(4);q(5);q(6)] stress3=Bar2D2Node_Stress(E,x2,y2,x3,y3,alpha1,u1) 其中求得个节点支反力为： P =（-27326N，0N，10004N，10001N，17322N，-10001N） 单元1受力：stress1 = 2732600pa； 单元2受力：stress2 = 4732980pa； 单元3受力：stress3 = 2366500pa ANSYS仿真过程及结果： ①进入ANSYS(设定工作目录和工作文件)程序→ANSYS→ANSYS Interactive→Working directory (设置工作目录)→Initialjobname(设置工作文件名): planetruss→Run→OK（也可以直接进入） ②设置计算类型： ANSYS Main Menu: Preferences…→Structural→OK ③选择单元类型： ANSYS Main Menu: Preprocessor→Element Type→Add/Edit/Delete…→Add…→Link：3D finit stn180→OK (返回到 Element Types 窗口)→Close ④定义材料参数： ANSYS Main Menu: Preprocessor→Material Props→Material Models→Structural→Linear →Elastic→Isotropic: EX:2.1e11 (弹性模量)， PRXY: 0.3(泊松比)→OK→鼠标点击该窗口右上角的“”来关闭该窗口 ⑤定义实常数以确定单元的截面积： ANSYS Main Menu: Preprocessor→Real Constants…→Add/Edit/Delete→Add→Type 1→OK→Real ConstantSetNo:1(第1号实常数),AREA:0.01(单元的截面积)→OK→Close ⑥生成单元以及划分网格： ANSYS Main Menu: Preprocessor→Modeling→Creat→Nodes→In Active CS→Node number 1→X:0,Y:0,Z:0→Apply→Node number 2→X:0.4,Y:0,Z:0→Apply→Node number 3→X:0.4,Y:0.3,Z:0→Apply→Node number 4→X:0,Y:0.3,Z:0→OK ANSYS Main Menu:Preprocessor→Modeling→Create→ Elements→Elem Attributes（接受默认值）→Usernumbered→Thrunodes→OK→选择 node1 和 node2→Apply→选择 node2 和 node3→Apply→选择 node1 和 node3