目录
.绪论 (2)
第一章.有限元课程设计 (4)
一.工程问题 (4)
二.简化模型 (4)
三.解析法求解 (5)
四.ANSYS求解 (8)
五.结果分析 (19)
第二章.机械优化设计说明 (20)
一.题目及解析 (20)
二.黄金分割法计算框图 (23)
三.C语言程序 (24)
四.运行结果 (27)
五.结果分析 (27)
第三章.设计感言 (28)
第四章.参考文献 (28)
前言
有限元法在解决圣维南扭转问题近似解时首先提出的。有限元在弹性力学平面问题的第一个成功应用是由美国学者于1956年解决飞机结构强度时提出的、经过几十年得发展,有限元一惊成为现代结构分析得有效方法和主要手段。它的应用已经从弹性力学的平面问题扩展到空间问题和板壳问题。对于有限元法,从选择基本未知量的角度来看,他可以分为三种方法:位移法,力法,混合法。从推导方法来看,它可以分为直线法,变分法,加权余数法。但随后随着计算机的发展,有限元法如虎添翼。国内外已有许多大型通用的有限元分析程序,并已经出现了将人工智能技术引入有限元分析软件,形成了比较完善得专家系统,逐步实现了有限元的智能化。
优化设计是现代设计方法的重要内容之一。它以数学规划为理论基础以电子计算机为工具,在充分考虑多种设计约束的前提下,寻求满足预订目标的最佳设计。优化设计理论于方法用于工程设计是在六十年代后期开始的,特别是今年来,随着有限元素法,可靠性设计,计算机辅助设计的理论与发展及优化设计方法的综合应用使整个工程设计过程逐步向自动化集成化智能化发展,其前景使令人鼓舞的。因而工程设计工作者必须适应这种发展变化,学习,掌握和应用优化设计理论与方法。
今年来随着计算机技术的普及和计算速度的不断提高,有限元分析在工程设计和分析中得到了越来越广泛的重视,已经成为解决复杂的工程分析计算问题的有效途径,现在从汽车到航天飞机几乎所有
的机械制造都已离不开有限元分析计算,其再机械制造,材料加工,航空航天,汽车,土木建筑,电子电器,国防军土,船舶,铁道,石化能源,科学研究等各个领域的广泛使用已使设计水平发生了质的飞跃,主要表现在以下几个方面:
增加产品和工程的可靠性
在产品的设计阶段发现潜在的问题
经过分析计算,采用优化设计方案,降低原材料成本
缩短产品投向市场的时间
模拟试验方案,减少试验次数,从而减少试验经费
ANSYS软件致力于耦合场的分析计算,能够进行结构,流体,热,电磁四种场的计算,已博得了世界上数千家用户的钟爱。ADINA 非线性有限元分析软件由著名的有限元专家,麻省理工学院的K.J.Bathe教授领导开发,其单一系统即可进行结构,流体,热的耦合计算。并同时具有隐式和显式两种时间积分算法。犹豫其在非线性求解,流固耦合分析等发面的强大功能,迅速成为有限元分析软件的后起之秀,现已成为非线性分析得首选软件。
第一章·有限元课程设计
一、工程问题
材料力学Ι辅导及习题精解(陈平第四版) P180
受力如图所示的梁,试根据 q,Fs,M之间的微分关系和积分关系画出Fs,M图。
二、简化模型
1、梁的参数:
长度 l=10m;
宽度 b=1m
厚度 h=0.5m
2、材料参数
梁选择线性、弹性、各向同性的材料。它的弹性模量EI=207e5 Pa.。
3、梁的边界条件
在节点A处梁受Y方向的约束;节点B受X方向和Y方向的约束。
4、梁的载荷
FB之间作用着均布载荷q=2kN/m,在节点E处作用着集中力偶Me=10kN·m。D 处作用沿Y负方向的集
中力 F=1 KN,在C处作用沿Y负方向的集中力F=3KN。
经上述分析,此外伸梁可简化成一个二维单元,其中一端固
支,一端外伸。另一个支撑点变化为一个绞支点,在梁的上部分布着均布载荷。如图所示(a)
三、解析法求解
(1)先求支座反力,得Fa=6KN Fb=6Kn
(2)做Fs图如图所示
从梁左端开始。由C截面有向下的集中力F1的作用,所以Fs图有零开始向下突变,其突变值为F1=-3KN。由于AC段内无分布载荷作用,所以C,A两截面的剪力增量为零。CA段得剪力图为一条水平直线,并从C点一直延伸到A点稍偏左的截面处。由于A截面处有向上的集中力Fa的作用,Fs图有一向上的突变,值为Fa=+6KN,所以A右截面的剪力值为-3+6=+3KN.在AD段,又因为没有均布载荷所以Fs图也为一条直线,从A点一直延伸到D点稍左的截面处。由于D 截面处有一向下的集中力F2的作用,所以Fs图有一向下的突变,值为F2=-1KN,所以D右截面的剪力值为3-1=+2KN。在DEF段,因为均布载荷也为零所以Fs图为一水平直线,从D点一直延伸到F截面处。由于在FB段有向下的均布载荷作用,q=-2KN/m,所以Fs图为一下降的斜直线。BF两截面上的剪力之差等于载荷q(x)图的面积,即-2乘4=-8KN,所以B左截面上的剪力值为2-8=-6KN.最后在B截面由于有向上的集中力Fb的作用,所以Fs图有一个向上的突变,值为Fb=+6KN 于是,Fs图画出如图所示(b)。
(3)做M图(如图所示(c))
因为C截面为自由端,又没有集中力偶作用,所以Mc=0;弯矩图从0开始。在CA段,因为没有均布载荷,又因为Fs=-3<0,所以M图为一下降的斜直线;而A C两截面上的弯矩之差,就是CA段的Fs图的面积,即-3乘2=-6KN·m ,所以,A截面的弯矩值为-6 KN·m.在AD段,因为均布载荷为0,Fs=3KN>0,所以M图为一上升的斜直线;而D A两截面上的弯矩之差,等于AD段的Fs图的面积,即3乘1=3 KN·m
所以D截面的弯矩值为-6+3=-3 KN·m。在DE段,均布载荷为0,Fs=2KN>0,所以M图也为一个上增的斜直线;而E D两截面上的弯矩之差等于DE段Fs 图的面积,即2乘1=2 KN·m,所以E左截面的弯矩值为-3+2=-1 KN·m。在E截面,因为有一个顺时针转向的集中力偶作用,所以M图有一向上的突变,值为Me==5 KN·m,所以E右截面的弯矩值为-1+5=+4 KN·m。在EF段,因为没有均布载荷,又因为Fs=2KN>0,所以M图为一上增的斜直线;F E两截面上的弯矩之差等于EF段的Fs图面积,即2乘2=4 KN·m,所以F截面的弯矩值为4+4=8 KN·m。在FB段,由于q(x)=-2 KN/m,所以M图为一上凸的曲线。又因为在F截面以右的1m的FG段内Fs>0,所以FG段得M图为一上增而向上凸的二次曲线。在G截面,因为Fs由正变负,所以,M图在这里有极大值。G F两截面上的弯矩之差等于FG段的Fs图的面积,即1/2乘2乘1=1 KN·m所以,G截面的弯矩值为8+1=9 KN·m。B G 两截面上的弯矩之差等于GB段的Fs图的面积即1/2乘(-6)乘3=-9 KN·m所以,B截面的弯矩值为9-9=0,M图回到0.全梁的M图绘于
(c)中。如图
四、ANSYS求解
将梁划分为20个单元,21个节点,用BEAM3来建立单元进行静力学分析。
交互式的求解过程
1.创建节点
1.1 创建梁的各个节点
1.Main Menu:Preprocessor→Modeling→Create→Node→In Active CS。
2.在创建节点窗口内,在NODE后的编辑框内输入节点号1,并在X,Y,Z后的编辑框内输入0,0,0作为节点
1的坐标值。
3.按下该窗口内的Apply按钮。
4.输入节点号21,并在X,Y,Z后的编辑框内输入10,0,0作为节点21的坐标值。
5.按下OK按钮。
6.Main Menu:Preprocessor→-Modeling-Create→Node→Fill between Nds。
7.在图形窗口内,用鼠标选择节点1和21。
8.按下Fill between Nds窗口内的Apply按钮。
9.按下OK按钮,完成在节点1到节点21之间节点的充1.2 显示各个节点
1.Utility Menu:Plotctrls→Numberings
2.将Node numbers项设置为On。
3.Utility Menu:Plot→Nodes
4.Utility Menu:List→Nodes
5.对出现的窗口不做任何操作,按下OK按钮。
6.浏览节点信息后,关闭该信息窗口。
2.定义单元类型和材料特性
2.1 定义单元类型
1.Main Menu:Preprocessor→Element Type→Add/Edit/Delete
2.按下Element Type窗口内的Add按钮。
3.在单元类型库中,选择左侧列表中的BEAM单元家族,及右侧列表中2D elastic 3类型。
4.按下OK按钮完成选择。
5.按下Close按钮关闭Element Type窗口。
2.2 定义材料特性
1.Main Menu:Preprocessor→Material Props→Material
Models。
1. 在材料定义窗口内选择:Structural→Linear→Elastic→
Isotropic。
2.在EX后的文本框内输入数值207e5作为弹性模量。
3.按下OK按钮完成定义。
2.3 定义几何参数
1.Main Menu:Preprocessor→Real Constants→Add/Edit/Delete。
2.按下Real Constants窗口内的Add按钮。
按下Real Constants for Element Type窗口内的OK按钮。
3.依次输入1,1,0.02088,0.5。
4.按下OK按钮完成定义。
5.按下Real Constants窗口内的Close按钮。
3.创建单元
3.1 创建单元
1.Main Menu:Preprocessor→Create→Elements→Auto-Numbered→Thru Nodes。
2.在图形窗口内,用鼠标点选节点1和2。
3.按下按下OK按钮完成单元1的定义。
4.Main Menu:Preprocessor→Model→Copy→Elements→Auto-Numbered。用光标选择单元1,然后点Apply。
5.在ITIME后的编辑框内输入20(包括被复制的单元1)作为要复制的单元总数。
6.按下按下OK按钮完成单元2到单元20的定义。
3.2 显示单元资料
1.Utility Menu:PlotCtrls→Numberings
2.在第一个下拉列表中,选择Elements numbers选项。
3.Utility Menu:Plot→Elements
4.Utility Menu:List→Elements→Nodes+Attributes
5.浏览单元信息后,关闭该窗口。
4.施加约束和载荷
4.1 节点自由度约束
1.Main Menu:Solution→Define Loads→Apply→Structural→Displacement→On nodes。
2.用鼠标在图形窗口内选择节点5。
3.按下选择窗口内的Apply按钮。
4.选择自由度UY,并在V ALUE后为其输入数值0。
5.按下Apply按钮。
6.用鼠标在图形窗口内选择节点21。
7.按下选择窗口内的Apply按钮。
8.选择自由度UX和UY,并在V ALUE后为其输入数值0。
9.按下OK按钮。
4.2施加载荷
4.2.1施加节点1处的集中载荷F1。
1.Main Menu:Solution→Define Loads→Apply→Structural→Force/Moment→On nodes。
2.用鼠标在图形窗口内选择节点1。
3.按下选择窗口内的Apply按钮。
4.在第一个下拉列表中选择FY,并在下面的文本框内输入其值-3(向上为Y轴正方向)。
5.按下Apply按钮。
4.2.2施加节点7处的集中载荷F2。
6.Main Menu:Solution→Define Loads→Apply→
Structural→Force/Moment→On nodes。
7.用鼠标在图形窗口内选择节点7。
8.按下选择窗口内的Apply按钮。
9.在第一个下拉列表中选择FY,并在下面的文本框内输入其值-1(向上为Y轴正方向)。
10.按下Apply按钮。
4.2.3施加节点3处的弯矩m。
1.Main Menu:Solution→Define Loads→Apply→Structural→Force/Moment→On nodes。
2.用鼠标在图形窗口内选择节点9。
3.按下选择窗口内的Apply按钮。
4.在第一个下拉列表中选择MZ,并在下面的文本框内输入其值-5(逆时针为正方向)(对照上面第4步)。
5.按下OK按钮。
4.2.4施加单元13到单元20上的的分布载荷q。
1.Main Menu:Solution→Define Loads→Apply→Structural →Pressure →On Beams。
2.用鼠标在图形窗口内选择单元13到单元20。
3.按下选择窗口内的Apply按钮。
4.在LKEY后的文本框内输入数值1。
5.在V ALI和V ALJ后的编辑框内分别输入10,
6.按下OK按钮。
5.求解
5.1 定义分析类型
1.Main Menu:Solution→Anslysis Type→New Analysis。
2.选中Static选项。
3.按下OK按钮。
5.2 求解
1.Main Menu:Solution→Solve→Current Ls。
2.按下OK按钮关闭Solve Current Load Step窗口。
3.按下Close按钮关闭求解结束后出现的Information窗口。
4.浏览/STATUS Command窗口内的信息后,将其关闭。
6.后处理
6.1 显示梁变形结果
1.Main Menu:General Postproc→Plot Results→Contour Plot Nodal Solu... →选择DOF Solution下的Displacement
vector sum
2.不改变对话框内的任何项,按下OK按钮。
6.2 建立单元结果表
6.2.1创建单元表,计算节点弯矩。
1.Main Menu:General Postproc→Element Table→Define Table。
2.按下Element Table Data窗口内的Add按钮。
3.在Lab后的文本框内输入IMOMENT。
4.在左侧列表中选择By sequence num项。
5.右侧列表中选择SMICS,项。
6.在右侧列表下的文本框内输入SMICS,6。
7.按下Apply按钮。在Lab后的文本框内输入JMOMENT。
8. 在Lab后的文本框内输入JMOMENT。
9.重复上面的步骤4和5。
10.右侧列表下的文本框内输入SMICS,12。
11.按下OK按钮。
6.2.2创建单元表,计算节点剪力。
1.Main Menu:General Postproc→Element Table→Define Table。
2.按下Element Table Data窗口内的Add按钮。
3.在Lab后的文本框内输入ISHEAR。
4.在左侧列表中选择By sequence num项。
5.右侧列表中选择SMICS,项。
6.右侧列表下的文本框内输入SMICS,2。
7.按下Apply按钮。
8.在Lab后的文本框内输入JSHEAR。
9.重复上面的步骤4和5。
10.右侧列表下的文本框内输入SMICS,8。
11.按下OK按钮。
6.3 列出所有表格资料
6.3.1列出资料
1.Main Menu:General Postproc→List Results→Element Table Data。
2.在List Element Table Data窗口内选择IMOMENT,JMOMENT,ISHEAR和JSHEAR。
3.按下OK按钮并在浏览资料窗口内的信息后,将其关闭。
6.3.2画剪力图
1.Main Menu:General Postproc→Plot Results→Line Elem Res
2.在第一个下拉列表中选择ISHEAR,在第二个下拉列表中选择JSHEAR。
3.按下OK按钮。
6.3.3画弯矩图
1.Main Menu:General Postproc→Plot Results→Line Elem Res
2.在第一个下拉列表中选择IMOMENT,在第二个下拉列表中选择JMOMENT。
3.按下OK按钮。
弯矩图及其相应数据如下(图形已经过反色处理):
7.退出程序
1.Toolbar:Quit。
2.选择Quit-No Save!
3.按下OK按钮。
五、结果分析
梁的弯矩在有集中力偶的地方会发生跳变,而剪力是在有集中力的地方会有跳变的
两种方法的求解结果一样,证明在运用正确的方法,选用正确的单元与节点进行有限元的分析,能得到与实际相符的结果,所以在。工程实际中将实际问题转化成数学与物理模型,能得到想要的结果
第二章·机械优化设计说明
一.题目及解析
用黄金分割法求f(x)=2X3+3X2-6X+1的最优解.设初始点ɑ0=0,初始步长
h=1,取迭代精度ε=0.01
解:a1=a0=0 , f1=f(a1)=1
a2=a1+h=1 , f2=f(a2)=0
∵f1>f2 ,
∴作前进运算
轴心受力构件的强度和刚度计算 1.轴心受力构件的强度计算 轴心受力构件的强度是以截面的平均应力达到钢材的屈服应力为承载力极限状态。轴心受力构件的强度计算公式为 f A N n ≤= σ (4-1) 式中: N ——构件的轴心拉力或压力设计值; n A ——构件的净截面面积; f ——钢材的抗拉强度设计值。 对于采用高强度螺栓摩擦型连接的构件,验算净截面强度时一部分剪力已由孔前接触面传递。因此,验算最外列螺栓处危险截面的强度时,应按下式计算: f A N n ≤= ' σ (4-2) 'N =)5 .01(1 n n N - (4-3) 式中: n ——连接一侧的高强度螺栓总数; 1n ——计算截面(最外列螺栓处)上的高强度螺栓数; ——孔前传力系数。 采用高强度螺栓摩擦型连接的拉杆,除按式(4-2)验算净截面强度外,还应按下式验算毛截面强度 f A N ≤= σ (4-4) 式中: A ——构件的毛截面面积。 2.轴心受力构件的刚度计算 为满足结构的正常使用要求,轴心受力构件应具有一定的刚度,以保证构件不会在运输和安装过程中产生弯曲或过大的变形,以及使用期间因自重产生明显下挠,还有在动力荷载作用下发生较大的振动。 轴心受力构件的刚度是以限制其长细比来保证的,即
][λλ≤ (4-5) 式中: λ——构件的最大长细比; [λ]——构件的容许长细比。 3. 轴心受压构件的整体稳定计算 《规范》对轴心受压构件的整体稳定计算采用下列形式: f A N ≤? (4-25) 式中:?——轴心受压构件的整体稳定系数,y cr f σ?= 。 整体稳定系数?值应根据构件的截面分类和构件的长细比查表得到。 构件长细比λ应按照下列规定确定: (1)截面为双轴对称或极对称的构件 ? ?? ==y y y x x x i l i l //00λλ (4-26) 式中:x l 0,y l 0——构件对主轴x 和y 的计算长度; x i ,y i ——构件截面对主轴x 和y 的回转半径。 双轴对称十字形截面构件,x λ或y λ取值不得小于t (其中b/t 为悬伸板件宽厚比)。 (2)截面为单轴对称的构件 以上讨论柱的整定稳定临界力时,假定构件失稳时只发生弯曲而没有扭转,即所谓弯曲屈曲。对于单轴对称截面,绕对称轴失稳时,在弯曲的同时总伴随着扭转,即形成弯扭屈曲。在相同情况下,弯扭失稳比弯曲失稳的临界应力要低。因此,对双板T 形和槽形等单轴对称截面进行弯扭分析后,认为绕对称轴(设为y 轴)的稳定应取计及扭转效应的下列换算长细比代替y λ [] 2 /122202022222)/1(4)()(2 1 z y z y z y yz i e λ λλλλλλ--+++= )/7.25//(2 202ωωλl I I A i t z +=
MATLAB: MATLAB是美国MathWorks公司出品的商业数学软件,用于数据分析、无线通信、深度学习、图像处理与计算机视觉、信号处理、量化金融与风险管理、机器人,控制系统等领域。 MATLAB是matrix&laboratory两个词的组合,意为矩阵工厂(矩阵实验室),软件主要面对科学计算、可视化以及交互式程序设计的高科技计算环境。它将数值分析、矩阵计算、科学数据可视化以及非线性动态系统的建模和仿真等诸多强大功能集成在一个易于使用的视窗环境中,为科学研究、工程设计以及必须进行有效数值计算的众多科学领域提供了一种全面的解决方案,并在很大程度上摆脱了传统非交互式程序设计语言(如C、Fortran)的编辑模式。 MATLAB和Mathematica、Maple并称为三大数学软件。它在数学类科技应用软件中在数值计算方面首屈一指。MATLAB可以进行矩阵运算、绘制函数和数据、实现算法、创建用户界面、连接其他编程语言的程序等。MATLAB的基本数据单位是矩阵,它的指令表达式与数学、工程中常用的形式十分相似,故用MATLAB来解算问题要比用C,FORTRAN等语言完成相同的事情简捷得多,并且MATLAB也吸收了像Maple等软件的优点,使MATLAB成为一个强大的数学软件。在新的版本中也加入了对C,FORTRAN,C++,JAVA的支持。 MATLAB有限元分析与应用:
《MATLAB有限元分析与应用》是2004年4月清华大学出版社出版的图书,作者是卡坦,译者是韩来彬。 内容简介: 《MATLAB有限元分析与应用》特别强调对MATLAB的交互应用,书中的每个示例都以交互的方式求解,使读者很容易就能把MATLAB用于有限分析和应用。另外,《MATLAB有限元分析与应用》还提供了大量免费资源。 《MATLAB有限元分析与应用》采用当今在工程和工程教育方面非常流行的数学软件MATLAB来进行有限元的分析和应用。《MATLAB有限元分析与应用》由简单到复杂,循序渐进地介绍了各种有限元及其分析与应用方法。书中提供了大量取自机械工程、土木工程、航空航天工程和材料科学的示例和习题,具有很高的工程应用价值。
Matlab 有限元分析20140226 为了用Matlab 进行有限元分析,首先要学会Matlab 基本操作,还要学会使用Matlab 进行有限元分析的基本操作。 1. 复习:上节课分析了弹簧系统 x 推导了系统刚度矩阵 11221 21200k k k k k k k k -????-????--+??
2. Matlab有限元分析的基本操作 (1)单元划分(选择何种单元,分成多少个单元,标号)(2)构造单元刚度矩阵(列出…) (3)组装系统刚度矩阵(集成整体刚度矩阵) (4)引入边界条件(消除冗余方程) (5)解方程 (6)后处理(扩展计算)
3. Matlab有限元分析实战【实例1】
分析: 步骤一:单元划分
步骤二:构造单元刚度矩阵 >>k1=SpringElementStiffness(100) >>…?
步骤三:构造系统刚度矩阵 a) 分析SpringAssemble库函数function y = SpringAssemble(K,k,i,j) % This function assembles the element stiffness % matrix k of the spring with nodes i and j into the % global stiffness matrix K. % function returns the global stiffness matrix K % after the element stiffness matrix k is assembled. K(i,i) = K(i,i) + k(1,1); K(i,j) = K(i,j) + k(1,2); K(j,i) = K(j,i) + k(2,1); K(j,j) = K(j,j) + k(2,2); y = K; b) K是多大矩阵? 今天的系统刚度矩阵是什么? 因为 11 22 1212 k k k k k k k k - ?? ?? - ????--+ ?? 所以 1000100 0200200 100200300 - ?? ?? - ?? ?? -- ?? ?
有限元分析软件 有限元分析是对于结构力学分析迅速发展起来的一种现代计算方法。它是50 年代首先在连续体力学领域--飞机结构静、动态特性分析中应用的一种有效的数值分析方法,随后很快广泛的应用于求解热传导、电磁场、流体力学等连续性问题。 有限元分析软件目前最流行的有:ANSYS、ADINA、ABAQUS、MSC 四个比较知名比较大的公司,其中ADINA、ABAQUS 在非线性分析方面有较强的能力目前是业内最认可的两款有限元分析软件,ANSYS、MSC 进入中国比较早所以在国内知名度高应用广泛。目前在多物理场耦合方面几大公司都可以做到结构、流体、热的耦合分析,但是除ADINA 以外其它三个必须与别的软件搭配进行迭代分析,唯一能做到真正流固耦合的软件只有ADINA。ANSYS是商业化比较早的一个软件,目前公司收购了很多其他软件在旗下。ABAQUS专注结构分析目前没有流体模块。MSC是比较老的一款软件目前更新速度比较慢。ADINA是在同一体系下开发有结构、流体、热分析的一款软件,功能强大但进入中国时间比较晚市场还没有完全铺开。 结构分析能力排名:ABAQUS、ADINA、MSC、ANSYS 流体分析能力排名:ANSYS、ADINA、MSC、ABAQUS 耦合分析能力排名:ADINA、ANSYS、MSC、ABAQUS 性价比排名:最好的是ADINA,其次ABAQUS、再次ANSYS、最后MSC ABAQUS 软件与ANSYS 软件的对比分析: 1.在世界范围内的知名度:两种软件同为国际知名的有限元分析软件,在世界范围内具有各自广泛的用户群。ANSYS 软件在致力于线性分析的用户中具有很好的声誉,它在计算机资源的利用,用户界面开发等方面也做出了较大的贡献。ABAQUS软件则致力于更复杂和深入的工程问题,其强大的非线性分析功能在设计和研究的高端用户群中得到了广泛的认可。由于ANSYS 产品进入中国市场早于ABAQUS,并且在五年前ANSYS 的界面是当时最好的界面之一,所以在中国,ANSYS 软件在用户数量和市场推广度方面要高于ABAQUS。但随着ABAQUS北京办事处的成立,ABAQUS软件的用户数目和市场占有率正在大幅度和稳步提高,并可望在今后的几年内赶上和超过ANSYS。 2.应用领域:ANSYS 软件注重应用领域的拓展,目前已覆盖流体、电磁场和多物理场耦合等十分广泛的研究领域。ABAQUS 则集中于结构力学和相关领域研究,致力于解决该领域的深层次实际问题。 3.性价比:ANSYS 软件由于价格政策灵活,具有多种销售方案,在解决常规的
《建筑结构》补修课导学三 2008年06月17日 王启平 第三章 受弯构件承载力计算 受弯构件的两种破坏形式:1.沿弯矩最大截面破坏,称为正截面破坏;2.是沿剪力最大或弯矩和剪力都较大的截面破坏,破坏截面与构件的轴线斜交,称为斜截面破坏。 (a )正截面破坏 (b )斜截面破坏 图3-1 受弯构件的两种破坏形式 3.1一般构造要求 3.1.1截面形式 在受弯构件中,仅在截面的受拉区配置纵向受力钢筋的截面,称为单筋截面。同时在截面的受拉区和受压区配置纵向受力钢筋的截面,称为双筋截面。 3.1.2梁的构造要求 梁中一般配置纵向受力钢筋、弯起钢筋、箍筋和架立钢筋,如图3-3所示。 图 梁的配筋 1. 截面尺寸 梁高与跨度之比l h /称为高跨比。对于肋形楼盖的主梁为1/8~1/14,次梁为1/12~1/18;独立梁不小于1/15(简支)和1/20(连续)。 矩形截面梁的高宽比b h /一般取2.0~3.0;T 形截面梁的b h /.一般取2.5~4.0 (此处b 为梁肋宽)。为便于统一模板尺寸,通常采用矩形截面梁的宽度或T 形截面梁的肋宽b = 100、120、150、(180)、200、(220)、250和300mm ,300mm 以上的级差为50mm ,括号中的数值仅用于木模;梁的高度h = 250、300、750、800、900、1000mm 等尺寸。当 有限元分析软件的比较(购买必看)-转贴 随着现代科学技术的发展,人们正在不断建造更为快速的交通工具、更大规模的建筑物、更大跨度的桥梁、更大功率的发电机组和更为精密的机械设备。这一切都要求工程师在设计阶段就能精确地预测出产品和工程的技术性能,需要对结构的静、动力强度以及温度场、流场、电磁场和渗流等技术参数进行分析计算。例如分析计算高层建筑和大跨度桥梁在地震时所受到的影响,看看是否会发生破坏性事故;分析计算核反应堆的温度场,确定传热和冷却系统是否合理;分析涡轮机叶片内的流体动力学参数,以提高其运转效率。这些都可归结为求解物理问题的控制偏微分方程式,这些问题的解析计算往往是不现实的。近年来在计算机技术和数值分析方法支持下发展起来的有限元分析(FEA,Finite Element A nalysis)方法则为解决这些复杂的工程分析计算问题提供了有效的途径。在工程实践中,有限元分析软件与CAD系统的集成应用使设计水平发生了质的飞跃,主要表现在以下几个方面: 增加设计功能,减少设计成本; 缩短设计和分析的循环周期; 增加产品和工程的可靠性; 采用优化设计,降低材料的消耗或成本; 在产品制造或工程施工前预先发现潜在的问题; 模拟各种试验方案,减少试验时间和经费; 进行机械事故分析,查找事故原因。 在大力推广CAD技术的今天,从自行车到航天飞机,所有的设计制造都离不开有限元分析计算,FEA在工程设计和分析中将得到越来越广泛的重视。国际上早20世纪在50年代末、60年代初就投入大量的人力和物力开发具有强大功能的有限元分析程序。其中最为著名的是由美国国家宇航局(NASA)在1965年委托美国计算科学公司和贝尔航空系统公司开发的NASTRAN有限元分析系统。该系统发展至今已有几十个版本,是目前世界上规模最大、功能最强的有限元分析系统。从那时到现在,世界各地的研究机构和大学也发展了一批规模较小但使用灵活、价格较低的专用或通用有限元分析软件,主要有德国的ASKA、英国的PA FEC、法国的SYSTUS、美国的ABQUS、ADINA、ANSYS、BERSAFE、BOSOR、COSMOS、ELAS、MARC和STARDYNE等公司的产品。 以下对一些常用的软件进行一些比较分析: 1. LSTC公司的LS-DYNA系列软件 轴心受压构件承载力计算 按照箍筋配置方式不同,钢筋混凝土轴心受压柱可分为两种:一种是配置纵向钢筋和普通箍筋的柱(图4.2.1a),称为普通箍筋 柱;一种是配置纵向钢筋和螺旋筋(图)或 焊接环筋(图4.2.1c)的柱,称为螺旋箍筋柱或 间接箍筋柱。 需要指出的是,在实际工程结构中,几 乎不存在真正的轴心受压构件。通常由于荷 载作用位置偏差、配筋不对称以及施工误差 等原因,总是或多或少存在初始偏心距。但 当这种偏心距很小时,如只承受节点荷载屋 架的受压弦杆和腹杆、以恒荷载为主的等跨 多层框架房屋的内柱等,为计算方便,可近 似按轴心受压构件计算。此外,偏心受压构件垂直于弯矩作用平面的承载力验算也按轴心受压构件计算。 一、轴心受压构件的破坏特征 按照长细比的大小,轴心受压柱可分为短柱和长柱两类。对方形和矩形柱,当≤8时属于短柱,否则为长柱。其中为柱的计算长度,为矩形截面的短边尺寸。 1.轴心受压短柱的破坏特征 配有普通箍筋的矩形截面短柱,在轴向压力N作用下整个截面的应变基本上是均匀分布的。N较小时,构件的压缩变形主要为弹性变形。随着荷载的增大,构件变形迅速增大。与此同时,混凝土塑性变形增加,弹性模量降低,应力增长逐渐变慢,而钢筋应力的增加则越来越快。对配置HPB235、HRB335、HRB400、RRB400级热轧钢筋的构件,钢筋将先达到其屈服强度,此后增加的荷载全部由混凝土来承受。在临近 破坏时,柱子表面出现纵向裂缝,混凝土保护层开始剥落,最后,箍筋之间的纵向钢筋压屈而向外凸出,混凝土被压碎崩裂而破坏(图4.2.2)。破坏时混凝土的应力达到棱柱体抗压强度。当短柱破坏时,混凝土达到极限压应变=,相应的纵向钢筋应力值=E s=2×105×mm2=400N/mm2。因此,当纵向钢筋为高强度钢筋时,构件破坏时纵向钢筋可能达不到屈服强度。设计中对于屈服强度超过400N/mm2的钢筋,其抗压强度设计值只能取400N/mm2。显然,在受压构件内配置高强度的钢筋不能充分发挥其作用,这是不经济的。 2.轴心受压长柱的破坏特征 对于长细比较大的长柱,由于各种偶然因素造成的初始偏心距的影响是不可忽略的,在轴心压力N作用下,由初始偏心距将产生附加弯矩,而这个附加弯矩产生的水平挠度又加大了原来的初始偏心距,这样相互影响的结果,促使了构件截面材料破坏较早到来,导致承截能力的降低。破坏时首先在凹边出现纵向裂缝,接着混凝土被压碎,纵向钢筋被压弯向外凸出,侧向挠度急速发展,最终柱子失去平衡并将凸边混凝土拉裂而破坏(图4.2.3)。试验表明,柱的长细比愈大,其承截力愈低,对于长细比很大的长柱,还有可能发生“失稳破坏”。 由上述试验可知,在同等条件下,即截面相同,配筋相同,材料相同的条件下,长柱承载力低于短柱承载力。在确定轴心受压构件承截力计算公式时,规范采用构件 有限元大作业程序设计 学校:天津大学 院系:建筑工程与力学学院 专业:01级工程力学 姓名:刘秀 学号:\\\\\\\\\\\ 指导老师: 连续体平面问题的有限元程序分析 [题目]: 如图所示的正方形薄板四周受均匀载荷的作用,该结构在边界 上受正向分布压力, m kN p 1=,同时在沿对角线y 轴上受一对集中压 力,载荷为2KN ,若取板厚1=t ,泊松比0=v 。 [分析过程]: 由于连续平板的对称性,只需要取其在第一象限的四分之一部分参加分析,然后人为作出一些辅助线将平板“分割”成若干部分,再为每个部分选择分析单元。采用将此模型化分为4个全等的直角三角型单元。利用其对称性,四分之一部分的边界约束,载荷可等效如图所示。 [程序原理及实现]: 用FORTRAN程序的实现。由节点信息文件NODE.IN和单元信息文件ELEMENT.IN,经过计算分析后输出一个一般性的文件DATA.OUT。模型基本信息由文件为BASIC.IN生成。 该程序的特点如下: 问题类型:可用于计算弹性力学平面问题和平面应变问题 单元类型:采用常应变三角形单元 位移模式:用用线性位移模式 载荷类型:节点载荷,非节点载荷应先换算为等效节点载荷 材料性质:弹性体由单一的均匀材料组成 约束方式:为“0”位移固定约束,为保证无刚体位移,弹性体至少应有对三个自由度的独立约束 方程求解:针对半带宽刚度方程的Gauss消元法 输入文件:由手工生成节点信息文件NODE.IN,和单元信息文件ELEMENT.IN 结果文件:输出一般的结果文件DATA.OUT 程序的原理如框图: 岩土工程有限元分析软件 PLAXIS 2D 2015? 案例教程 北京市古城西街19号研发主楼4层,100043 目录 新奥法(NATM)隧道开挖 (1) 1.1 输入 (2) 1.2生成网格 (5) 1.3计算 (6) 1.4 结果 (8) 新奥法(NATM)隧道开挖 本例利用PLAXIS分析NATM隧道施工过程。NATM是在地下开挖时,利用喷射混凝土作为临时支护,保证开挖稳定性的一种施工方法。 图1.1 项目几何尺寸 目标: ●模拟NATM隧道施工(β法)。 ●用重力加载生成初始应力。 PLAXIS 2D AE案例教程:新奥法(NATM)隧道开挖 1.1 输入 1.1.1一般设置 ●打开PLAXIS 2D AE软件,在出现的快速选择对话框中选择一个新的项目。 ●在工程属性窗口的工程标签下,键入一个合适标题。 ●在模型标签下,模型(平面应变)和单元(15-Node)保持默认选项。 ●保持单位和一般设置框为默认值。 ●在几何形状设定框中设定土层模型尺寸xmin=-50,xmax=50,ymin=0,ymax=35。 ●点击OK即关闭工程属性窗口,完成设定。 1.1.2土层定义 利用钻孔生成土层,模型中考虑11m厚的泥灰岩,这层的底部y min=0作为参考点,定义土层: 在x=-22处创建第一个钻孔。 ●修改土层窗口将出现。为钻孔添加三层土。钻孔Borehole_1第一层的深度为0.指 定第一层土的顶部和底部值为24。第二层土层的顶部=24和底部=11.第三层土层的 顶部=11和底部=0。 ●单击在修改土层窗口的底部钻孔按钮。 ●在出现的菜单中选择添加选项。添加钻孔窗口出现。 ●指定第二个钻孔的位置为x=-14. ●注意:钻孔Borehole_1的特性复制给了Borehole_2。 ●Borehole_2第一层的深度也是0。修改土层的顶部=30和底部=30.第二层土顶部=30 和底部=11。第三层土顶部=11和底部=0. ●指定第三个钻孔的位置为x=-7. ●Borehole_3第一层土顶部=35和底部=30.第二层土顶部=30和底部=11,。第三层土 的顶部=11和底部=0. ●所有钻孔设置水头高度为y=0m。土层分布如图1.2。 ●根据表1.1定义土层材料属性,并分别指定给相应土层(图1.2). ●关闭修改土层窗口,切换到结构模式定义结构单元。 图1.2 土层分布 受压构件承载力计算复习题 一、填空题: 1、小偏心受压构件的破坏都是由于 而造成 的。 【答案】混凝土被压碎 2、大偏心受压破坏属于 ,小偏心破坏属 于 。 【答案】延性 脆性 3、偏心受压构件在纵向弯曲影响下,其破坏特征有两 种类型,对长细比较小的短柱属于 破坏,对长细比较大的细长柱,属于 破坏。 【答案】强度破坏 失稳 4、在偏心受压构件中,用 考虑了纵向弯曲的 影响。 【答案】偏心距增大系数 5、大小偏心受压的分界限是 。 【答案】b ξξ= 6、在大偏心设计校核时,当 时,说明s A '不屈 服。 【答案】s a x '2 7、对于对称配筋的偏心受压构件,在进行截面设计时, 和 作为判别偏心受压类型的唯一依据。 【答案】b ξξ≤ b ξξ 8、偏心受压构件 对抗剪有利。 【答案】轴向压力N 9、在钢筋混凝土轴心受压柱中,螺旋钢筋的作用是使截面中间核心部分的混凝土形成约束混凝土,可以提高构件的______和______。 【答案】承载力 延性 10、偏心距较大,配筋率不高的受压构件属______受压情况,其承载力主要取决于______钢筋。 【答案】大偏心 受拉 11、受压构件的附加偏心距对______受压构件______受压构件影响比较大。 【答案】轴心 小偏心 12、在轴心受压构件的承载力计算公式中,当f y <400N /mm 2 时,取钢筋抗压强度设计值f y '=______;当f y ≥400N /mm 2时,取钢筋抗压强度设计值f y '=______N /mm 2。 【答案】f y 400 二、选择题: 1、大小偏心受压破坏特征的根本区别在于构件破坏时,( )。 A 受压混凝土是否破坏 B 受压钢筋是否屈服 C 混凝土是否全截面受压 D 远离作用力N 一侧钢筋是否屈服 Matlab有限元分析20140226 为了用Matlab进行有限元分析,首先要学会Matlab基本操作,还要学会使用Matlab进行有限元分析的基本操作。 1. 复习:上节课分析了弹簧系统 x 推导了系统刚度矩阵 2. Matlab有限元分析的基本操作 (1)单元划分(选择何种单元,分成多少个单元,标号)(2)构造单元刚度矩阵(列出…) (3)组装系统刚度矩阵(集成整体刚度矩阵) (4)引入边界条件(消除冗余方程) (5)解方程 (6)后处理(扩展计算) 3. Matlab有限元分析实战【实例1】 分析: 步骤一:单元划分 >>k1=SpringElementStiffness(100) a) 分析SpringAssemble库函数 function y = SpringAssemble(K,k,i,j) % This function assembles the element stiffness % matrix k of the spring with nodes i and j into the % global stiffness matrix K. % function returns the global stiffness matrix K % after the element stiffness matrix k is assembled. K(i,i) = K(i,i) + k(1,1); K(i,j) = K(i,j) + k(1,2); K(j,i) = K(j,i) + k(2,1); K(j,j) = K(j,j) + k(2,2); y = K; b) K是多大矩阵? 今天的系统刚度矩阵是什么? 因为 11 22 1212 k k k k k k k k - ?? ?? - ????--+ ?? 所以 1000100 0200200 100200300 - ?? ?? - ????-- ??? 有限元软件与多体动力学软件 数值分析技术与传统力学的结合在结构力学领域取得了辉煌的成就,出现了以ANSYS 、NASTRAN 等为代表的应用极为广泛的结构有限元分析软件。计算机技术在机构的静力学分析、运动学分析、动力学分析以及控制系统分析上的应用,则在二十世纪八十年代形成了计算多体系统动力学,并产生了以ADAMS 和DADS 为代表的动力学分析软件。两者共同构成计算机辅助工程(CAE )技术的重要内容。 商业通用软件的广泛应用给我们工程师带来了极大的便利,很多时候我们不需要精通工程问题中的力学原理,依然可以通过商业软件来解决问题,不过理论基础的缺失还是会给我们带来不少的困扰。随着动力有限元与柔性多体系统分析方法的成熟,有时候正确区分两者并不是很容易。 机械领域应用比较广泛的有两类软件,一类是有限元软件,代表的有:ANSYS, NASTRAN, ABAQUS, LS-DYNA, Dytran 等;另一类是多体动力学软件,代表的有ADAMS, Recurdyn , Simpack 等。在使用时,如何选用这两类软件并不难,但是如果深究这两类软件根本区别并不容易。例如,有限元软件可以分析静力学问题,也可以分析“动力学”问题,这里的“动力学”与多体动力学软件里面的动力学一样吗?有限元软件在分析动力学问题时,可以模拟物体的运动,它与多体动力学软件中模拟物体运动相同吗?多体动力学软件也可以分析柔性体的应力、应变等,这与有限元软件分析等价吗? 1 有限元软件 有限单元法是一种数学方法,不仅可以计算力学问题,还可以计算声学,热,磁等多种问题,我们这里只探讨有限元法在机械领域的应用。 计算结构应力、应变等的力学基础是弹性力学,弹性力学亦称为弹性理论,主要研究弹性体在外力作用或温度变化等外界因素下所产生的应力、应变和位移,从而为工程结构或构件的强度、刚度设计提供理论依据和计算方法。也就是说用有限元软件分析力学问题时,是用有限元法计算依据弹性力学列出的方程。 考虑下面这个问题,在()0t , 时间内给一个结构施加一个随时间变化的载荷()P t ,我们希望得到结构的应力分布,在刚刚施加载荷的时候,结构中的应力会有波动,应力场是变化的,但很久以后,应力场趋于稳定。 如果我们想得到载荷施加很久以后,稳定的应力场分布,那么应该用静力学分析方法分析 轴心受力构件习题及答案 一、选择题 的构件,在拉力N作用下的强度计算公1. 一根截面面积为A,净截面面积为A n 式为______。 2. 轴心受拉构件按强度极限状态是______。 净截面的平均应力达到钢材的抗拉强度 毛截面的平均应力达到钢材的抗拉强度 净截面的平均应力达到钢材的屈服强度 毛截面的平均应力达到钢材的屈服强度 3. 实腹式轴心受拉构件计算的容有______。 强度强度和整体稳定性强度、局部稳定和整体 稳定强度、刚度(长细比) 4. 轴心受力构件的强度计算,一般采用轴力除以净截面面积,这种计算方法对下列哪种连接方式是偏于保守的? 摩擦型高强度螺栓连接承压型高强度螺栓连 接普通螺栓连接铆钉连接 5. 工字型组合截面轴压杆局部稳定验算时,翼缘与腹板宽厚比限值是根据 ______导出的。 6. 图示单轴对称的理想轴心压杆,弹性失稳形式可能为______。 X轴弯曲及扭转失稳Y轴弯曲及扭转失稳 扭转失稳绕Y轴弯曲失稳 7. 用Q235号钢和16锰钢分别建造一轴心受压柱,其长细比相同,在弹性围屈曲时,前者的临界力______后者的临界力。 大于小于等于或接近无法 比较 8. 轴心受压格构式构件在验算其绕虚轴的整体稳定时采用换算长细比,是因为______。 格构构件的整体稳定承载力高于同截面的实腹构件 考虑强度降低的影响 考虑剪切变形的影响 考虑单支失稳对构件承载力的影响 9. 为防止钢构件中的板件失稳采取加劲措施,这一做法是为了______。 改变板件的宽厚比增大截面面积改变截面上 的应力分布状态增加截面的惯性矩 10. 轴心压杆构件采用冷弯薄壁型钢或普通型钢,其稳定性计算______。 完全相同 仅稳定系数取值不同 仅面积取值不同 完全不同 11. 工字型截面受压构件的腹板高度与厚度之比不能满足按全腹板进行计算的要求时,______。 姓名:刘刚学号:15 平面应力应变分析有限元法 Abstruct:本文通过对平面应力/应变问题的简要理论阐述,使读者对要分析的问题有大致的印象,然后结合两个实例,通过MATLAB软件的计算,将有限元分析平面应力/应变问题的过程形象的展示给读者,让人一目了然,快速了解有限元解决这类问题的方法和步骤! 一.基本理论 有限元法的基本思路和基本原则以结构力学中的位移法为基础,把复杂的结构或连续体看成有限个单元的组合,各单元彼此在节点出连接而组成整体。把连续体分成有限个单元和节点,称为离散化。先对单元进行特性分析,然后根据节点处的平衡和协调条件建立方程,综合后做整体分析。这样一分一合,先离散再综合的过程,就是把复杂结构或连续体的计算问题转化简单单元分析与综合问题。因此,一般的有限揭发包括三个主要步骤:离散化单元分析整体分析。 二.用到的函数 1. LinearTriangleElementStiffness(E,NU,t,xi,yi,xj,yj,xm,ym,p) (K k I f) (k u) (k u A) (E NU t) 三.实例 例1.考虑如图所示的受均布载荷作用的薄平板结构。将平板离散化成两个线性三角元,假定E=200GPa,v=,t=0.025m,w=3000kN/m. 1.离散化 2.写出单元刚度矩阵 通过matlab 的LinearTriangleElementStiffness 函数,得到两个单元刚度矩阵1k 和2k ,每个矩阵都是6 6的。 >> E=210e6 E = >> k1=LinearTriangleElementStiffness(E,NU,t,0,0,,,0,,1) k1 = +006 * Columns 1 through 5 0 0 0 0 0 0 0 0 Column 6 >> NU= NU = >> t= t = >> k2=LinearTriangleElementStiffness(E,NU,t,0,0,,0,,,1) 各种有限元分析软件比较 有限元分析(FEA,Finite Element Analysis)利用数学近似的方法对真实物理系统(几何和载荷工况)进行模拟。还利用简单而又相互作用的元素,即单元,就可以用有限数量的未知量去逼近无限未知量的真实系统。 有限元分析具有确保产品设计的安全合理性,同时采用优化设计,找出产品设计最佳方案,降低材料的消耗或成本; 在产品制造或工程施工前预先发现潜在的问题; 模拟各种试验方案,减少试验时间和经费等作用,越来越被应用,越来越的有限元分析也不断被开发出来,当我们在做有限元分析时,我们该选择什么样的软件?或者我们该学习什么软件?成了大多数人困惑的问题。看板网根据自己超过十年的有限元分析项目经验和培训经验,对各种有限元分析软件进行了一些比较,希望大家在选择时能够大家做参考。 有限元分析常用软件 国外软件 大型通用有限元商业软件:如ANSYS可以分析多学科的问题,例如:机械、电磁、热力学等;电机有限元分析软件NASTRAN等。还有三维结构设计方面的UG,CATIA,Proe等都是比较强大的。 国内软件 国产有限元软件:FEPG,SciFEA,JiFEX,KMAS等。 当然首先要明确你要用这个软件进行什么分析,一般会用到有限元分析的地方有:1.模流分析;2.结构强度分析;3.电磁场分析;4.谐响应分析(比如查找共振频率);5. 铸造分析。等等 ANSYS是商业化比较早的一个软件,目前公司收购了很多其他软件在旗下。ABAQUS专注结构分析目前没有流体模块。MSC是比较老的一款软件目前更新速度比较慢。ADINA是在同一体系下开发有结构、流体、热分析的一款软件,功能强大但进入中国时间比较晚市场还没有完全铺开。 workbench是一个综合性的有限元分析软件,几乎囊括了所有有限元分析领域,传统的优势领域有强度分析、谐响应分析和电磁分析。workbench是ansys 3受弯构件承载力计 算 1 、一般构造要求 受弯构件正截面承载力计算 1 、配筋率对构件破坏特征的影响及适筋受弯构件截面受力的几个阶段 受弯构件正截面破坏特征主要由纵向受拉钢筋的配筋率ρ大小确定。配筋率是指纵受受拉钢筋的截面面积与截面的有效面积之比。 (3-1) 式中As——纵向受力钢筋的截面面积,; b——截面的宽度,mm; ——截面的有效高度, ——受拉钢筋合力作用点到截面受拉边缘的距离。 根据梁纵向钢筋配筋率的不同,钢筋混凝土梁可分为适筋梁、超筋梁和少筋梁三种类型,不同类型梁的破坏特征不同。 (1)适筋梁 配置适量纵向受力钢筋的梁称为适筋梁。 适筋梁从开始加载到完全破坏,其应力 变化经历了三个阶段,如图3.8。 第I阶段(弹性工作阶段):荷载很小 时,混凝土的压应力及拉应力都很小, 梁截面上各个纤维的应变也很小,其应 力和应变几乎成直线关系,混凝土应力 分布图形接近三角形,如图3.8(a)。 当弯矩增大时,混凝土的拉应力、压应 力和钢筋的拉应力也随之增大。由于混 凝土抗拉强度较低,受拉区混凝土开始 表现出明显的塑性性质,应变较应力增 加快,故应力和应变不再是直线关系, 应力分布呈曲线, 当弯距增加到开裂弯距时,受拉边缘纤维的应变达到混凝土的极限拉应变,此时, 截面处于将裂未裂的极限状态,即第I阶段末,用Ia表示,如图3.13(b)所示。这时受压区塑性变形发展不明显,其应力图形仍接近三角形。Ia阶段的应力状态是抗裂验算的依据。 第Ⅱ阶段(带裂缝工作阶段):当弯矩继续增加时,受拉区混凝土的拉应变超过其极其拉应变,受拉区 出现裂缝,截面即进入第Ⅱ阶段。裂缝出现后,在裂缝截面处,受拉区混凝土大部分退出工作,未开裂部分混凝土虽可继续承担部分拉力,但因靠近中和轴很近,故其作用甚小,拉力几乎全部由受拉钢筋承担,在裂缝出现的瞬间,钢筋应力突然增加很大。随着弯矩的不断增加,裂缝逐渐向上扩展,中和轴逐渐上移。由于受压区应变不断增大,受压区混凝土呈现出一定的塑性特征,应力图形呈曲线形,如图3.8?所示。第Ⅱ阶段的应力状态代表了受弯构件在使用时的应力状态,故本阶段的应力状态作为裂缝宽度和变形验算的依据。 当弯矩继续增加,钢筋应力不断增大,直至达到屈服强度,这时截面所能承担的弯矩称为屈服弯矩。 它标志截面即将进入破坏阶段,即为第Ⅱ阶段极限状态,以Ⅱa表示,如图3.8(d)所示。 第Ⅲ阶段(破坏阶段):弯矩继续增加,截面进入第Ⅲ阶段。这时受拉钢筋的应力保持屈服强度不变,钢筋的应变迅速增大,促使受拉区混凝土的裂缝迅速向上扩展,中和轴继续上移,受压区混凝土高度缩小,混凝土压应力迅速增大,受压区混凝土的塑性特征表现得更加充分,压应力呈显著曲线分布[图3.8(e)]。到本阶段末(即Ⅲa阶段),受压边缘混凝土压应变达到极限应变,受压区混凝土产生近乎水平的裂缝,混凝土被压碎,甚至崩脱[图3.8(a)],截面宣告破坏,此时截面所承担的弯矩即为破坏弯矩Mu,这时的应力状态作为构件承载力计算的依据[图3.8(f)]。 板结构有限元分析实例详解1:带孔平板结构静力分析本节介绍带孔平板结构静力分析问题,同时介绍布尔操作的基本用法。 8.3.1 问题描述与分析 有孔的矩形平板,左侧边缘固定,长400mm,宽200 mm,厚度为10 mm,圆孔在板的正中心,半径为40 mm,左侧全约束,右侧边缘均布应力1MPa,如图8.7所示。求板的变形、位移及应力变化情况。(材料的材料属性为:弹性模量为300000 MPa,剪切模量为0.31。) 图8.7 带孔的矩形平板 由于小孔处边缘不规则,本文采用PLANE82高阶平面单元进行分析。 8.3.2 求解过程 8.3.2.1 定义工作目录及文件名 启动ANSYS Mechanical APDL Product Launcher窗口,如图8.8所示。在License下 拉选框中选择ANSYS Multiphysics产品,在Working Directory输入栏中输入工作目 录:C:\ANSYS12.0 Structural Finite Elements Analysis and Practice\Chapter 8\8-1,在Job Name一栏中输入工作文件名:Chapter8-1。以上参数设置完毕后,单 击Run按钮运行ANSYS。 图8.8 ANSYS设置窗口菜单 可以先在目标文件位置建立工作目录,然后单击Browse按钮选择工作目录;也 可以通过单击Browse按钮选择工作文件名。 8.3.2.2 定义单元类型和材料属性 选择Main Menu>Preferences命令,出现Preferences for GUI Filtering对话框, 如图8.9所示,在Individual discipline(s) to show in the GUI中勾选Structural,过滤掉ANSYS GUI菜单中与结构分析无关的选项,单击OK按钮关闭该对话框。 图8.9 Preferences for GUI Filtering对话框 1.物理现象:这个对工程师来说是直观的物理现象和物理量,温 度多少度,载荷是多大等等。通常来说,用户界面中呈现的、用户对工程问题进行设置时输入的都是此类信息。 2.数学方程:将物理现象翻译成相应的数学方程,例如流体对应 的是NS方程,传热对应的是传热方程等等;大部分描述这些现象的方程在空间上都是偏微分方程,偶尔也有ODE(如粒子轨迹、化学反应等)。在这个层面,软件把物理现象“翻译” 为以解析式表示的数学模型。 3.数值模型:在定义了数学模型,并执行了网格剖分后,商业软 件会将数学模型离散化,利用有限元方法、边界元法、有限差分法、不连续伽辽金法等方法生成数值模型。软件会组装并计算方程组雅可比矩阵,并利用求解器求解方程组。这个层面的计算通常是隐藏在后台的,用户只能通过一些求解器的参数来干预求解。 有限元是一种数值求解偏微分方程的方法。 基本过程大致是设置形函数,离散,形成求解矩阵,数值解矩阵,后处理之类的。 MATLAB要把这些过程均自己实现,不过在数值求解矩阵时可以调用已有函数。可以理解为MATLAB是一个通用的计算器,当然它的功能远不止如此。 而ANSYS之类的叫做通用有限元软件,针对不同行业已经将上述过程封装,前后处理也比较漂亮,甚至不太了解有限元理论的人也能算些简单的东西,当然结果可靠性又另说了。 比较两者,ANSYS之类的用起来容易得多,但灵活性不如MATLAB。MATLAB用起来很困难,也有人做了一些模块,但大多数只能解决一些相对简单的问题。 对于大多数工程问题,以及某些领域的物理问题,一般都用通用有限元软件,这些软件还能添加一些函数块,用以解决一些需要额外设置的东西。但是对于非常特殊的问题,以及一般性方程的有限元解,那只能用MATLAB或C,Fortran之类的了。 随着现代科学技术的发展,人们正在不断建造更为快速的交通工具、更大规模的建筑物、更大跨度的桥梁、更大功率的发电机组和更为精密的机械设备。这一切都要求工程师在设计阶段就能精确地预测出产品和工程的技术性能,需要对结构的静、动力强度以及温度场、流场、电磁场和渗流等技术参数进行分析计算。例如分析计算高层建筑和大跨度桥梁在地震时所受到的影响,看看是否会发生破坏性事故;分析计算核反应堆的温度场,确定传热和冷却系统是否合理;分析涡轮机叶片内的流体动力学参数,以提高其运转效率。这些都可归结为求解物理问题的控制偏微分方程式,这些问题的解析计算往往是不现实的。近年来在计算机技术和数值分析方法支持下发展起来的有限元分析(FEA,Finite Element Analysis)方法则为解决这些复杂的工程分析计算问题提供了有效的途径。在工程实践中,有限元分析软件与CAD系统的集成应用使设计水平发生了质的飞跃,主要表现在以下几个方面: 增加设计功能,减少设计成本; 缩短设计和分析的循环周期; 增加产品和工程的可靠性; 采用优化设计,降低材料的消耗或成本; 在产品制造或工程施工前预先发现潜在的问题; 模拟各种试验方案,减少试验时间和经费; 进行机械事故分析,查找事故原因。 在大力推广CAD技术的今天,从自行车到航天飞机,所有的设计制造都离不开有限元分析计算,FEA 在工程设计和分析中将得到越来越广泛的重视。国际上早20世纪在50年代末、60年代初就投入大量的人力和物力开发具有强大功能的有限元分析程序。其中最为著名的是由美国国家宇航局(NASA)在1965年委托美国计算科学公司和贝尔航空系统公司开发的NASTRAN有限元分析系统。该系统发展至今已有几十个版本,是目前世界上规模最大、功能最强的有限元分析系统。从那时到现在,世界各地的研究机构和大学也发展了一批规模较小但使用灵活、价格较低的专用或通用有限元分析软件,主要有德国的ASKA、英国的PAFEC、法国的SYSTUS、美国的ABQUS、ADINA、ANSYS、BERSAFE、BOSOR、COSMOS、ELAS、MARC和STARDYNE等公司的产品。 以下对一些常用的软件进行一些比较分析: 1. LSTC公司的LS-DYNA系列软件 LS-DYNA是一个通用显式非线性动力分析有限元程序,最初是1976年在美国劳伦斯利弗莫尔国家实验室(Lawrence Livermore National Lab.)由J.O.Hallquist 主持开发完成的,主要目的是为核武器的弹头设计提供分析工具,后经多次扩充和改进,计算功能更为强大。此软件受到美国能源部的大力资助以及世界十余家著名数值模拟软件公司(如ANSYS、MSC.software、ETA等)的加盟,极大地加强了其的前后处理能力和通用性,在全世界范围内得到了广泛的使用。在软件的广告中声称可以求解各种三维非线性结构的高速碰撞、爆炸和金属成型等接触非线性、冲击载荷非线性和材料非线性问题。即使是这样一个被人们所称道的数值模拟软件,实际上仍在诸多不足,特别是在爆炸冲击方面,功能相对较弱,其欧拉混合单元中目前最多只能容许三种物质,边界处理很粗糙,在拉格朗日——欧拉结合方面不如DYTRAN灵活。虽然提供了十余种岩土介质模型,但每种模型都有不足,缺少基本材料数据和依据,让用户难于选择和使用。2. MSC.software公司的DYTRAN软件 当前另一个可以计算侵彻与爆炸的商业通用软件是MSC.Software Corporation ( MSC公司) 的MSC.DYTR AN程序。该程序在是在LS-DYNA3D的框架下,在程序中增加荷兰PISCES INTERNATIONAL公司开发的PICSES的高级流体动力学和流体——结构相互作用功能,还在PISCES的欧拉模式算法基础上,开发了物质流动算法和流固耦合算法。在同类软件中,其高度非线性、流—固耦合方面有独特之处。MSC.DYTR AN的算法基本上可以概况为:MSC.DYTRAN采用基于Lagrange格式的有限单元方法(FEM)模拟结构的变形和应力,用基于纯Euler格式的有限体积方法(FVM)描述材料(包括气体和液体)流动,对通过流体与固体界面传递相互作用的流体—结构耦合分析,采用基于混合的Lagrange格式和纯Euler 格式的有限单元与有限体积技术,完成全耦合的流体-结构相互作用模拟。MSC.DYTRAN用有限体积法跟踪(完整word版)有限元分析软件的比较
4.2 轴心受压构件承载力计算
(完整版)有限元大作业matlab---课程设计例子
岩土工程有限元分析软件
受压构件承载力计算复习题(答案)详解
Matlab有限元分析操作基础共11页
多体动力学软件和有限元软件的区别(优.选)
轴心受力构件习题及问题详解
基于matlab的有限元法分析平面应力应变问题刘刚
各种有限元分析软件比较
最新3受弯构件承载力计算汇总
板结构有限元分析实例详解
matlab有限元分析实例
主流CAE有限元分析软件的比较