《【2017年整理】ansys加载与求解》由会员分享，可在线阅读，更多相关《【2017年整理】ansys加载与求解（11页珍藏版）》请在金锄头文库上搜索。
1、广州工程仿真科技有限公司 工程仿真网CAE加载与求解施加载荷是有限元分析中关键的一步，可以对网格划分之后的有限元模型施加载荷，也可以直接对实体模型施加载荷。当对模型进行了划分网格和施加载荷之后，就可以选择适当的求解器对问题进行求解。第一节基础知识一、载荷的分类ANSYS 中载荷(Loads)包括边界条件和模型内部或外部的作用力。在不同的学科中，载荷的定义如下。结构分析：位移、力、压力、弯矩、温度和重力。热分析：温度、热流率、对流、内部热生成、无限远面。磁场分析：磁势、磁流通、磁电流段、源电密度、无限远面。电场分析；电势(电压)、电流、电荷、电荷密度、无限远面。流场分析：速度、压力。在 ANSY
2、S 中，载荷主要分为六大类：DOF 约束(自由度约束)、力(集中载荷)、表面载荷、体载荷、惯性力及耦合场载荷，它们的含义为如下。DOF 约束(DOF constraint)：用户指定某个自由度为已知值。在结构分析中约束是位移和对称边界条件：在热力学分析中约束是温度和热流量等。力(集中载荷)(Fome)：施加于模型节点的集中载荷。如结构分析中的力和力矩，热分析中的热流率。表面载荷(SurfaceLoad)：作用在某个表面上的分布载荷。如结构分析中的压力，热分析中的对流和热流量。体载荷(Body loads)：作用在体积或场域内。如结构分析中的温度和重力，热分析中的热生成率。惯性载荷(Inerti
3、a loads)：结构质量或惯性引起的载荷。如重力加速度、角速度和角加速度，主要在结构分析中使用。耦合场载荷(Coupled-field loads)：它是一种特殊的情况，从一种分析中得到的结果用作另一种分析的载荷；如热分析中得到的节点温度可作为结构分析中的体载荷施加到每一个节点。二、载荷步、子步和平衡迭代1载荷步载荷步是指分步施加的载荷，在线性静态或稳态分析中，可以使用不同的载荷步施加不第 82 页广州工程仿真科技有限公司 工程仿真网CAE同的载荷组合。如图 51 所示，第一个载荷步用于线性载荷，第二个载荷步用于常数载荷部分，第二个载荷步用于卸载。2子步子步是指在一个特定的载荷步中每一次增加
4、的步长，也称为时间步，代表一段时间。对于不同的分析类型，子步的作用不同：在非线性静态分析或稳态分析中，使用子步逐渐施加载荷以便能获得精确解。在线性或非线性瞬态分析或稳态分析中，使用于步满足瞬态时间积分法则(为获得精确解，通常规定一个最小的时间步长)。在谐波分析中，使用于步可获得谐波频率范围内多个频率处的解。3阶跃载荷和坡度载荷在一个载荷步中，有两个或者两个以上的载荷步子步时，就必须选择所施加的载荷应该为阶跃载荷还是为坡度载荷。所谓阶跃载荷，就是指在第一个子步全部施加上去了，载荷在以后的每个子步中保持不变。坡度载荷就是指在每一个载荷步子步，载荷值都是递增的，直到最后一个载荷步子步，全部的载荷才施
5、加上去。4平衡迭代平衡迭代是指在给定子步下为了收敛而计算的附加解。平衡迭代仅应用于收敛起着很重要作用的非线性分析(静态或瞬态)中的迭代修正。三、通用选项通用选项包括瞬态或者静态分析当中载荷步结束的时间、子步步数或者说时间步大小、阶跃载荷、热应力当中的参考温度。选择 Main MenuSolutionLoad Step Opts 命令展开载荷步选项菜单。选择 Main MenuSolutionLoad Step OptsTimeFrequentTime-Time Step 命令，弹出如图 5-2 所示的对话框。注意：如果展开的载荷步选项菜单不完全，选择 Main MenuSolutionUnab
6、ridgad Menu 命令即可TIME 命令在与速率有关的问题当中是指实际的时间，要求指定一个时间值；在与速率无关的问题里面，时间是一个用作跟踪载荷的参数。显然，无论哪一种情况，都不能将时间设置为 0。DELTIM 命令是给 ANSYS 程序分析指定时间步的大小，在通用选项的另外一个窗口 Time and Substep Options 当中，ANSYS 程序是要通过 NSUBST 命令来指定分析过程当中的子步的大小。广州工程仿真科技有限公司 工程仿真网CAEKBC 命令是指定载荷的施加是采用阶跃式还是采用坡度式(线性方式)。概述迄今为止, 我们已经知道了如何施加以下类型的载荷:位移(DOF
7、 约束)压力和对流载荷(表面载荷)重力(惯性载荷)“结构”温度(体载荷)这些载荷占了五种载荷类型中的 4 种。本文将讲述剩下的一种载荷集中载荷, 比如应力分析中的节点载荷。将就以下问题进行讨论:A. 集中载荷B. 节点坐标C. 求解器D. 多重载荷步A. 力载荷一个力就是可以在一个节点或关键点处施加的集中载荷(也可以叫“点载荷”)和力一样，点载荷适合于线状模型，如梁，桁架，弹簧等。在实体单元或壳单元中, 点载荷往往引起应力奇异，但当您忽略了附近的应力时，它仍然是可接受的。记住，您可以通过选择来忽略附近施加了点载荷的单元。广州工程仿真科技有限公司 工程仿真网CAE在左下角展示的二维实体单元中，我
8、们注意到在加力位置出现最大应力 SMAX (23,854)。当在力附近的节点和单元不被选中时，SMAX (12,755) 就会移到底部角点处，这是由于在该角点处约束引起的另一处应力奇异。广州工程仿真科技有限公司 工程仿真网CAE通过不选底部角点附近的节点和单元，您就可以在上孔附近得到预期的应力 SMAX (8,098)。注意，对于轴对称模型:在全部 360范围内输入力的值。同样在全部 360范围内输出力的值(反力)。例如, 设想一个半径为 r 的圆柱形壳体边缘施加有 P lb/in 的载荷。把这个载荷施加在二维轴对称壳体模形上(比如 SHELL51 单元), 您就要施加一 2rP 的力。广州工
9、程仿真科技有限公司 工程仿真网CAE施加一个力需要有以下信息:节点号(您可以通过施取确定)力的大小(单位应与您正在使用的单位系统保持一致)力的方向FX, FY, 或 FZ使用:Solution -Loads-Apply Force/Moment或命令 FK 或 F问题: 在哪一个坐标系中 FX, FY, 和 FZ 有说明？B.节点坐标系所有的力，位移，和其它与方向有关的节点量都可以在节点坐标中说明。输入量:力和力矩 FX, FY, FZ, MX, MY, MZ位移约束 UX, UY, UZ, ROTX, ROTY, ROTZ耦合和约束方程其它输出量:计算出的位移 UX, UY, UZ, ROT
10、X, ROTY, ROTZ反力 FX, FY, FZ, MX, MY, MZ其它节点坐标系和模型中的每个节点有关。缺省时, 节点坐标系与总体笛卡尔坐标系一致，例如，所有施加的力和位移约束缺省时都是在笛卡尔坐标中。广州工程仿真科技有限公司 工程仿真网CAE必要时，您可以转换节点坐标的方向。例如:模拟一个斜的滚动支座。施加径向力.施加径向约束(比如模拟一个承受压力的刚性销)。旋转节点坐标分为四步:1.选择需要旋转的节点。2.激活您要旋转的节点所在的坐标系(或者生成一个局部坐标系)例如：CSYS,1。3.选择 Preprocessor Move/Modify -Rotate Node CS-To A
11、ctive CS, 然后在拾广州工程仿真科技有限公司 工程仿真网CAE取器中按Pick All或者使用 NROTAT,ALL 命令。4.激活所有节点。注意: 当您在反对称边界条件下施加对称约束时，ANSYS 自动旋转边界上的所有节点。演示:恢复 rib.db 文件。把工作平面移至底圆圆心处(使用关键点的中间位置)。在工作平面上激活柱坐标系（或建立一个局部坐标系）。选择半径在 r = 0.35 上的点并画出它们。在当前激活的坐标系中旋转所有被选择节点的坐标系。在所有选择的节点上施加一个 UX 位移约束(或者施加一个 FX 的力)。注意极径方向。激活笛卡尔坐标系(CSYS,0).在当前激活的坐标系
12、中旋转所有选择节点的坐标系。重新画出节点，注意新载荷的方向。C. 求解器求解器的功能是求解代表结构自由度个数的线性联立方程。求解的速度主要取决于模型的大小和您计算机的速度，所用时间可以是几秒，也可以是几小时。广州工程仿真科技有限公司 工程仿真网CAE只有一个载荷步的线性静态分析只需一次求解，而非线性或瞬态分析可能需要几十个，几百个甚至几千次求解。因此，选择求解器的类型是很重要的。ANSYS 中可用的求解器可以分为三类:直接消去求解器波前求解器稀疏求解器(缺省)迭代求解器PCG (预制条件共轭梯度求解器)ICCG (不完全乔利斯基共轭梯度求解器)JCG (雅可比共轭梯度求解器)并行求解器(需要特
13、殊的授权文件)AMG (Algebraic Multigrid)DDS (分布区域求解器)直接消去求解器求解以下内容:1.计算单元刚度矩阵。2.读取第一个单元的自由度。3.删除所有已知自由度或通过其它自由度可以表示的自由度, 然后把方程写入.tri 文件。保留的自由度构成波前。4.对所有单元重复 2，3 步骤直到所有的自由度都被消去。现在的.tri 就包含了一个三角化的矩阵。5.回代求解自由度,然后使用单元矩阵计算单元解。广州工程仿真科技有限公司 工程仿真网CAE由于自由度仍不能被排除，所以波前是对角化过程中保留在求解器中的自由度数。它随着求解的进行增大或缩小，当最后自由度消去完成后，波前数变为零。波前直接影响到计算的速度：值越大，速度越慢。重新进行单元编号为求解器选择一个合适的单元编号顺序可以减小波前数。ANSYS 在一开始就自动对单元重新编号。迭代求解器通过以下步骤求解:1.计算单元刚度矩阵。2.组集总体刚度矩阵。3.开始时把所有自由度的值设为零，然后一直迭代到收敛(基于输入的残余力的容许值)。4.用单元刚度矩阵计算单元解。在 ANSYS 中迭代求解器和 PCG, JCG, ICCG 的主要区别是所使用的预条件控制不同。广州工程仿真科技有限公司 工程仿真网CAE
