ANSYS-中使用函数加载的一个简单例子
ansys加载
表面效应单元
• 怎样施加如下的压力荷载:
– 象剪切荷载一样与表面相切的荷载? – 象螺栓荷载一样在表面上变化的荷载? – 象屋顶上风载荷一样与面成一定角度的载荷?
• 表面效应单元为处理一些问题提供了有效的方法。
ANSYS TRAINING
...表面效应单元
• 特点:
– 象“皮肤”一样覆盖在网格表面 – 如同面载荷的管道 – 很容易创建:
ANSYS TRAINING
自由度约束
自由度约束就是给某个自由度(DOF)指定一已知数 值 (值不一定是零)。
定义
ANSYS TRAINING
集中载荷
集中载荷 就是作用在模型的一个点上的载荷。
定义
ANSYS TRAINING
集中载荷 (续)
举例
• 结构分析中的力和弯矩。 • 热分析中热流率。 • 集中载荷可以添加到节点和关键点上。(添加到关键点
载。如下图位置1轮廓线所示。
ANSYS TRAINING
...表面效应单元
– KEYOPT(11)=2 对整个面施加压力。
• 对斜面有用(如屋顶)或风载荷。 • 如: sfe,eslope,5,pres,,1000,0,-1,0 对整个斜面定义了完全相同
的压力。如下图所示
ANSYS TRAINING
惯性载荷
惯性载荷 是由物体的惯性(质量矩阵)引起的载荷, 例如重力加速度,加速度,以及角加速度。
定义
ANSYS TRAINING
惯性载荷(续)
• 惯性载荷只有结构分析中有。 • 惯性载荷是对整个结构定义的,
是独立于实体模型和有限元模型 的。
绕Y轴的角速度
• 考虑惯性载荷就必须定义材料密 度 (材料特性DENS)。
ansys-一些函数命令
在ANSYS帮助系统中关于*SET命令的注释下列出了ANSYS中可以使用的数学函数。
所有这些数学函数均可以在ANSYS环境中使用,这些数学函数包括:ABS(X) 求绝对值ACOS(X) 反余弦ASIN(X) 反正弦ATAN(X) 反正切ATAN2(X,Y) 反正切, ArcTangent of (Y/X) , 可以考虑变量X,Y 的符号COS(X) 求余弦COSH(X) 双曲余弦EXP(X) 指数函数GDIS(X,Y) 求以X为均值,Y为标准差的高斯分布,在使用蒙地卡罗法研究随机荷载和随机材料参数时,可以用该函数处理计算结果LOG(X) 自然对数LOG10(X) 常用对数(以10为基)MOD(X,Y) 求 X/Y的余数. 如果 Y=0, 函数值为 0NINT(X) 求最近的整数RAND(X,Y) 取随机数,其中X 是下限, Y是上限SIGN(X,Y) 取 X的绝对值并赋予Y的符号. Y>=0, 函数值为|X|, Y<0, 函数值为-|X|,.SIN(X) 正弦SINH(X) 双曲正弦SQRT(X) 平方根TAN(X) 正切TANH(X) 双曲正切2004-4-15在ANSYS帮助系统中关于*SET命令的注释下列出了ANSYS中可以使用的数学函数。
所有这些数学函数均可以在ANSYS环境中使用,这些数学函数包括:ABS(X) 求绝对值ACOS(X) 反余弦ASIN(X) 反正弦ATAN(X) 反正切ATAN2(X,Y) 反正切, ArcTangent of (Y/X) , 可以考虑变量X,Y 的符号COS(X) 求余弦COSH(X) 双曲余弦EXP(X) 指数函数GDIS(X,Y) 求以X为均值,Y为标准差的高斯分布,在使用蒙地卡罗法研究随机荷载和随机材料参数时,可以用该函数处理计算结果LOG(X) 自然对数LOG10(X) 常用对数(以10为基)MOD(X,Y) 求 X/Y的余数. 如果 Y=0, 函数值为 0NINT(X) 求最近的整数RAND(X,Y) 取随机数,其中X 是下限, Y是上限SIGN(X,Y) 取 X的绝对值并赋予Y的符号. Y>=0, 函数值为|X|, Y<0, 函数值为-|X|,.SIN(X) 正弦SINH(X) 双曲正弦SQRT(X) 平方根TAN(X) 正切TANH(X) 双曲正切有限元分析软件ANSYS命令流中文说明CommandVSBV, NV1, NV2, SEPO, KEEP1, KEEP2 — Subtracts volumes from volumes,用于2个solid相减操作,最终目的是要nv1-nv2=?通过后面的参数设置,可以得到很多种情况:sepo 项是2个体的边界情况,当缺省的时候,是表示2个体相减后,其边界是公用的,当为sepo 的时候,表示相减后,2个体有各自的独立边界。
ANSYS110基础与实例教程第4章加载与后处理(1)精品PPT课件
场分析中计算的磁力作为结构分析中的力载荷。
学 系
《力学分析应用软件基础》多媒体课件 2020/10/15
11/12第二学期 4
ANSYS11.0基础与实例教程
载荷步和子步
河
第 四 1)载荷步(load step)
南 科 技
章
载荷步就是我们平时讲的分步施加载荷,以模拟真实的载荷配
大 学
加 载 与
置。下图所示显示了一个需要三个载荷步的载荷历程曲线:第一个 载荷步用于线性载荷,第二个载荷步用于不变载荷,第三个载荷步 用于卸载。载荷值在载荷部的结束点达到全值。
程 学 院 力
步。
学
系
《力学分析应用软件基础》多媒体课件 2020/10/15
11/12第二学期 2
ANSYS11.0基础与实例教程
关于载荷的一些概念
河
在ANSYS的术语中,载荷(Loads)包括边界条件和外部或者内部 南
第四 作用力。在不同的学科中,载荷的具体含义也不尽相同,下面为不同
科 技
章 学科中所指的载荷术语。
理 面上的载荷需要逐个单元来拾取,非常麻烦。
学 院
力
学
系
《力学分析应用软件基础》多媒体课件 2020/10/15
11/12第二学期 8
ANSYS11.0基础与实例教程
第 实体模型载荷 四 • 易于施加
有限元模型载荷
河 南
科
• 不便施加
技
章 • 与有限元模型独立
加 载
• 坐标系与节点和单
• 跟随有限元模型 • 载荷显示比较直观
程
理 •缺点:
学 院
不能显示所有的实体模型载荷;关键点上的载荷很难施加。
ansys 加载
题目:加载2.1 载荷概述有限元分析的主要目的是检查结构或构件对一定载荷条件的响应。
因此,在分析中指定合适的载荷条件是关键的一步。
在ANSYS程序中,可以用各种方式对模型加载,而且借助于载荷步选项,可以控制在求解中载荷如何使用。
2.2 什么是载荷在ANSYS术语中,载荷(loads)包括边界条件和外部或内部作用力函数,如图2-1所示。
不同学科中的载荷实例为:结构分析:位移,力,压力,温度(热应变),重力热分析:温度,热流速率,对流,内部热生成,无限表面磁场分析:磁势,磁通量,磁场段,源流密度,无限表面电场分析:电势(电压),电流,电荷,电荷密度,无限表面流体分析:速度,压力图2-1 “载荷”包括边界条件以及其它类型的载荷载荷分为六类:DOF约束,力(集中载荷),表面载荷,体积载荷、惯性力及耦合场载荷。
·DOF constraint(DOF约束)将用一已知值给定某个自由度。
例如,在结构分析中约束被指定为位移和对称边界条件;在热力分析中指定为温度和热通量平行的边界条件。
·Force(力)为施加于模型节点的集中载荷。
例如,在结构分析中被指定为力和力矩;在热力分析中为热流速率;在磁场分析中为电流段。
·Surface load(表面载荷)为施加于某个表面上的分布载荷。
例如,在结构分析中为压力;在热力分析中为对流和热通量。
·Body load(体积载荷)为体积的或场载荷。
例如,在结构分析中为温度和fluences;在热力分析中为热生成速率;在磁场分析中为流密度。
·Inertia loads(惯性载荷)由物体惯性引起的载荷,如重力加速度,角速度和角加速度。
主要在结构分析中使用。
·Coupled-field loads(耦合场载荷)为以上载荷的一种特殊情况,从一种分析得到的结果用作为另一分析的载荷。
例如,可施加磁场分析中计算出的磁力作为结构分析中的力载荷。
ANSYS-中使用函数加载的一个简单例子
ANSYS-中使用函数加载的一个简单例子ANSYS 中使用函数加载的一个简单例子本文将通过一个具体实例说明在ANSYS 中如何使用函数加载,后续将通过该实例在分析过程中遇到的一个问题提出自己的一点看法。
实例的具体说明:一个1/4 圆柱,内半径30 mm,外半径42 mm,长度100mm,如图 1 所示:所用材料为双线性弹塑性材料,其机械性能为:P (y) = 8e7 + 7E7 * (Y/42)即:X = 0 ,Y = 42 (最高点) 时,P = 15E7;X = 42,Y = 0 (最低点)时,P = 8E7。
我们采用函数方式来施加这一压力载荷,首先定义函数:在Solution 模块中,点击菜单路径:Solution > Define Loads > Apply > Functions >Define/Edit将会弹出一个函数编辑器,可以在其中定义所需的函数。
在函数编辑器中,函数类型选择为Single equation,即单值函数;计算函数值时使用的插值坐标系( (x,y,z) interpreted in CSYS) 选择0,即总体直角坐标系,如图3 所示:然后,在函数编辑器中间位置的“Result = “ 小窗口中输入要定义的函数表达式,如果表达式中有x, y, z, time 等变量(供定义函数时使用的“自变量”),可以用{X},{Y},{Z},{TIME} 等的形式输入;或者点击下面一个小窗口右边的小箭头,会出现一个下拉列表,列出可以选择的变量,然后从该列表中选择某个自变量,则该自变量会按照上述格式写入函数中,如图 5 所示:接下来最好检查一下函数定义是否正确。
做法是:点击函数编辑器中的Graph 按钮,弹出一个“Plot information” 对话框,其中显示了所定义的函数方程,列出了方程中所有的变量。
然后可以激活要变化的变量(对于本函数,只有一个变量– X,会自动激活),并在下面的输入窗口中输入该变量的最小和最大值,以及在该变量范围中的计算点数,再点击Graph 按钮,则会绘制函数的曲线;或者点击List 按钮,则列出各计算点上的函数值。
ANSYS瞬态动力学时间函数加载
•施加已有表载荷
•函数载荷步控制
•基本变量:程序中使用的独立变量。方程式中允许的基本 变量有TIME, X, Y, Z,和TEMP等。 •状态变量:定义用来管理分段函数的变量。
函数边界条件加载过程: •定义函数关系式 •转换为表格载荷 •施加已有表载荷 •函数载荷步控制
Force
22.5 10 0.5 1.0 1.5 t
顶面压力载荷P=1000*sin(2πt) (a)单段函数
顶面压力载荷 (b)分段函数
二、单段函数载荷步加载
命令流
顶面压力载荷P=1000*sin(2πt)
•定义函数关系式 •转换为表格载荷
•施加已有表载荷
•函数载荷步控制
顶面位移-时间曲线
•定义函数关系式
在函数编辑器 中使用一组基本变 量,方程变量,和 数学函数去建构函 数方程。
1
2
•定义函数的详细过程: 1.打开函数编辑器: 2.选择函数类型 3.定义函数方程 4.绘制函数图形
•定义函数关系式 •转换为表格载荷
•施加已有表载荷
•函数载荷步控制
顶面位移-时间曲线
•定义函数关系式
•定义函数的详细过程:
1.打开函数编辑器: 2.选择函数类型 3.定义自变量
1
2
4.定义各函数段
5.保存函数文件
3
5
4
•转换为表格载荷
在分析中,你准备为方程变量指定值,指 定表格型参数名并使用函数时。 1.载入函数: 2.定义表格: 3.保存表格:43来自5.保存函数文件5
•转换为表格载荷
在分析中,你准备为方程变量指定值,指 定表格型参数名并使用函数时。 1.载入函数: 2.定义表格: 3.保存表格:
ansys加载函数载荷func
1.在经典ANSYS里面Apply/Functions/Define/Edit打开函数编辑器●Functions Type:选择函数类型。
选择单个方程或多值函数。
如果选择后者,必须键入状态变量名,也就是管理函数中方程的变量。
当选择一个多值函数时,六状态表格将被激活。
●Degr ees/Radians:选择度或弧度,这一选择仅决定方程如何被运算,而不会影响*AFUN设置。
●使用初始变量方程和键区定义结果方程(单个方程)或描述状态变量的方程(多值函数),出如果定义单方程函数,保存方程。
如果是定义多值函数,则继续下面的步骤。
●单击Regime1,键入在函数表格下定义的状态变量的相应的最大最小值限制。
●定义这个状态的方程。
●单击Regime2,注意状态变量的最小值限制已被定义并且不可更改,这一特征确保状态保持连续而无间隙。
定义这个状态的最高值限制。
●定义这个状态的方程。
●在六个状态中连续如上操作。
在每个状态里,不必储存或保存单个方程,除非想在另一状态中重用某个方程。
●输入一个注释描述函数(可选)选择File/Comments。
●计算器区域使用计算器,你可以在输入表达式时,加入标准的数学操作符和函数调用,你只需点击序列数字,运算符或者函数等按钮,就可把函数加入表达式中,点击INV按钮,可轮流改变部分按钮的函数功能。
▪按钮“(”与“)”按钮,成对使用圆括号强制改变表达式中的运算顺序。
▪MAX/MIN按钮:查找变量中最大值/变量中最小值。
▪COMPLEX/CONJUGATE按钮:形成一个复变量/对一个复数变量执行共轭运算,利用INV按钮进行函数功能切换。
▪LN/e^X按钮:求一个变量的自然对数/求变量的e次幂,利用INV按钮进行函数功能切换。
▪STO/RCL按钮:将表达式区域信息存储在内存中/从内存中恢复重复使用的表达式,利用INV 按钮进行功能切换。
▪CVAR按钮:计算两个变量之间的协方差(covariance),只适用于PSD求解。
用ANSYS计算拱形橡胶减震垫的极限荷载
第23卷第3期海南大学学报自然科学版V ol.23N o.3 2005年9月NATURA L SCIENCE JOURNA L OF H AINAN UNIVERSIT Y Sep.2005 文章编号:1004-1729(2005)03-0238-04用ANSYS计算拱形橡胶减震垫的极限荷载李正中1,陈玉骥2(1.怀化铁路总公司张家界工务段,湖南张家界427001;2.中南大学铁道校区土建学院,湖南长沙410075)摘 要:通过实验确定了橡胶材料的M ooney2Rivlin常数,再以拱形橡胶减震垫为对象,建立了用于计算的空间有限元模型.考虑材料非线性和几何非线性,用工程软件ANSY S得到了拱形橡胶减震垫的极限荷载以及与之相对应的拱形橡胶减震垫的压缩量.关键词:拱形橡胶减震垫;超弹性材料;极限荷载中图分类号:O344.5 文献标识码:A在铁路线路上,钢轨和混凝土轨枕之间设置了减震垫,目的在于减少火车在运行时对混凝土轨枕面的冲击,改善行车的舒适度.目前,主要应用的减震垫有板式橡胶垫和圆柱式橡胶垫,其各有优缺点.板式橡胶减震垫的承载力较大,但是减震效果欠佳;而圆柱式橡胶减震垫,因为采用空心的构造,具有较好的减震效果,但由于构造上的原因,在圆柱与上部块体的结合处易产生应力集中现象,因而极限荷载偏低,疲劳寿命也较短.拱形橡胶减震垫是一种新型减震产品(见图1),与板式橡胶减震垫和圆柱式橡胶减震垫相比,拱形橡胶减震垫具有独特的优点.由于拱形橡胶减震垫采用半球形空腔的构造,减少了局部应力集中所带来的危害,从而提高了极限荷载和使用寿命,而且空心结构在工作时所产生的空气压力可以起到气体弹簧的作用,有助于减少振动的传递和冲击.可见,拱形橡胶减震垫是一种较好的减震装置.对于这种新型减震垫的力学性能研究,还未见公开报道.为了考察拱形橡胶减震垫的强度储备和变形情况,本文采用工程软件ANSY S对其极限荷载进行了分析,以下分4个方面进行讨论.1 拱形橡胶减震垫有限元模型的建立1.1 有关橡胶材料的超弹性理论 拱形橡胶垫的橡胶是主要的受力材料.橡胶不同于一般的弹性材料,它属于超弹性材料,即它的应力应变关系不是线性规律,也不同于应力应变关系曲线的切线斜率逐渐下降的塑性材料.而是随着应力的增长,应变的增长速度逐渐减小的超弹性材料.橡胶具有超弹性,在较小的外力作用下就能显示出高度的变形能力,由此决定了橡胶在变形过程中会呈现较强的几何非线性和材料非线性[1].对于橡胶超弹性性质的描述有不同的方式[2~4],其中,有从统计理论入手分析的,也有从分子之间的连接力入手分析的.随着计算机技术的发展以及有限元理论的广泛应用,一种被称之为唯象理论的概念逐渐引入了有关橡胶材料的理论计算中,此理论的基本思想是直接从橡胶材料的大变形现象出发,不考虑具体的分子结构,推导出可以完全或比较近似表达橡胶材料的应力应变关系的一种数学格式,从而同经典的弹性理论一样,不用考虑分子和结构等细节.这种方收稿日期:2004-11-26作者简介:李正中(1967-),男,四川广元人,怀化铁路总公司张家界工务段助理工程师.法由M ooney (1940)创始,而后Rivlin 将其推广到了一种更通用更合理的形式.图1 拱形橡胶减震垫的尺寸(单位:mm )Rivlin (1956)假设橡胶材料是不可压缩的,在小应变状态下是各向同性的,而各向同性不可压缩材料的应变能函数W 可以用下式表达W =∑∞i =0,j =0C ij (I 1-3)i (I 2-3)j ,(1)式中I 1和I 2为应变不变量.(1)式中每一项都是(I 1-3)和(I 2-3)的幂函数,此处之所以选用(I 1-3)和(I 2-3)而不用I 1和I 2是为了在零应变时(即I 1=I 3=3),W 将自动消失.1.2 拱形橡胶减震垫的材料常数的确定 M ooney 2Rivlin 常数是一种适用于不可压缩橡胶类材料的材料常数,可使用2、5、9等常数来描述.在有关手册中通常查不到超弹材料的M ooney 2Rivlin常数,但可通过对材料的标准试样的实验所得的实测数据来确定这些常数[5].为得到橡胶材料的物理特性,选取了部分样品做单向压缩实验,测得材料的应力应变关系,并将所测得的数据通过工程软件ANSY S 自带的超弹性分析程序进行材料特性的拟合计算[2],得出应变能函数所需要的常数,从而确定所用橡胶材料的材料属性.按上述办法得到的由M ooney 2Rivlin 模型所描述的应力应变关系如图2所示. 图2 在ANSY S 中定义的应力应变曲线 图3 1Π4拱形橡胶减震垫有限元模型1.3 拱形橡胶垫有限元建模 由于拱形橡胶减震垫具有对称性,故采用1Π4拱形橡胶减震垫建立有限元模型.建模时,采用H yper158单元,材料属性由以上实验所得的数据以及在ANSY S 中计算所得的M ooney 2Rivlin 超弹性模型确定.拱形橡胶减震垫的网格划分有限元模型如图3所示.实际的橡胶减震垫是支撑在混凝土轨枕上的,但为了简化计算,不考虑橡胶减震垫与轨枕的接触摩932 第3期 李正中等:用ANSY S 计算拱形橡胶减震垫的极限荷载擦,故在模型底面只约束竖向位移,而在对称面上则约束在对称面法线方向上的位移.2 荷载施加形式的确定作用于拱形橡胶减震垫的荷载有顶面铁轨所传递的列车轴重荷载以及半球内部的气体压力荷载.轴重荷载采用标准特种荷载,即火车的轴重为25t ,每个拱形橡胶减震垫上的荷载为25Π2=12.5t.半球内部的气体压力荷载是由顶面荷载的作用使半球体积发生变化所引起,它与半球体积的变化有关,故事先无法知道该气体压力的值.为了得到半球内气体的荷载,需进行必要的试算.为此,先做如下假设:1)在施加荷载的过程中,半球内气体不会漏掉;2)半球内的气体变化为等温变化;3)气体的压力随着轴重的增加而线性增加;4)半球内气体的初始压力为标准大气压.试算时,首先忽略半球内气体的存在,对橡胶减震垫施加顶面的列车轴重荷载,加载采用分级加载方式,每级加载10%(即1.25t ),求解过程考虑橡胶材料的材料非线性和由于大变形引起的几何非线性.然后根据橡胶减震垫在12.5t 荷载作用下半球的变形情况,计算半球在变形后的体积,通过与变形前的体积做比较,得到此时的半球内气体的压力.计算表明,所有半球的变化形状都很接近,几乎都是椭球形,故为了简化计算,假设所有半球变化后的尺寸都一样,并取平均尺寸进行计算.根据假设2计算出在顶面荷载12.5t 作用下,有如下关系:P 0V 0=P 1V 1,(2)式中P 0为标准大气压,取为1.01×105Pa ,P 1为变形后的半球内的气体压力,V 0、V 1分别为半球变形之前、后的体积.由于半球在变形后成为一个椭球形的空腔,所以V 0=43πR 3,V 1=43πab 2.(3)式中R 为半球在变形之前的半径;a ,b 分别为半球在变形之后椭球长半轴和短半轴的长度.由(2)、(3)式得P 1=V 0V 1P 0=R 3ab 2P 0,(4)通过上述公式,得到受压变形后半球内的气体压力P 1.在拱形橡胶减震垫未受到荷载时,整个橡胶减震垫边界(包括半球内)都受到一个标准大气压作用,该标准大气压为一平衡荷载,故这时的气压不计入荷载,即以后加载时所有气体压力荷载均应以标准大气压为起点开始计算,所以在变形后半球内的气体压力应该用以下气体压力作为荷载P =P 1-P 0.(5)由于该气体压力的数值是在没有考虑气体压力的情况下得到的,与实际情况的气体压力有所不同.因此,应该将气体压力按假定3施加于半球表面上,再进行求解计算,求得的结果再与上一次的计算结果比较,依此循环计算,直到相邻两次计算结果相差较小时方可认为这时的压力为顶面荷载12.5t 荷载时的气体压力.经过反复试算,最终得到在12.5t 荷载时的平均气体压力为0.35个大气压.3 拱形橡胶减震垫的极限荷载为简化计算,先引入以下假设:042海南大学学报自然科学版 2005年 1)计算时,只考虑到材料的超弹性性质,不考虑到其它的因素,如橡胶材料的疲劳以及其工作环境对材料的影响.图4 顶面位移随荷载变化2)半球内气体压力在加载过程中随顶面荷载同比例逐渐加大,到正常使用荷载时,假设气体压力达到最大值,然后随着荷载的增加,认为气体出现了泄漏现象,气体压力逐渐减小,当达到极限荷载时,气体正好全部漏完,即此时半球内无气体压力.采用上一节的有限元模型,考虑材料非线性,以及橡胶材料受力过程中产生的大变形,通过迭代计算,得到拱形橡胶垫的极限荷载.具体计算结果见图4.从以上计算结果可以看出,拱形橡胶减震垫顶面的位移(压缩量)随着荷载的增加而增加,在荷载均匀增加的情况下,压缩量的增加率是先快后慢,表现出超弹性材料的性质.该拱形橡胶减震垫的极限荷载为37.5t ,相对应的顶面位移为6.67m m.4 结 论1)拱形橡胶减震垫的橡胶材料属于超弹性材料,在计算过程中必须考虑材料的非线性特征,同时还要考虑超弹性材料的大变形能力;2)拱形橡胶减震垫顶面的位移随着荷载的增加逐渐增加,但增量则逐渐减小;3)该拱形橡胶减震垫的极限荷载为37.5t ,安全系数37.5Π12.5=3,顶面的极限位移大小为6.67mm ,约为其厚度的一半.参考文献:[1]史守峡.平面应力不可压缩橡胶薄片的非线性有限元分析[J ].哈尔滨工程大学学报,1998,19(3):11-15.[2]特雷劳尔L G R.橡胶弹性物理学[M].王梦蛟,等译.北京:化学工业出版社,1982.[3]弗雷克利P K,佩恩A R.橡胶在工程中应用的理论与实践[M].杜承泽,等译.北京:化学工业出版社,1985.[4]朱敏.橡胶化学与物理[M].北京:化学工业出版社,1984.[5]易太连,翁雪涛,朱石坚.不可压缩橡胶体的静态性能分析[J ].海军工程大学学报,2002,1(14):76-80.C alculation on the U ltim ate Load of the Arch Rubber Shock P adby E ngineering Softw are ANSYSLI Zheng 2zhong 1,CHE N Y u 2ji 2(1.Zhangjiajie Engineering Workgroup ,Huaihua Railway Parent 2company ,Zhangjiajie 427001,China ;2.Institute of Civil Engineering ,Railway Campus of Mid 2South University ,Changsha 410075,China )Abstract :This paper determined the M ooney 2Rivlin constants for rubber materials by test ,and then taking the arch rubber shock pad as an object ,established a spatial finite element m odel used for calculation.When their material non 2linearity and geometrical non 2linearity were taken into account ,the ultimate load and the corresponding com pressive load of the arch rubber shock pad could be obtained by engineering s oft 2ware ANSY S.K ey w ords :arch rubber shock pad ;super 2elastic materials ;ultimate load 142 第3期 李正中等:用ANSY S 计算拱形橡胶减震垫的极限荷载。
ANSYS基础教程—加载&求解
ANSYS基础教程—加载&求解关键字:ANSYS ANSYS常用命令力载荷求解器多重载荷步信息化调查找茬投稿收藏评论好文推荐打印社区分享本文主要讲述五种载荷类型中剩下一种载荷—集中载荷, 比如应力分析中的节点载荷,包括以下内容:集中载荷、节点坐标、求解器、多重载荷步。
概述·迄今为止, 我们已经知道了如何施加以下类型的载荷:–位移(DOF 约束)–压力和对流载荷(表面载荷)–重力(惯性载荷)–“结构”温度(体载荷)这些载荷占了五种载荷类型中的4种。
本文将讲述剩下的一种载荷—集中载荷, 比如应力分析中的节点载荷。
·将就以下问题进行讨论:A. 集中载荷B. 节点坐标C. 求解器D. 多重载荷步A. 力载荷·一个力就是可以在一个节点或关键点处施加的集中载荷(也可以叫“点载荷”)·和力一样,点载荷适合于线状模型,如梁,桁架,弹簧等。
在实体单元或壳单元中, 点载荷往往引起应力奇异,但当您忽略了附近的应力时,它仍然是可接受的。
记住,您可以通过选择来忽略附近施加了点载荷的单元。
·在左下角展示的二维实体单元中,我们注意到在加力位置出现最大应力SMAX (23,854)。
当在力附近的节点和单元不被选中时,SMAX (12,755)就会移到底部角点处,这是由于在该角点处约束引起的另一处应力奇异。
通过不选底部角点附近的节点和单元,您就可以在上孔附近得到预期的应力SMAX (8,098)。
注意,对于轴对称模型:·在全部360°范围内输入力的值。
·同样在全部360°范围内输出力的值(反力)。
·例如, 设想一个半径为r的圆柱形壳体边缘施加有P lb/in 的载荷。
把这个载荷施加在二维轴对称壳体模形上(比如SHELL51单元), 您就要施加一2πrP的力。
·施加一个力需要有以下信息:–节点号(您可以通过施取确定)–力的大小(单位应与您正在使用的单位系统保持一致)–力的方向—FX, FY, 或FZ使用:–Solution > -Loads-Apply > Force/Moment–或命令FK或F·问题:在哪一个坐标系中FX, FY, 和FZ 有说明?B.节点坐标系·所有的力,位移,和其它与方向有关的节点量都可以在节点坐标中说明。
ANSYS 函数
一、ANSYS 中可以使用的数学函数在ANSYS 帮助系统中关于*SET 命令的注释下列出了ANSYS 中可以使用的数学函数。
所有这些数学函数均可以在ANSYS 环境中使用,这些数学函数包括,其中三角函数的单位默认是弧度:ABS(X)——xACOS(X)——)arccos(xASIN(X)——)arcsin(xATAN(X)——)(x arctgATAN2(X,Y)——)(xy arc COS(X) ——)cos(xCOSH(X)——)(x cohEXP(X)——x eGDIS(X,Y) 求以X 为均值,Y 为标准差的高斯分布,在使用蒙地卡罗法研究随机荷载和随机材料参数时,可以用该函数处理计算结果LOG(X)——)ln(xLOG10(X)——)lg(xMOD(X,Y) 求 X/Y 的余数. 如果 Y=0, 函数值为 0NINT(X) 求最近的整数RAND(X,Y) 取随机数,其中X 是下限, Y 是上限SIGN(X,Y) 取 X 的绝对值并赋予Y 的符号. Y>=0, 函数值为|X|, Y<0, 函数值为-|X|,.SIN(X) 正弦SINH(X) 双曲正弦SQRT(X)——xTAN(X)——)(x tg TANH(X) 双曲正切二、ansys 内嵌函数1. distnd( i,j) — I,j 两点的距离2. node(x,y,z)—提取距离位置(x,y,z)最近的节点号3. kp(x,y,z)—提取距离位置(x,y,z)最近的关键点号4.基本函数GDIS(x,y) Random sample of a Gaussian (normal) distribution with mean x and standard deviation y.VALCHR(CPARM) Numerical value of CPARM (if CPARM is non-numeric, returns 0.0).CHRVAL(PARM) Character value of numerical parameter PARM. Number of decimal places depends on magnitude.UPCASE(CPARM) Upper case equivalent of CPARM.LWCASE(CPARM) Lower case equivalent of CPARM.5. kx(i) 表示关键点i 的x坐标值;同理 ky(i);kz(i)6. nx(i)表示节点i 的x坐标值;同理 ny(i);nz(i)7. nsel(k)是节点k在就是1,不在就是0.8. NDNEXT(N) Next higher node number above N in selected set (or zero if none found).9. NELEM(ENUM,NPOS) returns the node number in position NPOS for element ENUM. Node number at position 1,2,... or 20 of elementN, where npos is 1,2,...20.10. UX(N), UY(N), UZ(N) X, Y, or Z structural displacement or vector sum.11.ROTX(N), ROTY(N), ROTZ(N) X, Y, or Z structural rotation or vector sum.12. TEMP(N) Temperature.13. PRES(N) Pressure.14. VX(N), VY(N), VZ(N) X, Y, or Z fluid velocity or vector sum.15. ENKE(N) Turbulent kinetic energy (FLOTRAN).16. ENDS(N) Turbulent energy dissipation (FLOTRAN).17. CENTRX(N), CENTRY(N), and CENTRZ(N) always retrieve the element centroid in global Cartesian coordinates, and are determined from the selected nodes on the elements.18. ELADJ(N,face) Element number adjacent to face 1,2,...6.。
ANSYS中函数边界条件加载
应用实例
一个空心圆壳结构,其侧视图如图所示,半径为 7850mm,壁厚44mm,承受风速60m/s的风激励,底部圆边 为简支约束。在分析中,根据风洞实验可将风激励简化为 加载函数: p 0.441 | sin t | ( 0.5 cos2 cos sin 2 )
在ANSYS的全局球坐标系中,加载函数表示为:
第一种,也是最常用的方法,可通过ANSYS主菜单中的函 数编辑器和加载器来实现,可通过GUI方式来实现整个过程, 操作界面简明易懂,方便查错和调试。 第二种方法是通过表参数的方法来实现。表参数的最大特点 就是提供按行、列和面的下标进行线性插值的功能,可以用 于定义随时间变化的边界条件或者载荷、响应谱曲线、压力 曲线、材料-温度曲线、磁性材料的B-H曲线等。这是一种非 常直接和有效的方法,就方法本质而言,第一种方法中用函 数编辑器定义的函数文件,最终还是要用函数加载器调入并 定义成表参数供用户使用。 第三种方法是使用ANSYS提供的APDL语言,编程实现复杂 边界条件的施加。这种方法适用于中、高级用户,利用 APDL的强大功能可以实现很多GUI操作无法实现的边界条 件的施加。
p 0.441* abs(sin({TIME})) * (0.5 * cos({Z })^2 * cos({ Y }) (sin({Z })^2)
设定风激励的频率OMEGA为2,即W=2.
fini /clear,nostart *afun,deg /prep7 et,1,93 r,1,44 mp,ex,1,2.094e5 mp,prxy,1,0.262 mp,dens,1,7.83e-6 csys,1 k,1,7850,90 k,2,7850,0 k,3,7850,-30 k,100, l,1,2 l,2,3 arotat,1,2,,,,,1,100,360 aesize,all,1000 csys mshape,0 mshkey,1 amesh,all
在ANSYS中施加任意面、方向、变化载荷的方法
在任意面施加任意方向任意变化的压力在某些特殊的应用场合,可能需要在结构件的某个面上施加某个坐标方向的随坐标位置变化的压力载荷,当然,这在一定程度上可以通过ANSYS表面效应单元实现。
如果利用ANSYS的参数化设计语言,也可以非常完美地实现此功能,下面通过一个小例子描述此方法。
!!!在执行如下加载命令之前,请务必用选择命令asel将需要加载的几何面选择出来!!!finish/prep7et,500,shell63press=100e6amesh, allesla, snsla,s,1! 如果载荷的反向是一个特殊坐标系的方向,可在此建立局部坐标系,并将! 所有节点坐标系旋转到局部坐标系下.*get,enmax,elem,,num,maxdofsel,s,fx,fy,fzfcum,add !!!将力的施加方式设置为"累加",而不是缺省的"替代"*do,i,1,enmax*if,esel,eq,1,then*get,ae,elem,i,area !此命令用单元真实面积,如用投影面积,请用下几条命令! *get,ae,elem,i,aproj,x !此命令用单元X投影面积,如用真实面积,请用上一条命令! *get,ae,elem,i,aproj,y !此命令用单元Y投影面积! *get,ae,elem,i,aproj,z !此命令用单元Z投影面积xe=centrx !单元中心X坐标(用于求解压力值)ye=centry !单元中心Y坐标(用于求解压力值)ze=centrz !单元中心Z坐标(用于求解压力值)! 下面输入压力随坐标变化的公式,本例的压力随X和Y坐标线性变化.p_e=(xe-10)*press+(ye-5)*pressf_tot=p_e*aeesel,s,elem,,insle,s,corner*get,nn,node,,countf_n=f_tot/nn*do,j,1,nnf,nelem(i,j),fx,f_n !压力的作用方向为X方向! f,nelem(i,j),fy,f_n !压力的作用方向为Y方向! f,nelem(i,j),fz,f_n !压力的作用方向为Z方向*enddo*endifesla,s*enddoaclear,allfcum,repl !!!将力的施加方式还原为缺省的"替代"dofsel,allallsel说明:本信息在任意面施加任意方向任意变化的压力在某些特殊的应用场合,可能需要在结构件的某个面上施加某个坐标方向的随坐标位置变化的压力载荷,当然,这在一定程度上可以通过ANSYS表面效应单元实现。
ANSYS加载与求解
1、对实体模型划分网格,(Mapped, Global setting为4) 2、选择Main Menu> Solution> Define Loads> Operate >Transfer to FE> Forces命令,弹出图形选取对话框,然后 单击[OK]按钮,
有限元及ANSYS
A 加载-集中载荷
FEM 加载
有问题
网格的修改均使载荷无效,将自
动删除先前的载荷,须在新网格 上重新施加载荷 不便使用图形拾取施加载荷
有限元及ANSYS
A 加载- DOF约束
DOF约束:给某个自由度(DOF)指定一已知数值 (值不一定是零)。
1. 结构分析中约束被指定为位移约束(平移、旋转或对称、反 对称边界条件);标识的方向均在节点坐标系中。 2. 热分析中约束被指定为温度。 3. 位移约束可施加于节点、关键点、线和面上,用来限制对象 某一方向上的自由度。
注意:只有将载荷直接加到 节点上或将载荷转换后,累 加、忽略方式才起作用。
有限元及ANSYS
A 加载-集中载荷
4、缩放 选择Main Menu> Solution> Define Loads> Operate >Scale FE Loads> Forces命令,弹出对话框。
注意:只有将载荷直接加到节点上或将载荷转换后, 比例缩放才起作用。
1000 500 L3 VALI = 500 VALJ = 1000
线的IJ方向显示,可选择Utility menu/plotctrl/symbols,弹出对话框 设置ldir=on 如果需要表面压力呈多箭头显示,可 选择Utility menu/plotctrl/symbols, 弹出对话框设置如下
ansys 内嵌获取函数的用法
内嵌获取函数的用法
对于某些项目,可以用内嵌的获取函数来代替*GET命令。
获取函数返回项目的值并直接用于当前运行之中。
这样就不必先把值赋给参数,然后再在运行中调用该参数,从而可以省去起中间作用的参数。
例如,要计算两个节点的X坐标的平均值,可以采用*GET函数:1.使用下面的命令来把节点1的X坐标值赋给参数L1:
*GET,L1,NODE,1,LOC,X
2. 再使用*GET命令来把节点2的X坐标值赋给参数L2;
3. 计算中间值MID=(L1+L2)/2.。
更简便的方法是使用节点坐标的获取函数NX(N),该函数返回节点N的X坐标值。
这样就可以不用中间参数L1 和L2。
如下所示:
MID=(NX(1)+NX(2))/2
获取函数的参数可以是参数也可以是其它的获取函数。
例如,获取函数NELEM(ENUM,NPOS)返回在单元ENUM上NPOS处的节点编号,则联合函数NX(NELEM(ENUM,NPOS))返回该节点的X坐标值。
如何实现在ANSYS中加载徐变
如何实现在ANSYS 中加载徐变
ANSYS 中是可以实现徐变的加载的,具体就是采用第6个徐变模型3241/13/(1)cr c c c T c t e c εσ+-=+
现在对各个参数赋值,c2=1,c3=0,c4=0
这样公式就简化为1cr c t εσ=,因此我们只要定义1c t =(,)c t τ
在命令流中写出来就可以。
下面是我的一段命令流,在第i 个载荷步,使用第i -1个载荷步的时间为龄期,持荷时间为第i 与第I-1步之间的时间差。
rate,1
TB,1,,,6
TBDA TA,1,((cr(j,1)+cr(j,2)*(t-0.5)**cr(j,3))*(1-exp(cr(j,4)*0.5))+(cr(j,5)+cr(j,6)*(t-0.5)**cr(j,7))*(1-exp(cr(j,8)*0.5)))*1e-12/(t-temarray(j,4)),1,0,0
当然大家对徐变公式可以有不同的理解,如定义龄期是一个定值:如7天,也可以编程实现的,不过这个时候加载的就是徐变增量了。
ANSYS应用技巧(一)
1.将扩散管的入口速度测量结果用作CFX的入口边界条件对一个扩散管的入口速度进行了实验测量,并打算把测量所得到的结果用作CFX 计算的入口边界条件,怎么做?怎么样才能够在CFX 中引入实验数据而非公式做为入口边界条件呢?只需要创建并使用一个边界条件分布文件就可以了,文件的格式非常简单,就是一个有文件头的CSV 文件。
一个简单的办法是用CFX5.7 把你需要的截面变量分布文件输出,并以此为基础创建自己的变量分布文件。
具体要输出入口边界的操作是这样的:1) 点击“文件”->“输出”命令,在类型中选择“BC”变量分布而不是“通用”或者别的;2) 选择你要输出的变量;3) 点击“OK”;然后打开“profilefile.scv”文件,保持原有格式和单位。
以此文件做为模板,你就可以创建自己的“profilefile.scv”文件了。
2.如何将扩散管的入口速度测量结果用作CFX的入口边界条件对一个扩散管的入口速度进行了实验测量,并打算把测量所得到的结果用作CFX 计算的入口边界条件,怎么做?怎么样才能够在CFX 中引入实验数据而非公式做为入口边界条件呢?只需要创建并使用一个边界条件分布文件就可以了,文件的格式非常简单,就是一个有文件头的CSV 文件。
一个简单的办法是用CFX5.7 把你需要的截面变量分布文件输出,并以此为基础创建自己的变量分布文件。
具体要输出入口边界的操作是这样的:1) 点击“文件”->“输出”命令,在类型中选择“BC”变量分布而不是“通用”或者别的;2) 选择你要输出的变量;3) 点击“OK”;然后打开“profilefile.scv”文件,保持原有格式和单位。
以此文件做为模板,你就可以创建自己的“profilefile.scv”文件了。
3.ANSYSFLOTRAN分析的几个处理技巧Flotran是CFD分析模块,与ANSYS的其它模块一起可以方便的进行流固耦合分析。
在流体分析以及流固耦合分析中,常见的几个问题的处理技巧如下:1)小、负主元问题:可以采用修正的惯性松弛因子的方法处理,惯性松弛因子可以设置为1.0。
手把手教你ansys函数加载
ANSYS函数加载详细步骤
——caemaster
20##12月27日星期四最近,很多同学在ANSYS分析中需要使用函数加载,有的载荷是时间的函数,有的载荷是位置坐标的函数,也有的载荷是温度的函数.
今天小哥我〔cae_mastersina 〕手把手教大家使用函数方式加载.
我举例使用的是一个1×1的面,面上进行自定义函数施加压力载荷Pressure,目的是让大家掌握函数加载方法的详细步骤:
打开ANSYS,先输入以下命令
/PREP7
ET,1,185
BLOCK,,1,,0.2,,1
SMRT,3
VSWEEP,ALL
1,打开函数编辑器
2,坐标系的选择和函数的输入〔在读取函数的时候也可以选择坐标系〕
猜猜这是一个什么载荷?〔中间大,四端小的斜坡载荷〕3,保存函数
4,读取函数
选择刚才保存的fp.func
表格名称自己填写,不要和ANSYS
命令使用相同的名字.
5,加载
输入命令〔显示载荷箭头,将载荷传递至FEA模型〕
/PSF,PRES,NORM,2,0,1
SBCT
/REPLOT
总结:
是不是很简单,五步搞定,把自己的函数写入函数编辑器,可以自定义任意的函数.如果载荷是时间的函数或者是温度的函数,方法也是一样的哦,赶快动手试试吧!。
使用函数边界条件加载及应用实例(Word)
使用函数边界条件加载及应用实例源自:《ANSYS9.0经典产品高级分析技术与实例讲解》1.定义材料特性,创建矩形,指定单元并网格划分命令流:finish/clear/prep7rect,1,10,,0.5et,1,55!定义流体特性mp,KXX,1,10mp,dens,1,1mp,C,1,10mp,VISC,1,0.01!定义平板特性mp,KXX,2,10mp,dens,2,10mp,c,2,5mat,2esize,,25amesh,all2.定义函数定义状态控制变量和函数方程3.定义雷诺数Re表达式并存储到M0。
单击sto键再单击M0 如下。
即将该方程表达式存入0号缓存单元中4.存储prandt1数Clear清楚结果文本框中的信息,并输入如下,基本变量为KXX。
储存到M15.定义Regime1的热传导系数方程6.定义Regime2的热传导系数方程7.添加注释file>comments8.file>save输入文件名htc.func 然后OK 9.曲线或者列表显示定义的函数10.关闭函数编辑对话框11.加载函数并为方程变量指定值。
Function>read from file 打开 htc.func 12.定义各变量的值Regime1Regime213.施加荷载和边界条件,完成计算过程。
nsel,s,loc,y,0d,all,temp,25nsel,s,loc,y,0.5SF,all,CONV,%heatcf%,100 !施加函数边界条件finish/solu %求解time,1deltim,0.1outres,all,allallselsolvefinish/post1set,last/psf,conv,hcoe,1,0.e+00,1/replotFinish50SET TIME/FREQ LOAD STEP SUBSTEP CUMULATIVE1 0.91910E-01 1 1 12 0.97141E-01 1 2 23 0.11306 1 3 34 0.13387 1 4 45 0.13727 1 5 56 0.18748 1 6 67 0.19939 1 7 78 0.21332 1 8 89 0.21337 1 9 910 0.21968 1 10 1011 0.23014 1 11 1112 0.23862 1 12 1213 0.23883 1 13 1314 0.23889 1 14 1415 0.24366 1 15 1516 0.26509 1 16 1617 0.27274 1 17 1718 0.28106 1 18 1819 0.28823 1 19 1920 0.32155 1 20 20友情提示:方案范本是经验性极强的领域,本范文无法思考和涵盖全面,供参考!最好找专业人士起草或审核后使用。
ansys加载与求解
三.检查计算结果,观察收敛情况,决定是否修改网格模型
1.节点最大应力检查
2. 修改网格模型:这 里选择将孔周围的单
元进行网格细化。运 用MeshTool的Refine 完成,细化级别=1。
孔周围要细化的单元
细化后的网格
3.重新计算:这时检查孔顶部最大等效应力(编号可能会改变)= 3.39N/mm2
2.施加约束:
3.施加载荷:板右侧边缘上有一个背离平板的20N/mm的均布线载荷,则均布压力=线 载荷除以板厚20mm=1N/mm2。对模型右侧边施加-1的均布表面压力。
Байду номын сангаас
Solution>Define Loads>Apply>Structural>Pressure>On Lines 4. 求解:Solution>Solve>Current LS
简化分析模型
3. 几何边界、载荷、网格模型以及求解过程的有限元模型
边界、载荷、网格模型(映射网格)
有限元模型
4. 求解结果及其分析
(1)一般性分析
查计算结果可知,平板的最右侧中点位移最大,最大位移=0.519E-06mm;
孔顶部或底部的应力最大,最大等效应力=0.2889M/mm2
其变形图及应力云图如下页所示。
(4)Body load(体积载荷):为施加于模型上的体积载荷或者场 载荷。在结构分析中为温度。
(5)Inertia load(惯性载荷):由物体的惯性引起的载荷,如 重力加速度、角加速度。主要在结构分析中使用。
(6)Coupled-field loads(耦合场载荷):为以上载荷的一种 特殊情况,是从一种分析得到的结果作为另一种分析的载荷。
ANSYS函数
ANSYS函数ANSYS 函数1. distnd( i,j) — I,j 两点的距离2. node(x,y,z) —提取距离位置(x,y,z)最近的节点号3. kp(x,y,z) —提取距离位置(x,y,z)最近的关键点号4. 基本函数ABS(x):求x的绝对值SIGN(x,y):反回值大小等于︳X︱,符号同y ,(y=0时也返回正数)EXP(x) :指数函数x eln(xLOG(x) :自然对数)logLOG10(x) :常用对数()x10SQRT(x) :求x的平方根NINT(x) :求最接近x的整数x/的余数(y=0时返回0)MOD(x,y) :求yx~间平均分布RAND(x,y)::随机数发生器,随机数在yGDIS(x,y):正态分布函数,平均值为x,方差为ySIN(x), COS(x), TAN(x) :三角函数,默认情况下x的单位为弧度,可用命令*AFUN转换为角度。
SINH(x), COSH(x), TANH(x) :双曲函数ASIN(x), ACOS(x), A TAN(x):反三角函数,默认情况下返回值单位为弧度ATAN2(y,x) :反双曲函数,默认情况下返回值单位为弧度V ALCHR(CPARM) :字符参数CPARM 数值,(如果CPARM含有非数字字符则返回0) CHRV AL(PARM) :数值参数PARM的字符串,小数位取决于量级。
UPCASE(CPARM):将字符串CPARM转换为大写LWCASE(CPARM) :将字符串CPARM 转换为小写5. kx(i):表示关键点i 的x坐标值;同理ky(i);kz(i)6. nx(i):表示节点i 的x坐标值;同理ny(i);nz(i)7. nsel(k):是节点k在就是1,不在就是0.8. NDNEXT(N)Next higher node number above N in selected set (or zero if none found).9. NELEM(ENUM,NPOS) returns the node number in position NPOS for element ENUM. Node number at position 1,2,... or 20 of elementN, where npos is 1,2,...20.10. UX(N), UY(N), UZ(N) X, Y, or Z structural displacement or vector sum.11. ROTX(N), ROTY(N), ROTZ(N) X, Y, or Z structural rotation or vector sum.12. TEMP(N) Temperature.13. PRES(N) Pressure.14. VX(N), VY(N), VZ(N) X, Y, or Z fluid velocity or vector sum.15. ENKE(N) Turbulent kinetic energy (FLOTRAN).16. ENDS(N) Turbulent energy dissipation (FLOTRAN).17. CENTRX(N), CENTRY(N), and CENTRZ(N)always retrieve the element centroid in global Cartesian coordinates, and are determined from the selected nodes on the elements.18. ELADJ(N,face)Element number adjacent to face 1,2,...6.。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
ANSYS 中使用函数加载的一个简单例子
本文将通过一个具体实例说明在ANSYS 中如何使用函数加载,后续将通过该实例在分析过程中遇到的一个问题提出自己的一点看法。
实例的具体说明:
一个1/4 圆柱,半径30 mm,外半径42 mm,长度100mm,如图1 所示:
所用材料为双线性弹塑性材料,其机械性能为:
弹性模量 E = 201000 Mpa;泊松比μ=0.3
屈服应力σ= 200 Mpa;切线模量Et = 2010使用单元类型solid185 (8 节点六面体单元)。
取整体单元边长4 mm,然后可以直接对该几何模型划分MAP 网格,划分网格结果如图2:
约束条件为:
轴向两个截面为对称边界条件;一个端面约束轴向位移Uz。
载荷条件为:
在外表面施加变化的压力载荷,载荷函数为:
P (y) = 8e7 + 7E7 * (Y/42)
即:
X = 0 ,Y = 42 (最高点) 时,P = 15E7;
X = 42,Y = 0 (最低点)时,P = 8E7。
我们采用函数方式来施加这一压力载荷,首先定义函数:
在Solution 模块中,点击菜单路径:
Solution > Define Loads > Apply > Functions >
Define/Edit
将会弹出一个函数编辑器,可以在其中定义所需的函数。
在函数编辑器中,函数类型选择为Single equation,即单值函数;计算函数值时使用的插值坐标系( (x,y,z) interpreted in CSYS) 选择0,即总体直角坐标系,如图3 所示:
然后,在函数编辑器中间位置的“Result = “ 小窗口中输入要定义的函数表达式,如果表达式中有x, y, z, time 等变量(供定义函数时使用的“自变量”),可以用{X},{Y},{Z},{TIME} 等的形式输入;或者点击下面一个小窗口右边的小箭头,会出现一个下拉列表,列出可以选择的变量,然后从该列表中选择某个自变量,则该自变量会按照上述格式写入函数中,如图5 所示:
接下来最好检查一下函数定义是否正确。
做法是:点击函数编辑器中的Graph 按钮,弹出一个“Plot information” 对话框,其中显示了所定义的函数方程,列出了方程中所有的变量。
然后可以激活要变化的变量(对于本函数,只有一个变量– X,会自动激活),并在下面的输入窗口中输入该变量的最小和最大值,以及在该变量围中的计算
点数,再点击Graph 按钮,则会绘制函数的曲线;或者点击List 按钮,则列出各计算点上的函数值。
对于本函数,点击Graph 后所显示的函数曲线如图7 所示,可以看到它是一个线性函数:
点击List 后列出的函数值:
如果我们不希望外表面的压力函数是一个线性函数,而是一个三角函数,使得压力从下向上以余弦方式变化,可以将当前坐标系转换为总体圆柱坐标系(csys,1),然后压力函数定义为:
P = 8E7 + 7E7 * cos(Y)
在总体圆柱坐标系中,Y 为角度,这里需要确定Y 的单位是角度(Degree) 还是弧度(Radians)。
为了方便,在函数编辑器中,将角度单位设置为Degree。
下面使用函数编辑器来定义这一函数:
然后同样对所定义的函数进行检查。
在函数编辑器中点击GRAPH 按钮,然后在弹出的对话框中输入X-Axis Range (X 轴围)为0 到90,再点击Graph 按钮以显示压力载荷曲线,点击List 按钮列表压力载荷值,结果如图10 所示:
可以看到,压力不再是线性变化的。
如果我们接受了这一函数定义,需要将它保存起来,以便后续定义载荷时使用。
回到函数编辑器,点击上部菜单中的File > Save:
图11 保存函数
然后弹出一个“另存为” 对话框,在其中,输入文件名(这里输入文件名为press_test),(如有需要还可以修改保存文件的目录) 再点击“保存” 按钮,就可以把所定义的函数保存为以 .func 为扩展名的函数文件。
这里将函数保持为文件press.func:。