当前位置:文档之家› 有限元-三杆桁架地优化设计

有限元-三杆桁架地优化设计

有限元作业作业题目:三杆桁架的优化设计学生姓名:苏宏斌学号:200900403064专业:数控技术班级:机自Y094指导老师:王虎奇三杆桁架的优化设计题目描述::下图所示为一个有3根杆组成的桁架,承受纵向和横向载荷,杆件的横截面面积和基本尺寸B在指定范围内变化,要求桁架的每根杆件承受的最大应力小于(800+学号最后两位数)MPa,试对该结构进行优化设计,使得桁架重量最少。

弹性模量E=220GPa;泊松比:0.3;密度ρ=7800kg/m3材料最大许用应力:σ= 864 MPa横截面面积变化范围:0.01~10cm2(初始值为10)基本尺寸B变化范围:1~2m(初始值为2)前处理:(1)定义工作文件名:utility menu-file-change jobname,在弹出的change jobname对话框中输入文件名为: suhongbin 单击ok按钮。

(2)定义工作标题:utility menu-file-change tile,在弹出的change tile对话框中输入suhongbin ,单击ok按钮。

(3)定义参数的初始值:utility menu-parameters-scalar parameters命令,弹出对话框,在selection下的文本框中输入B=2,按下enter键;A1=0.001, 按下enter键;A2=0.001, 按下enter键;A3=0.001, 单击Close按钮。

参数将在菜单中显示。

(4)设置材料属性:main menu-preprocessor-material props-material model命令,设置EX=2.2e11 , PRXY=0.3 , DENS=7800。

(5)定义单元类型:main menu-preprocessor-element type-add/edit/delete命令,弹出element type对话框。

单击add按钮,弹出library of element type 对话框,在左边列中选择structural link,在右边列中选择2D spar 1,单击ok。

(6)定义实常数:main menu-preprocessor-real constant命令,弹出real constants对话框。

单击add按钮,又单击ok按钮,弹出real constant set number 1,for link1对话框。

在real constant set no后面的文本框中输入1,在cross-section area后面的文本框中输入A1,单击apply按钮,在弹出同样对话框的同样位置分别输入2,A2, 单击apply按钮,在弹出同样对话框的同样位置分别输入3,A3单击ok按钮,单击close按钮,完成实常数的设置。

(self:此步最好按写的步骤做,不要用0.001替代A1,即使A1=0.001,因为A1也许在后面是变量,就不是0.001了)(7)生成有限元节点:main menu-preprocessor-moding-create-nodes-in active cs命令,弹出的create nodes in active coordinate system对话框,在node number后面的文本中输入1,其坐标值为:-B,0,0,单击apply按钮,在node number后面的文本中输入2,其坐标值为:0,0,0,单击apply按钮,在node number后面的文本中输入3,其坐标值为:B,0,0,单击apply按钮,在node number后面的文本中输入4,其坐标值为:0,-2,0,单击ok按钮.(8)关闭坐标符号显示:utility menu-plotctrls-window options,弹出一个对话框,在location of triad 后面的下拉式选择栏中选择no shown ,单击ok 。

(9)打开节点编号显示:utility menu-plotctrls-numbering命令,弹出plot numbering controls对话框。

选择node number复选框,单击按ok钮。

(10)生成第一个单元:main menu-preprocessor-moding-create-elements-auto numbered-thru nodes命令,弹出一个拾取框,拾取编号为1和4(按照该顺序)的节点,单击ok按钮。

(11)改变第二个单元属性:mainmenu-preprocessor-moding-create-elements-elem attributes命令,弹出element attributes对话框。

在real constant set number下拉列表中选择2,单击ok按钮。

(12)生成第二个单元:main menu-preprocessor-moding-create-elements-auto numbered-thru nodes命令,弹出一个拾取框,拾取编号为2和4(按照该顺序)的节点,单击ok按钮。

(13)改变第二个单元属性:main menu-preprocessor-moding-create-elements-elem attributes命令,弹出element attributes对话框。

在real constant set number下拉列表中选择3,单击ok按钮。

(14)生成第二个单元:main menu-preprocessor-moding-create-elements-auto numbered-thru nodes命令,弹出一个拾取框,拾取编号为3和4(按照该顺序)的节点,单击ok按钮。

求解:(1)施加边界条件:main menu-solution-loadsapply-structural-displacement-on nodes命令,给编号为1,2,3的节点进行全约束。

(2)施加集中载荷:main menu-solution-loadsapply-structural force/moment-on nodes命令,给编号为4的节点施加Fx方向-20000,Fy方向-20000的力。

(3)保存数据:单击工具栏中的save_db按钮。

(4)求解运算:main menu-solution-solve-current ls命令。

(5)保存优化结果到文件:utility menu-file-save as,在弹出的对话框中输入文件名为truss_resu,单击ok按钮。

优化设置:(1)定义单元表:main menu-general postproc-element table-define table命令,弹出element table date对话框,单击add按钮,弹出define additional element table items对话框。

在user label for item后面的文本框中输入EVOLUME。

又在item,comp results date item的左栏中选择geometry,在右栏中选择elem volume volu.单击ok按钮,单击close按钮。

(2)计算单元体积的总和:main menu-general postproc-element table-sum of each item命令,弹出tabular sum of each element table item对话框,单击ok按钮,弹出一个信息窗口,在窗口中显示体积总和为76.5685。

选择该窗口菜单栏上的file-close,关闭该窗口。

(3)取出体积的值:utility menu-parameters-get scalar date命令,弹出get scalar data对话框。

在type of data to be retrived的左栏中选择results data,在右栏中选择elem table sums,单击ok按钮,弹出get element table sum results 对话框,在name of parameter to be defined后面的文本中输入VTOT,单击ok按钮关闭该对话框。

(4)计算初始重量:utility menu-parameters-scalar parameters命令,弹出scalarparameters对话框,在selection下面的文本中输入DENS=7800,并按下enter 键,又在selection下面的文本中输入WT=DENS*VTOT,并按下enter键,总的体积将计算并在该对话框中显示,总的重量应为59.7234631.单击close按钮关闭该对话框。

(5)设置单元表:main menu-general postproc-element table-define table命令,弹出element table data对话框,单击add按钮,弹出define additional element table items对话框,在user lable for item后面的文本框中输入SIGMA,又在item,comp results data item的左栏中选择by sequence num,在右栏中选择LS, 在其下面出现的LS, 后面的文本框中输入1(表示序列号为1),单击ok按钮,又单击element table data对话框的close按钮,关闭该对话框。

(6)得到第一杆的轴向应力:utility menu-parameters-get scalar data命令,弹出get scalar data对话框,在the type of data to be retrived的左栏中选中results data,在右栏中选择elem table data,单击ok按钮,弹出get element table data对话框,在name of parameter to be defined后面的文本中输入sig1,在element number n后面的文本中输入1,在elem table data to be retrived后面的下拉列表中选取SIGMA,单击apply按钮。

(7)得到第二杆的轴向应力:又弹出get scalar data对话框,在the type of data tobe retrived的左栏中选中results data,在右栏中选择elem table data,单击ok按钮,弹出get element table data对话框,在name of parameter to be defined后面的文本中输入sig2,在element number n后面的文本中输入2,在elem table data to be retrived后面的下拉列表中选取SIGMA,单击apply 按钮。

相关主题