workbench APDL

workbench APDL
workbench APDL

随着V12的发布,ANSYS公司把ANSYS,也就是通常称之为ANSYS Classic的,重命名为Mechanical APDL,并将其作为一个模块,置入Workbench界面下。或许这就在暗示经典ANSYS以后的位置。

对于很多习惯用Mechanical APDL的用户来说,一时半会还不会放弃它,更何况目前WB并不能完全取代Mechanical APDL。值得庆幸的是,ANSYS加强了Mechanical APDL 跟其他各个模块的连接。由此我们可以兼得鱼跟熊掌。

在一个分析模块中,我们可以在四处连上Mechanical APDL:Geomerty,Model,Setup

跟Solution。

连接方法有两种,一是从左边的component systems下将Mechanical APDL扯到合适的位置后放开

二是右键在相应处点击,选择Transfer data to new: Mechanical APDL

选择edit in Mechanical APDL就可以查看导入的模型

通过gplot,我们可以很明显地看到,在Geometry下只有几何模型导入

在Model下则只有有限元模型导入。值得注意的是,这时候的网格只是mesh200单元

当然事先定义的NamedSelection也可以通过componet的形式导入

在Setup下可以导入除结构网格外,还可以导入约束荷载。但在WB中荷载是通过table定

义,直接显示是其值为零。

在Solution可以模型结果导入Mechanical APDL。update后analysis仍然显示问号,但并不影响查看结果。怀疑是个bug。

现在假设我们在WB中导入几何模型,划分网格后,想在Mechanical APDL里施加荷载,定义边界,计算后进行后处理。

在Setup跟Solution中插入Mechanical APDL,Workbench总是要更新原有模型,有时候很骂人,干脆插入一个/eof以阻止更新

这是在Mechanical APDL中所需要的命令流,tapdl.inp。其中dnode,fforce是事先定义的Named Selection。并提出最大应力,以进行下一步的参数化分析。

在将tapdl.inp导入后,点击edit

将kraft跟seqvmax分别定义为input跟output参数

打开paramenter set,可以查看我们定义的参数。另外几个是在DM中定义的几何模型参数。

我们可以定义多个design point,一次点击便可完成所有计算。由input file定义的参数,同样可以参与DX优化,这以后再说。

【线性代数】之矩阵的乘法运算

Born T o Win 考研数学线性代数之矩阵的乘法运算 任意两个矩阵不一定能够相乘,即两个矩阵要相乘必须满足的条件是:只有当第一个矩阵A 的列数与第二个矩阵B 的行数相等时A ×B 才有意义。一个m ×n 的矩阵A 左乘一个n ×p 的矩阵B ,会得到一个m ×p 的矩阵C 。左乘:又称前乘,就是乘在左边(即乘号前),比如说,A 左乘E 即AE 。 一个m 行n 列的矩阵与一个n 行p 列的矩阵可以相乘,得到的结果是一个m 行p 列的矩阵,其中的第i 行第j 列位置上的数为第一个矩阵第i 行上的n 个数与第二个矩阵第j 列上的n 个数对应相乘后所得的n 个乘积之和。比如,下面的算式表示一个2行2列的矩阵乘以2行3列的矩阵,其结果是一个2行3列的矩阵。其中,结果矩阵的那个4(结果矩阵中第二(i )行第二(j)列)= 2(第一个矩阵第二(i)行第一列)*2(第二个矩阵中第一行第二(j)列) + 0(第一个矩阵第二(i)行第二列)*1(第二个矩阵中第二行第二(j)列): 矩阵乘法的两个重要性质:一,矩阵乘法满足结合律; 二,矩阵乘法不满足交换律。为什么矩阵乘法不满足交换律呢?这是由矩阵乘法定义决定的。因为矩阵AB=C ,C 的结果是由A 的行与B 的列相乘和的结果;而BA=D ,D 的结果是由B 的行与A 的列相乘和的结果。显然,得到的结果C 和D 不一定相等。同时,交换后两个矩阵有可能不能相乘。 因为矩阵乘法不满足交换律,所以矩阵乘法也不满足消去律。即由AB=AC 是得不到B=C 的,这是因为()AB AC A B C O =?-=是得不到A=O 或B-C=O 即B=C.例 111000010A B ????=≠=≠ ? ?-????0, 但0000AB O ??== ??? 那么由AB=O 一定得不到A=O 或B=O 吗?回答是否定的。比如A 是m ×n 阶矩阵,B 是n ×s 阶矩阵,若A 的秩为n ,则AB=O ,得B=O ;若B 的秩为m ,则AO ,得A=O.为什么吗?原因会在有关齐次线性方程组的文章里进行讲解.

ANSYS命令流解释大全

A N S Y S命令流解释大 全 Document serial number【LGGKGB-LGG98YT-LGGT8CB-LGUT-

一、定义材料号及特性 mp,lab, mat, co, c1,…….c4 lab: 待定义的特性项目(ex,alpx,reft,prxy,nuxy,gxy,mu,dens) ex: 弹性模量 nuxy: 小泊松比 alpx: 热膨胀系数 reft: 参考温度 reft: 参考温度 prxy: 主泊松比 gxy: 剪切模量 mu: 摩擦系数 dens: 质量密度 mat: 材料编号(缺省为当前材料号) c 材料特性值,或材料之特性,温度曲线中的常数项 c1-c4: 材料的特性-温度曲线中1次项,2次项,3次项,4次项的系数二、定义DP材料: 首先要定义EX和泊松比:MP,EX,MAT,…… MP,NUXY,MAT,…… 定义DP材料单元表(这里不考虑温度):TB,DP,MAT 进入单元表并编辑添加单元表:TBDATA,1,C TBDATA,2,ψ TBDATA,3,…… 如定义:EX=1E8,NUXY=,C=27,ψ=45的命令如下:

MP,EX,1,1E8 MP,NUXY,1, TB,DP,1 TBDATA,1,27 TBDATA,2,45这里要注意的是,在前处理的最初,要将角度单位转化到“度”,即命令:*afun,deg 三、单元生死载荷步 !第一个载荷步 TIME,... !设定时间值(静力分析选项) NLGEOM,ON !打开大位移效果 NROPT,FULL !设定牛顿-拉夫森选项 ESTIF,... !设定非缺省缩减因子(可选) ESEL,... !选择在本载荷步中将不激活的单元 EKILL,... !不激活选择的单元 ESEL,S,LIVE !选择所有活动单元 NSLE,S !选择所有活动结点 NSEL,INVE !选择所有非活动结点(不与活动单 元相连的结点) D,ALL,ALL,0 !约束所有不活动的结点自由度(可 选) NSEL,ALL !选择所有结点 ESEL,ALL !选择所有单元

ANSYS网格划分技巧

【分享】复杂几何模型的系列网格划分技术 众所周知,对于有限元分析来说,网格划分是其中最关键的一个步骤,网格划分的好坏直接影响到解算的精度和速度。在ANSYS中,大家知道,网格划分有三个步骤:定义单元属性(包括实常数)、在几何模型上定义网格属性、划分网格。在这里,我们仅对网格划分这个步骤所涉及到的一些问题,尤其是与复杂模型相关的一些问题作简要阐述。 一、自由网格划分 自由网格划分是自动化程度最高的网格划分技术之一,它在面上(平面、曲面)可以自动生成三角形或四边形网格,在体上自动生成四面体网格。通常情况下,可利用ANSYS的智能尺寸控制技术(SMARTSIZE命令)来自动控制网格的大小和疏密分布,也可进行人工设置网格的大小(AESIZE、LESIZE、KESIZE、ESIZE等系列命令)并控制疏密分布以及选择分网算法等(MOPT命令)。对于复杂几何模型而言,这种分网方法省时省力,但缺点是单元数量通常会很大,计算效率降低。同时,由于这种方法对于三维复杂模型只能生成四面体单元,为了获得较好的计算精度,建议采用二次四面体单元(92号单元)。如果选用的是六面体单元,则此方法自动将六面体单元退化为阶次一致的四面体单元,因此,最好不要选用线性的六面体单元(没有中间节点,比如45号单元),因为该单元退化后为线性的四面体单元,具有过刚的刚度,计算精度较差;如果选用二

次的六面体单元(比如95号单元),由于其是退化形式,节点数与其六面体原型单元一致,只是有多个节点在同一位置而已,因此,可以利用TCHG命令将模型中的退化形式的四面体单元变化为非退化的四面体单元,减少每个单元的节点数量,提高求解效率。在有些情况下,必须要用六面体单元的退化形式来进行自由网格划分,比如,在进行混合网格划分(后面详述)时,只有用六面体单元才能形成金字塔过渡单元。对于计算流体力学和考虑集肤效应的电磁场分析而言,自由网格划分中的层网格功能(由LESIZE命令的LAYER1和LAYER2域控制)是非常有用的。 二、映射网格划分     映射网格划分是对规整模型的一种规整网格划分方法,其原始概念是:对于面,只能是四边形面,网格划分数需在对边上保持一致,形成的单元全部为四边形;对于体,只能是六面体,对应线和面的网格划分数保持一致;形成的单元全部为六面体。在ANSYS中,这些条件有了很大的放宽,包括: 1 面可以是三角形、四边形、或其它任意多边形。对于四边以上的多边形,必须用LCCAT命令将某些边联成一条边,以使得对于网格划分而言,仍然是三角形或四边形;或者用AMAP命令定义3到4个顶点(程序自动将两个顶点之间的所有线段联成一条)来进行映射划分。 2 面上对边的网格划分数可以不同,但有一些限制条件。

矩阵相乘的快速算法

矩阵相乘的快速算法 算法介绍 矩阵相乘在进行3D变换的时候是经常用到的。在应用中常用矩阵相乘的定义算法对其进行计算。这个算法用到了大量的循环和相乘运算,这使得算法效率不高。而矩阵相乘的计算效率很大程度上的影响了整个程序的运行速度,所以对矩阵相乘算法进行一些改进是必要的。 这里要介绍的矩阵算法称为斯特拉森方法,它是由v.斯特拉森在1969年提出的一个方法。 我们先讨论二阶矩阵的计算方法。 对于二阶矩阵 A= a11a12 B= b11b12 a21a22b21 b22 先计算下面7个量(1) x1 = (a11 + a22) * (b11 + b22); x2 = (a21 + a22) * b11; x3 = a11 * (b12 - b22); x4 = a22 * (b21 - b11); x5 = (a11 + a12) * b22; x6 = (a21 - a11) * (b11 + b12); x7 = (a12 - a22) * (b21 + b22); 再设C = AB。根据矩阵相乘的规则,C的各元素为(2) c11 = a11 * b11 + a12 * b21 c12 = a11 * b12 + a12 * b22 c21 = a21 * b11 + a22 * b21 c22 = a21 * b12 + a22 * b22 比较(1)(2),C的各元素可以表示为(3) c11 = x1 + x4 - x5 + x7 c12 = x3 + x5 c21 = x2 + x4 c22 = x1 + x3 - x2 + x6 根据以上的方法,我们就可以计算4阶矩阵了,先将4阶矩阵A和B划分成四块 2阶矩阵,分别利用公式计算它们的乘积,再使用(1)(3)来计算出最后结果。

ANSYS中的APDL命令总结

在ANSYS中,命令流是由一条条ANSYS的命令组成的一个命令组合,这些命令按照一定顺序排布,能够完成一定的ANSYS功能,这些功能一般来说通过菜单操作也能够实现(而那些命令流能够实现,菜单操作实现不了的单个命令比较少见)。以下命令是结合我自身经验,和前辈们的一些经验而总结出来的,希望对大家有帮助。 (1).Lsel, type, item, comp, vmin, vmax, vinc, kswp选择线 type: s 从全部线中选一组线 r 从当前选中线中选一组线 a 再选一部线附加给当前选中组 au none u(unselect) inve: 反向选择 item: line 线号 loc坐标 length 线长 comp: x,y,z kswp: 0 只选线 1 选择线及相关关键点、节点和单元 (2).Nsel, type, item, comp, vmin, vmax, vinc, kabs选择一组节点 type: S: 选择一组新节点(缺省) R: 在当前组中再选择 A: 再选一组附加于当前组 U: 在当前组中不选一部分 All: 恢复为选中所有 None: 全不选 Inve: 反向选择 Stat: 显示当前选择状态 Item: loc: 坐标 node: 节点号 Comp: 分量 Vmin,vmax,vinc: ITEM范围 Kabs: “0”使用正负号 “1”仅用绝对值 (3).Esel, type, item, comp, vmin, vmax, vinc, kabs选择一组单元 type: S: 选择一组单元(缺省) R: 在当前组中再选一部分作为一组 A: 为当前组附加单元 U: 在当前组中不选一部分单元 All: 选所有单元 None: 全不选 Inve: 反向选择当前组 Stat: 显示当前选择状态 Item:Elem: 单元号 Type: 单元类型号

ANSYS 13.0 Workbench 网格划分及操作案例

第 3章 ANSYS 13.0 Workbench网格划分及操作案例 网格是计算机辅助工程(CAE)模拟过程中不可分割的一部分。网格直接影响到求解精 度、求解收敛性和求解速度。此外,建立网格模型所花费的时间往往是取得 CAE 解决方案所 耗费时间中的一个重要部分。因此,一个越好的自动化网格工具,越能得到好的解决方案。 3.1 ANSYS 13.0 Workbench 网格划分概述 ANSYS 13.0 提供了强大的自动化能力,通过实用智能的默认设置简化一个新几何体的网 格初始化,从而使得网格在第一次使用时就能生成。此外,变化参数可以得到即时更新的网 格。ANSYS 13.0 的网格技术提供了生成网格的灵活性,可以把正确的网格用于正确的地方, 并确保在物理模型上进行精确有效的数值模拟。 网格的节点和单元参与有限元求解,ANSYS 13.0在求解开始时会自动生成默认的网格。 可以通过预览网格,检查有限元模型是否满足要求,细化网格可以使结果更精确,但是会增 加 CPU 计算时间和需要更大的存储空间,因此需要权衡计算成本和细化网格之间的矛盾。在 理想情况下,我们所需要的网格密度是结果随着网格细化而收敛,但要注意:细化网格不能 弥补不准确的假设和错误的输入条件。 ANSYS 13.0 的网格技术通过 ANSYS Workbench的【Mesh】组件实现。作为下一代网格 划分平台, ANSYS 13.0 的网格技术集成 ANSYS 强大的前处理功能, 集成 ICEM CFD、 TGRID、 CFX-MESH、GAMBIT网格划分功能,并计划在 ANSYS 15.0 中完全整合。【Mesh】中可以根 据不同的物理场和求解器生成网格,物理场有流场、结构场和电磁场,流场求解可采用 【Fluent】、【CFX】、【POLYFLOW】,结构场求解可以采用显式动力算法和隐式算法。不同的 物理场对网格的要求不一样,通常流场的网格比结构场要细密得多,因此选择不同的物理场, 也会有不同的网格划分。【Mesh】组件在项目流程图中直接与其他 Workbench分析系统集成。 3.2 ANSYS 13.0 Workbench 网格划分 ANSYS 网格划分不能单独启动,只能在 Workbench 中调用分析系统或【Mesh】组件启 动,如图 3-1 所示。 图3-1 调入分析系统及网格划分组件

strassen矩阵相乘算法C++代码

Strassen 矩阵相乘算法代码 #include #include #include #include usingnamespace std; template class Strassen_class { public: void ADD(T** MatrixA, T** MatrixB, T** MatrixResult, int MatrixSize); void SUB(T** MatrixA, T** MatrixB, T** MatrixResult, int MatrixSize); void MUL(T** MatrixA, T** MatrixB, T** MatrixResult, int MatrixSize);//朴素算法实现void FillMatrix(T** MatrixA, T** MatrixB, int length);//A,B矩阵赋值 void PrintMatrix(T **MatrixA, int MatrixSize);//打印矩阵 void Strassen(int N, T **MatrixA, T **MatrixB, T **MatrixC);//Strassen算法实现 }; template void Strassen_class::ADD(T** MatrixA, T** MatrixB, T** MatrixResult, int MatrixSize) { for (int i = 0; i void Strassen_class::SUB(T** MatrixA, T** MatrixB, T** MatrixResult, int MatrixSize) { for (int i = 0; i void Strassen_class::MUL(T** MatrixA, T** MatrixB, T** MatrixResult, int MatrixSize) {

ansysworkbench中划分网格的几种方法

转自宋博士的博客 如何在ANSYS WORKBENCH中划分网格经常有朋友问到这个问题。我整理了一下,先给出第一个入门篇,说明最基本的划分思路。以后再对某些专题问题进行细致阐述。 ANSYS WORKBENCH中提供了对于网格划分的几种方法,为了便于说明问题,我们首先创建一个简单的模型,然后分别使用几种网格划分方法对之划分网格,从而考察各种划分方法的特点。 1. 创建一个网格划分系统。 2. 创建一个变截面轴。 先把一个直径为20mm的圆拉伸30mm成为一个圆柱体 再以上述圆柱体的右端面为基础,创建一个直径为26mm的圆,拉伸30mm得到第二个圆柱体。 对小圆柱的端面倒角2mm。

退出DM. 3.进入网格划分程序,并设定网格划分方法。双击mesh进入到网格划分程序。 下面分别考察各种网格划分方法的特点。(1)用扫掠网格划分。 对整个构件使用sweep方式划分网格。

结果失败。 该方法只能针对规则的形体(只有单一的源面和目标面)进行网格划分。(2)使用多域扫掠型网格划分。 结果如下

可见ANSYS把该构件自动分成了多个规则区域,而对每一个区域使用扫略网格划分,得到了很规则的六面体网格。这是最合适的网格划分方法。 (3)使用四面体网格划分方法。 使用四面体网格划分,且使用patch conforming算法。 可见,该方式得到的网格都是四面体网格。且在倒角处网格比较细密。 其内部单元如下图(这里剖开了一个截面) 使用四面体网格划分,但是使用patch independent算法。忽略细节。

?、网格划分结果如下图 此时得到的仍旧是四面体网格,但是倒角处并没有特别处理。(4)使用自动网格划分方法。 得到的结果如下图

矩阵连乘最佳加括号方式动态规划算法

矩阵连乘最佳加括号方式-动态规划算法 一、问题描述 给定n个矩阵{A1,A2,…,A n},其中A i与A i+1是可乘的,i=1,2,…,n-1。要算出这n个矩阵的连乘积A1A2…A n。由于矩阵乘法满足结合律,故计算矩阵的连乘积可以有许多不同的计算次序。这种计算次序可以用加括号的方式来确定。若一个矩阵连乘积的计算次序完全确定,也就是说该连乘积已完全加括号,则可以依此次序反复调用2个矩阵相乘的标准算法计算出矩阵连乘积。完全加括号的矩阵连乘积可递归地定义为: (1)单个矩阵是完全加括号的; (2)矩阵连乘积A是完全加括号的,则A可表示为2个完全加括号的矩阵连乘积B和C 的乘积并加括号,即A=(BC)。 例如,矩阵连乘积A1A2A3A4有5种不同的完全加括号的方式:(A1(A2(A3A4))),(A1((A2A3)A4)),((A1A2)(A3A4)),((A1(A2A3))A4),(((A1A2)A3)A4)。每一种完全加括号的方式对应于一个矩阵连乘积的计算次序,这决定着作乘积所需要的计算量。若A是一个p×q矩阵,B是一个q×r矩阵,则计算其乘积C=AB的标准算法中,需要进行pqr次数乘。 为了说明在计算矩阵连乘积时,加括号方式对整个计算量的影响,先考察3个矩阵 {A1,A2,A3}连乘的情况。设这三个矩阵的维数分别为10×100,100×5,5×50。加括号的方式只有两种:((A1A2)A3),(A1(A2A3)),第一种方式需要的数乘次数为10×100×5+10×5×50=7500,第二种方式需要的数乘次数为100×5×50+10×100×50=75000。第二种加括号方式的计算量时第一种方式计算量的10倍。由此可见,在计算矩阵连乘积时,加括号方式,即计算次序对计算量有很大的影响。于是,自然提出矩阵连乘积的最优计算次序问题,即对于给定的相继n个矩阵{A1,A2,…,A n}(其中矩阵A i的维数为p i-1×p i,i=1,2,…,n),如何确定计算矩阵连乘积A1A2…A n的计算次序(完全加括号方式),使得依此次序计算矩阵连乘积需要的数乘次数最少。 穷举搜索法的计算量太大,它不是一个有效的算法,本实验采用动态规划算法解矩阵连乘积的最优计算次序问题。 二、算法思路

ANSYS命令流学习笔记workbench中命令流的一些应用

!ANSYS命令流学习笔记5 workbench中命令流的一些应用 学习重点: 1. 定义单元类型 2. 使用各向异性材料时,定义其单元为圆柱坐标系 3. 有角度吊装时,定义吊装约束 4. workbench、APDL的联合仿真 案例如下: 如下图模型,四个顶点通过杆件连接,进行吊装时的有限元分析。 1. 建立模型,设定必要的坐标系。 分成两个solid,内部的圆柱,剩余的矩形部分。建立两个坐标系,分别用于指定各向异性材料的属性、吊装的固定点。 下图,建立圆柱坐标系,编号100,用于指定各向异性材料。 下图,建立直角坐标系,编号12,用于指定吊装固定点。 2. 建立named selection,方便在命令流中选择必要的元素。 下图,将四个吊装点,中间的圆柱,分别定义为任何名称,必须是英文才能用于APDL命令中。 3. 定义边界条件,施加重力加速度,在static structural 下插入command(APDL) ,内容如下 /prep7 alls cmsel,s,c1,elem !选择c1单元所有节点,既圆柱体的所有单元 emodif,all,esys,100 !其坐标系转换为100坐标系,因为缠绕的各向异性材料必须在圆柱坐标系下定义单元的坐标。 !完成对各向异性材料的坐标系设定。 et,10,10 !定义编号为10的,link10单元 r,10,0.01 !定义编号为10的实常数0.01,用于定义link10单元的截面积0.01mm^2 *get,nmax,node,,num,max !获取node的最大数值,储存在nmax名称的变量里 csys,12 !调用csys12坐标系 n,nmax+1, !csys坐标原点建立node,后面会将其固定,既吊装的固定点 mat,1 type,10 !选取编号10的单元类型 real,10 !选取编号10的实常数 cmsel,s,kk1,node !选择kk1点,kk1已经定义为named selection *get,k1,node,,num,max !获取已选节点中的节点数最大值,既kk1的节点编号,取值为k1 cmsel,s,kk2,node ! k2,k3,k4方法类型k1 *get,k2,node,,num,max cmsel,s,kk3,node *get,k3,node,,num,max cmsel,s,kk4,node *get,k4,node,,num,max

Ansys命令流大全(整理)

1、A,P1,P2,P3,P4,P5,P6,P7,P8,P9 此命令用已知的一组关键点点(P1~P9)来定义面(Area), 最少使用三个点才能围成面,同时产生转围绕些面的线。 点要依次序输入,输入的顺序会决定面的法线方向。 如果超过四个点,则这些点必须在同一个平面上。 Menu Paths:Main Menu>Preprocessor>Create>Arbitrary>Through KPs 2、*ABBR,Abbr,String--定义一个缩略语. Abbr:用来表示字符串"String"的缩略语,长度不超过8个字符. String:将由"Abbr"表示的字符串,长度不超过60个字符.3、ABBRES,Lab,Fname,Ext-从一个编码文件中读出缩略语. Lab:指定读操作的标题, NEW:用这些读出的缩略语重新取代当前的缩略语(默认)CHANGE:将读出的缩略语添加到当前缩略语阵列,并替代现存同名的缩略语. Ext:如果"Fname"是空的,则缺省的扩展命是"ABBR".4、ABBSA V,Lab,Fname,Ext-将当前的缩略语写入一个文本文件里 Lab:指定写操作的标题,若为ALL,表示将所有的缩略语都写入文件(默认) 5、add, ir, ia,ib,ic,name,--,--,facta, factb, factc 将ia,ib,ic变量相加赋给ir变量 ir, ia,ib,ic:变量号 name: 变量的名称 6、Adele,na1,na2,ninc,kswp !kswp=0时只删除掉面积本身,=1时低单元点一并删除。 7、Adrag, nl1,nl2,nl3,nl4,nl5,nl6, nlp1,nlp2,nlp3,nlp4,nlp5,nlp6 !面积的建立,沿某组线段路径,拉伸而成。 8、Afillt,na1,na2,rad !建立圆角面积,在两相交平面间产生曲面,rad为半径。 9、*AFUN,Lab 在参数表达式中,为角度函数指定单位. Lab:指定将要使用的角度单位.有3个选项. RAD:在角度函数的输入与输出中使用弧度单位(默认)DEG:在角度函数的输入与输出中使用度单位. STAT:显示该命令当前的设置(即是度还是弧度). 10、Agen, itime,na1,na2,ninc,dx,dy,dz,kinc,noelem,imove !面积复制命令。itime包含本身所复制的次数;na1,na2,ninc为现有的坐标系统下复制到其他位置(dx,dy,dz);kinc为每次复制时面积号码的增加量。 11、AINV, NA, NV 面与体相交生成一个相交面. NA, NV :分别为指定面,指定体的编号.其中NA可以为P.说明:面与体相交生成新面.如果相交的区域是线,则生成新线. 指定源实体的单元属性和边界条件不会转换到新生成的实 体上. 12、AL,L1,L2,L3,L4,L5,L6,L7,L8,L9,L10 此命令由已知的一组直线(L1,…L10)围绕成面(Area), 至少须要3条线才能形成面,线段的号码没有严格的顺序限制,只要它们能完成封闭的面积即可。 同时若使用超过4条线去定义面时,所有的线必须在同一平面上,以右手定则来决定面积的方向。如果L1为负号,则反向。Menu Paths:Main Menu>Preprocessor>Create>Arbitrary>By Lines 13、ALLSEL, LABT, ENTITY 选中所有项目 LABT: ALL: 选所有项目及其低级项目 BELOW: 选指定项目的直接下属及更低级项目 ENTITY: ALL: 所有项目(缺省) VOLU:体高级 AREA:面 LINE :线 KP:关键点 ELEM:单元 NODE:节点低级 14、Amesh, nA1,nA2,ninc 划分面单元网格nA1,nA2,ninc 待划分的面号,nA1如果是All,则对所有选中面划分 15、ANORM, ANUM, NOEFLIP 修改面的正法线方向. ANUM:面的编号,改变面的正法线方向与面的法线方向相同.NOEFLIP:确定是否要改变重定向面上单元的正法线方向,这样可以使他们与面的正法线方向一致 若为0,改变单元的正法线方向; 若为1,不改变已存在单元的正法线方向; 说明:重新改变面的方向使得他们与指定的正法线方向相同. 不能用"ANORM"命令改变具体或面载荷的任何单元的正法线方向. 16、数学函数 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, 函数值为0 NINT(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) 双曲正切 17、antype, status, ldstep, substep, action 声明分析类型,即欲进行哪种分析,系统默认为静力学分析。 antype: static or 1 静力分析 buckle or 2 屈曲分析 modal or 3 模态分析 trans or 4 瞬态分析 status: new 重新分析(缺省),以后各项将忽略 rest 再分析,仅对static,full transion 有效 ldstep: 指定从哪个荷载步开始继续分析,缺省为最大的,runn数(指分析点的最后一步) substep: 指定从哪个子步开始继续分析。缺省为本目录中,runn文件中最高的子步数 action, continue: 继续分析指定的ldstep,substep 说明:继续以前的分析(因某种原因中断)有两种类型singleframe restart: 从停止点继续 需要文件:jobname.db 必须在初始求解后马上存盘 jobname.emat 单元矩阵 jobname.esav 或.osav : 如果.esav坏了,将.osav 改为.esav results file: 不必要,但如果有,后继分析的结果也将很好地附加到它后面 注意:如果初始分析生成了.rdb, .ldhi, 或rnnn 文件。必须删除再做后继分析 步骤:(1)进入anasys 以同样工作名 (2)进入求解器,并恢复数据库 (3)antype, rest (4)指定附加的荷载 (5)指定是否使用现有的矩阵(jobname.trl)(缺省重新生成) kuse: 1 用现有矩阵

ansys_网格网格划分教程

血管模型网格划分 网格划分即将所用的模型划分为有限体积或单元,这里我们使用Ansys自身的网格划分器对提取出来的血管模型(STL格式的三维模型)进行网格划分。具体步骤为: 一、软件启动 单击开始---所有程序--Ansys14.0---Meshing---ICEM CFD 14.0。 二、模型导入 1、单击主菜单栏中的File---Import Geometry---STL,如下图: 主菜单烂

2、在下拉菜单对话框中选取血管模型并确定后,血管模型导入完成,如下图; 三、图形参数设置(封闭模型) 1、单击Geometry 工具栏中的Repair Geometry 图标,后单击Build Diagnostic 按钮单击Apply 按钮运行,如下图: 单击此处可调整显示边框与实体 Geometry 工具栏 Repair Geometry 按钮 单击此按钮后单击Apply 按钮 运行后端口显示封闭黄线

2、运行完成后,单击左侧工具栏中的Close Holes图标,然后单击鼠标图形按钮,而后单击模型端口处黄线并单击Apply运行,从而使模型端口封闭,模型有几个端口则反复操作几次,如下图: Close Holes按 钮运行结束 后端口封 闭 四、网格化分参数设置 1、单击Mesh工具栏的Compute Mesh按钮,选取Volume mesh按钮,后单击Compute按钮,而后单击YES按钮,进行初步的电脑网格划分; Mesh工具栏,Compute Mesh 按钮 单击此按钮,而后在下面的Mesh方法中选择,一般选择默 认的方法

运行结束后显示 网格 2、完成上步操作后,单击Mesh工具栏中的Globe mesh setup按钮,一般不更改默认设置,单击Apply按钮运行完成; 3、再次重复单击Compute Mesh按钮,后单击Apply按钮,而后单击YES按钮,进行初步的电脑网格划分,注:(1)可用多种方式进行划分而后单击融合操作,直至满意;(2)若模型有其他漏洞,程序会提示是否修复,一般选择不修复; 五、设置模型边界---共包括:出口端、入口端及墙壁 1、右键单击屏幕左侧控制树中的Part,在下拉菜单中选取Part create,而后在下面的窗口栏中Part部分对端口命名(如:input),然后在Creat Part by Selection部分单击鼠标箭头图形,再在浏览界面中左键单击所需设定的端口(如:输入端),注:只单击一次后进行下步操作,是否选中可能显示不明显,单击Apply运行完成此步操作; 2、依次完成输入端、输出端以及墙的设定。如下图:

AnsysWorkbench划分网格

Ansys Workbench 划分网格 (张栋zd0561@https://www.360docs.net/doc/dc3385188.html,) 1、对于三维几何体(对于三维几何体(3D 3D 3D) )有几种不同的网格化分方法。如图1下部所示。 图1网格划分的种类 1.1、Automatic(自动划分法) 1.2、Tetrahedron(四面体划分法) 它包括两种划分方法:Patch Conforming(A W 自带功能),Patch Independent(依靠ICEM CFD Tetra Algorithm 软件包来实现)。

步骤:Mesh(右键)——Insert——Method (操作区上方)Meshcontrl——Method (左下角)Scope——Geometry Method——Tetrahedrons(四面体网格) Algorithm——Patch Conforming (补充:Patch Independent该算法是基于Icem CFD Tetra的,Tetra部分具有膨胀应用,其对CAD许多面的修补均有用,包括碎面、短边、较差的面参数等。在没有载荷或命名选项的情况下,面和边无需考虑。) 图2四面体网格分两类

图3四面体划分法的参数设置 1.3、Hex Dominant(六面体主导法) 1.4、Sweep(扫掠划分法) 1.5、MultiZone(多区划分法) 2、对于面体或者壳二维几何 对于面体或壳二维(2D),A W有一下: Quad Dominant(四边形单元主导) Triangles(三角形单元) Uniform Quad/Tri(均匀四面体/三角形单元) Uniform Quad(均匀四边形单元) 3、网格参数设置 下图为缺省设置(Defaults)下的物理环境(Physics Preferance)

ANSYS软件APDL命令流建模的体会

ANSYS软件APDL命令流建模的体会ANSYS软件APDL命令流建模的体会首先申明,本人学习ANSYS基本上是靠自己一点一点琢磨出来的,由于本人喜欢用APDL命令流,故总结出来的几点经验也就比较适合用APDL命令的朋友。 1、多看help,ANSYS的help为我们提供了很强大的功能,我最喜欢的是其中对各个命令有关参数的说明和解释部分,不管是建模、加载、后处理等,都可以通过apdl命令来实现。只要你知道命令,如“aatt ”,在help搜索栏输入“aatt”,回车,弹出aatt的有关页码,一般其中有一个只有“aatt”的一项,确认,即可看到你要查询的aatt命令的有关参数意义,本人常用的命令有: et---定义单元类型 mp---定义材料属性 k----建关键点, l----建线条 a---由关键点建立面 al---由线建立面 v----由关键点建立体 vl---由线建立体 va--由面建立体 lsel---在很多很多线中选择你需要的目标线,数量可以无限多…… asel---在很多很多面中选择你需要的目标面,数量也可以无限多…… vsel---在很多很多体中选择你需要的目标体,数量也可以无限多…… latt----给选中的线按材料编号赋属性(前提是首先已定义好材料) aatt---给选中的面按材料编号赋属性

vatt-----给选中的体按材料编号赋属性 acel---按坐标轴赋体积力, lmesh,amesh,vmesh---对线、面、体进行剖分 d---在节点上加约束边界 dl---在线上加载约束边界 da----在面上加载约束边界 2、以上只是列出了常见的几个命令,但是ansys提供的命令是很多的,我们不可能都记得,计算记得,也不知道其有关参数是如何定义的,那不要紧,我们可以与界面操作结合起来学习。我们先利用界面操作实现,然后在保存路径里面找到文件“file.log”,在该文件里有该操作等价的apdl命令,那以后我们就可以使用了。 3、复合命令,很多命令是复合命令,通过几个命令的组合以实现一定的目标,如FITEM、FLST等。这里不予以详述,大家可在学习中慢慢体会。 4、ansys提供的apdl语言可像fortain、c语言一样,可以编程,有条件语句、逻辑语句、文件读写等,但是这些语句语法有个特点,就是在相应的语句前要加“*”,以示其与以上apdl命令的区别。 以上只是一点小小的总结,希望对大家有帮助。 K, NPT, X, Y, Z Defines a keypoint. Npt: Reference number for keypoint. If zero, the lowest available number is assigned X,y,z: Keypoint location in the active coordinate system (may be R, θ, Z or R, θ, Φ). If X = P, gra phical picking is enabled and all other fields (including NPT) are ignored (valid only in the GUI).

ANSYS最常用命令流+中文注释(超级大全)

ANSYS最常用命令流+中文注释 VSBV, NV1, NV2, SEPO, KEEP1, KEEP2 —Subtracts volumes from volumes,用于2个solid相减操作,最终目的是要nv1-nv2=?通过后面的参数设置,可以得到很多种情况:sepo项是2个体的边界情况,当缺省的时候,是表示2个体相减后,其边界是公用的,当为sepo的时候,表示相减后,2个体有各自的独立边界。keep1与keep2是询问相减后,保留哪个体?当第一个为keep时,保留nv1,都缺省的时候,操作结果最终只有一个体,比如:vsbv,1,2,sepo,,keep,表示执行1-2的操作,结果是保留体2,体1被删除,还有一个1-2的结果体,现在一共是2个体(即1-2与2),且都各自有自己的边界。如vsbv,1,2,,keep,,则为1-2后,剩下体1和体1-2,且2个体在边界处公用。同理,将v换成a 及l是对面和线进行减操作! mp,lab, mat, co, c1,…….c4 定义材料号及特性 lab: 待定义的特性项目(ex,alpx,reft,prxy,nuxy,gxy,mu,dens) ex: 弹性模量 nuxy: 小泊松比 alpx: 热膨胀系数 reft: 参考温度 reft: 参考温度 prxy: 主泊松比 gxy: 剪切模量 mu: 摩擦系数 dens: 质量密度 mat: 材料编号(缺省为当前材料号) co: 材料特性值,或材料之特性,温度曲线中的常数项 c1-c4: 材料的特性-温度曲线中1次项,2次项,3次项,4次项的系数 定义DP材料: 首先要定义EX和泊松比:MP,EX,MA T,…… MP,NUXY,MAT,…… 定义DP材料单元表(这里不考虑温度):TB,DP,MA T 进入单元表并编辑添加单元表:TBDATA,1,C TBDATA,2,ψ TBDATA,3,…… 如定义:EX=1E8,NUXY=0.3,C=27,ψ=45的命令如下:MP,EX,1,1E8 MP,NUXY,1,0.3 TB,DP,1 TBDATA,1,27 TBDATA,2,45这里要注意的是,在前处理的最初,要将角度单位转化到“度”,即命令:*afun,deg VSEL, Type, Item, Comp, VMIN, VMAX, VINC, KSWP Type,是选择的方式,有选择(s),补选(a),不选(u),全选(all)、反选(inv)等,其余方式不常用 Item, Comp 是选取的原则以及下面的子项 如volu 就是根据实体编号选择, loc 就是根据坐标选取,它的comp就可以是实体的某方向坐标! 其余还有材料类型、实常数等 MIN, VMAX, VINC,这个就不必说了吧! ,例:vsel,s,volu,,14 vsel,a,volu,,17,23,2 上面的命令选中了实体编号为14,17,19,21,23的五个实体 VDELE, NV1, NV2, NINC, KSWP: 删除未分网格的体 nv1:初始体号 nv2:最终的体号 ninc:体号之间的间隔 kswp=0:只删除体 kswp=1:删除体及组成关键点,线面 如果nv1=all,则nv2,ninc不起作用 其后面常常跟着一条显示命令VPLO,或aplo,nplo,这个湿没有参数的命令,输入后直接回车,就可以显示刚刚选择了的体、面或节点,很实用的哦! Nsel, type, item, comp, vmin, vmax, vinc, kabs 选择一组节点为下一步做准备 Type: S: 选择一组新节点(缺省) R: 在当前组中再选择 A: 再选一组附加于当前组 U: 在当前组中不选一部分 All: 恢复为选中所有 None: 全不选 Inve: 反向选择 Stat: 显示当前选择状态 Item: loc: 坐标 node: 节点号

ANSYS网格划分浅谈

ANSYS网格划分浅谈 在本学期,我们学习了CAX这门课程。通过对这一门课程8周的学习使我对本模块的认识和了解有了一种新的看法。在老师的认真教育和带领下把我们引入了一个新的领域。在CAX这个领域中包括CAD CAM CAE CAPP等的各项技术,这些技术都是将理论知识和计算机辅助集合在一起的新兴工业工程技术,是要将理论和实践的学科。在下面我主要将我这段学习期间对于ANSYS软件的学习中关于有限元网格划分的一些认识和经验做个报告总结。 1、ANSYS网格划分简述 ANSYS软件是融结构、流体、电场、磁场、声场分析于一体的大型通用有限元分析软件。由世界上最大的有限元分析软件公司之一的美国ANSYS开发,它能与多数CAD软件接口,实现数据的共享和交换,如Pro/Engineer, NASTRAN, Alogor, I-DEAS, AutoCAD等,是现代产品设计中的高级CAD工具之一。 ANSYS有限元网格划分是进行数值模拟分析至关重要的一步,它直接影响着后续数值计算分析结果的精确性。在划分网格前,用户首先需要对模型中将要用到的单元属性进行定义。在ANSYS中,大家知道,网格划分有三个步骤:定义单元属性(包括实常数)、在几何模型上定义网格属性、划分网格。单元属性主要包括:单元类型、实常数、材料常数。典型的实常数包括:厚度、横截面面积、高度、梁的惯性矩等。材料属性包括:弹性模量、泊松比、密度、热膨胀系数等。 ANSYS网格划分的指导思想是首先进行总体模型规划,包括物理模型的构造、单元类型的选择、网格密度的确定等多方面的内容。在网格划分和初步求解时,做到先简单后复杂,先粗后精,2D单元和3D单元合理搭配使用。为提高求解的效率要充分利用重复与对称等特征,由于工程结构一般具有重复对称或轴对称、镜象对称等特点,采用子结构或对称模型可以提高求解的效率和精度。利用轴对称或子结构时要注意场合,在进行模态分析、屈曲分析整体求解时,则应采用整体模型,同时选择合理的起点并设置合理的坐标系,可以提高求解的精度和效率,例如,轴对称场合多采用柱坐标系。有限元分析的精度和效率与单元的密度和几何形状有着密切的关系,按照相应的误差准则和网格疏密程度,避免网格的畸形。在网格重划分过程中常采用曲率控制、单元尺寸与数量控制、穿透控制等控制准则。在选用单元时要注意剪力自锁、沙漏和网格扭曲、不可压缩材料的体积自锁等问题。 ANSYS软件平台提供了网格映射划分和自由适应划分的策略。映射划分用于曲线、曲面、实体的网格划分方法,可使用三角形、四边形、四面体、五面体和六面体,通过指定单元边长、网格数量等参数对网格进行严格控制,映射划分只用于规则的几何图素,对于裁剪曲面或者空间自由曲面等复杂几何体则难以控制。自由网格划分用于空间自由曲面和复杂实体,采用三角形、四边形、四面体进行划分,采用网格数量、边长及曲率来控制网格的质量。

相关文档
最新文档