接触分析实例

实例1：带孔方板的接触分析实例

静力分析分为线性静力分析和非线性静力分析。线性分析，即外载荷与系统的响应之间为线性关系。而在真实的物理结构中，结构刚度随变形而发生变化，即所谓的“非线性分析”。

非线性问题的分类：

材料非线性（material nonlinearity）:材料的应力—应变关系为非线性，即本构关系为非线性。几何非线性：指位移的大小对结构响应的影响，包括大位移、大转动、初始应力、几何刚性化和突然翻转等，实例：过盈装配

边界条件非线性：即边界条件在分析过程中不断的发生变化，实例：接触问题

后处理：

(1)延展平面应力单元来构造三维视图

View—odb display option—sweep&extrude

(2)显示接触压强(CPRESS)和接触状态(COPEN)

Result—field output

接触分析中的主要问题及解决收敛问题的常用方法：

1.接触面的分类：

由单元构成的柔性或刚性接触面；节点构成的接触面；解析刚体接触面

2.一个接触对最多有一个由节点构成的接触面；如果只有一个接触面，称为“自接触”

3.模拟接触问题的两种算法：

通用接触算法和接触对算法。其中后一种算法，对接触面的类型有较多的限制，需要指定相互接触的面

4.定义接触对：

主面必须连续，且由节点组成；如果是有限滑移，主面发生接触的部位必须光滑，不能有尖角，不要让从面节点落在主面之外（尤其不能落在主面的背面）；接触面在发生接触的部位有较大尖角或凹角时，应该将其分别定义为两个面

5.有限滑移和小滑移

在有限滑移过程中，需不断的判断从面节点和主面的哪部分接触；而在小滑移中，在分析开始时就定义了接触的部位，而在分析过程中不会发生变化

小滑移也可以用于几何非线性问题，有两种算法：点对面和面对面算法

实例2：过盈装配过程模拟

边界条件给出的位移值是相对于模型初始状态的绝对位移值，而不是当前分析步中的增量值；可以控制在每个增量步中的位移均与；轴对称模型将三维模型转化为二维平面问题，只有三个自由度。

后处理：

1.延展轴对称单元来构造等效的三维视图；显示接触压强；显示压头节点的反作用力RF

接触分析中用到的输出变量：

CFS: Total force due to frictional stress (CFS n, n = 1, 2, 3).，摩擦应力的合力

CFN:接触压力的合力

CAREA:接触面积

CMN:接触压力相对于原点的合力矩

CMS:摩擦压力相对于原点的合力矩

CFT:接触压力和摩擦压力的合力

CMT:接触压力和摩擦压力相对于原点的合力矩

面节点输出变量：

STATUS:接触状态

CPRESS:接触压强

CSHEAR1：局部方向1上的摩擦剪应力

COPEN:从面节点与主面的距离

CSLIP1:在局部方向1上的相对切向滑移

Flac3D命令--完整经典版

实例分析命令： 1. X ，Y ，Z 旋转 Shift+ X ，Y ，Z 反向旋转 Gen zone ……；model ……；prop ……（材料参数）；set grav 0,0,-9.81（重力加速度） plot add block group red yellow 把在group 中的部分染成红色和黄色 plot add axes black 坐标轴线为黑色；print zone stress% K 单元应力结果输出 ini dens 2000 ran z a b （设置初始密度，有时不同层密度不同）；ini ……（设置初始条件）；fix ……（固定界面） set plot jpg ；set plot quality 100 ；plot hard file 1.jpg 图像输出（格式、像素、名称） plot set magf 1.0视图的放大倍数为1.0；plo con szz z 方向应力云图 2. ini z add -1 range group one 群one 的所有单元，在z 方向上向下移动1m ；然后合并 命令 gen merge 1e-5 range z 0此命令是接触面单元合并成一个整体，1e-5是容差 3. (基坑开挖步骤)：Step 1: create initial model state （建立初始模型）Step 2: excavate trench （开挖隧道） 4. group Top range group Base not 定义（群组Base 以外的为）群组Top 5. plot blo gro 使得各个群组不同颜色显示 6. （两个部分间设置界面；切割法）：gen separate Top 使两部分的接触网格分离 为两部分；interface 1 wrap Base Top 在（Base 和Top ）这两部分之间添加接触单元；plot create view_int 显示，并创建标题view_int ；plot add surface 显示表面；plot add interface red 界面颜色红色 7. （简单的定义函数及运行函数）new ；def setup 定义函数setup ；numy = 8定义常 量numy 为8；depth = 10.0 定义depth 为10；end 结束对函数的定义；setup 运行函数setup 8. （隧道生成）上部圆形放射性圆柱及下部块体单元体的建立，然后镜像。 9. 模拟模型的材料问题时为什么要去定义某个方向上的初始速度？— 10. 渐变应力施加：apply nstress -1e6 gradient 0,0,1e5 range z 3.464,0 plane dip 60 dd 270 origin .1 0 0；施加法向应力：apply nstress -1e6 range plane dip 60 dd 270 origin .1 0 0 11. d ip dd 确定平面位置使用：（纠结） 12. p rint gp position range id=14647 输出节点坐标 13. a pply sxx -10e6 gradient 0 , 0, 1e5 range z -100 , 0在这个求解方程中，z 为变量，所以xx σ为：65=-1010+10xx z σ?? ；原点（0,0,0） 14. f ree x range x -.1 .1 z 6.9 10.1放松x=0 平面上，z=7，10 这一部分在x 方向的约 束（可以在此处产生破坏） 15. 体积模量K 和剪切模量G 与杨氏模量及泊松比v 之间的转换关系如下： =3(1-2v)E K G=2(1+v) E 16. 一般而言，大多数问题可以采用FLAC 3D 默认的收敛标准（或称相对收敛标准），即当体 系最大不平衡力与典型内力的比率R 小于定值10-5；（也可由用户自定义该值，命令：

abaqus屈曲分析实例

整个计算过程包括2个分析步，第1步做屈曲分析，笫2步做极限强度分析。 第1步：屈曲分析 载荷步定义如下： Step 1-Initial Step 2- Buckle

? Re Mbs M^nce C^wvoini live 2oc*$ *l^*?4 tjdp V ：i.Jsa&# 录 +r A AJIu fffiC? fe3 Ha ? ；r????y fa-t n>rr ?: OfEYcm v Se?今 gh 3, gqcvKeiry C*p*?9r ? ? O?lec? ■ %?no?v C5 廉 H5Wr> MM fa Tin* Forti Sv Al€ *dep6?? ve^ tbjUx9)lo t JeiWA Tc?D -^lQZlll?hQ we' E ejewwiw b>w* biE Glcte 」r?>w* 69D eJe*MKi r?jw* bee CWfcr*?9*^ s￡ Zac? “ Iraftet H U 匕“rb ? 2 更 K?4dCu^u!R? 虫 Hntwr GUput b 伽》ezi5 &■心 AcUxv? V H H?*?ctnr? 易 htecMtlar. hra, 日 CcrtadCcrtra 0 C?Wl >?wt K Ccctect sub lx 權 CwMoarSt Hj fiUdi _n ,.. ? ?! ? MCg WtW Swtfc lk2 pe**j". liwar p?nwbia?ko ▼ freque." 拯 sufAuun The 11?-51>^ )L>4ldH9jjn-2 “9 wioZ S *0 Sxe U>* oil^ 51 “ed S iU* TO . 0 . -ISO -MO mtb rew ：t no 心 &逐Ply OCCOIIMV * 巧恪tc ?：?L -5Moe>?* bw tZfft to ?D7cp 炉、？ZlHWr? Me" “乡“r?x HMldrann ?2 vd 乡 tygeJa* 400 0 0 with x*w ：? ?o tfi* oc

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).查看结果 模态分析的结果包括结构的频率、振型、相对应力和力等

FLAC3D 实例命令流1

第1部分命令流按照顺序进行 2-1定义一个FISH函数 new def abc abc = 25 * 3 + 5 End print abc 2-2使用一个变量 new def abc hh = 25 abc = hh * 3 + 5 End Print hh Print abc 2-3对变量和函数的理解 new def abc hh = 25 abc = hh * 3 + 5 End set abc=0 hh=0 print hh print abc print hh new def abc abc = hh * 3 + 5 end set hh=25 print abc set abc=0 hh=0 print hh print abc print hh 2-4获取变量的历史记录 new gen zone brick size 1 2 1 model mohr prop shear=1e8 bulk=2e8 cohes=1e5 tens=1e10 fix x y z range y -0.1 0.1 apply yvel -1e-5 range y 1.9 2.1 plot set rotation 0 0 45

plot block group def get_ad ad1 = gp_near(0,2,0) ad2 = gp_near(1,2,0) ad3 = gp_near(0,2,1) ad4 = gp_near(1,2,1) end get_ad def load load=gp_yfunbal(ad1)+gp_yfunbal(ad2)+gp_yfunbal(ad3)+gp_yfunbal(ad4) end hist load hist gp ydis 0,2,0 step 1000 plot his 1 vs -2 2-5用FISH函数计算体积模量和剪砌模量 new def derive s_mod = y_mod / (2.0 * (1.0 + p_ratio)) b_mod = y_mod / (3.0 * (1.0 - 2.0 * p_ratio)) end set y_mod = 5e8 p_ratio = 0.25 derive print b_mod print s_mod 2-6 在FLAC输入中使用符号变量 New def derive s_mod = y_mod / (2.0 * (1.0 + p_ratio)) b_mod = y_mod / (3.0 * (1.0 - 2.0 * p_ratio)) end set y_mod = 5e8 p_ratio = 0.25 derive gen zone brick size 2,2,2 model elastic prop bulk=b_mod shear=s_mod print zone prop bulk print zone prop shear 2-7 控制循环 New def xxx sum = 0 prod = 1 loop n (1,10)

ANSYS高级分析-子结构

1 引言 在ANSYS平台上，所谓子结构技术就是将一组单元用矩阵凝聚为一个单元过程的技术，切吧这个单一的矩阵单元称为超单元。在ANSYS分析中，超单元可以象其他单元类型一样使用。唯一的区别就是必须先进行结构生成分析以生成能够利用的超单元。但子结构并非在所有ANSYS模块中都能利用，目前ANSYS子结构技术可以在ANSYS/Mutiphysics，ANSYS/Mechanical和ANSYS/Structural中使用。 在ANSYS平台上，使用子结构的目的主要是为了节省机时，并且允许在比较有限的计算机设备资源的基础上求解超大规模的问题。比如进行非线性分析和带有大量重复几何结构的分析。在非线性分析中，可以将模型线性部分作成子结构，这部分的单元矩阵就不用在非线性迭代过程中重复计算。而在有重复几何结构的模型中（如有四条腿的桌子），可以对于重复的部分生成超单元，然后将它拷贝到不同的位置，这样做可以节省大量的计算时间和计算机资源。 子结构还用于模型有大转动的情况下。对于这些模型，ANSYS假定每个结构都是围绕其质心转动的。在三维情况下，子结构有三个转动自由度和三个平动自由度。在大转动模型中，用户在使用部分之前无须对子结构施加约束，因为每个子结构都是作为一个单元进行处理，是允许刚体位移的。 对于大型三维问题的分析而言，需用磁盘空间相对于一个普通计算机系统来说太庞大了，在这种情况下，用户可以通过子结构将问题分块进行分析，从而使得每一块对于计算机系统来说都是可以计算和承受的。 2 ANSYS子结构使用步骤 ANSYS子结构使用过程分为以下三个步骤： 1）ANSYS子结构生成部分 生成部分就是将普通的有限元单元凝聚为一个超单元。凝聚是通过定义一组主自由度来实现的。主自由度用于定义超单元与模型中其他单元的边界，提取模型的动力学特性。图1是一个板状构件用接触单元分析的示意。由于接触单元需要迭代计算，将板状构件形成子结构将显著地节省机时。本例中，主自由度是板与接触单元相连的自由度。 图1 子结构使用示例 2）ANSYS子结构使用部分 用部分就是将超单元与模型整体相连进行分析的部分。整个模型可以是一个超单元，也可以象上例一样是超单元与非超单元相连的。使用部分的计算只是超单元的凝聚（自由度计算仅限于主自由度）和非超单元的全部计算。

(完整word版)abaqus6.12-典型实例分析

1.应用背景概述 随着科学技术的发展，汽车已经成为人们生活中必不可少的交通工具。但当今由于交通事故造成的损失日益剧增，研究汽车的碰撞安全性能，提高其耐撞性成为各国汽车行业研究的重要课题。目前国内外许多著名大学、研究机构以及汽车生产厂商都在大力研究节省成本的汽车安全检测方法，而汽车碰撞理论以及模拟技术随之迅速发展，其中运用有限元方法来研究车辆碰撞模拟得到了相当的重视。而本案例就是取材于汽车碰撞模拟分析中的一个小案例―――保险杠撞击刚性墙。 2.问题描述 该案例选取的几何模型是通过导入已有的*.IGS文件来生成的（已经通过Solidworks软件建好模型的），共包括刚性墙（PART-wall）、保险杠（PART-bumper）、平板（PART-plane）以及横梁（PART-rail）四个部件，该分析案例的关注要点就是主要吸能部件（保险杠）的变形模拟，即发生车体碰撞时其是否能够对车体有足够的保护能力？这里根据具体车体模型建立了保险杠撞击刚性墙的有限元分析模型，为了节省计算资源和时间成本这里也对保险杠的对称模型进行了简化，详细的撞击模型请参照图1所示，撞击时保险杠分析模型以2000mm/s的速度撞击刚性墙，其中分析模型中的保险杠与平板之间、平板与横梁之间不定义接触，采用焊接进行连接，对于保险杠和刚性墙之间的接触采用接触对算法来定义。 1.横梁（rail） 2.平板（plane） 3.保险杠（bumper） 4.刚性墙（wall） 图2.1 碰撞模型的SolidWorks图

为了使模拟结果尽可能真实，通过查阅相关资料，定义了在碰撞过程中相关的数据以及各部件的材料属性。其中，刚性墙的材料密度为7.83×10-9，弹性模量为2.07×105，泊松比为0.28；保险杠、平板以及横梁的材料密度为7.83×10-9，弹性模量为2.07×105，泊松比为0.28，塑形应力-应变数据如表2.1所示。 表2.1 应力-应变数据表 应力210 300 314 325 390 438 505 527 应变0.0000 0.0309 0.0409 0.0500 0.1510 0.3010 0.7010 0.9010 注：本例中的单位制为：ton，mm，s。 3.案例详细求解过程 本案例使用软件为版本为abaqus6.12，各详细截图及分析以该版本为准。3.1 创建部件 （1）启动ABAQUS/CAE，创建一个新的模型数据库，重命名为The crash simulation，保存模型为The crash simulation.cae。 （2）通过导入已有的*.IGS文件来创建各个部件，在主菜单中执行【File】→【Import】→【Part】命令，选择刚刚创建保存的的bumper_asm.igs文件，弹出【Create Part From IGS File】对话框如图3.1所示，根据图3.1所示设定【Repair Options】的相关选项，其它参数默认，单击【Ok】按钮，可以看到在模型树中显示了导入的部件bumper_asm。 图3.1 Create Part From IGS File对话框

ABAQUS时程分析实例

ABAQUS时程分析法计算地震反应得简单实例ABAQUS时程分析法计算地震反应得简单实例（在原反应谱模型上 修改） 问题描述： 悬臂柱高12m，工字型截面（图1），密度7800kg/m3，EX=2、1e11Pa，泊松比0、3，所有振型得阻尼比为2%，在3m高处有一集中质量160kg，在6m、9m、12m处分别有120kg 得集中质量。反应谱按7度多遇地震，取地震影响系数为0、08，第一组，III类场地，卓越周期Tg=0、45s。 图1 计算对象 第一部分：反应谱法 几点说明： λ本例建模过程使用CAE； λ添加反应谱必须在inp中加关键词实现，CAE不支持反应谱； λ *Spectrum不可以在keyword editor中添加，keyword editor不支持此关键词读入。 λ ABAQUS得反应谱法计算过程以及后处理要比ANSYS方便得多。 操作过程为： （1）打开ABAQUS/CAE，点击create model database。

（2）进入Part模块，点击create part，命名为column，3D、deformation、wire。continue （3）Create lines，在 分别输入0，0回车；0，3回车；0，6回车；0，9回车；0，12回车。

（4）进入property模块，create material，name：steel，general-->>density，mass density：7800 mechanical-->>elasticity-->>elastic，young‘s modulus：2、1e11，poisson’s ratio：0、3、

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面与面接触分析实例

面与面接触实例：插销拨拉问题分析 定义单元类型 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

FLAC3D 实例命令流1

第1部分命令流按照顺序进行2-1定义一个FISH函数 new def abc abc = 25 * 3 + 5 End print abc 2-2使用一个变量 new def abc hh = 25 abc = hh * 3 + 5 End Print hh Print abc 2-3对变量和函数的理解 new def abc hh = 25 abc = hh * 3 + 5 End set abc=0 hh=0 print hh print abc print hh new def abc abc = hh * 3 + 5 end set hh=25 print abc set abc=0 hh=0 print hh print abc print hh 2-4获取变量的历史记录 new gen zone brick size 1 2 1 model mohr prop shear=1e8 bulk=2e8 cohes=1e5 tens=1e10

fix x y z range y -0.1 0.1 apply yvel -1e-5 range y 1.9 2.1 plot set rotation 0 0 45 plot block group def get_ad ad1 = gp_near(0,2,0) ad2 = gp_near(1,2,0) ad3 = gp_near(0,2,1) ad4 = gp_near(1,2,1) end get_ad def load load=gp_yfunbal(ad1)+gp_yfunbal(ad2)+gp_yfunbal(ad3)+gp_yfunbal(ad4) end hist load hist gp ydis 0,2,0 step 1000 plot his 1 vs -2 2-5用FISH函数计算体积模量和剪砌模量 new def derive s_mod = y_mod / (2.0 * (1.0 + p_ratio)) b_mod = y_mod / (3.0 * (1.0 - 2.0 * p_ratio)) end set y_mod = 5e8 p_ratio = 0.25 derive print b_mod print s_mod 2-6 在FLAC输入中使用符号变量 New def derive s_mod = y_mod / (2.0 * (1.0 + p_ratio)) b_mod = y_mod / (3.0 * (1.0 - 2.0 * p_ratio)) end set y_mod = 5e8 p_ratio = 0.25 derive gen zone brick size 2,2,2 model elastic prop bulk=b_mod shear=s_mod print zone prop bulk print zone prop shear

(完整word版)ABAQUS实例分析

《现代机械设计方法》课程结业论文 （ 2011 级） 题目：ABAQUS实例分析 学生姓名 XXXX 学号 XXXXX 专业机械工程 学院名称机电工程与自动化学院 指导老师 XX 2013年 5 月8 日

目录 第一章Abaqus简介 (1) 一、Abaqus总体介绍 (1) 二、Abaqus基本使用方法 (2) 1.2.1 Abaqus分析步骤 (2) 1.2.2 Abaqus/CAE界面 (3) 1.2.3 Abaqus/CAE的功能模块 (3) 第二章基于Abaqus的通孔端盖分析实例 (4) 一、工作任务的明确 (6) 二、具体步骤 (6) 2.2.1 启动Abaqus/CAE (4) 2.2.2 导入零件 (5) 2.2.3 创建材料和截面属性 (6) 2.2.4 定义装配件 (7) 2.2.5 定义接触和绑定约束（tie） (10) 2.2.6 定义分析步 (14) 2.2.7 划分网格 (15) 2.2.8 施加载荷 (19) 2.2.9 定义边界条件 (20) 2.2.10 提交分析作业 (21) 2.2.11 后处理 (22) 第三章课程学习心得与作业体会 (23)

第一章： Abaqus简介 一、Abaqus总体介绍 Abaqus是功能强大的有限元分析软件，可以分析复杂的固体力学和结构力学系统，模拟非常庞大的模型，处理高度非线性问题。Abaqus不但可以做单一零件的力学和多物理场的分析，同时还可以完成系统级的分析和研究。 Abaqus使用起来十分简便，可以很容易的为复杂问题建立模型。Abaqus具备十分丰富的单元库，可以模拟任意几何形状，其丰富的材料模型库可以模拟大多数典型工程材料的性能，包括金属、橡胶、聚合物、复合材料、钢筋混泥土、可压缩的弹性泡沫以及地质材料（例如土壤、岩石）等。 Abaqus主要具有以下分析功能： 1.静态应力/位移分析 2.动态分析 3.非线性动态应力/位移分析 4.粘弹性/粘塑性响应分析 5.热传导分析 6.退火成形过程分析 7.质量扩散分析 8.准静态分析 9.耦合分析 10.海洋工程结构分析 11.瞬态温度/位移耦合分析 12.疲劳分析 13.水下冲击分析 14.设计灵敏度分析 二、Abaqus基本使用方法 1.2.1 Abaqus分析步骤 有限元分析包括以下三个步骤： 1.前处理（Abaqus/CAE）:在前期处理阶段需要定义物理问题的模型，并生 成一个Abaqus输入文件。提交给Abaqus/Standard或 Abaqus/Explicit。 2.分析计算（Abaqus/Standard或Abaqus/Explicit）:在分析计算阶段， 使用Abaqus/Standard或Abaqus/Explicit求解输入文件中所定义的

第20章 接触分析实例

第20章接触分析实例 在这个实例中，将对一个盘轴紧配合结构进行接触分析。第一个载荷步分析轴和盘在过盈配合时的应力，第二个载荷步分析将该轴从盘心拔出时轴和盘的接触应力情况。 20.1 问题描述： 在旋转机械中通常会遇到轴与轴承、轴与齿轮、轴与盘连接的问题，根据各自的不同情况可能有不同的连接形式。但大多数连接形式中存在过盈配合，也就是涉及到接触问题的分析。这里我们以某转子中轴和盘的连接为例，分析轴和盘的配合应力以及将轴从盘中拔处时盘轴连接处的应力情况。 本实例的轴为一等直径空心轴，盘为等厚度圆盘，其结构及尺寸如图20.1所示。由于模型和载荷都是轴对称的，可以用轴对称方法进行分析。这里为了后处理时观察结果更直观，我们采用整个模型的四分之一进行建模分析，最后将其进行扩展，来观察整个结构的变形及应力分布、变化情况。盘和轴用同一种材料，其性质如下： 弹性模量：EX=2.1E5 泊松比：NUXY=0.3 接触摩擦系数：MU=0.2 20.1 盘轴结构图

20.2 建立有限元模型 在ANSYS6.1中，首先我们通过完成如下工作来建立本实例的有限元模型，需要完成的工作有：指定分析标题，定义单元类型，定义材料性能，建立结构几何模型、进行网格划分等。根据本实例的结构特点，我们将首先建立代表盘和轴的两个1/4圆环面，然后对其进行网格划分，得到有限元模型。 20.2.1设置分析标题 本实例为进行如图20.1所示的盘轴结构的接触分析，属于非线性结构分析范畴。跟前面实例一样，为了在后面进行菜单方式操作时的方便，需要在开始分析时就指定本实例分析范畴为“Structural”。本实例的标题可以命名为：“Analysis of a Axis Contacting a hole in a Disc”，具体的操作过程如下： 1．选取菜单路径Utility Menu | File | Change Jobname，将弹出Change Jobname (修改文件名)对话框，如图20.2所示。在Enter new jobname (输入新文件名)文本框中输入文字“CH20”，为本分析实例的数据库文件名。并单击New log and error files (新的日志和错误 文件)单选框，使其变为“Yes”，为本实例的分析过程创建新的日志。单击按钮关闭对话框，完成文件名的修改。 图20.2 修改文件名对话框 2．选取菜单路径Utility Menu | File | Change Title，将弹出Change Title (修改标题)对话框，如图20.3所示。在Enter new title (输入新标题)文本框中输入文字“Analysis of a Axis Contacting a hole in a Disc”，为本分析实例的标题名。单击按钮，完成对标题名的指定。 图20.3 修改标题对话框 3．选取菜单路径Utility Menu | Plot | Replot，指定的标题“Analysis of a Axis Contacting a hole in a Disc”将显示在图形窗口的左下角(图略)。 4．选取菜单路径Main Menu | Preference，将弹出Preference of GUI Filtering (菜单过滤参数选择)对话框。单击Structual(结构)选项使之被选中，以将菜单设置为与结构分析相关

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 所示：

abaqus实例

一．创建部件 1.打开abaqus; 开始/程序/Abaqus6.10-1/Abaque CAE 2.Model/Rename/Model-1,并输入名字link4

3.单击Create part弹出Create part对话框， Name输入link-4; Modeling Space 选择2D Planar Type 选择Deformable Base Feature 选择Wire Approximate size 输入800；然后单击continue 4.单击(Create Lines:connected)通过点(0,0)、(400,0)、(400,300)、(0,300)单击(Create Lines:connected)连接(400,300)和(0,0)两点，单击提示区中的Done按钮(或者单击鼠标滚轮，也叫中键)，形成四杆桁架结构

5.单击工具栏中的(Save Model Database),保存模型为link4.cae 二.定义材料属性 6.双击模型树中的Materials(或者将Module切换到Property,单击Create Material -ε) 弹出Edit Material对话框后。 执行对话框中Mechanical/Elasticity/Elastic命令， 在对话框底部出现的Data栏中输入Young’s Module为29.5e4， 单击OK.完成材料设定。

7.单击“Create Section ”,弹出Create Section对话框， Category中选择Beam; Type中选择Truss; 单击continue按钮 弹出Edit Section对话框， 材料选择默认的Material-1,输入截面积(Cross-sectional area)为100，单击ok按钮。

FLAC3D快速入门及简单实例

FLAC3D快速入门 及简单实例 李佳宇编 LJY指南针教程

前言 FLAC及FLAC3D是由国际著名学者、英国皇家工程院院士、离散元的 发明人Peter Cundall博士在70年代中期开始研究的，主要面对岩土工程的通 用软件系统，目前已经在全球70多个国家得到广泛应用，在岩土工程学术界 和工业界赢得了广泛的赞誉。前国际岩石力学会主席 C.Fairhurst(1994)对 FLAC程序的评价是：“现在它是国际上广泛应用的可靠程序。” 我从研二(2010年)开始接触FLAC3D，最初的原因是导师要求每一个人至 少学会一个数值计算软件，而他嘴里每天念叨最多的就是FLAC，自己当时对数 值计算一无所知，便答应老师要学会FLAC3D。第一次打开软件界面，我心里 就凉了大半截，面对着一个操作界面跟记事本无异的所谓“功能强大”的岩土工 程专业软件，半点兴趣也提不起来。年底，从项目工地回到学校准备论文开题， 老师对我的开题报告非常不满意，当着全教研室师生的面，劈头盖脸大批一顿， 第二天又找谈话。在巨大的压力和强烈的自尊心驱使下，我硬着头皮开始啃 FLAC3D，一个半月之后，终于有了初步的计算结果，对老师有个交代，我也能 回家过年了。 前面这一段过程可能是大多数FLAC3D初学者的必经阶段，或者是即将 开始软件学习的人惧怕的事情。毫无疑问，FLAC3D极其不友好的界面是阻碍 初学者前进的很大障碍，当然还包括它是一个全英文的软件。但是当你费尽周折 的走进FLAC3D的世界，你就会发现它独特的魅力，比如简洁的界面，快捷的 命令流操作，高效的计算方法，不易报错等等。另外一个拿不上台面的优点就是 它非常小巧，包括Manual在内一共才几十兆大小，而且已经被破解成绿色版， 只要把它和命令流装进U盘，你就可以随便找一个身边功能最强大的电脑开始 计算了，如果你有过ANSYS、ABAQUS等大型软件痛苦的安装经历，你便能 毕业之后，本以为不用再接触数值计算，但工作需要使得我又一次开始与理解“绿色版”的含义，当然还请大家尊重知识产权，支持正版。 FLAC3D进行亲密接触。我的领导给了我很多新思路和很大的支持，如今我的 水平比研究生时有了不少提高。于是，我想把我的经验总结成文，希望对初学 者起

ANSYS高级接触问题

ANSYS高级接触问题 第一章接触问题概述 在工程中会遇到大量的接触问题，如齿轮的啮合、法兰联接、机电轴承接触、卡头与卡座、密封、板成形、冲击等等。接触是典型的状态非线性问题，它是一种高度非线性行为。接触例子如图1： 分析中常常需要确定两个或多个相互接触物体的位移、接触区域的大小和接触面上的应力分布。 接触分析存在两大难点： 在求解之前，你不知道接触区域、表面之间是接触或分开是未知的，表面之间突然接触或突然不接触会导致系统刚度的突然变化。 大多数接触问题需要计算摩擦。摩擦是与路径有关的现象，摩擦响应还可能是杂乱的，使问题求解难以收敛。 1.1 接触分类 1.1.1 刚－柔 一个表面是完全刚性的—除刚体运动外无应变、应力和变形，另一表面为软材料构成是可变形的。 只在一个表面特别刚硬并且不关心刚硬物体的应力时有效。 1.1.2 柔－柔

两个接触体都可以变形。 1.2 接触单元 ANSYS采用接触单元来模拟接触问题：跟踪接触位置；保证接触协调性（防止接触表面相互穿透）；在接触表面之间传递接触应力（正压力和摩擦）。 接触单元就是覆盖在分析模型接触面上的一层单元。在ANSYS中可以采用三种不同的单元来模拟接触：面—面接触单元；点—面接触单元；点—点接触单元。 不同的单元类型具有完全不同的单元特性和分析过程。 1.2.1 面—面接触单元 用于任意形状的两个表面接触，不必事先知道接触的准确位置；两个面可以具有不同的网格；支持大的相对滑动；支持大应变和大转动。例如：面一面接触可以模拟金属成型，如轧制过程。 1.2.2 点—面接触单元 用于某一点和任意形状的面的接触，可使用多个点－面接触单元模拟棱边和面的接触；不必事先知道接触的准确位置；两个面可以具有不同的网格；支持大的相对滑动；支持大应变和大转动。例：点面接触可以模拟棱边和面之间的接触。 1.2.3 点－点接触单元 用于模拟单点和另一个确定点之间的接触。建立模型时必须事先知道确切的接触位置；多个点－点接触单元可以模拟两个具有多个单元表面间的接触；每个

ABAQUS有限元接触分析的基本概念

ABAQUS有限元接触分析的基本概念2009-11-24 00:06:28 作者：jiangnanxue 来源：智造网—助力中国制造业创新—https://www.wendangku.net/doc/7f393963.html, CAE(计算机辅助工程)是一门复杂的工程科学，涉及仿真技术、软件、产品设计和力学等众多领域。世界上几大CAE公司各自以其独到的技术占领着相应的市场。ABAQUS有限元分析软件拥有世界上最大的非线性力学用户群，是国际上公认的最先进的大型通用非线性有限元分析软件之一。它广泛应用于机械制造、石油化工、航空航天、汽车交通、土木工程、国防军工、水利水电、生物医学、电子工程、能源、地矿、造船以及日用家电等工业和科学研究领域。ABAQUS在技术、品质和可靠性等方面具有卓越的声誉，可以对工程中各种复杂的线性和非线性问题进行分析计算。 《ABAQUS有限元分析常见问题解答》以问答的形式，详细介绍了使用ABAQUS建模分析过程中的各种常见问题，并以实例的形式教给读者如何分析问题、查找错误原因和尝试解决办法，帮助读者提高解决问题的能力。 《ABAQUS有限元分析常见问题解答》一书由机械工业出版社出版。 16.1.1 点对面离散与面对面离散 【常见问题16-1】 在ABAQUS/Standard分析中定义接触时，可以选择点对面离散方法(node-to-surface-dis - cre-tization)和面对面离散方法(surface-to-surface discretization)，二者有何差别? 『解答』 在点对面离散方法中，从面(slave surface)上的每个节点与该节点在主面(master surface)上的投影点建立接触关系，每个接触条件都包含一个从面节点和它的投影点附近的一组主面节点。 使用点对面离散方法时，从面节点不会穿透(penetrate)主面，但是主面节点可以穿透从面。 面对面离散方法会为整个从面(而不是单个节点)建立接触条件，在接触分析过程中同时考虑主面和从面的形状变化。可能在某些节点上出现穿透现象，但是穿透的程度不会很严重。 在如图16-l和图16-2所示的实例中，比较了两种情况。