ansys实例分析2

  • 格式:pdf
  • 大小:480.08 KB
  • 文档页数:6

ansys实例分析2

实验四 压杆的最优化设计

一、问题描述

图45所示的空心压杆两端受轴向外载荷P。轴的内径为1d,外径为2d,支承间距尺

寸为l。试确定压杆的结构尺寸1d、2d和l,以保证在压杆不产生屈服并且不破坏压杆稳定

性条件下,压杆的体积和重量最小。性条件下,压杆的体积和重量最小。 该问题的分析过程如下:该问题的分析过程如下:

压杆为细长直杆,承受轴向压力,会因轴向压力达到临界值时突然弯曲而失去稳定性。

设计压杆,除应使其压力不超过材料的弹性极限外,还必须使其承受的轴向压力小于压杆的

临界载荷。临界载荷。

压杆在机械装置中应用的例子较多,例如在液压机构中当活塞的行程足够大时,会导致

活塞杆足够长,这种细长的活塞杆便是压杆。活塞杆足够长,这种细长的活塞杆便是压杆。

根据欧拉压杆公式,对与两端均为铰支的压杆,其临界载荷为根据欧拉压杆公式,对与两端均为铰支的压杆,其临界载荷为

2

2cPEJlp

=

式中E为压杆材料的弹性模量;J为压杆横截面的

最小惯性矩,EJ为抗弯刚度;L为压杆长度。为压杆长度。

将欧拉公式推广到端部不同约束的压杆,则

上式变为上式变为

()2

2cPEJlp

m=

式中m为长度折算系数,其值将随压杆两端约束

形式的不同而异。形式的不同而异。当两端铰支时取当两端铰支时取1m=;一端固

定,另一端自由时取2m=;一端固定,一端固定,另一端铰另一端铰

支时取0.75m=;两端均固定时,取0.5m=。

由欧拉公式可知,cP与

J成正比。合理的设计压杆截面形状,使其材料尽量远离形心图45 压杆机构简图压杆机构简图

分布,就能使J增大而提高压更的抗弯刚度EJ,增大临界载荷cP。所以在相同截面面积的

条件下,管状压杆比实心压杆有更大的临界载荷。条件下,管状压杆比实心压杆有更大的临界载荷。

以管状压杆的内径1d,外径2d,长度l作为设计变量,以其体积或重量作为目标函数,

以压杆不产生屈服和不破坏轴向稳定性以及其尺寸约束条件,则管状压杆最优化设计的数学

模型为模型为

()()

()()[]

()()()2221221

1122221442212122

1min11max

2min22max

minmaxmin,,4

4,0

,0,64fddlddl

Pgdddd

ddgddPEJJl

ddd

ddd

lllp

asp

pp

bm=-

=-£-

-=-£=

££

££

££对圆管:

式中,ab分别为大于1的安全系数;P为设计给定的外载荷。为设计给定的外载荷。

二、实训目的

理解“最优化设计”的概念和方法;了解“最优化设计”的建模方法;了解“最优化设计”问题的解决过程;了解LINGO的编程、方法和求解过程。的编程、方法和求解过程。

三、实验步骤

1.运行LINGO程序。LINGO程序主界面如图3-3所示。所示。

2. 点击下拉菜单,选择或者单击工具栏中的按钮,新建一个优化实例模型LINGO 

Model(*.lg4)。

3.在“LINGO Model”窗口建立模型,输入以下程序语句:”窗口建立模型,输入以下程序语句:”窗口建立模型,输入以下程序语句:

Model: 

min=3.14*(d2^2-d1^2)*L/4;=3.14*(d2^2-d1^2)*L/4;! ! 目标函数;

!内径约束条件如下;

d1min=0.02; 

d1max=0.03; 

d1>=dimin; 

d1<=d1max; 

!外径约束条件如下; d2min=0.02; 

d2max=0.03; d2>=d2min; 

d2<=d2max; 

!长度约束条件如下; Lmin=0.1; 

Lmax=1; 

L>=Lmin; 

L<=Lmax; !定义安全系数;

a=1.2; 

b=1.3; 

!外加载荷20kN; P=20*10^3; 

E=200*10^9; !弹性模量; C=350*10^6; !材料抗压极限;

(4*a*P)/(3.14*(d2^2-d1^2))<=C; !不产生屈服的约束条件;

u=0.5; !长度折算系数;

J=3.14*(d2^4-d1^4)/64; !刚度系数;

P<=(3.14^2*E*J)/(

b*(u*L)^2); 

END

4.模型建立完

成后要对模型进行

求解:单击下拉菜

单“LINGO”选择

“Sovle”,或者单

击工具栏中的模型

求解命令按钮。

LINGO程序会弹

出“LINGO LINGO Solve Solve 

Status”窗口和

“Solution Report”

窗口。在“Solution 

Report”窗口中会

列出完整的模型求

解过程

,在图46 求解状态窗口求解状态窗口

“LINGO Solve Status”窗口中列出了求解器的有关信息,如图”窗口中列出了求解器的有关信息,如图46所示。所示。

实验五 悬臂梁的可靠度分析

一、问题描述

如图47所示,长度L=2m的悬臂梁,自由端受集中力F=50KN作用。弹性模量作用。弹性模量

E=3.0e7kPa,泊松比v=0.17。悬臂梁截面:面积240.04,0.003,0.2AmImHm===

二、实验步骤

1)建立有限元模型建立有限元模型

(1)定义分析标题:Utility 

Menu->File->Chan

ge Title,输

入”Probabilistic Design of a beam”

,单击OK。

(2)定义参数初始值:Utility 

Menu->Parameters->Scalar Parameters,打开数值参数对话打开数值参数对话打开数值参数对话

框。在Selection编辑栏中输入F=50,按Enter键:输入E=3.0E7,按Enter键;单击C1ose

关闭数值参数对话框。关闭数值参数对话框。

(3)定义单元类型:Main Menu->Preprocessor->Element Type->Add/Edit/Delete,出现Element Types列表框。单击Add出现单元类型库对话框,左侧列表中选择左侧列表中选择 Structural Beam,

右侧列表中选择2DElastic3,单击OK。

(4)定义实常数:Main Main Menu->Preprocessor->Real Menu->Preprocessor->Real Menu->Preprocessor->Real Constants->Add/Edit/Delete,Constants->Add/Edit/Delete,出现Real 

Constants列表框。单击Add出现E1ement Type for Real Constants对话框,单击OK,出现

Real Constant Set Number 1, For Beam3对话框。在AREA中输入0.04,IZZ中输入0.0003,

HEIGHT中输入0.2,单击OK 。单击Close关闭Real Constants列表框. 

(5)定义材料特性:Main Menu->Preprocessor->Material Props->Material Models,在Define 

Material Model Behavior窗口中,双击Structural->Linear->Elastic->Isotropic。在出现的对话。在出现的对话

框中,在Ex中输入3.0E7,在PRXY中输入0.17,单击OK, Material Model l出现在Material 

Models Defined窗口左侧列表中,选择菜单Material ->Exit退出。退出。

(6)建立模型:Main Main Menu->Preprocessor->Create->Keypoints->In Menu->Preprocessor->Create->Keypoints->In Menu->Preprocessor->Create->Keypoints->In Active Active Active CS,CS,出现Create 

Keypoints in Active Coordinate System对话框。在对话框。在Location in Active CS中输入中输入0,0,

0,单图47 自由端受集中力F=50KN作用作用

击Apply,输入2,0,0,单击OK,图形窗口出现2个关键点。执行Main 

Menu->Preprocessor->Create->Lines->Straight Line,出现Create Straight Line选择框。在图形

窗口中依次选择l、2关键点,单击OK生成一直线段。生成一直线段。 (7)定义线的属性: Main Menu->Preprocessor->Define->All Lines,出现Line attributes对话框对话框

。将Pick Orientation Keypoint激活为YES,单击OK,出现Line Line attributesattributes选择框,在选择

框单击OK。

(8)划分网格:Menu->Preprocessor->->Mesh Tool,出现Mesh Tool对话框。单击Lines右

侧的Set,出现Element Size on Picked选择框,选择框,单击Pick A11,出现Element Size on Picked lines 

对话框。在NDIV中输入20,单击OK,在图形窗口中每段弧线以20段均分线段表示。在段均分线段表示。在

Mesh Tool对话框中单击Mesh,出现Mesh Lines选择框,单击Pick All生成单元。生成单元。 (9)定义分析类型:Menu->Solution->Unabridged Menu->New Analysis.出现New Analysis

对话框,选样Static,单击OK。

(10)施加约束与荷载:Menu->Solution->Loads->Apply->Displacement->On Menu->Solution->Loads->Apply->Displacement->On NodesNodes,出现Apply UROT on Nodes选择框。选择框。在图形窗口中选择左端节点,在选择框中单击Ok出现Apply 

UROT UROT on on on NodesNodes对话框。选中A11 A11 DOFDOF单击OK。约束标识出现在图形窗口。执行

Menu->Solution->Loads->Apply->Force/Moment->On Nodes,出现Apply F/M on Nodes选择

框,选择右端节点,单击OK,出现Apply F/M on Nodes对话框。对话框。在Lab中选样FY,在Value中输入“F”,单击OK,集中力标识出现在梁的自由端。,集中力标识出现在梁的自由端。

(11)保存数据库文件:Utility Menu->File->Save as Jobname.db。

(12)求解:Menu->Solution->Loads->Solve->Current Menu->Solution->Loads->Solve->Current LSLS,查看求解信息,关闭求解状态