非线性基础知识以及ansys例题讲解
弹塑性分析
在这一册中,我们将详细地介绍由于塑性变性引起的非线性问题--弹塑性分析,我们的介绍人为以下几个方面:
?什么是塑性
?塑性理论简介
?ANSYS程序中所用的性选项
?怎样使用塑性
?塑性分析练习题
什么是塑性
塑性是一种在某种给定载荷下,材料产生永久变形的材料特性,对大多的工程材料来说,当其应力低于比例极限时,应力一应变关系是线性的。另外,大多数材料在其应力低于屈服点时，表现为弹性行为，也就是说，当移走载荷时，其应变也完全消失。
由于屈服点和比例极限相差很小，因此在ANSYS程序中，假定它们相同。在应力一应变的曲线中，低于屈服点的叫作弹性部分，超过屈服点的叫作塑性部分，也叫作应变强化部分。塑性分析中考虑了塑性区域的材料特性。
路径相关性：
即然塑性是不可恢复的，那么这种问题的就与加载历史有关，这类非线性问题叫作与路径相关的或非保守的非线性。
路径相关性是指对一种给定的边界条件，可能有多个正确的解—内部的应力，应变分布—存在，为了得到真正正确的结果，我们必须按照系统真正经历的加载过程加载。
率相关性：
塑性应变的大小可能是加载速度快慢的函数，如果塑性应变的大小与时间有关，这种塑性叫作率无关性塑性，相反，与应变率有关的性叫作率相关的塑性。
大多的材料都有某种程度上的率相关性，但在大多数静力分析所经历的应变率范围，两者的应力－应变曲线差别不大，所以在一般的分析中，我们变为是与率无关的。
工程应力，应变与真实的应力、应变：
塑性材料的数据一般以拉伸的应力—应变曲线形式给出。材料数据可能是工程应力
）。（P A0）与工程应变（?l l0），也可能是真实应力（P/A）与真实应变（n L l l()
0大应变的塑性分析一般采用真实的应力，应变数据而小应变分析一般采用工程的应力、应变数据。
什么时候激活塑性：
当材料中的应力超过屈服点时，塑性被激活（也就是说，有塑性应变发生）。而屈服应力本身可能是下列某个参数的函数。
?温度
?应变率
?以前的应变历史
?侧限压力
?其它参数
塑性理论介绍
在这一章中，我们将依次介绍塑性的三个主要方面：
?屈服准则
?流动准则
?强化准则
屈服准则：
对单向受拉试件，我们可以通过简单的比较轴向应力与材料的屈服应力来决定是否有塑性变形发生，然而，对于一般的应力状态，是否到达屈服点并不是明显的。
屈服准则是一个可以用来与单轴测试的屈服应力相比较的应力状态的标量表示。因此，知道了应力状态和屈服准则，程序就能确定是否有塑性应变产生。
屈服准则的值有时候也叫作等效应力，一个通用的屈服准则是Von Mises 屈服准则，当等效应力超过材料的屈服应力时，将会发生塑性变形。
可以在主应力空间中画出Mises 屈服准则，见 图3－1。
在3－D 中，屈服面是一个以
123σσσ==为轴的圆柱面，在2－D 中，屈服面是一个椭圆，在屈服面内部的任何应力状态，都是弹性的，屈服面外部的任何应力状态都会引起屈服。注意：静水压应力状态（123
σσσ==）不会导致屈服：屈服与静水压应力无关，而只与偏差应力有关，因此，1180σ=，230σσ==的应力状态比123
180σσσ==的应力状态接近屈 服。Mises 屈服准则是一种除了土壤和脆性材料外典型使用的屈服准则，在土壤和脆性材料中，屈服应力是与静水压应力（侧限压力）有关的，侧限压力越高，发生屈服所需要的剪应力越大。
流动准则：
流动准则描述了发生屈服时，塑性应变的方向，也就是说，流动准则定义了单个塑性应变分量（x pl
ε，y pl
ε 等）随着屈服是怎样发展的。
一般来说，流动方程是塑性应变在垂直于屈服面的方向发展的屈服准则中推导出来的。这种流动准则叫作相关流动准则，如果不用其它的流动准 则（从其它不同的函数推导出来）。则叫作不相关的流动准则。
强化准则：
强化准则描述了初始屈服准则随着塑性应变的增加是怎样发展的。
一般来说，屈服面的变化是以前应变历史的函数，在ANSYS 程序中，使用了两种强化准则。
等向强化是指屈服面以材料中所作塑性功的大小为基础在尺寸上扩张。对M ises 屈服准则来说，屈服面在所有方向均匀扩张。见图3-2。
图3-2 等向强化时的屈服面变化图
由于等向强化，在受压方向的屈服应力等于受拉过程中所达到的最高应力。
随动强化假定屈服面的大小保持不变而仅在屈服的方向上移 动，当某个方向的屈服应力升高时，其相反方向的屈服应力应该降低。见图3-3。
图3-3 随动强化时的屈服面变化图
在随动强化中，由于拉伸方向屈服应力的增加导致压缩方向屈服应力的降低，所以在对应的两个屈服应力之间总存一个2
y σ的差值，初始各向同性的材料在屈服后将不再是向同性的。
塑性选项
ANSYS 程序提供了多种塑性材料选项，在此主要介绍四种典型的材料选项可以通过激活一个数据表来选择这些选项。
? 经典双线性随动强化 BKIN
? 双线性等向强化 BISO
? 多线性随动强化 MKIN
? 多线性等向强化 MISO
经典的双线性随动强化（BKIN ）使用一个双线性来表示应力应变曲线，所以有两个斜率，弹性斜率和塑性斜率，由于随动强化的Vonmises 屈服准 则被使用，所以包含有鲍辛格效应，此选项适用于遵守Von Mises 屈服准则，初始为各向同性材料的小应变问题，这包括大多数的金属。
需要输入的常数是屈服应力
y σ和切向斜率T E ，可以定义高达六条不同温度下的曲线。
注意：
?使用MP命令来定义弹性模量
?弹性模量也可以是与温度相关的
?切向斜率Et不可以是负数，也不能大于弹性模量
在使用经典的双线性随动强化时，可以分下面三步来定义材料特性。
1、定义弹性模量
2、激活双线性随动强化选项
3、使用数据表来定义非线性特性
双线性等向强化（BIS0），也是使用双线性来表示应力－应变曲线，在此选项中，等向强化的Von Mises 屈服准则被使用，这个选项一般用于初始各向同性材料的大应变问题。需要输入的常数与BKIN选项相同。
多线性随动强化（MKIN）使用多线性来表示应力－应变曲线，模拟随动强化效应，这个选项使用Von Mises 屈服准则，对使用双线性选项（BKIN）不能足够表示应力－应变曲线的小应变分析是有用的。
需要的输入包括最多五个应力－应变数据点（用数据表输入），可以定义五条不同温度下的曲线。
在使用多线性随动强化时，可以使用与BKIN相同的步骤来定义材料特性，所不同的是在数据表中输入的常数不同，下面是一个用命令流定义多线性随动强化的标准输入。
MPTEMP，，10，70
MPDATA，EX，3，，30ES，25ES
TB，MK2N，3
TBTEMP，，STRA2N
TBDA TA，，0.01，0.05，0.1
TBTEMP，10
TBDA TA，，30000，37000，38000
TBTEMP，70
TBDA TA，，225000，31000，33000
多线性等向强化（MISO）使用多线性来表示使用Von Mises屈服准则的等向强化的应力－应变曲线，它适用于比例加载的情况和大应变分析。
需要输入最多100个应力－应变曲线，最多可以定义20条不同温度下的曲线。
其材料特性的定义步骤如下：
1、定义弹性模量
2、定义MISO数据表
3、为输入的应力－应变数据指定温度值
4、输入应力－应变数据
5、画材料的应力－应变曲线
与MKIN 数据表不同的是，MISO的数据表对不同的温度可以有不同的应变值，因此，每条温度曲线有它自己的输入表。
怎样使用塑性
在这一章中，我们将介绍在程序中怎样使用塑性，重点介绍以下几个方面
?可用的ANSYS 输入
?ANSYS 输出量
?使用塑性的一些原则
?加强收敛性的方法
?查看塑性分析的结果
ANSYS 输入：
当使用T B命令选择塑性选项和输入所需常数时，应该考虑到：
?常数应该是塑性选项所期望的形式，例如，我们总是需要应力和总的应变，而不是应力与塑性应变。
?如果还在进行大应变分析，应力－应变曲线数据应该是真实应力－真实应变。
对双线性选项（BK I N ，BISO ），输入常数
y σ和T E 可以按下述方法来决定，如果材料没有明显的屈服应力y σ，通常以产生0.2%的塑性应变所对应的应力作为屈服应力，而T E 可以通过在分析中所预期的应变范围内来拟合实验曲线得到。
其它有用的载荷步选项：
? 使用的子步数（使用的时间步长），既然塑性是一种与路径相关的非线性，因此需要使用许多载荷增量来加载
? 激活自动时间步长
? 如果在分析所经历的应变范围内，应力－应变曲线是光滑的，使用预测器选项，这能够极大的降低塑性分析中的总体迭代数。
输出量
在塑性分析中，对每个节点都可以输出下列量：
EPPL －塑性应变分量x pl
ε, y pl
ε等等
EPEQ －累加的等效塑性应变
SEPL －根据输入的应力－应变曲线估算出的对于EPEQ 的等效应 力
HPRES －静水压应力
PSV －塑性状态变量
PLWK －单位体积内累加的塑性功
上面所列节点的塑性输出量实际上是离节点最近的那个积分点的值。
如果一个单元的所有积分点都是弹性的（EPEQ ＝0），那么节点的弹性应变和应力从积分点外插得到，如果任一积分点是塑性的（EPEQ>0），那么节点的弹性应变和应力实际上是积分点的值，这是程序的缺省情况，但可 以人为的改变它。
程序使用中的一些基本原则：
下面的这些原则应该有助于可执行一个精确的塑性分析
1、 所需要的塑性材料常数必须能够足以描述所经历的应力或应变范围内的材料特性。
2、 缓慢加载，应该保证在一个时间步内，最大的塑性应变增量小于5%，一 般 来说，如果Fy 是系统刚开始屈服时的载荷，那么在塑性范围内的载荷增量应近似为：
? 0.05*Fy － 对用面力或集中力加载的情况
? Fy － 对用位移加载的情况
3、 当模拟类似梁或壳的几何体时，必须有足够的网格密度，为了能够足够的模拟弯曲反应，在厚度方向必须至少有二个单元。
4、 除非那个区域的单元足够大，应该避免应力奇异，由于建模而导致的应力奇异有：
? 单点加载或单点约束
? 凹角
? 模型之间采用单点连接
? 单点耦合或接触条件
5、 如果模型的大部分区域都保持在弹性区内，那么可以采用下列方法来降低计算时间：
? 在弹性区内仅仅使用线性材料特性（ 不 使 用TB 命 令）
? 在线性部分使用子结构
加强收敛性的方法：
如果不收敛是由于数值计算导致的，可以采用下述方法来加强问题的收敛性：
1、使用小的时间步长
2、 如果自适应下降因子是关闭的，打开它，相反，如果它是打开的 ，且割线刚度正在被连续地使用，那么关闭它。
3、使用线性搜索，特别是当大变形或大应变被激活时
1、预测器选项有助于加速缓慢收敛的问题，但也可能使其它的问题变得不稳定。
5、可以将缺省的牛顿－拉普森选项转换成修正的（MODI ）或初始刚度（INIT ）牛顿－拉普森选项，这两个选项比全牛顿－拉普森选项更稳定（ 需要更的迭代），但这两个选项仅在小挠度和小应变塑性分析中有效。
查看结果
1、感兴趣的输出项（例如应力，变形，支反力等等）对加载历史的响应应该是光滑的，一个不光滑的曲线可能表明使用了太大的时间步长或太粗的网格。
2、每个时间步长内的塑性应变增量应该小于5％，这个值在输出文件中以“Max plastic Strain Step”输出，也可以使用POST26来显示这个值（Main Menu:Time Hist Postpro Define Variables）。
3、塑性应变等值线应该是光滑的，通过任一单元的梯度不应该太大。
4、画出某点的应力—应变图，应力是指输出量SEQV（Mises 等效应力），总应变由累加的塑性应变EPEQ和弹性应变得来。
塑性分析实例（GUI方法）
在这个实例分析中，我们将进行一个圆盘在周期载荷作用下的塑性分析。
问题描述：
一个周边简支的圆盘，在其中心受到一个冲杆的周期作用。由于冲杆被假定是刚性的，因此在建模时不考虑冲杆，而将圆盘上和冲杆接触的结点的Y方向上的位移耦合起来。
由于模型和载荷都是轴对称的，因此用轴对称模型来进行计算。求解通过四个载荷步实现。
问题详细说明：
材料性质：
EX=70000 (杨氏模量）
NUXY=0.325（泊松比）
塑性时的应力—应变关系如下：
应变应力
0.0007857 55
0.00575 112
0.02925 172
0.1 241
加载历史：
时间载荷
0 0
1 -6000
2 750
3 -6000
问题描述图：
步骤一：建立计算所需要的模型。
在这一步中，建立计算分析所需要的模型，包括定义单元类型，划分网格，给定边界条件。并将数据库文件保存为“exercise2.db”。在此，对这一步的过程不作详细叙述。
步骤二：恢复数据库文件“exercise.db”
Utility Menu>File>Resume from
步骤三：定义材料性质
1、选择菜单路径Main Menu>Preprocessor>Matersal Props>-Constant-Isotropic. Isotropic Matersal Properties (各向同性材料性质）对话框出现。
2、单击OK来指定材料号为1。另一个I sotropic Material Properties对话框出现。
3、对杨氏模量（EX）键入EXX 。
4、对泊松比（NUXY）键入0.325。
5、单击OK。
步骤四：定义和填充多线性随动强化数据表（MKIN）
1、选择菜单路径Main Menu>Preprocessor>Material Props>DataTables>Define/Activate. Define/Activate Data Table(激活数据表）对话框出现。
2、在关于type of data table(数据表类型）的卷动框中，卷动到“Multi kinem MKIN”
且选中它。
3、在material refersuce number(材料参考号）中，健入1。
4、对number of temperatures(温度数）键入1，单击OK。
5、选择菜单路径Main Menu>Preprocessor>Material Props>Data Tables>Edit Active.。. Data Table MKIN对话框出现。
6、在“Strain”一行中，从第二列起分别输入STN1,STN2,STN3,STN4。
7、在“Curve 1”一行中，从第二列起分别输入STS1,STS2,STS3,STS4。
8、选择File>Apply & Quit。
9、选择菜单路径Main Menu>Preprosessor>Material Porps>Data Tables>Graph. Graph Data Tables(图形表示数据表）对话框出现。
1、单击OK接受绘制MKIN表的缺省。一个MKIN表的标绘图出现在ANSYS 图形窗口中。
步骤五：进入求解器
选择菜单路径Main Menu>Solution。
步骤六：定义分析类型和选项
1、选择菜单路径Main Menu>Solution>-Analysis Type-New Analysis.
2、单击“Static”来选中它然后单击OK。
步骤七：打开预测器，设置输出控制。
1、选择菜单路径Main menu>solution-Load Set Opts-Nonlinear>Predictor。
2、将predictor的状态设置为“ON”。
3、选择菜单路径Main Menu>Solution>-Load Step Options- Output Ctrls> DB/Results File. Coutrols for Database and Results File Writing (对数据库和结果文件写入的控制）对话框出现。
3、单击“Every substep”且选中它。
步骤八：设置载荷步选项
1、选择菜单路径Main Menu>Solution>-Load Step Options-Time/Frequenc > time&Substep。 Time&Substep Option(时间和子步数选项）对话框出现。
2、对time at end of Load Step(载荷步终止时间）键入 1e-6
3、对Number of substeps (子步数）键入1。
步骤九：对第一个载荷步加载
在结点3的Y方向施加一大小为 0的集中力载荷。
步骤十：将第一个载荷步写入载荷步文件。
1、选择菜单路径Main Menu>Solution>-Write Ls File，出现对话框。
2、在“LSNUM”的输入框中键入 1
步骤十一：对第二个载荷步加载，并写入载荷步文件。
1、选择菜单路径Main Menu>Solution>-Load Step Options-Time/Frequenc>time&Substep。 Time&Substep Option(时间和时间步选项）对话
框出现。
2、对time at end of Load Step(载荷步终止时间）键入1
3、对Number of substeps (子步数）键入10。
4、单击automatic time stepping option（自动时间步长选项）使之为ON，然后单击OK。
5、在结点3的Y方向施加一大小为 -6000的集中力载荷。
6、选择菜单路径Main Menu>Solution>-Write Ls File，出现对话框。
7、在“LSNUM”的输入框中键入 2
步骤十二：对第三个载荷步加载，并写入载荷步文件。
1、选择菜单路径Main Menu>Solution>-Load Step Options-Time/Frequenc>time&Substep。 Time&Substep Option(时间和时间步选项）对话框出现。
2、对time at end of Load Step(载荷步终止时间）键入2
3、在结点3的Y方向施加一大小为 750的集中力载荷。
4、选择菜单路径Main Menu>Solution>-Write Ls File，出现对话框。
5、在“LSNUM”的输入框中键入3
步骤十三：对第四个载荷步加载，并写入载荷步文件。
1、选择菜单路径Main Menu>Solution>-Load Step Options-
Time/Frequenc>time&Substep。 Time&Substep Option(时间和
时间步选项）对话框出现。
2、对time at end of Load Step(载荷步终止时间）键入3
3、在结点3的Y方向施加一大小为 -6000的集中力载荷。
4、选择菜单路径Main Menu>Solution>-Write Ls File，出现对话框。
5、在“LSNUM”的输入框中键入4
步骤十三：求解问题
1、选择菜单路径Main Menu>Solution>-Solve-From Ls Files，对话框出现。
2、对“LSMIN”键入1，对“LSMAX”键入4。
3、单击对话框中的OK开始求解。
步骤十四：进行后处理。
在这一步中，可以进行所想要的后处理，在此不进行详述。
非线性静态实例分析（命令流方式）
你可以用下面显示的ANSYS命令替代GUI选择，进行上面这个例题的塑性分析。
fini
/cle
/title,circular plate loaded by a circular punch - kinematic hardening
rpl=65
rpu=5
h=6.5
exx=70000
sts1=55 !yield stress
stn1=sts1/exx
sts2=112
stn2=0.00575
sts3=172
stn3=0.02925
sts4=241
stn4=0.1
nex=15 ! No. of elements along the radius
net=2 ! No. of elements in the plate's thickness
/prep7
et,1,42,,,1
! define node for convenient postprocessing
n,1,rpl,0
n,2,0,0
n,3,0,h/2
!define geometry
k,1,,-(h/2)
k,2,rpu,-(h/2)
k,3,rpl,-(h/2)
kgen,2,1,3,1,,h,,3
nex1=nint(0.8*net)
nex2=nex-nex1
l,1,2
l,4,5
l,2,3
l,5,6
a,1,2,5,4
a,2,3,6,5
lesize,1,,,nex1
lesize,2,,,nex1
lesize,3,,,nex2,2.5 lesize,4,,,nex2,2.5 esize,(h/net)
amesh,all
nummrg,node
nodes
nsel,s,loc,x,0,rpu
nsel,r,loc,y,(h/2)
cp,1,uy,all
nsel,all
esel,all
fini
/solu
nsel,s,loc,x,rpl
nsel,r,loc,y,0
d,all,uy
nsel,s,loc,x,0
d,all,ux
nsel,all
outres,all,all
fini
/prep7
mp,ex,1,exx
mp,nuxy,1,0.325
tb,mkin,1
tbtemp,,strain tbdata,,stn1,stn2,stn3,stn4 tbtemp,,
tbdata,,sts1,sts2,sts3,sts4 fini
/solu
pred,on
outres,all,all
nsubst,1
time,1e-6
f,3,fy,0
lswrite
ansys非线性分析指南
ANSYS 非线性分析指南(1) 基本过程 第一章结构静力分析 1. 1 结构分析概述 结构分析的定义： 结构分析是有限元分析方法最常用的一个应用领域。结构这个术语是一个广义的概念，它包括土木工程结构，如桥梁和建筑物；汽车结构，如车身、骨架；海洋结构，如船舶结构；航空结构，如飞机机身、机翼等，同时还包括机械零部件，如活塞传动轴等等。 在ANSYS 产品家族中有七种结构分析的类型，结构分析中计算得出的基 本未知量- 节点自由度，是位移；其他的一些未知量，如应变、应力和反力， 可通过节点位移导出。 七种结构分析的类型分别是： a. 静力分析- 用于求解静力载荷作用下结构的位移和应力等。静力分析 包括线性和非线性分析。而非线性分析涉及塑性、应力刚化、大变形、大应变、超弹性、接触面和蠕变，等。 b. 模态分析- 用于计算结构的固有频率和模态。 c. 谐波分析- 用于确定结构在随时间正弦变化的载荷作用下的响应。 d. 瞬态动力分析- 用于计算结构在随时间任意变化的载荷作用下的响应，并且可计及上述提到的静力分析中所有的非线性性质。 e. 谱分析- 是模态分析的应用拓广，用于计算由于响应谱或PSD 输入 随机振动引起的应力和应变。 f. 屈曲分析- 用于计算屈曲载荷和确定屈曲模态，ANSYS 可进行线性特征值和非线性屈曲分析。 g. 显式动力分析- ANSYS/LS-DYNA可用于计算高度非线性动力学和复 杂的接触问题。 除了前面提到的七种分析类型，还有如下特殊的分析应用： ? 断裂力学 ? 复合材料 ? 疲劳分析
? p-Method 结构分析所用的单元：绝大多数的ANSYS 单元类型可用于结构分析。单元类型从简单的杆单元和梁单元一直到较为复杂的层合壳单元和大应变实体单元 1.2 结构线性静力分析 静力分析的定义： 静力分析计算在固定不变的载荷作用下结构的响应。它不考虑惯性和阻尼的影响，如结构受随时间变化载荷的情况。可是静力分析可以计算那些固定不变的惯性载荷对结构的影响，如重力和离心力；以及那些可以近似为等价静力作用的随时间变化载荷，如通常在许多建筑规范中所定义的等价静力风载和地震载荷。 静力分析中的载荷： 静力分析用于计算由那些不包括惯性和阻尼效应的载荷作用于结构或部件上引起的位移、应力、应变和力。固定不变的载荷和响应是一种假定，即假定载荷和结构的响应随时间的变化非常缓慢，静力分析所施加的载荷包括： ? - 外部施加的作用力和压力 ? - 稳态的惯性力如中力和离心力 ? - 位移载荷 ? - 温度载荷 线性静力分析和非线性静力分析 静力分析既可以是线性的也可以是非线性的。非线性静力分析包括所有的非线性类型：大变形、塑性、蠕变、应力刚化、接触、间隙单元、超弹性单元等，本节主要讨论线性静力分析，非线性静力分析在下一节中介绍。 线性静力分析的求解步骤 1 建模 2 施加载荷和边界条件求解 3 结果评价和分析
ansys经典例题步骤
Project1 梁的有限元建模与变形分析 计算分析模型如图1-1 所示, 习题文件名: beam。 NOTE:要求选择不同形状的截面分别进行计算。 梁承受均布载荷：1.0e5 Pa 图1-1梁的计算分析模型 梁截面分别采用以下三种截面（单位：m）： 矩形截面：圆截面：工字形截面： B=0.1, H=0.15 R=0.1 w1=0.1，w2=0.1，w3=0.2, t1=0.0114，t2=0.0114，t3=0.007 1.1进入ANSYS 程序→ANSYSED 6.1 →Interactive →change the working directory into yours →input Initial jobname: beam→Run 1.2设置计算类型 ANSYS Main Menu: Preferences →select Structural →OK 1.3选择单元类型 ANSYS Main Menu: Preprocessor →Element Type→Add/Edit/Delete… →Add… →select Beam 2 node 188 →OK (back to Element Types window)→Close (the Element Type window) 1.4定义材料参数 ANSYS Main Menu: Preprocessor →Material Props →Material Models →Structural→Linear→Elastic→Isotropic→input EX:2.1e11, PRXY:0.3→OK 1.5定义截面 ANSYS Main Menu: Preprocessor →Sections →Beam →Common Sectns→分别定义矩形截面、圆截面和工字形截面：矩形截面:ID=1,B=0.1，H=0.15 →Apply →圆截面：ID=2,R=0.1 →Apply →工字形截面：ID=3,w1=0.1，w2=0.1，w3=0.2，t1=0.0114，t2=0.0114，t3=0.007→OK
ANSYS例题
1.ANSYS中混凝土模式预应力模拟的算例 !简支梁实体与预应力钢筋分析 /PREP7 egjx=2e5 !Ey agjx=140 !单根钢绞线面积 ehnt=4e4 !Eh xzxs=1.0e-5 !线胀系数 yjl=200000 !定义预加力 et,1,link8 !定义link8单元 et,2,solid95 !定义solid95单元 r,1,agjx !定义link8单元的面积 r,2 !定义第2种实常数 mp,ex,1,egjx !定义link8单元的弹性模量mp,prxy,1,0.3 !定义link8单元的泊松系数mp,alpx,1,1.0e-5 !定义线膨胀系数 mp,ex,2,ehnt !定义solid95单元的弹性模量mp,prxy,2,0.3 !定义solid95单元的泊松系数blc4, , ,100,200,3000 !定义梁体 /view,1,1,1,1 !定义ISO查看 /ang,1 vplot !绘制梁体 kwpave,6 !工作平面移动到关键点6
wpoff,-30 !工作平面移动-30mm(X) wprot,0,0,90 !工作平面旋转 vsbw,1 !分割梁体 wpoff,0,0,-40 !工作平面移动-40mm(Z) vsbw,2 !分割梁体 wpoff,0,40 !工作平面移动40mm(Y) wprot,0,90 !工作平面旋转 vsbw,all !分割梁体 wpstyl !关闭工作平面显示 nummrg,all,,,,low !整理 numcmp,all !压缩编号 esize,30 !定义网分时边长控制 lsel,s,,,28,38,10 !定义line28和38为新的选择集latt,1,1,1 !定义选择集的属性 lmesh,all !对线划分单元 allsel,all !新的选择集为所有的实体 gplot !绘制所有的实体 vsel,s,,,all !定义所有体为选择集 vatt,2,2,2 !定义选择集的属性 mshape,0,3d !将体划分单元的形状定位HEX mshkey,1 !采用MAPPED划分器 vmesh,all !对体进行划分单元
ansys学习-非线性静态分析实例
ansys学习－非线性静态分析实例 问题描述 一个子弹以给定的速度射向壁面。壁面假定是刚性的和无摩擦的。将研究子弹和壁面接触后达80微秒长的现象。目的是确定子弹的整个变形，速度历程，以及最大等效Von Mises应变。求解使用SI单位。 用轴对称单元模拟棒。求解最好能通过单一载荷步实现。在这个载荷步中，将同时施加初始速度和约束。将圆柱体末端的节点Y方向约束住以模拟一固壁面。打开自动时间分步来允许ANSYS 确定时间步长。定义分析结束的时间为8E-5秒，以确保有足够长的时间来扑捉整个变形过程。 问题详细说明 下列材料性质应用于这个问题： EX=117.0E09 (杨氏模量） DENS=8930.0 （密度） NUXY=0.35（泊松比） Yield Strength=400.0OE06（屈服强度） Tangent Modulus （剪切模量） 下列尺寸应用于这个问题： 长=32.4E-3m 直径=6.4E-3m 对于这个问题的初始速度是227.0。 图1铜圆柱体图解 求解步骤： 步骤一：设置分析标题 1、选择菜单路径：Utility Menn>File>ChangeTitle。
2、键入文字“Coppery Cylinder Impacting a Rigid Wall” 3、单击OK。 步骤二：定义单元类型 1、选择菜单路径Mail Menu>Preprocessor>Element Type>All/Edit/Delete。 2、单击Add。Library of Element Types(单元类型库）对话框出现。 3、在靠近左边的列表中，单击“Visio Solid”仅一次。 4、选靠近右边的列表中，单击“4node Plas 106”仅一次。 5、单击OK。Library of Element Types 对话框关闭。 6、单击Options (选项）。VISCO106 element type Options(visco106单元类型选项）对话框出现。 7、在关于element behavior（单元特性）的卷动柜中，卷动到“Axisymmetric” 且选中它。 8、单击OK。 9、单击Element Types (单元类型）对话框中的Close。 步骤三：定义材料性质 1、选择菜单路径Main Menu>Preprocessor>Material Props>-Constant-Isotropic. Isotropic Matersal Properties (各向同性材料性质）对话框出现。 2、单击OK来指定材料号为1。另一个I sotropic Material Properties对话框出现。 3、对杨氏模量（EX）键入117.0E09 4、对密度（DENS）键入8930。 5、对泊松比（NUXY）键入0.35。 6、单击OK。 步骤四：定义双线性各向同性强化数据表（BISO） 1、选择菜单路径Main Menu>Preprocessor>Matersal Props>Data Tables> Define/Activate . Define/Activate Data Table(定义数据表）对话柜出现。 2、在关于type of data table(数据表类型）的卷动框中，卷动到“Bilin isotr BISO”且选中它。 3、对material reference number(材料参考号）健入1。 4、对number of temperatures(温度数）键入1和单击OK。 5、选择菜单路径Main Menu>Preprocessor>Material Props>Data Tables>Edit Active. Data Table BISO对话框出现。 6、对YLD Strs(屈服应力）键入400.0e06。 7、对 Tang Mod(剪切模量）键入100.0e06。 8、选择File>Apply & Quit。 9、选择菜单路径Main Menu>Preprosessor>Material Porps>Data Tables>Graph. Graph Data Tables(图形表示数据表）对话框出现。 10、单击OK接受绘制BISO表的缺省。一个BISO表的标绘图出现在ANSYS图形窗口中。 11、在ANSYS TooLbar上单击SAVE_DB。 步骤五、产生矩形 在这一步中，你产生一个代表柱体半横截面积的矩形。
[整理]《ANSYS120宝典》习题.
第1章 习题 1.ANSYS软件程序包括几大功能模块？分别有什么作用？ 2.如何启动和退出ANSYS程序？ 3.ANSYS程序有哪几种文件类型？ 4.ANSYS结构有限元分析的基本过程是什么？ 5.两杆平面桁架尺寸及角度如习题图1.1所示，杆件材料的弹性模量为2.1×1011Pa，泊松 比为0.3，截面面积为10cm2，所受集中力载荷F=1000N。试采用二维杆单元LINK1计算集中力位置节点的位移和约束节点的约束反力。 习题图1.1 两杆平面桁架 第2章 习题 1.建立有限元模型有几种方法？ 2.ANSYS程序提供了哪几种坐标系供用户选择？ 3.ANSYS程序中如何平移和旋转工作平面？ 4.试分别采用自底向上的建模方法和自顶向下的建模方法建立如习题图2.1所示的平面图 形，其中没有尺寸标注的图形读者可自行假定，并试着采用布尔运算的拉伸操作将平面图形沿法向拉伸为立体图形。
习题图2.1 平面图形 5.试分别利用布尔运算建立如习题图2.2所示的立体图形，其中没有尺寸标注的图形读者 可自行假定。 习题图2.2 立体图形 6.试对习题图2.3所示的图形进行映射网格划分，并任意控制其网格尺寸，图形尺寸读者 可自行假定。 习题图2.3 映射网格划分
第3章 习题 1.试阐述ANSYS载荷类型及其加载方式。 2.试阐述ANSYS主要求解器类型及其适用范围。 3.如何进行多载荷步的创建，并进行求解？ 4.试建立如习题图3.1所示的矩形梁，并按照图形所示施加约束和载荷，矩形梁尺寸及载 荷位置大小读者可自行假定。 习题图3.1 矩形梁约束与载荷 5.试建立如习题图3.2所示的平面图形，并按照图形所示施加约束和载荷，平面图形的尺 寸及载荷大小读者可自行假定。 习题图3.2 平面图形约束与载荷 第4章 习题
(完整版)ansys内部例题详解
郑重申明：本人能力有限，文中不可避免会有错误，欢迎朋友们批评指正，希望大家相互提高，呵呵，谢谢啦！ 12.8. Sample Rigid Body Dynamic Analysis 刚体动力学分析实例 This sample analysis demonstrates how to model a flexible component in ANSYS and export the flexible body information to a file for use in ADAMS. The example also provides brief instructions on how to perform the rigid body dynamic analysis in ADAMS, and details on how to transfer the loads from ADAMS to ANSYS in order to perform a stress analysis. 该实例演示了如果在ANSYS中制作柔性部件及输出可在ADAMS中使用的柔性体信息文件。同样该例子也提供了有关于在ADAMS中进行动力学分析的简单介绍，和如何将载荷信息从ADAMS转换到ANSYS中进行应力分析的详细介绍。 12.8.1. Problem Description 问题描述 In the linkage assembly shown below, Link3 is a flexible component. Link3 is modeled as a rectangular rod in ANSYS using SOLID45elements. The joints in ADAMS will be attached to interface points (nodes) at the middle of the holes at either end of Link3. These middle points are connected to the cylindrical joint surfaces by a spider web of BEAM4 elements. 联动装置装配如下图所示，连杆3是一个柔性部件，为矩形杆件在ANSYS中采用SOLID45单元构造。ADAMS中连接铰将连接在位于杆两端的孔中心接触节点上。这些节点会通过BEAM4单元构造的蜘蛛网格与圆柱铰表面连接。 Figure 12.5: Linkage Assembly联动装置转配图
ANSYS命令流学习笔记10-利用APDL在WorkBench中进行非线性屈曲分析
!ANSYS命令流学习笔记10-利用APDL在WorkBench中进行非线性屈曲分析 !学习重点： !1、强化非线性屈曲知识 首先了解屈曲问题。在理想化情况下,当F < Fcr时, 结构处于稳定平衡状态，若引入一个小的侧向扰动力，然后卸载, 结构将返回到它的初始位置。当F > Fcr时, 结构处于不稳定平衡状态, 任何扰动力将引起坍塌。当F = Fcr时,结构处于中性平衡状态，把这个力定义为临界载荷。在实际结构中, 几何缺陷的存在或力的扰动将决定载荷路径的方向。在实际结构中, 很难达到临界载荷，因为扰动和非线性行为, 低于临界载荷时结构通常变得不稳定。 要理解非线性屈曲分析，首先要了解特征值屈曲。特征值屈曲分析预测一个理想线弹性结构的理论屈曲强度，缺陷和非线性行为阻止大多数实际结构达到理想的弹性屈曲强度，特征值屈曲一般产生非保守解, 使用时应谨慎。 !理论解，根据Euler公式。其中μ取决于固定方式。 !有限元方法， 已知在特征值屈曲问题： 求解，即可得到临界载荷 而非线性屈曲问题： 其中为结构初始刚度，为有缺陷的结构刚度，为位移矩阵，为载荷矩阵。 非线性屈曲分析时考虑结构平衡受扰动（初始缺陷、载荷扰动）的非线性静力分析，该分析时一直加载到结构极限承载状态的全过程分析，分析中可以综合考虑材料塑性、几何非线性、接触、大变形。非线性屈曲比特征值屈曲更精确，因此推荐用于设计或结构的评价。 !2、熟悉WB中非线性屈曲分析流程 (1) 前处理，施加单元载荷，进行预应力静力分析。 (2) 基于预应力静力分析，指定分析类型为特征值屈曲分析，完成特征值屈曲分析。 (3) 在APDL模块将一阶特征屈曲模态位移乘以适当系数，将此变形后的形状当做非线性分析的初始模型。
土木工程专业ANSYS初学者经典例题
巷道开挖过程的有限元模型与力学分析 学院： 班级： 学号： 姓名： 指导老师：XXX老师 XX大学 2009－2010年度上学期
目录 巷道开挖过程的有限元模型与力学分析 (3) 有限元模型概述 (3) 1 建立有限元模型 (4) 1.1 定义工作文件和工作标题 (4) 1.2 定义单元类型、实常数和材料 (4) 1.3 建立几何模型 (5) 2 网格划分 (7) 3 加载与初始地应力模拟 (8) 3.1 设置分析类型 (8) 3.2施加边界条件 (9) 3.3 施加上部面压力： (10) 3.4 施加重力加速度： (10) 3.4 设置加载步骤： (10) 3.5 求解初始地应力： (11) 3.6 保存分析结果： (11) 4 浏览初始地应力的计算结果（后处理） (12) 4.1显示变形形状： (12) 4.2查看节点结果等值线图： (12) 5 开挖巷道求解 (16) 5.1杀死巷道对应的单元： (16) 5.2输入加载步文件： (16) 5.3查看杀死巷道单元后的受力情况： (16) 5.4求解开挖后的有限元模型： (16) 6 浏览查看开挖后的计算结果（后处理） (17) 6.1显示变形形状： (17) 6.2查看节点结果等值线图： (17) 7 应力集中分析 (22) 8 总结体会 (22)
巷道开挖过程的有限元模型与力学分析 某半圆形拱巷道断面，参数见图，其所处地质条件为IV级围岩，上覆盖层厚度为100米后，各材料的力学参数见表。 有限元模型概述 本题采用ANSYS有限元分析软件模拟巷道开挖过程。由于地下巷道属于细长结构物，即巷道的横断面相对于纵向的长度来说很小，可且假定在围着荷载作用下，在其纵向没有位移，只有横向发生位移. 所以，巷道的力学分析可以采用弹性力学理论中的平面应变模型进行，这是一个较复杂的非线性力学问题。采用ANSYS有限元分析软件对巷道开挖进行模拟时，应首先根据地质条件建立合适的地下有限元分析模型，由于巷道对整个地下空间来说是属于“小孔口问题”，巷道周围出现孔口应力集中，并且应力集中区域影响范围约大于1.5倍的孔口尺寸，因此建几何模型时要选择合适的尺寸，其次采用PLANE42单元类型来分析平面应变问题。接下来设定单元尺寸大小划分网格。然后分步求解载荷，第一步求解初始地应力，第二步杀死巷道对应的单元后求解开挖后的载荷，最后分析巷道周围岩石的、位移、应力、应变的变化。
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非线性分析的几点忠告
关于非线性分析的几点忠告 了解程序的运作方式和结构的表现行为 如果你以前没有使用过某一种特别的非线性特性，在将它用于大的，复杂的模型前，构造一个非常简单的 模型（也就是，仅包含少量单元），以及确保你理解了如何处理这种特性。 通过首先分析一个简化模型，以便使你对结构的特性有一个初步了解。对于非线性静态模型，一个初步的 线性静态分析可以使你知道模型的哪一个区域将首先经历非线性响应，以及在什么载荷范围这些非线性将 开始起作用。对于非线性瞬态分析，一个对梁，质量块及弹簧的初步模拟可以使你用最小的代价对结构的 动态有一个深入了解。在你着手最终的非线性瞬时动态分析前，初步非线性静态，线性瞬时动态，和/或模 态分析同样地可以有助于你理解你结构的非线性动态响应的不同的方面。 阅读和理解程序的输出信息和警告。至少，在你尝试后处理你的结果前，确保你的问题收敛。对于与路程 相关的问题，打印输出的平衡迭代记录在帮助你确定你的结果是有效还是无效方面是特别重的。 简化 尽可能简化最终模型。如果可以将3─D结构表示为2─D平面应力，平面应变或轴对称模型，那么这样做， 如果可以通过对称或反对称表面的使用缩减你的模型尺寸，那么这样做。（然而，如果你的模型非对称加 载，通常你不可以利用反对称来缩减非线性模型的大小。由于大位移，反对称变成不可用的。）如果你可 以忽略某个非线性细节而不影响你模型的关键区域的结果，那么这样做。 只要有可能就依照静态等效载荷模拟瞬时动态加载。 考虑对模型的线性部分建立子结构以降低中间载荷或时间增量及平衡迭代所需要的计算时间。 采用足够的网格密度 考虑到经受塑性变形的区域要求一个合理的积分点密度。每个低阶单元将提供和高阶单元所能提供的一样
ansys实例解析SOLID45例题 手把手教你学ansys
2.1 三维实体元SOLID45 2.1.1 单元简介 1．单元形状 SOLID45单元是三维实体单元，该单元在每个节点上有三个自由度（X、Y、Z方向的位移）。 SOLID45的几何形状 2．单元参数输入 对于SOLID45单元而言，其必须或可选的参数输入包括： ●节点：I、J、K、L、M、N、O、P ●自由度：UX, UY, UZ ●材料参数 ?EX（弹性模量）、EY、EZ ?PRXY, PRYZ, PRXZ (or NUXY, NUYZ, NUXZ)（泊松比） ?ALPX, ALPY, ALPZ (or CTEX, CTEY, CTEZ or THSX, THSY, THSZ)（热膨胀系数） ?DENS（密度） ?其它参数 3．结果的输出
SOLID45单元的结果输出 对于三维实体单元SOLID45而言，工程上主要关注其在荷载作用下的变形（挠度）及应力各分量等。对于实体结构的变形，可以直接通过ANSYS中的位移来获得，而应力各分量，包括S：X, Y, Z, XY, YZ, XZ，主应力分量S：1, 2, 3，等效应力S：EQV，都可以在后处理中绘制等值线图，便于观察。 2.1.2 实例分析 考虑一根10m高的混凝土柱，EX＝2e10Pa，Prxy=0.167；Dens=2400；截面尺寸为B*H=0.6m×0.8m，在自由端作用100KN的水平推力，分析其变形和内力。 ANSYS的分析过程： ●进入ANSYS工作界面；改变工作目录Utility Menu>Change Directory… ●进入前处理模块：Main Menu>Preprocessor ●定义单元类型； 点击Element Type > Add/Edit/Delete，在弹出的对话框中添加；
几个ansys经典实例(长见识)
平面问题斜支座的处理 如图5-7所示，为一个带斜支座的平面应力结构，其中位置2及3处为固定约束，位置4处为一个45o的斜支座，试用一个4节点矩形单元分析该结构的位移场。 (a)平面结构(b)有限元分析模型 图5-7 带斜支座的平面结构 基于ANSYS平台，分别采用约束方程以及局部坐标系的斜支座约束这两种方式来进行处理。 (7) 模型加约束 左边施加X,Y方向的位移约束 ANSYS Main Menu: Solution →Define Loads →Apply →-Structural→Displacement On Nodes →选取2，3号节点→OK →Lab2: All DOF(施加X,Y方向的位移约束) →OK 以下提供两种方法处理斜支座问题，使用时选择一种方法。 ?采用约束方程来处理斜支座 ANSYS Main Menu：Preprocessor →Coupling/ Ceqn →Constraint Eqn ：Const :0, NODE1:4, Lab1: UX,C1:1,NODE2:4,Lab2:UY,C2:1→OK 或者?采用斜支座的局部坐标来施加位移约束 ANSYS Utility Menu：WorkPlane →Local Coordinate System →Create local system →At specified LOC + →单击图形中的任意一点→OK →XC、YC、ZC分别设定为2，0，0,THXY：45 →OK ANSYS Main Menu：Preprocessor →modeling →Move / Modify →Rotate Node CS →To active CS → 选择4号节点 ANSYS Main Menu：Solution →Define Loads →Apply →Structural →Displacement On Nodes →选取4号节点→OK →选择Lab2:UY(施加Y方向的位移约束) →OK 命令流； !---方法1 begin----以下的一条命令为采用约束方程的方式对斜支座进行处理 CE,1,0,4,UX,1,4,UY,-1 !建立约束方程(No.1): 0=node4_UX*1+node_UY*(-1) !---方法1 end --- !--- 方法2 begin --以下三条命令为定义局部坐标系，进行旋转，施加位移约束 !local,11,0,2,0,0,45 !在4号节点建立局部坐标系 !nrotat, 4 !将4号节点坐标系旋转为与局部坐标系相同 !D,4,UY !在局部坐标下添加位移约束 !--- 方法2 end
Ansys第25例非线性分析综合应用实例
第25例非线性分析综合应用实例----钢板卷制成圆筒 本例介绍了综合利用ANSYS非线性分析功能模拟将钢板卷制成圆筒的方法和步骤。25.1问题描述 将钢板卷制成圆筒一般要使用卷板机。图25-1所示为对称式三辊卷板机， 该机器将钢板卷制成圆筒时分为三个步骤：首先，上辊下降使钢板发生挠曲，钢板挠曲线的最低点首先发生屈服；然后，下辊转动驱动钢板向前移动，使钢板各点发生同样的屈服形成圆筒；最后，圆筒卷制完成，上辊上升卸下筒体。 图25-1对称式三辊卷板机 用ANSYS模拟将钢板卷制成圆筒，相应地也分为三个步骤。由于第二个步骤需要模拟上、下辊转动，而ANSYS的SOLIDn单元不支持大转动，位移边界条件不能施加大的转动角度，所以上、下辊需要用壳单元建立有限元模型。上、下辊与钢板的作用需要用接触模拟，钢板卷制成圆筒材料发生屈服，产生大变形, 所以钢板卷制成圆筒包括状态非线性、材料非线性和结构非线性三种非线性。 用ANSYS模拟将钢板卷制成圆筒，计算结果可以得到圆筒直径与上辊下压量的关系，上、下辊受力大小，上、下辊的变形，下辊驱动力矩及卸载回弹等重
25.2 命令流 /CLEAR /FILNAM, EXAMPLE25 /CONFIG, NRES, 2000 /PREP7 /PNUM, VOLU, ON ET, 1, SHELL181 ET, 2, SOLID186 MP, EX, 1, 2E11 MP, DENS, 1, 7800 MP, NUXY, 1, 0.3 MP, EX, 2, 2E11 MP, DENS, 2, 7800 MP, NUXY, 2, 0.3 TB, BKIN, 2, 1 TBTEMP, 0 TBDATA,, 240E6, 0 SECTYPE, 1, SHELL SECDATA, 0.02 CYLIND, 0.38/2, 0, 0.2, 1.7, 0, 360 要数据。因为分析过程复杂，步骤较多，所以本例只采用命令流法执行命令。 !清除数据库，新建文件 ! 指定任务名为?EXAMPLE25 “ !设置最大子步数 !前处理 !进入前处理器 !打开体号 !选择单元类型，壳单元用于划分上、下辊 !实体单元用于划分钢板 !定义材料模型 1 的弹性模量 ! 定义材料模型 1 的密度 !定义材料模型 1 的泊松比 !定义材料模型 2 的弹性模量 ! 定义材料模型 2 的密度 !定义材料模型 2 的泊松比 ! 定义材料模型 2 的屈服极限、切向模量 !定义截面 !壳厚度
ansys学习非线性静态分析实例
a n s y s学习非线性静态分 析实例 Newly compiled on November 23, 2020
ansys学习－非线性静态分析实例 问题描述 一个子弹以给定的速度射向壁面。壁面假定是刚性的和无摩擦的。将研究子弹和壁面接触后达80微秒长的现象。目的是确定子弹的整个变形，速度历程，以及最大等效Von Mises应变。求解使用SI单位。 用轴对称单元模拟棒。求解最好能通过单一载荷步实现。在这个载荷步中，将同时施加初始速度和约束。将圆柱体末端的节点Y方向约束住以模拟一固壁面。打开自动时间分步来允许ANSYS确定时间步长。定义分析结束的时间为8E-5秒，以确保有足够长的时间来扑捉整个变形过程。 问题详细说明 下列材料性质应用于这个问题： EX= (杨氏模量） DENS= （密度） NUXY=（泊松比） Yield Strength=（屈服强度） Tangent Modulus （剪切模量） 下列尺寸应用于这个问题： 长=-3m 直径=-3m 对于这个问题的初始速度是。 图1铜圆柱体图解 求解步骤： 步骤一：设置分析标题 1、选择菜单路径：Utility Menn>File>ChangeTitle。 2、键入文字“Coppery Cylinder Impacting a Rigid Wall” 3、单击OK。 步骤二：定义单元类型
1、选择菜单路径Mail Menu>Preprocessor>Element Type>All/Edit/Delete。 2、单击Add。Library of Element Types(单元类型库）对话框出现。 3、在靠近左边的列表中，单击“Visio Solid”仅一次。 4、选靠近右边的列表中，单击“4node Plas 106”仅一次。 5、单击OK。Library of Element Types 对话框关闭。 6、单击Options (选项）。VISCO106 element type Options(visco106单元类型选项）对话框出现。 7、在关于element behavior（单元特性）的卷动柜中，卷动到“Axisymmetric” 且选中它。 8、单击OK。 9、单击Element Types (单元类型）对话框中的Close。 步骤三：定义材料性质 1、选择菜单路径Main Menu>Preprocessor>Material Props>-Constant-Isotropic. Isotropic Matersal Properties (各向同性材料性质）对话框出现。 2、单击OK来指定材料号为1。另一个I sotropic Material Properties对话框出现。 3、对杨氏模量（EX）键入 4、对密度（DENS）键入8930。 5、对泊松比（NUXY）键入。 6、单击OK。 步骤四：定义双线性各向同性强化数据表（BISO） 1、选择菜单路径Main Menu>Preprocessor>Matersal Props>Data Tables> Define/Activate . Define/Activate Data Table(定义数据表）对话柜出现。 2、在关于type of data table(数据表类型）的卷动框中，卷动到“Bilin isotr BISO”且选中它。 3、对material reference number(材料参考号）健入1。 4、对number of temperatures(温度数）键入1和单击OK。 5、选择菜单路径Main Menu>Preprocessor>Material Props>Data Tables>Edit Active. Data Table BISO对话框出现。 6、对YLD Strs(屈服应力）键入。 7、对 Tang Mod(剪切模量）键入。
一个经典的ansys热分析实例(流程序)
/PREP7 /TITLE,Steady-state thermal analysis of pipe junction /UNITS,BIN ! 英制单位；Use U. S. Customary system of units (inches) ! /SHOW, ! Specify graphics driver for interactive run ET,1,90 ! Define 20-node, 3-D thermal solid element MP,DENS,1,.285 ! Density = .285 lbf/in^3 MPTEMP,,70,200,300,400,500 ! Create temperature table MPDATA,KXX,1,,8.35/12,8.90/12,9.35/12,9.80/12,10.23/12 ! 指定与温度相对应的数据材料属性；导热系数；Define conductivity values MPDATA,C,1,,.113,.117,.119,.122,.125 ! Define specific heat values（比热） MPDATA,HF,2,,426/144,405/144,352/144,275/144,221/144 ! Define film coefficient;除144是单位问题，上面的除12也是单元问题 ! Define parameters for model generation RI1=1.3 ! Inside radius of cylindrical tank RO1=1.5 ! Outside radius Z1=2 ! Length RI2=.4 ! Inside radius of pipe RO2=.5 ! Outside pipe radius Z2=2 ! Pipe length CYLIND,RI1,RO1,,Z1,,90 ! 90 degree cylindrical volume for tank WPROTA,0,-90 ! 旋转当前工作的平面；从Y到Z旋转-90度；；Rotate working plane to pipe axis CYLIND,RI2,RO2,,Z2,-90 ! 角度选择在了第四象限；90 degree cylindrical volume for pipe WPSTYL,DEFA ! 重新安排工作平面的设置；另外WPSTYL,STAT to list the status of the working plane；；Return working plane to default setting BOPT,NUMB,OFF ! 关掉布尔操作的数字警告信息；Turn off Boolean numbering warning VOVLAP,1,2 ! 交迭体；Overlap the two cylinders /PNUM,VOLU,1 ! 体编号打开；Turn volume numbers on /VIEW,,-3,-1,1
ANSYS温度场例题分析
短圆柱体的热传导过程 问题：一短圆柱体，直径和高度均为1m，现在其上端面施加大小为100℃的均匀温度载荷，圆柱体下端面及侧面的温度均为0℃，试求圆柱体内部的温度场分布（假设圆柱体不与外界发生热交换）。圆柱体材料的热传导系数为30W/(m·℃)。 求解： 第一步：建立工作文件名和工作标题 在ANSYS软件中建立相应的文件夹，并选择Thermal复选框。 第二部：定义单元类型 在单元类型（element type）中选择thermal solid和quad 4node 55，在单元类型选择数字（element type reference number）输入框中输入1，在单元类型选择框里选择Axisymmetric，其余默认即可。 第三步：定义材料性能参数 在材料性能参数对话框中输入圆柱体的导热系数30. 第四步：创建几何模型、划分网格 创建数据点，输入点坐标。在第一个输入框中输入关键点编号1，并输入第一个关键点坐标0、0、0，重复输入第二个、第三个、第四个关键点，相应的坐标分别为2（0.5,0,0）；3（0.5,1,0）；4（0,1,0）。结果如下图1所示：
在模型中创建直线，选择编号为1、2的关键点生成一条直线，在选取2、3生成一条直线，同样选择编号为3、4和编号为4、1的关键点生成另外两条直线。 结果如下图2所示： 之后在plot numbering controls对话框，分别打开KP Keypoint numbers、LINE line numbers、AREA Area numbers，建立直线L1、L2、L3、L4线段。生成几 何模型，如下图所示：
ansys学习非线性静态分析实例
a n s y s学习非线性静态 分析实例 Document serial number【NL89WT-NY98YT-NC8CB-NNUUT-NUT108】
ansys学习－非线性静态分析实例 问题描述 ??? 一个子弹以给定的速度射向壁面。壁面假定是刚性的和无摩擦的。将研究子弹和壁面接触后达80微秒长的现象。目的是确定子弹的整个变形，速度历程，以及最大等效Von Mises应变。求解使用SI单位。 用轴对称单元模拟棒。求解最好能通过单一载荷步实现。在这个载荷步中，将同时施加初始速度和约束。将圆柱体末端的节点Y方向约束住以模拟一固壁面。打开自动时间分步来允许ANSYS确定时间步长。定义分析结束的时间为8E-5秒，以确保有足够长的时间来扑捉整个变形过程。 问题详细说明 下列材料性质应用于这个问题： EX=? (杨氏模量） DENS=?? （密度） NUXY=（泊松比） Yield Strength=（屈服强度） Tangent Modulus??? （剪切模量） 下列尺寸应用于这个问题： 长=-3m 直径=-3m 对于这个问题的初始速度是。 ? ?????????????????????????????????????????????????????????? 图1铜圆柱体图解 求解步骤： 步骤一：设置分析标题 1、选择菜单路径：Utility Menn>File>ChangeTitle。 2、键入文字“Coppery Cylinder Impacting a Rigid Wall”
3、单击OK。 步骤二：定义单元类型 1、选择菜单路径Mail Menu>Preprocessor>Element Type>All/Edit/Delete。 2、单击Add。Library of Element Types(单元类型库）对话框出现。 3、在靠近左边的列表中，单击“Visio Solid”仅一次。 4、选靠近右边的列表中，单击“4node Plas 106”仅一次。 5、单击OK。Library of Element Types 对话框关闭。 6、单击Options (选项）。VISCO106 element type Options(visco106单元类型选项）对话框出现。 7、在关于element behavior（单元特性）的卷动柜中，卷动到“Axisymmetric” 且选中它。 8、单击OK。 9、单击Element Types (单元类型）对话框中的Close。 步骤三：定义材料性质 1、选择菜单路径Main Menu>Preprocessor>Material Props>-Constant-Isotropic. Isotropic Matersal Properties (各向同性材料性质）对话框出现。 2、单击OK来指定材料号为1。另一个I sotropic Material Properties对话框出现。 3、对杨氏模量（EX）键入 4、对密度（DENS）键入8930。 5、对泊松比（NUXY）键入。 6、单击OK。 步骤四：定义双线性各向同性强化数据表（BISO） 1、选择菜单路径Main Menu>Preprocessor>Matersal Props>Data Tables>????? Define/Activate . Define/Activate Data Table(定义数据表）对话柜出现。 2、在关于type of data table(数据表类型）的卷动框中，卷动到“Bilin isotr BISO”且选中它。 3、对material reference number(材料参考号）健入1。 4、对number of temperatures(温度数）键入1和单击OK。