ANSYS中单位统一的误区分析(转)

ANSYS中单位统一的误区分析(转)
ANSYS中单位统一的误区分析(转)

ANSYS中单位统一的误区分析:

在ANSYS中没有规定单位,需要用户自己去定义自己的单位制,这就会涉及到单位统一的问题。下边的误区可能是多数初学者经常范的:

EXAMPLE:

计算一个圆柱体的固有频率(为分析简便,采用最简单的形状作为例子),其尺寸如下:

圆柱体长:L=1m;

圆柱体半径:R=0.1m;

材料特性:

弹性模量:2.06e11 Pa;

材料密度:7800kg/m^3;

泊松比:0.3

计算结果如下:

***** INDEX OF DATA SETS ON RESULTS FILE *****

SET TIME/FREQ LOAD STEP SUBSTEP CUMULATIVE

1 0.0000 1 1 1

2 0.0000 1 2 2

3 0.0000 1 3 3

4 0.0000 1 4 4

5 0.0000 1 5 5

6 0.29698E-03 1 6 6

7 834.79 1 7 7

8 834.79 1 8 8

9 1593.7 1 9 9

10 2022.4 1 10 10

如果在建模时采用毫米为单位(在解决实际工程问题时,经常需要从其他CAD 软件导入实体模型,而这些模型常常以毫末为单位),则必须修改材料特性参数,已达到单位统一。如果将材料参数修改为如下:

弹性模量:2.06e5 Pa;(理由:Pa=N/m^2,m=10^3mm,所以要将原值缩小10^6倍)

材料密度:7800e-9kg/m^3;(理由:Pa=N/m^2,m=10^3mm,所以要将原值缩小10^9倍)

泊松比:0.3

计算结果:

***** INDEX OF DATA SETS ON RESULTS FILE *****

SET TIME/FREQ LOAD STEP SUBSTEP CUMULATIVE

1 0.0000 1 1 1

2 0.0000 1 2 2

3 0.0000 1 3 3

4 0.0000 1 4 4

5 0.0000 1 5 5

6 0.46152E-05 1 6 6

7 26.398 1 7 7

8 26.398 1 8 8

9 50.397 1 9 9

10 63.954 1 10 10

二次计算为何出现如此大的差别?????

其实高手可以发现我们在单位转换中有一个转换是错误的,对!你猜对了——弹性模量。

弹性模量的单位是Pa=N/m^2,我们在转换的时候只考虑了m——>mm,而忽略了N 这个单位,他其实与m也有关系。不知道你现在是否有点晕。如果晕了,下面该是清醒的时候了:

看下边这个表:

国际单位制中常用的单位

名称长度力时间质量压力(压强)速度加速度密度 Stress 杨氏模量

单位m N s kg Pa

m/s m/s^2 kg/m^3 Pa Pa

量纲m kgm/s^2 s kg kg/ms^2 m/s m/s^2 kg/m^3 kg/ms^2 kg/ms^2

这下看清楚了吗,注意力的量纲——kgm/s^2,哈哈,原来这里也有一个m,好了,知道为什么前边两次计算结果相差那么大了吧。

接下来改总结一下了:

在ANSYS单位统一变换时,要将单位转换到量刚上(也就是将单位用kg——m——s表示)!!!!!!

观察其特性,然后再做相应的变换。

最后看看如何对前边变换做修改:(在以毫米为单位的模型中)

弹性模量输入值:2.06e8 ;(想想看弹性模量的量纲 kg/ms^2只含一个m,看来只能将原来的数值缩小10^3倍啦)

材料密度输入值:7800e-9;

泊松比:0.3

用上边的数再算算吧看看会有什么结果等着你!!!

弹性模量,密度是材料的特性,其大小是不能变的,在上述分析中变的是在不同的单位系统下,输入ANSYS的数据值.

对于kg, mm,

密度为7.8e-6 kg/mm^3 (假设为钢)

杨氏模量为:2.06e8 kg*mm/s^2

我也觉得十楼的对,对MMKS即mm,kg,s单位制来说就是密度为7.8e-6 kg/mm^3 (假设为钢)

杨氏模量为:2.06e8 kg*mm/s^2

而十一楼说的应该是MMNS单位制,

这些单位制在Pro/E讲的相当清楚。

主系统'毫米千克秒(mmKs)'的单位信息

基本数量:

长度 mm

质量 kg

力 mm kg / sec^2

时间 sec

温度 C

重力 9806.65 mm / sec^2

衍生数量:

面积< mm^2

体积 mm^3

速度 mm / sec

加速度 mm / sec^2

角速度 rad / sec

角度加速度 rad / sec^2

频率 1 / sec

密度 kg / mm^3

压力 kg / (mm sec^2)

应力 kg / (mm sec^2)

杨氏模量 kg / (mm sec^2)

主系统'毫米牛顿秒(mmNs)'的单位信息

基本数量:

长度 mm

质量 tonne

力 N

时间 sec

温度 C

重力 9806.65 mm / sec^2

衍生数量:

面积< mm^2

体积 mm^3

速度 mm / sec

加速度 mm / sec^2

角速度 rad / sec

角度加速度 rad / sec^2

频率 1 / sec

密度 tonne / mm^3

压力 N / mm^2

应力 N / mm^2

杨氏模量 N / mm^2

单位制在惯性、动力、热分析中很重要,单位不正确计算就是错误的。

楼主的说法是完全错误的,虽然最后的计算结果是对的。

(楼主原文)

最后看看如何对前边变换做修改:(在以毫米为单位的模型中)

弹性模量:2.06e8 Pa;(想想看弹性模量的量纲kg/ms^2只含一个m,看来只能将原来的数值缩小10^3倍啦)

材料密度:7800e-9kg/m^3;

(1)、上面文字中2.06e8 Pa,7800e-9kg/m^3的说法不严谨。

(2)、最主要的错误是楼主对弹性模量和材料密度的换算值。

长度为m时用SI国际单位制(m·Kg·s),长度为mm时一般用MPA单位制(mm·T·s)基本物理量及其量纲:

质量:m

长度:L

时间:t

温度:T

(首先记住的是基本物理量:质量、长度、时间、温度)

国际单位制中长度m,质量kg,时间s

MPA单位制中长度mm,质量T(吨),时间s

导出物理量量纲:

速度:

v=L/t

加速度:

a=L/t^2

密度:

ρ=m/L^3

力:

f=m*a=m* L/t^2

压力、应力、弹性模量:P=f/L^2=m/(L*t^2)

因此在MPA单位制下

力的单位为:f=m*a=m* L/t^2=吨*毫米/秒^2=10^3Kg*10^-3米/秒^2 =Kg*米/秒^2=牛顿

弹性模量P=f/L^2=牛顿/毫米^2=MPa

密度ρ=m/L^3=吨/毫米^3

重力加速度:毫米/秒^2

楼主用到的材料特性:

弹性模量:2.06e11 Pa

材料密度:7800kg/m^3

泊松比:0.3

在MPA下分别为:

弹性模量:2.06e5 MPa

材料密度:7800kg/m^3=7.8吨/(10^9mm^3)=7.8e-9吨/mm^3

泊松比:0.3

长度为毫米时楼主的换算:

弹性模量:2.06e8 Pa

材料密度:7800e-9kg/m^3

泊松比:0.3

姑且不论楼主后面单位名称混乱,楼主的换算在计算频率时数值是

对的,因为固有频率的量纲式为(E*L/m)^1/2,在MPA单位制下弹性模量和密度分别乘上1000(楼主的单位,只是数值上相等)不影响频率值的大小。但楼主考虑了没有,你的单位制计算静力问题,

计算结果应力的单位是什么?约束反力的单位是什么?质量的单位

是什么?

我一般采用mm-N-s的单位

长度:mm

弹性模量:2.06e^5 N/mm^2(也就是MPa)

密度:7.85x10^-9

重力加速度:9815

ansys错误汇总大全-史上最全

ANSYS分析出现问题 NSYS error message 错误信息汇总 2011-10-19 12:57:12| 分类:ANSYS | 标签:ansys 错误error |举报|字号订阅以前很多的心得全丢了,现在把新遇到的error message及解决方法逐一添加如下: 1\ Too many expressions. 表达式太长,ansys要求一个表达式不要超过6个分段,比如以下不对 A22=y1*z2-y1*z3-y2*z1+y2*z3+y3*z1-y3*z2 有7个段 改为 A22=y1*z2-y1*z3-y2*z1+y2*z3+y3*z1 A22=A22-y3*z2 就行了 2\ No *DO trips needed, enter *ENDDO . 循环执行次数为0,说明下标的变化范围越界,就是形如下面的循环 *do,i,0,-1 .....

..... *enddo 3\ *** NOTE *** CP = 227.688 TIME= 12:30:54 One or more elements have become highly distorted. Excessive distortion of elements is usually a symptom indicating the need for corrective action elsewhere. Try incrementing the load more slowly (increase the number of substeps or decrease the time step size). You may need to improve your mesh to obtain elements with better aspect ratios. Also consider the behavior of materials, contact pairs, and/or constraint equations. If this message appears in the first iteration of first substep, be sure to perform element shape checking. 为什么上面的错误信息用深色底纹标出呢?原因很简单,上面的错误出现在非线性计算中意味着致命错误,说明计算无法收敛,遇到这个错误是非常头疼的,下面重点讨论这个问题的由来和解决办法。 1、错误信息的内容。这段英文的意思是:一个或多个单元出现严重扭曲。单元的过度扭曲通常意味着需要一些改进措施,比如:减缓载荷的施加速度(增加子步数或者减少时间步长),改进网格质量,同时考虑材料、接触和/或约束方程。

ANSYS常见警告信息相关解释

ANSYS常见警告信息相关解释 NO.0001、ESYS is not valid for line element. 原因:是因为我使用LATT的时候,把“--”的那个不小心填成了“1”。经过ANSYS的命令手册里说那是没有用的项目,但是根据我的理解,这些所谓的没有用的项目实际上都是ANSYS在为后续的版本留接口。对于LATT,实际上那个项目可能就是单元坐标系的设置。当我发现原因后,把1改成0——即使用全局直角坐标系,就没有WARNING了。当然,直接空白也没有问题。 NO.0002、使用*TREAD的时候,有的时候明明看文件好好的,可是却出现*TREAD end-of-file in data read.后来仔细检查,发现我TXT的数据文件里,分隔是采用TAB键分隔的。但是在最后一列后面,如果把鼠标点上去,发现数据后面还有一个空格键。于是,我把每个列最后多的空格键删除,,然后发现上面的信息就没有了。 NO.0003、Coefficient ratio exceeds1.0e8-Check results. 这个大概是跟收敛有关,但是我找不到具体的原因。我建立的一个桥梁分析模型,尽管我分析的结果完全符合我的力学概念判断,规律完全符合基本规律,数据也基本符合实际观测,但是却还是不断出现这个警告信息。有人知道这个信息是什么意思,怎么调试能消除吗? NO.0004、*TREAD end-of-file in data read txt中的表格数据不完整! NO.0005、No*CREATE for*END.The*END command is ignored 忘了写*END了吧。 NO.0006、Keypoint1is referenced by only one line.Improperly connected line set for AL command 两条线不共点,尝试nummrg命令 NO.0007、L1is not a recognized PREP7command,abbreviation,or macro.This command will be ignored 还没有进入prep7,先:/prep7 NO.0008、Keypoint2belongs to line4and cannot be moved 同一位置点2已经存在了,尝试对同位置的生成新点换个编号,比如1002 NO.0009、Shape testing revealed that32of the640new or modified elements violate shape warning limits.To review test results,please see the output file or issue the CHECK command. 单元形状奇异,在我的模型中6面体单元的三个边长差距较大,可忽略该错误。 NO.0010、用命令流建模的时候遇到的

ansys模态分析及详细过程

压电变换器的自振频率分析及详细过程 1.模态分析的定义及其应用 模态分析用于确定设计结构或机器部件的振动特性(固有频率和振型),即结构的固有频率和振型,它们是承受动态载荷结构设计中的重要参数。同时,也可以作为其它动力学分析问题的起点,例如瞬态动力学分析、谐响应分析和谱分析,其中模态分析也是进行谱分析或模态叠加法谐响应分析或瞬态动力学分析所必需的前期分析过程。 ANSYS的模态分析可以对有预应力的结构进行模态分析和循环对称结构模态分析。前者有旋转的涡轮叶片等的模态分析,后者则允许在建立一部分循环对称结构的模型来完成对整个结构的模态分析。 ANSYS提供的模态提取方法有:子空间法(subspace)、分块法(block lancets),缩减法(reduced/householder)、动态提取法(power dynamics)、非对称法(unsymmetric),阻尼法(damped), QR阻尼法(QR damped)等,大多数分析都可使用子空间法、分块法、缩减法。 ANSYS的模态分析是线形分析,任何非线性特性,例如塑性、接触单元等,即使被定义了也将被忽略。 2.模态分析操作过程 一个典型的模态分析过程主要包括建模、模态求解、扩展模态以及观察结果四个步骤。 (1).建模 模态分析的建模过程与其他分析类型的建模过程是类似的,主要包括定义单元类型、单元实常数、材料性质、建立几何模型以及划分有限元网格等基本步骤。 (2).施加载荷和求解 包括指定分析类型、指定分析选项、施加约束、设置载荷选项,并进行固有频率的求解等。 指定分析类型,Main Menu- Solution-Analysis Type-New Analysis,选择Modal。 指定分析选项,Main Menu-Solution-Analysis Type-Analysis Options,选择MODOPT(模态提取方法〕,设置模态提取数量MXPAND. 定义主自由度,仅缩减法使用。 施加约束,Main Menu-Solution-Define Loads-Apply-Structural-Displacement。 求解,Main Menu-Solution-Solve-Current LS。 (3).扩展模态 如果要在POSTI中观察结果,必须先扩展模态,即将振型写入结果文件。过程包括重新进入求解器、激话扩展处理及其选项、指定载荷步选项、扩展处理等。 激活扩展处理及其选项,Main Menu-Solution-Load Step Opts-Expansionpass-Single Expand-Expand modes。 指定载荷步选项。 扩展处理,Main Menu-solution-Solve-Current LS。 注意:扩展模态可以如前述办法单独进行,也可以在施加载荷和求解阶段同时进行。本例即采用了后面的方法 (4).查看结果 模态分析的结果包括结构的频率、振型、相对应力和力等

ANSYS框架结构分析

有限元分析大作业报告 一、结构形式及参数 1、结构基本参数 某框架结构如下图所示,为两榀、三跨七层框架。结构由梁板柱组成,梁板柱之间刚结。材料为C35混凝土,弹性模量为3.15e10N/m2,泊松比取0.25,质量密度为2500kg/m3,梁截面为300mm×700 mm,柱截面为500mm×500mm,楼板厚度为120mm。梁和柱采用beam44 单元,板采用shell 63单元。单位采用国际单位制。 二、静力分析及结果 1、荷载详情 荷载包括自重荷载,采用命令acel,0,0,9.8施加;以及垂直板面向下的均布恒荷载0.35 kN/m2和活荷载0.15 kN/m,两者合并后采用命令*do,mm,204,245,1 sfe,mm,2,pres,,500,500,500,500 *end do施加。 2、结构变形:最大变形发生在91号节点,数值为1.573mm,方向竖直向下(-Z方向)。

3、位移云图 4、等效应力云图:最大等效应力发生在78号节点,数值为175064Pa。

5、支座反力(保留两位小数,单位如表中所示) 节点编码FX(kN) FY(kN) FZ(kN) MX(kN﹒m) MY(kN﹒m) MZ(kN﹒m) 1 -3.87 5.33 514.15 -5.19 -3.74 0.00 2 -6.36 0.09 774.5 3 -0.12 -6.13 0.00 3 -6.36 -0.09 774.53 0.12 -6.13 0.00 4 -3.87 -5.33 514.1 5 5.19 -3.74 0.00 5 0.00 8.2 6 693.8 7 -8.00 0.00 0.00 6 0.00 0.06 107.28 -0.08 0.00 0.00 7 0.00 -0.06 107.28 0.08 0.00 0.00 8 0.00 -8.26 693.87 8.00 0.00 0.00 9 3.87 5.33 514.15 -5.19 3.74 0.00 10 6.36 0.09 774.53 -0.12 6.13 0.00 11 6.36 -0.09 774.53 0.12 6.13 0.00 12 3.87 -5.33 514.15 5.19 3.74 0.00 三、模态分析结果 1、各阶振型频率及类型 振型阶次自振频率(Hz)振动形式 1 1.838 2 弯曲振型 2 1.8627 弯曲振型 3 2.2773 扭转振型 4 5.6636 弯曲振型 5 5.7097 弯曲振型

ANSYS错误提示及其含义

1 在Ansys中出现“Shape testing revealed that 450 of the 1500 new or modified elements violate shape warning limits.”,是什么原因造成的呢? 单元网格质量不够好,尽量用规则化网格,或者再较为细密一点。 2 在Ansys中,用Area Fillet对两空间曲面进行倒角时出现以下错误:Area 6 offset could not fully converge to offset distance 10. Maximum error between the two surfaces is 1% of offset distance.请问这是什么错误?怎么解决?其中一个是圆柱接管表面,一个是碟形封头表面。ansys的布尔操作能力比较弱。如果一定要在ansys里面做的话,那么你试试看先对线进行倒角,然后由倒角后的线形成倒角的面。建议最好用UG、PRO/E这类软件生成实体模型然后导入到ansys。 3在Ansys中,出现错误“There are 21 small equation solver pivot t erms。”,是否是在建立接触contact时出现的错误? 不是建立接触对的错误,一般是单元形状质量太差(例如有接近零度的锐角或者接近180度的钝角)造成small equation solver pivot terms 4在Ansys中,出现警告“SOLID45 wedges are recommended only in regions of relatively low stress gradients.”,是什么意思? "这只是一个警告,它告诉你:推荐SOLID45单元只用在应力梯度较低的区域。 它只是告诉你注意这个问题,如果应力梯度较高,则可能计算结果不可信。" 5ansys向adams导的过程中,出现如下问题“There is not enough memory for the Sparse Matrix Solver to proceed.Please shut down other applications that may be running or increase the virtual memory on your system and return ANSYS.Memory currently allocated for the Sparse Matrix Solver=50MB.Memory currently required for the Sparse Matrix Solver to continue=25MB”,是什么原因造成的? 不清楚你ansys导入adams过程中怎么还需要使用Sparse Matrix Solver(稀疏矩阵求解器)。估计是scrachmemery太低了,从ansys product launcher 进入设置内存,total workspace和dataspace的差就是scrachmemery。如:total workspace 1150MB,dataspace200MB,scrachmemery就是1150-200=950MB。 6在Ansys中,出现错误“error:element type 1 is PLANE42,which can't be used with the VMES command, meshing of volume 3 aborted.”,是什么意思? 意思是:单元类型1是PLANE42,不能使用划分体网格的命令VMES,划分体3中止。 改进办法:1修改单元类型为适合体网格的单元类型。如solid,或shell。2不使用VMES 命令,使用AMESH。 7在Ansys中,出现错误“error: key point 10 is referenced by only one line. Improperly connected line set for AL command.”,是什么意思?该怎么解决? 意思是:关键点10只在一条线上。不适合使用AL命令连接线。 AL命令是用线来定义面,而选择两条线可能只有两个关键点,因为关键点10不在线上,而定义面至少有三个点。 改进办法:再选一条线

用ANSYS进行桥梁结构分析

用ANSYS进行桥梁结构分析 谢宝来华龙海 引言:我院现在进行桥梁结构分析主要用桥梁博士和BSACS,这两种软件均以平面杆系为计算内核,多用来解决平面问题。近来偶然接触到ANSYS,发现其结构分析功能强大,现将一些研究心得写出来,并用一个很好的学习例子(空间钢管拱斜拉桥)作为引玉之砖,和同事们共同研究讨论,共同提高我院的桥梁结构分析水平而努力。 【摘要】本文从有限元的一些基本概念出发,重点介绍了有限元软件ANSYS平台的特点、使用方法和利用APDL语言快速进行桥梁的结构分析,最后通过工程实例来更近一步的介绍ANSYS进行结构分析的一般方法,同时进行归纳总结了各种单元类型的适用范围和桥梁结构分析最合适的单元类型。 【关键词】ANSYS有限元APDL结构桥梁工程单元类型 一、基本概念 有限元分析(FEA)是利用数学近似的方法对真实物理系统(几何和载荷工况)进行模拟。还利用简单而又相互作用的元素,即单元,就可以用有限数量的未知量去逼近无限未知量的真实系统。 有限元模型是真实系统理想化的数学抽象。 真实系统有限元模型 自由度(DOFs)用于描述一个物理场的响应特性。

节点和单元 荷载 1、每个单元的特性是通过一些线性方程式来描述的。 2、作为一个整体,单元形成了整体结构的数学模型。 3、信息是通过单元之间的公共节点传递的。 4、节点自由度是随连接该节点单元类型变化的。 单元形函数 1、FEA仅仅求解节点处的DOF值。 2、单元形函数是一种数学函数,规定了从节点DOF值到单元内所有点处DOF值的计算方法。 3、因此,单元形函数提供出一种描述单元内部结果的“形状”。 4、单元形函数描述的是给定单元的一种假定的特性。 5、单元形函数与真实工作特性吻合好坏程度直接影响求解精度。 6、DOF值可以精确或不太精确地等于在节点处的真实解,但单元内的平均值与实际情况吻合得很好。 7、这些平均意义上的典型解是从单元DOFs推导出来的(如,结构应力,热梯度)。 8、如果单元形函数不能精确描述单元内部的DOFs,就不能很好地得到导出数据,因为这些导出数

ANSYS工程分析 基础与观念Chapter04

第4章 ANSYS结构分析的基本观念Basic Concepts for ANSYS Structural Analysis 这一章要介绍关于ANSYS结构分析的基本观念,熟悉这些基本观念有助于让你很快地区分你的工程问题的类别,然后依此选择适当的ANSYS分析工具。在第1节中我们会对分析领域(analysis fields)做一个介绍,如结构分析、热传分析等。第2节则对分析类别(analysis types)作一介绍,如静力分析、模态分析、或是瞬时分析等。第3节解释何谓线性分析,何谓非线性分析。第4节要对结构材料模式(material models)作一个讨论并作有系统的分类。第5节讨论结构材料破坏准则。第6、7节分别举两个实例,一个是结构动力分析,一个是非线性分析来总合前面的讨论。这两个例子再加上第3章介绍过的静力分析例子,这三个例子可以说是用来做为正式介绍ANSYS命令(第5、6、7章)之前的准备工作。最后(第8节)我们以两个简单的练习题做本章的结束。

第4.1节学科领域与元素类型 Disciplines and Element Types 4.1.1 学科领域(Disciplines) 我们之前提过,ANSYS提供了五大学科领域的分析能力:结傋分析、热传分析、流场分析、电场分析、磁场分析(电场分析及磁场分析可统称为电磁场分析),此外ANSYS也提供了偶合场分析(coupled-field analysis)的能力。为了能分析横跨多学科领域的偶合场,ANSYS提供了一些偶合场元素(coupled-field elements),但是这些元素还是无法涵盖所有偶合的可能性(举例来说,ANSYS 并没有流场与结构的偶合场元素)。但是在ANSYS的操作环境下,再加上利用APDL [Ref. 20],理论上可以进行各种偶合场分析(但是计算时间及收敛性常是问题所在)。下一小节将举几个例子来解说偶合场分析的含义,更详细的偶合场分析步骤你必须参阅Ref. 15。 4.1.2 偶合场分析 以下我们举三个例子来说明何谓偶合场分析。 第一个例子是热应力的计算,这是最常会遇到的问题之一。当你进行热应力分析时,通常分成两个阶段:先做热传分析解出温度分布后,再以温度分布作为结构负载来进行结构分析,而解出应力值。在第一个阶段,热边界条件(thermal boundary conditions)是热传分析的负载,我们希望知道在此热边界条件之下,温度是怎么分布的。因为不均匀的温度分布会造成结构的翘曲变形,所以第二个阶段是希望知道在这些温度分布下结构的变形及应力。这是一个很典型的偶合场分析问题,因为结构怎么变形是依温度怎么分布而定,而温度如何分布则与结构如何变形(变形量很大时,几何形状会改变)有关,这种相依的关系就称为偶合(coupling)。严格来说,前述的分析程序(先做热传分析再做结构分析)观念上不是很正确的,较正确的做法应该是热传与结构分析必须同时进行,也就是说温

ANSYS中文翻译官方手册_接触分析

一般的接触分类 (2) ANSYS接触能力 (2) 点─点接触单元 (2) 点─面接触单元 (2) 面─面的接触单元 (3) 执行接触分析 (4) 面─面的接触分析 (4) 接触分析的步骤: (4) 步骤1:建立模型,并划分网格 (4) 步骤二:识别接触对 (4) 步骤三:定义刚性目标面 (5) 步骤4:定义柔性体的接触面 (8) 步骤5:设置实常数和单元关键字 (10) 步骤六: (21) 步骤7:给变形体单元加必要的边界条件 (21) 步骤8:定义求解和载步选项 (22) 第十步:检查结果 (23) 点─面接触分析 (25) 点─面接触分析的步骤 (26) 点-点的接触 (35) 接触分析实例(GUI方法) (38) 非线性静态实例分析(命令流方式) (42) 接触分析 接触问题是一种高度非线性行为,需要较大的计算资源,为了进行实为有效的计算,理解问题的特性和建立合理的模型是很重要的。 接触问题存在两个较大的难点:其一,在你求解问题之前,你不知道接触区域,表面之间是接触或分开是未知的,突然变化的,这随载荷、材料、边界条件和其它因素而定;其二,大多的接触问题需要计算摩擦,有几种摩擦和模型供你挑选,它们都是非线性的,摩擦使问题的收敛性变得困难。

一般的接触分类 接触问题分为两种基本类型:刚体─柔体的接触,半柔体─柔体的接触,在刚体─柔体的接触问题中,接触面的一个或多个被当作刚体,(与它接触的变形体相比,有大得多的刚度),一般情况下,一种软材料和一种硬材料接触时,问题可以被假定为刚体─柔体的接触,许多金属成形问题归为此类接触,另一类,柔体─柔体的接触,是一种更普遍的类型,在这种情况下,两个接触体都是变形体(有近似的刚度)。 ANSYS接触能力 ANSYS支持三种接触方式:点─点,点─面,平面─面,每种接触方式使用的接触单元适用于某类问题。 为了给接触问题建模,首先必须认识到模型中的哪些部分可能会相互接触,如果相互作用的其中之一是一点,模型的对立应组元是一个结点。如果相互作用的其中之一是一个面,模型的对应组元是单元,例如梁单元,壳单元或实体单元,有限元模型通过指定的接触单元来识别可能的接触匹对,接触单元是覆盖在分析模型接触面之上的一层单元,至于ANSTS使用的接触单元和使用它们的过程,下面分类详述。 点─点接触单元 点─点接触单元主要用于模拟点─点的接触行为,为了使用点─点的接触单元,你需要预先知道接触位置,这类接触问题只能适用于接触面之间有较小相对滑动的情况(即使在几何非线性情况下) 如果两个面上的结点一一对应,相对滑动又以忽略不计,两个面挠度(转动)保持小量,那么可以用点─点的接触单元来求解面─面的接触问题,过盈装配问题是一个用点─点的接触单元来模拟面─与的接触问题的典型例子。 点─面接触单元 点─面接触单元主要用于给点─面的接触行为建模,例如两根梁的相互接触。 如果通过一组结点来定义接触面,生成多个单元,那么可以通过点─面的接触单元来模拟面─面的接触问题,面即可以是刚性体也可以是柔性体,这类接触问题的一个典型例子是插头到插座里。

ANSYS错误集锦 李

ansys分析出现问题 NO.0052 some contact elements overlap with the other contact element which can cause over constraint.这是由于在同一实体上,即有绑定接触(MPC)的定义,又有刚性区或远场载荷(MPC)的定义,操作中注意在定义刚性区或远场载荷时 避免选择不必要的DOF自由度,以消除过约束 NO.0053 Shape testing revealed that 450 of the 1500 new or modified elements violate shape warning limits. 是什么原因造成的呢? 单元网格质量不够好 尽量,用规则化网格,或者再较为细密一点 NO.0054在用Area Fillet对两空间曲面进行倒角时出现以下错误:Area 6 offset could not fully converge to offset distance 10. Maximum error between the two surfaces is 1% of offset distance.请问这是什么错误?怎么解决?其 中一个是圆柱接管表面,一个是碟形封头表面。 ansys的布尔操作能力比较弱。如果一定要在ansys里面做的话,那么你试试看先对线进行倒角,然 后由倒角后的线形成倒角的面。 建议最好用UG、PRO/E这类软件生成实体模型然后导入到ansys

NO.0055 There are 21 small equation solver pivot terms.; SOLID45 wedges are recommended only in regions of relatively low stress gradients.第一个问题我自己觉得是在建立contact时出现的错误,但自己还没有 改正过来;第二个也不知道是什么原因。 还有一个:initial penetration 4.44089×10E-6 was detacted between contact element 53928 and target element 53616;也是建立接触是出现 的,也还没有接近。 第一个问题:There are 21 small equation solver pivot terms.;不是建立接触对的错误,一般是单元形状质量太差(例如有i接近零度的锐角或者接近180度的钝角)造成small equation solver pivot terms 第二个问题:SOLID45 wedges are recommended only in regions of relatively low stress gradients.这只是一个警告,它告诉你:推荐SOLID45单元只用在应力梯度较低 的区域。它只是告诉你注意这个问题,如果应力梯度较高,则可能计算结果不 可信。 NO.0056 ansys向adams导的过程中,出现如下问题There is not enough memory for the Sparse Matrix Solver to

ansys面与面接触分析实例

面与面接触实例:插销拨拉问题分析 定义单元类型 Element/add/edit/delete 定义材料属性 Material Props/Material Models Structural/Linear/Elastic/Isotropic 定义材料的摩擦系数 … 建立几何模型 Modeling/Create/Volumes/Block/By Dimensions X1=Y1=0,X2=Y2=2,Z1=,Z2=

Modeling/Create/Volumes/Cylinder/By Dimensions Modeling/Operate/Booleans/Subtract/Volumes 先拾取长方体,再拾取圆柱体。 Modeling/Create/Volumes/Cylinder/By Dimensions 、 划分掠扫网格 Meshing/Size Cntrls/ManualSize/Lines/Picked Lines 拾取插销前端的水平和垂直直线,输入NDIV=3再拾取插座前端的曲线,输入NDIV=4

PlotCtrls/Style/Size and Shape,在Facets/element edge列表中选择2 facets/edge 建立接触单元 : Modeling/Create/Contact pair,弹出Contact Manager对话框,如图所示。 单击最左边的按钮,启动Contact Wizard(接触向导),如图所示。

单击Pick Target,选择目标面。 选择接触面 定义位移约束 施加对称约束,Define Loads/Apply/Structural/Displacement/Symmetric On Areas,选择对称面。 再固定插座的左侧面。 ) 设置求解选项 Analysis Type/Sol’s Control

ANSYS 非线性_结构分析

目录 非线性结构分析的定义 (1) 非线性行为的原因 (1) 非线性分析的重要信息 (3) 非线性分析中使用的命令 (8) 非线性分析步骤综述 (8) 第一步:建模 (9) 第二步:加载且得到解 (9) 第三步:考察结果 (16) 非线性分析例题(GUI方法) (20) 第一步:设置分析标题 (21) 第二步:定义单元类型 (21) 第三步:定义材料性质 (22) 第四步:定义双线性各向同性强化数据表 (22) 第五步:产生矩形 (22) 1

第六步:设置单元尺寸 (23) 第七步:划分网格 (23) 第八步:定义分析类型和选项 (23) 第九步:定义初始速度 (24) 第十步:施加约束 (24) 第十一步:设置载荷步选项 (24) 第十二步:求解 (25) 第十三步:确定柱体的应变 (25) 第十四步:画等值线 (26) 第十五步:用Post26定义变量 (26) 第十六步:计算随时间变化的速度 (26) 非线性分析例题(命令流方法) (27) 非线性结构分析 非线性结构的定义 在日常生活中,会经常遇到结构非线性。例如,无论何时用钉书针钉书,金 2

属钉书钉将永久地弯曲成一个不同的形状。(看图1─1(a))如果你在一个木架上放置重物,随着时间的迁移它将越来越下垂。(看图1─1(b))。当在 汽车或卡车上装货时,它的轮胎和下面路面间接触将随货物重量的啬而变化。(看图1─1(c))如果将上面例子所载荷变形曲线画出来,你将发现它们都显示了非线性结构的基本特征--变化的结构刚性. 图1─1 非线性结构行为的普通例子 3

非线性行为的原因 引起结构非线性的原因很多,它可以被分成三种主要类型: 状态变化(包括接触) 许多普通结构的表现出一种与状态相关的非线性行为,例如,一根只能拉伸的电缆可能是松散的,也可能是绷紧的。轴承套可能是接触的,也可能是不接触的, 冻土可能是冻结的,也可能是融化的。这些系统的刚度由于系统状态的改变在不同的值之间突然变化。状态改变也许和载荷直接有关(如在电缆情况中),也可能由某种外部原因引起(如在冻土中的紊乱热力学条件)。ANSYS程序中单元的激活与杀死选项用来给这种状态的变化建模。 接触是一种很普遍的非线性行为,接触是状态变化非线性类型形中一个特殊而重要的子集。 几何非线性 如果结构经受大变形,它变化的几何形状可能会引起结构的非线性地响应。一个例的垂向刚性)。随着垂向载荷的增加,杆不断弯曲以致于动力臂明显地减少,导致杆端显示出在较高载荷下不断增长的刚性。 4

ansys错误提示汇总

1把体用面分割的时候出现的错误提示: B ool e a n op er a t i on fa i l e d. t r y a d ju st i n g t h e t ol er a n c e v a l u e on t h e BT O L c om m m a n d t o s om e f r a ct i on of t h e m i n i m u m k ey p oi n t di st a n c e. M od e l S i z e(c u r r e n t p r obl em)1.183933e+000,BT O L s et t i n g 1.000 00e-005,m i n m u m K P T di st a n c e 4.308365e-006 先在要分割的地方设置一个工作平面,用布尔运算“d i v i d ed--v ol u m e by w or k i n g pl a n e”进行分割的时候,出现上述错误,主要愿意可能是设置的公差太小, 当时试了几次都么有成功,最后干脆把体重新建立了一个,又画了一个很大的面,终于成功了。 2.一个常见的代表性错误! 原来我的虚拟内存设置为“无分页文件”,现在改为“系统管理”,就不在出现计算内存不够的情况了。 Error! Element type 1 is Solid95,whichcan not be used with the AMES command, meshing of are a 2 aborted. 刚开始学习的人经常出这种错误,这是因为不同单元类型对应不同的划分网格操作。 上面的错误是说单元类型为Solid95(实体类型),不能用AMES命令划分面网格。 3 Meshing of volume 5 has been aborted because of a lack of memory. Closed down othe r processes and/or choose a larger element size, then try the VMESH command again. Min imum additional memory required=853MB(by kitty_zoe) 说你的内存空间不够,可能因为你的计算单元太多,增加mesh尺寸,减少数量或者增加最小内存设定(a nsys10中在customization preferences菜单存储栏可以修改) 你划分的网格太细了,内存不足。建议将模型划分为几个部分,分部分进行划分,可以减少内存使用,试一下! 4.The input volumes do not meet the conditions required for the VGLU operation . No new entities were created . The VOVLAP operation is a possible alternative VGLU是将两个或多个体粘到一块,体之间的交集应该是面,帮助里的说法,This operation is only v alid if the intersections of the input volumes are areas along the boundaries of those volu mes。你粘结glue的体可能有重叠,所以后面提示了一个VOVLAP命令,该命令是将两个或多个体的重叠部分拿出来作为结果 VMESH划分时,精度不同,单元数量差别太大了,如果是自由网格划分,那么尝试几个SMRT等级看看。还有就是单元形状不同,产生的网格质量也差别很大,我前几天才重新划了一次网格,印象很深。shape,

ansys结构分析基本原理

1 应力-应变关系 本文将介绍结构分析中材料线性理论,在后续再介绍材料非线性的理论。在线弹性理论中应力-应变关系: (1) 其中: {σ}:应力分量,即在ANSYS软件里以S代替σ形式出现。 [D]:弹性矩阵或弹性刚度矩阵或应力-应变矩阵。利用(14)~(19)给出了其具体表达式。(4)给出了其逆矩阵的表达式。通过给出完整的[D]可以定义少数的各向异性单元。在ANSYS中利用命令:TB,ANEL来输入具体数值。 :弹性应变矢量。在ANSY中以EPEL形式输出。 {ε}:总的应变矢量,即 {εth}:热应变矢量,(3)给出了其定义式,在ANSYS中以EPTH形式给出。 注意: {εel}:是由应力引起的应变。 软件中的剪切应变( εxy、εyz和εxz)是工程应变,他们是拉伸应变的两倍。ε通常用来表示拉伸应变,但为了简化输出而采用此表示。将在材料的非线性分析中说明总应变的分量,以EPTO形式输出。 图1 单元的应力矢量图 如图1给出了单元应力矢量图。ANSYS程序中规定正应力和正应变拉伸是为正,压缩时为负。 (1)式还可以被写作以下形式:

(2) 三维情况下,热应变矢量为: (3) 其中: :方向的正割热膨胀系数。 ΔT=T-T ref T:问题中节点当前温度。 :参考温度也就是应变自由时的温度。用TREF或MP命令输入。 T ref 柔度矩阵的定义: (4) 其中: E x: 方向上的杨氏模量,在MP命令中用EX输入。 v xy:主泊松比,在MP命令中用PRXY输入。 :次泊松比,在MP命令中用NUXY输入。 v yx G : 平面上的剪切模量,在MP命令中用GXY输入。 xy 此外,[D]-1是对称矩阵,因此 (5)

ANSYS 提示错误解决办法

ESYS is not valid for line element. 原因:是因为我使用LATT的时候,把“--”的那个不小心填成了“1”。经过ANSYS的命令手册里说那是没有用的项目,但是根据我的理解,这些所谓的没有用的项目实际上都是ANSYS 在为后续的版本留接口。对于LATT,实际上那个项目可能就是单元坐标系的设置。当我发现原因后,把1改成0——即使用全局直角坐标系,就没有W ARNING了。当然,直接空白也没有问题。 NO.0002 使用*TREAD的时候,有的时候明明看文件好好的,可是却出现*TREAD end-of-file in data read. 后来仔细检查,发现我TXT的数据文件里,分隔是采用TAB键分隔的。但是在最后一列后面,如果把鼠标点上去,发现数据后面还有一个空格键。于是,我把每个列最后多的空格键删除,然后发现上面的信息就没有了。 NO.0003 Coefficient ratio exceeds 1.0e8 - Check results. 这个大概是跟收敛有关,但是我找不到具体的原因。我建立的一个桥梁分析模型,尽管我分析的结果完全符合我的力学概念判断,规律完全符合基本规律,数据也基本符合实际观测,但是却还是不断出现这个警告信息。 NO.0004 *TREAD end-of-file in data read txt中的表格数据不完整! NO.0005 No *CREATE for *END. The *END command is ignored 忘了写*END了吧,呵呵 NO.0006 Keypoint 1 is referenced by only one line. Improperly connected line set for AL command 两条线不共点,尝试nummrg命令。 NO.0007 L1 is not a recognized PREP7 command, abbreviation, or macro. This command will be ignored 还没有进入prep7,先:/prep7 NO.0008 Keypoint 2 belongs to line 4 and cannot be moved 关键点2属于线4,移动低级体素时先移动高级体素! NO.0009 Shape testing revealed that 32 of the 640 new or modified elements violate shape warning limits. To review test results, please see the output file or issue the CHECK command. 单元形状奇异,在我的模型中6面体单元的三个边长差距较大,可忽略该错误 NO.0010 用命令流建模的时候遇到的 The drag direction (from the keypoint on drag line 27 that is closest to a keypoint KP of the given area 95) is orthogonal to the area normal at that KP. Area cannot be dragged by the VDRAG command. 意思是拉伸源面的法向与拉伸路径垂直,不能使用VDRAG命令

ANSYS 中使用接触向导定义多个接触对详细实例(图文)

ANSYS 中如何使用接触向导定义接触对 在ANSYS 中定义接触通常有两种方法: 1. 用户自己手工创建接触单元和目标单元。这种方法,在定义接触和目标单元时还比较简单,但是在设置或修改单元属性和定义实常数时却比较复杂。需要用户对接触有较深刻的理解和通过实践积累丰富的经验。 2. 使用接触管理器中的接触向导定义接触对:使用接触管理器 (接触向导) 定义接触对(即接触单元和目标单元) 时,可以定义除了点-点接触以外的各种接触类型;它可以自动生成接触单元和目标单元,并提供了一组默认的单元属性和实常数值。使用这些默认的设置,加上适当的求解设置,对于多数接触问题都能够获得收敛的结果。而且,如果使用默认设置时,计算不收敛或对结果不太满意,也可以通过接触管理器(接触向导) 对单元属性和实常数方便的进行修改和调整。 因此,我们推荐,在可能的情况下,尽量使用接触管理器(接触向导) 来定义接触。本文将通过一个实例介绍接触管理器的基本使用方法。 所使用的例子如下: 两块平板,中间夹一个圆球。上面平板的上表面承受压力,分析模型的变形和应力随压力的变化。 两块平板,尺寸都是(100*100*20),相距100。中间夹一个半径50 的圆球。两个平板分别与圆球的上下边缘接触。尺寸单位为mm。几何模型如图1。

图 1 中,为了能够划分映射网格,分别对体积进行了切割材料属性为:两块平板: E = 201000 Mpa;μ= 0.3 圆球: E = 70100 Mpa;μ= 0.33 接下来对各个Volumes 划分网格,单元类型采用solid186 (20 节点六面体),单元边长统一取 6 mm。网格划分结果如图 2 所示:

ANSYS 常见错误

1、把体用面分割的时候出现的错误提示: Boolean operation failed. try adjusting the tolerance value on the BTOL commmand to some fraction of the minimum keypoint distance. Model Size (current problem)1.183933e+000,BTOL setting 1.00000e-005,minmum KP T distance 4.308365e-006 先在要分割的地方设置一个工作平面,用布尔运算“divided --volume by working plane”进行分割的时候,出现上述错误,主要原因可能是设置的公差太小, 当时试了几次都么有成功,最后干脆把体重新建立了一个,又画了一个很大的面,终于成功了。 2、一个常见的代表性错误! 原来我的虚拟内存设置为“无分页文件”,现在改为“系统管理”,就不在出现计算内存不够的情况了。 Error! Element type 1 is Solid95,which can not be used with the AMES command, meshi ng of area 2 aborted. 刚开始学习的人经常出这种错误,这是因为不同单元类型对应不同的划分网格操作。 上面的错误是说单元类型为Solid95(实体类型),不能用AMES命令划分面网格。 3、Meshing of volume 5 has been aborted because of a lack of memory. Closed d own other processes and/or choose a larger element size, then try the VMESH co mmand again. Minimum additional memory required=853MB(by kitty_zoe ) 说你的内存空间不够,可能因为你的计算单元太多,增加mesh尺寸,减少数量或者增加最小内存设定(ansys10中在customization preferences菜单存储栏可以修改) 你划分的网格太细了,内存不足。建议将模型划分为几个部分,分部分进行划分,可以减少内存使用,试一下! 4、The input volumes do not meet the conditions required for the VGLU operation . No new entities were created . The VOVLAP operation is a possible alternative VGLU 是将两个或多个体粘到一块,体之间的交集应该是面,帮助里的说法,This operati on is only valid if the intersections of the input volumes are areas along the bound aries of those volumes。你粘结glue的体可能有重叠,所以后面提示了一个VOVLAP命令,该命令是将两个或多个体的重叠部分拿出来作为结果 VMESH划分时,精度不同,单元数量差别太大了,如果是自由网格划分,那么尝试几个S MRT等级看看。还有就是单元形状不同,产生的网格质量也差别很大,我前几天才重新划了一次网格,印象很深。shape,0,3d和shape,1,3d就是划分体时控制单元形状的

相关文档
最新文档