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

htbbzzg的博客

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

 
 
 

日志

 
 

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

2012-02-02 07:40:36|  分类: Nastran 超单元 |  标签: |举报 |字号 订阅

  下载LOFTER 我的照片书  |

表 5-3   输入文件 inrel3.dat

      $

      $ file - inrel3.dat

      $

      $ all 7 s.e. in main bulk data section

      $

      $ inertia relief - note PARAM,INREL,-1

      $ GRDSET removed and K6ROT is set in case control

      $ constraints removed on grids 1 & 2

      $

      $ ----------------------------------------------------------------

      $

      SOL 101

      TIME 15

      CEND

      TITLE = File inrel1.dat - inertia relief using main bulk data

      SUBTITLE = reference point is point 13 in residual

      DISP = ALL

      stress = all

      SET 999 = 0,1,2,3,4,5,6,7

      SUPER = 999 $ ALL CASE CONTROL IS FOR ALL SUPERELEMENTS

      PARAM,K6ROT,1000.

      SUBCASE 101

      LABEL = PRESSURE LOAD

      LOAD = 101

      $

      SUBCASE 201

      LABEL = 2# NORMAL LOADS

      LOAD = 201

      $

      SUBCASE 301

      LABEL = OPPOSING LOADS

      LOAD = 301

      $

      include ’plot.blk’

      $

      BEGIN BULK

      $

      $ unique entries for inertia relief

      $

      SUPORT,13,123456     

      PARAM,INREL,-1

      $

      param,post,-1

      $

      include ’modeli.dat’

      include ’seset.dat’

      include ’load1.dat’

      $

      Enddata


从 inrel2.f06 中选择的输出数据

以下为从使用零件超单元惯性释放的求解结果中选择的输出结果。同样,所有三种方法提供了相同的结果。

单个超单元的输出:

每个超单元的输出都包含 OLOAD RESULTANT 的内容,将在第 8 章“超单元静力分析中的输出”中说明。除了各子情况的合成结果外,输出结果中还包含 6 个附加的载荷矢量,代表在参考点 (PARAM,GRDPNT) 或基本坐标系原点上分别施加 6 个方向的单位加速度载荷。这些载荷不是直接施加的,而是用来计算惯性载荷以与施加的外载荷平衡。这一计算作为残余结构运算的一部分 (实际矩阵运算的说明见 NX Nastran User’s Guide  291 页上的“Solution Sequences” – 求解序列)。

      SUBCASE 101

      USER INFORMATION MESSAGE

      ORIGIN OF SUPERELEMENT BASIC COORDINATE SYSTEM WILL BE USED AS REFERENCE LOCATION.

      RESULTANTS ABOUT ORIGIN OF SUPERELEMENT BASIC COORDINATE SYSTEM IN SUPERELEMENT BASIC SYSTEM

      COORDINATES.0 OLOAD RESULTANT

                    T1                        T2                       T3                         R1                      R2                     R3

      101 0.0000000E+00 0.0000000E+00 -1.6000000E+01 -1.2800000E+02 -5.1200001E+01 0.0000000E+00

      201 0.0000000E+00 0.0000000E+00 2.0000000E+00 2.0000000E+01 1.0400000E+01 0.0000000E+00

      301 0.0000000E+00 0.0000000E+00 2.0000000E+00 2.0000000E+01 1.0400000E+01 0.0000000E+00

      302 -5.8637600E-04 0.0000000E+00 0.0000000E+00 0.0000000E+00 0.0000000E+00 4.6910080E-03

      303 0.0000000E+00 -5.8637600E-04 0.0000000E+00 0.0000000E+00 0.0000000E+00 1.8764031E-03

      304 0.0000000E+00 0.0000000E+00 -5.8637600E-04 -4.6910080E-03 -1.8764031E-03 0.0000000E+00

      305 0.0000000E+00 0.0000000E+00 -4.6910080E-03 -3.8372442E-02 -1.5011225E-02 0.0000000E+00

      306 0.0000000E+00 0.0000000E+00 -1.8764031E-03 -1.5011225E-02 -6.8488712E-03 0.0000000E+00

      307 4.6910080E-03 1.8764031E-03 0.0000000E+00 0.0000000E+00 0.0000000E+00 -4.5221314E-02

 

上面结果中的 OLOAD RESULTANT 为执行 inrel2.dat 文件的合成结果,包含三种外载荷的输出 (子情况 101, 201 和 301),随后是对于 6 个单位加速度的输出。注意:各子情况的合成结果用子情况的编号识别,而单位加速度载荷的合成结果使用整数编号, 其值从最大子情况编号加 1 开始。同样对模型中的每一个超单元进行输出,对于主模型数据段定义的超单元,将对施加到内部单元和节点上的载荷求和 (对于 PARTs,将对该零件的模型数据中定义的所有载荷求和)。

最后 6 个载荷是施加到当前超单元的物理质量上的单位加速度载荷的合成 (不包括上游超单元的质量). 这些载荷不包括在主模型数据段中定义的超单元的外部节点上附着的集中质量单元,但是包括对于零件超单元的分区模型数据中定义的所有质量。

在处理超单元时,对每个超单元都输出 OLOAD RESULTANT。

在残余结构中,程序计算使结构达到平衡状态所需的惯性载荷。

在残余结构中,程序也求解矩阵 DM,其中每一列代表一个 Suport 自由度移动一个单位而其它 Suport 自由度被约束时,模型的静力解。这个矩阵不输出,但该解的合应变能以表格形式在数据块 KLL 中输出,标题为 User Information Message 4158。

上表中输出的值应当是数值零 (即本质上是零)。但程序不做检查,因为 NX Nastran 是无量纲的。在使用一种单位制时,一个值可以代表数值零,而在另一种单位制中却可能代表有限的值。在使用惯性释放时,总要检查这个表。

SUPPORT PT. NO. 是参考点集 (r-set) 的自由度编号,由 SUPORT 或 SUPORT1 卡定义。如果任何 (参考自由度载荷矢量) 的应变能值不是数值零,或如果 EPSILON 值不是数值零,则表示模型受到过度的约束 (在关联的方向上原本没有刚体运动的可能) 或有其它问题。

 

      *** USER INFORMATION MESSAGE 4158 (DFMSA)

      ----STATISTICS FOR SPARSE DECOMPOSITION OF DATA BLOCK KLL FOLLOW

      NUMBER OF NEGATIVE TERMS ON FACTOR DIAGONAL = 0

      MAXIMUM RATIO OF MATRIX DIAGONAL TO FACTOR DIAGONAL = 3.0E+03 AT ROW NUMBER 38

      0*** USER INFORMATION MESSAGE 3035 FOR DATA BLOCK KLR

      SUPPORT PT.NO. EPSILON STRAIN ENERGY EPSILONS LARGER THAN .001 ARE FLAGGED

      WITH ASTERISKS

      1 1.4363500E-13 7.3050614E-09

      2 1.4363500E-13 2.5611371E-09

      3 1.4363500E-13 2.5465852E-11

      4 1.4363500E-13 1.4478019E-09

      5 1.4363500E-13 2.4965630E-10

      6 1.4363500E-13 1.3865602E-07

       

接下来输出的是几个矩阵。矩阵 QRR 表示整个结构的刚体质量矩阵,包括所有超单元。计算是相对参考节点进行的,矩阵 QRR 是一个 6 x 6 矩阵,使用总体坐标系。计算时忽略标量点上的质量:

 

      INTERMEDIATE MATRIX ... QRR

      COLUMN 1

      1.583215E-03 -9.274533E-14 0.000000E+00 0.000000E+00 0.000000E+00 -5.204673E-03

      COLUMN 2

      2.745478E-14 1.583215E-03 0.000000E+00 0.000000E+00 0.000000E+00 6.332861E-04

      COLUMN 3

      0.000000E+00 0.000000E+00 1.583215E-03 5.204673E-03 -6.332861E-04 0.000000E+00

      COLUMN 4

      0.000000E+00 0.000000E+00 1.090425E-02 4.393803E-02 -4.361699E-03 0.000000E+00

      COLUMN 5

      0.000000E+00 0.000000E+00 2.091180E-14 1.040331E-13 1.562481E-02 0.000000E+00

      COLUMN 6

      -1.090425E-02 1.056872E-12 0.000000E+00 0.000000E+00 0.000000E+00 5.956284E-02

       

矩阵 QRL 是外载荷之和,相对参考点进行计算,每个载荷条件占用一列:

 

      0 SUBCASE 101

      INTERMEDIATE MATRIX ... QRL

      COLUMN 1

      0.000000E+00 0.000000E+00 3.200000E+01 1.408000E+02 -1.280000E+01 0.000000E+00

      COLUMN 2

      0.000000E+00 0.000000E+00 -4.000000E+00 -2.560000E+01 1.600000E+00 0.000000E+00

      COLUMN 3

      0.000000E+00 0.000000E+00 -3.758416E-11 -1.821405E-10 -2.08000E+01 0.000000E+00

 

矩阵 URACCEL 为用户施加的刚体加速度,以卡片DMlG,UACCEL 定义。如果不存在该卡片,则不出现这一输出 (本例中未提供)。

矩阵 URA 为刚体加速度矩阵,由施加的载荷计算。在算出这些加速度之后,将计算使结构达到静力平衡状态所需的惯性载荷:

      INTERMEDIATE MATRIX ... URA

      COLUMN 1

      0.000000E+00 0.000000E+00 1.009351E+04 -4.400139E+03 -6.921340E-10 0.000000E+00

      COLUMN 2

      0.000000E+00 0.000000E+00 -8.071358E+03 1.538730E+03 5.292301E-09 0.000000E+00

      COLUMN 3

      0.000000E+00 0.000000E+00 -3.727374E-09 1.434969E-09 1.331216E+03 0.000000E+00

 

在求出使结构达到平衡状态所需的惯性载荷之后,将施加这些惯性载荷并与外载荷组合到一起。

然后程序将约束 SUPORT (参考点) 自由度并进行静力求解。此后的输出是常规的超单元输出。

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

历史上的今天

在LOFTER的更多文章

评论

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

页脚

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