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

htbbzzg的博客

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

 
 
 

日志

 
 

NX Nastran 超单元用户指南 一份翻译资料 (34)  

2012-05-27 08:03:36|  分类: Nastran 超单元 |  标签: |举报 |字号 订阅

  下载LOFTER 我的照片书  |
第 10 章 动力减缩的输入和输出

■ 动力减缩的情况控制

■ 单级动力减缩

■ 多级动力减缩

 

本章说明动力超单元的处理过程,以及与动力减缩有关的输出。

10.1 动力减缩的情况控制

    上一章说明了动力减缩的处理过程。现在定义控制这一过程的输入。本章的第一部分处理情况控制,第二部分处理模型数据。为了清楚,模型数据段分为两类 – 一类是在主模型数据段定义超单元,另一类是使用 PARTs 超单元。对于 PARTs 或主模型数据超单元,情况控制命令是一样的。 

            作者注:只在必要时进行动力减缩。提供情况控制命令和模型数据卡。

    当存在必须的情况控制命令和模型数据卡时,调用动力减缩。如果缺少任何一致信息,将执行静力减缩。

 

动力减缩的情况控制:

    指南的这一部分讨论动力减缩的情况控制和对于超单元情况控制是如何工作的。因此,如果你还没有看过第 4 章 “静力分析的载荷、约束和情况控制”,建议你先看一下。

    动力分析的情况控制是对静力分析的情况控制的简单扩展。事实上,如果对所有超单元执行静力减缩 (默认),对动力问题不必考虑任何新的情况控制卡。但是,如果要进行动力减缩,对几乎所有情况,必须使用扩展的情况控制,只有一个例外 – 对所有超单元和残余结构使用同一特征值 (或动力减缩) 求解控制。这一例外很少发生,除了一些小的简单问题。

             作者注:控制超单元动力减缩的命令必须出现在代表该超单元的第一个 SUBCASE 卡中。

 

    各超单元的动力减缩方法由该超单元的第一个 SUBCASE 卡中的命令定义。METHOD 命令控制动力过程。

    如果要对超单元进行部件模态综合,必须有一个 METHOD 命令指向模型数据部分的 EIGR 或 EIGRL 卡。如果求解装配模态 (只需要 SOL 103 或模态解),则在情况控制段,每一个残余结构子情况中必须有一个 METHOD 命令。由于在 SOL 103 中允许引入多种边界条件,也可以求解多个多种残余结构的特征值问题。如果你利用这一特性,在每个残余结构的 SUBCASE 中,BC 必须有不同的值 (对数据库中的结果进行区分)。

下面是情况控制段的一个例子,它说明对简单模型的不同动力减缩方法:

 

TITLE = S.E. SAMPLE PROBLEM - NORMAL MODES

SUBTITLE = S.E. DYNAMICS - USE DIFFERENT REDUCTION TECHNIQUES

DISP = ALL

PARAM,GRDPNT,0

$

SUBCASE 1

SET 12 = 1,2

SUPER = 12

LABEL = COMPONENT MODE REDUCTION FOR S.E. 1 & 2

METHOD = 1

$

SUBCASE 2

SET 13 = 3,4

SUPER = 13

LABEL = STATIC REDUCTION FOR S.E. 3 & 4

$

SUBCASE 3

SET 15 = 5,6

SUPER = 15

LABEL = COMPONENT MODE REDUCTION FOR S.E. 5 & 6

METHOD = 1

$

SUBCASE 4

SUPER = 7

LABEL = STATIC REDUCTION FOR S.E. 7

$

SUBCASE 5

LABEL = RES STR - SOLVE PROBLEM FOR SYSTEM MODES

SUPER = 0

METHOD = 2

$

 

如果输入文件有正确的模型数据卡,则以上情况控制段需要如下内容:

    SUBCASE 1  (用于超单元 1 和 2)。这一 SUBCASE 要求采用特征值 METHOD 1.进行部件模态综合;

    SUBCASE 2  (用于超单元 3 和 4) 。由于没有要求动力减缩的命令,对这些超单元只进行静力减缩

    SUBCASE 3  (用于超单元 5 和 6)。这一 SUBCASE 告诉程序对这些超单元执行部件模态综合。

    SUBCASE 4  (用于超单元 7)。这一 SUBCASE 是执行部件模态综合的最常用格式。对这一子情况,要求程序根据 EIGRL (或 EIGE) 卡的说明计算超单元 7 的模态。

    SUBCASE 5  (用于残余结构)。这一 SUBCASE 告诉程序使用特种咱 METHOD 2 求解残余结构的装配模态。

 

10.2 单级动力减缩

    单级和多级动力减缩的模型数据卡是相似的。这些卡显示在不同的部分以强调方法的差别。如以前提到的,情况控制与所用的树的类型 (单级或多级) 无关。

    在单级动力分析中,所有超单元直接附着在残余结构上。因此,在系统求解时,部件模态直接传送给残余结构。对多级动力分析,可能要创建装配,在传送给残余结构之前,上游超单元的部件模态传送给其下游的超单元,以得到下游超单元的装配模态。

    在单级树中,代表部件模态的自由度必须是残余结构的内部自由度,对超单元是外部的,使用如下方法进行。

 

使用主模型数据超单元的单级动力减缩的模型数据卡:

为了在一个超单元上进行动力减缩,除了要求的情况控制命令,需要定义代表动力变换矢量的自由度 (部件模态)。还需要定义在求解特征值时如何处理外部自由度 (默认为固定边界)。

在进行超单元动力减缩时,在主模型数据卡中使用如下模型数据卡:

SEQSET 和/或 SEQSET1:可以使用其中一种或两种卡,定义在后续处理中代表超单元动力变换矢量 (部件模态) 的自由度。这些卡指向 SPOINTs 和/或节点,它们是残余结构的内部自由度。如果只用主模型数据卡 定义模型,可以使用节点。如果存在任何 PARTs,只能用 SPOINTs。

     作者注:对于列在 SESUP 卡上的每一个自由度,舍弃一个动力变换矢量 (从最低的自然频率开始)。由于程序不检查

                      舍弃的矢量是否静力变换矢量的线形组合,不推荐这种方法。

 SEBSET 和/或 SEBSET1:置于超单元 B 集中的外部自由度表,在计算动力变换时,这些自由度是约束的。(默认 = 所有外部自由度)。

SECSET 和/或 SECSET1:置于超单元 C 集中的外部自由度表,在计算动力变换时,这些自由度不约束。(默认 = 没有自由度)。

SESUP:置于超单元 R 集中的外部自由度表,在计算动力变换时,这些自由度不约束。但是,对应 SESUP 卡上的每一个自由度,将舍弃一个动力变换矢量 (从第一个模态开始),但不作任何检查。

SENQSET:对于所选的或全部超单元,自动生成 Q 集自由度的 (超单元) 编号表。这张卡只能用于模型中游一个或多个 PARTs 的情况。如果使用这张卡,会对所选的超单元,自动生成 Q 集自由度。这些自由度传送给残余结构,不与下游超单元耦合,直到处理残余结构。因此,对于多级超单元动力分析,这张卡可能不是最好的选择。

GRID 和/或 SPOINT:必须为NX Nastran 定义在部件模态中使用的自由度 (注意:如果使用了 SENQSET,这一步不需要)。

     注:对主模型数据超单元,部件模态必须附着在节点或标量点上 (如存在 PARTs,则节点不能用;SENQSET 卡可以用于主模型数据超单元和 PARTs 两种情况)。

    为将部件模态附着到节点或标量点上,使用 SEQSET 和/或 SEQSET1 卡。使用这些卡时,需要指明对应的超单元并选择代表部件模态的节点自由度或标量点。确保所选择的自由度数足以代表所计算的所有部件模态,否则,无代表的模态将不包含在解中。在单级过程中,代表部件模态的 GRID 自由度和 SPOINTs 必须是时残余结构的内部自由度。

    如果用 GRID 点代表超单元的部件模态,这些节点必须不附着在结构单元上。如果他们附着在结构单元上,将会遇到致命错误信息:User Fatal Message 5290,指明有一个单元附着在 Q 集自由度上:

*** USER FATAL MESSAGE 5290 (GPSP)

EXTERNAL ID = 1001 DOF = 1 IS IN THE Q-SET, BUT HAS A STIFFNESS TERM BEFORE REDUCTION.

^^^ USER FATAL MESSAGE 9032 (ERRPH1) - RUN IS TERMINATED DUE TO MESSAGE(S) ABOVE. TO CONTINUE

PROCESSING ALL SUPERELEMENTS,

INSERT PARAM,ERROR,0 INTO BULK DATA.

上一信息说明:节点 1001 处于当前超单元的 Q 集中,在计算部件模态综合前有一个刚度项,计算终止。注意:此信息提到的 PARAM,ERROR 可以用于继续计算。对 PARAM,ERROR 给以 0 值,告诉程序当出现致命错时仅终止当前超单元的计算。程序将对其它超单元继续计算,直到不能进一步计算为止。如果使用这一参数,必须仔细查看输出文件,以确信没有遇到致命错。

         作者注:如果一个模型仅定义了主模型数据,则 SPOINTs (按照定义) 是残余结构的内部自由度。

 

  如果一个模型仅定义了主模型数据,则很容易用 SPOINTs 代表部件模态。当一个模型仅定义了主模型数据时,按照定义,所有 SPOINTs 都是残余结构的内部自由度,而不能是任何其它超单元的内部自由度,这与有 PARTs 的模型的行为是不同的。在有 PARTs 的模型中,SPOINTs 可以是任何超单元的内部自由度。

    下两行显示如何用 SPOINTs 定义一个超单元的 Q 集自由度:

               SPOINT,1001,thru,1010

               SEQSET1,1,0,1001,thru,1010

上述代码建立 10 个 SPOINTs (1001 到 1010),并告诉程序用这些 SPOINTs 表示超单元 10 的头 10  个模态。

 

使用主模型数据的单级超单元的例子:

    首先用一个悬臂梁的例子显示不同方法的输入。这一模型也用于第9 章 “在动力分析中引入超单元” 中生成比较不同减缩方法结果的表格。在第9 章 “在动力分析中引入超单元” 中,在超单元级采用不同的计算模态数,对每一种动力减缩方法进行了几次计算以进行比较。所有这些计算都是用 SOL 103 (SEMODES) 进行的。

这里将使用的方法是:

* 只用静力计算 (seg10_a.dat)

* 固定边界 CMS (seg10_b1.dat)

* 自由-自由 CMS (seg10_c1.dat)

* 混合边界 CMS (seg10_d1.dat)

以下列出本例的模型数据,这一模型是图 9-1 所示的悬臂梁。用 SESET 卡将此模型分为两个超单元 (1 和 2)。超单元 1 有内部点 7  到 11,超单元 2 有内部点 1 到 5。两个超单元都以节点 6 作为它们的外部点。

        (文件 cantbeam.dat  等,见下一日志)

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

历史上的今天

在LOFTER的更多文章

评论

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

页脚

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