注册 登录  
 加关注
   显示下一条  |  关闭
温馨提示!由于新浪微博认证机制调整,您的新浪微博帐号绑定已过期,请重新绑定!立即重新绑定新浪微博》  |  关闭

htbbzzg的博客

与朋友们分享 CAE 资料和经验

 
 
 

日志

 
 

ANSYS 入门教程 (34) - 分析类型与求解控制选项 (c)  

2010-08-25 07:00:50|  分类: ANSYS 入门基础 |  标签: |举报 |字号 订阅

  下载LOFTER 我的照片书  |

4.5   重启动分析

    在完成一个初始分析过程之后,可能需要再次运行并继续分析,重启动则允许修改并继续一个分析,而不是从开始进行分析。
    重启动分析的目的有如下几个:
        ① 在线性静力分析中增加更多的荷载步继续分析,如另外施加荷载等。
        ② 在非线性分析中从收敛失败处恢复,改进收敛策略后继续分析。
        ③ 在瞬态分析中加入另外的时间-历程曲线,以研究不同的加载序列等。
    重启动分析的两种类型:
        ① 单点重启动分析:只能从初始分析的终止点开始继续计算。
        ② 多点重启动分析:能够在初始分析的任何点开始继续分析。
    重启动分析的前提条件:
        ① 分析类型必须是静态(稳态)或完全法瞬态分析。
        ② 重启动前的初始分析至少完成了一次迭代。
        ③ 初始分析的运行不能是系统中断、系统崩溃或强行杀死进程而终止的。
        ④ 初始分析和所生成的重启动分析文件必须是在同一版本下完成的。

 

一、  单点重启动分析

    单点重启动分析除要满足上述前提条件外,还必须存在初始分析所产生的下列文件:
        ① Jobname.DB:初始分析后保存的数据库文件。该文件在初始求解完成之后应立即保存,如果在以后的任何点保存数据库文件,边界条件和其它变量的初始值有可能被改变,从而导致重启动不能正常运行。
        ② Jobname.EMAT:单元矩阵文件。
        ③ Jobname.ESAV 或 Jobname.OSAV:单元数据(ESAV)或旧的单元数据(OSAV)。只有在 ESAV 文件丢失、损坏、或求解发散、位移超限、负主元等时才需要 OSAV 文件。在重启动分析前应将 OSAV 文件改为 ESAV 文件。
        ④ Jobname.RST:结果文件。该文件不是必须的,但如果存在该文件,将把重启动分析的结果用适当的荷载步和子步号追加到该文件中,以方便用户在后处理中的操作。

    同时应注意下列问题:
        ① 对因收敛失败、时间超限、异常中断等引起初始分析终止的情况,数据库会自动保存。
        ② 如果初始分析运行产生了 RDB、LDHI、Rnnn 文件,在单点重启动分析前应删除。
        ③ 交互式运行中会产生 DBB 文件,即 DB 文件的备份文件,但批处理方式不产生备份。
        ④ 如果子步数超过 1000(缺省),则程序将在该点终止,即刚刚重启动分析但子步数超限而终止;为解决这个问题可在重启动分析前,将初始分析的结果文件改名。

    在完成一个初始分析过程后,单点重启动分析的过程如下:
        ⑴  进入 ANSYS,用 /filname 命令将工作名指定为初始分析工作名。
        ⑵  进入求解层,用 RESUME 命令恢复数据库。
        ⑶  用 ANTYPE,,REST 设置为重启动分析。
        ⑷  按需要修改或施加另外的荷载,修改的斜坡荷载从以前的值开始,新施加的斜坡荷载从零开始,新施加的体荷载从初始值开始。删除荷载后重新施加按新施加处理,而不是修改荷载。
        ⑸  指定三角化矩阵是否重新使用,用 KUSE 命令设置。
        ⑹  用 SOLVE 命令进行求解。
        ⑺  对另外的荷载步重复第 ⑷ ⑸ ⑹ 步。当然也可采用荷载步文件进行求解。

    当数据库前后不一致时如何进行重启动分析呢?
    如在重启动分析之前进入了后处理,并执行了 SET 命令和 SAVE 命令,就可能导致数据库中的边界条件与重启动所需要的不一致。
    缺省情况下,ANSYS 在求解完成后,最后一个荷载步的边界条件是自动保存在内存中的。为在数据库不一致情况下进行重启动分析,可以首先进行一个“虚假”荷载步,以保证边界条件等的正确性。其方法是在上述第 ⑷ 步中,重新定义与初始分析的最后一子步相同的边界条件即可。

 1.   线性/非线性静力分析的重启动
    以图所示的悬臂梁为例,采用单点重启动分析的过程如下:

        !  EX4.17  线性静力分析的单点重启动分析
        finish  $ /clear  $ /filname,rtest1                                                ! 定义工作文件名为 RTEST1
        l0=1000  $ b0=10$h0=20  $ /prep7  $ k,1  $ k,2,l0  $ l,1,2       ! 定义参数,创建几何模型
        et,1,beam3  $ mp,ex,1,2.0e5  $ mp,prxy,1,0.3                          ! 定义单元类型及材料特性
        r,1,b0*h0,b0*h0*h0*h0/12,h0                                                    ! 定义实常数
        lesize,all,,,10  $ lmesh,all  $ d,1,all                                            ! 划分单元,施加约束
        /solu  $ antype,0                                                                       ! 进入求解层,定义静力分析类型
        !NLGEOM,ON                                                                           ! 关闭则为线性分析,打开则为非线性
        outres,all,all  $ autots,off                                                           ! 定义输出选项,关闭自动时间步
        nsubst,10  $ f,2,fy,2000  $ solve                                               ! 定义子步数、施加荷载,求解
        save                                                                                          ! 保存数据库文件,即 RTEST1.DB

                                ANSYS 入门教程 (34) - 分析类型与求解控制选项 (c) - htbbzzg - htbbzzg的博客

 

                                ANSYS 入门教程 (34) - 分析类型与求解控制选项 (c) - htbbzzg - htbbzzg的博客

 

        finish  $ /clear  $ /filname,rtest1                 ! 重新开始一个工作,并指定同名工作名
        /solu  $ resume                                          ! 进入求解层,恢复数据库
        antype,,rest                                               ! 定义分析类型为重启动分析
        f,7,fy,4440  $ nsubst,20                             ! 施加荷载,定义子步数
        solve                                                          ! 求解
        finish  $ /post26                                         ! 进入时间后处理,查看结果
        nsol,2,2,u,y  $ rforce,3,1,f,y
        prod,4,3,,,,,,-1  $ xvar,2  $ plvar,4  $ prvar,2,4
    从上述分析可以看出,重启动分析是“后续荷载步在前步的基础上计算”。

                ANSYS 入门教程 (34) - 分析类型与求解控制选项 (c) - htbbzzg - htbbzzg的博客

 

                ANSYS 入门教程 (34) - 分析类型与求解控制选项 (c) - htbbzzg - htbbzzg的博客

 

2. 重启动分析中的生死单元示例
    以平面框架为例,说明有生死单元时使用重启动的过程.
        !  ex4.18 A 有生死单元时的重启动分析
        finish  $ /clear  $ /filname,frame1
        !  初始分析:建模、划分单元、求解等
        /prep7  $ l0=1000  $ d0=10  $ et,1,beam3  $ mp,ex,1,2.0e5  $ mp,prxy,1,0.3
        a0=acos(-1)*d0**2/4  $ i0=acos(-1)*d0**4/64  $ r,1,a0,i0,d0  $ k,1  $ k,2,,l0  $ k,3,l0,l0
        k,4,l0$l,1,2  $ l,2,3  $ l,3,4  $ lesize,all,,,10  $ lmesh,all  $ dk,1,all  $ dk,4,all
        /solu  $ antype,0  $ outres,all,all  $ autots,off  $ nropt,full  $ time,1 $ nsubst,10
        f,17,fy,-20  $ lsel,s,loc,x,l0  $ esll,s  $ cm,e1cm,elem  $ ekill,all  $ allsel
        !  求解并立即保存初始分析的数据库
        solve   $ save

 

                ANSYS 入门教程 (34) - 分析类型与求解控制选项 (c) - htbbzzg - htbbzzg的博客

 

                ANSYS 入门教程 (34) - 分析类型与求解控制选项 (c) - htbbzzg - htbbzzg的博客

 

        !  重新开始一个工作,以说明重启动分析过程
        finish  $ /clear  $ /filname,frame1                         ! 注意工作名与初始分析相同
        /solu  $ resume  $ antype,,rest                            ! 进入求解层、恢复数据库、设置重启动分析
        time,2  $ cmsel,s,e1cm  $ ealive,all                     ! 激活原来杀死的单元
        allsel  $ f,16,fy,-44  $ nsubst,20                          ! 施加新的荷载及求解选项
        solve  $ finish  $ /post26                                     ! 求解完毕后进入时间后处理,考察结果
        nsol,2,17,u,y  $ rforce,3,1,f,y  $ rforce,4,22,f,y  $ add,5,3,4
        prod,6,2,,,,,,-1  $ xvar,6  $ plvar,5  $ prvar,2,3,4,5,6

                ANSYS 入门教程 (34) - 分析类型与求解控制选项 (c) - htbbzzg - htbbzzg的博客

 

                ANSYS 入门教程 (34) - 分析类型与求解控制选项 (c) - htbbzzg - htbbzzg的博客

 

                ANSYS 入门教程 (34) - 分析类型与求解控制选项 (c) - htbbzzg - htbbzzg的博客

 

        !  ex4.18B  正常生死单元分析
        finish  $ /clear  $ /prep7  $ l0=1000  $ d0=10$et,1,beam3  $ mp,ex,1,2.0e5
        mp,prxy,1,0.3
        a0=acos(-1)*d0**2/4  $ i0=acos(-1)*d0**4/64  $ r,1,a0,i0,d0
        k,1  $ k,2,,l0  $ k,3,l0,l0  $ k,4,l0  $ l,1,2  $ l,2,3  $ l,3,4  $ lesize,all,,,10  $ lmesh,all
        dk,1,all  $ dk,4,all
        /solu  $ antype,0  $ outres,all,all  $ autots,off  $ nropt,full  $ time,1  $ nsubst,10
        f,17,fy,-20  $ lsel,s,loc,x,l0  $ esll,s  $ cm,e1cm,elem  $ ekill,all  $ allsel
        solve  $ time,2  $ cmsel,s,e1cm  $ ealive,all  $ allsel
        f,16,fy,-44  $ nsubst,20  $ solve  $ finish
        /post26  $ nsol,2,17,u,y  $ rforce,3,1,f,y  $ rforce,4,22,f,y  $ add,5,3,4
        prod,6,2,,,,,,-1
        xvar,6  $ plvar,5  $ prvar,2,3,4,5,6
    通过分析可以看出,重启动分析实质上是在原来分析的某个基础上继续计算,与常规连续计算是一样的,不应赋予重启动分析太多的其它作用(如下文的时变结构及其处理)。当然,重启动分析的其它目的(如非线性分析中继续收敛计算)等还是可以达到的。

  评论这张
 
阅读(2023)| 评论(0)
推荐 转载

历史上的今天

在LOFTER的更多文章

评论

<#--最新日志,群博日志--> <#--推荐日志--> <#--引用记录--> <#--博主推荐--> <#--随机阅读--> <#--首页推荐--> <#--历史上的今天--> <#--被推荐日志--> <#--上一篇,下一篇--> <#-- 热度 --> <#-- 网易新闻广告 --> <#--右边模块结构--> <#--评论模块结构--> <#--引用模块结构--> <#--博主发起的投票-->
 
 
 
 
 
 
 
 
 
 
 
 
 
 

页脚

网易公司版权所有 ©1997-2016