培训体系MARC培训手册

培训体系MARC培训手册

MARC

培训手册

美国曼科分析研究公司北京代表处

录

1.MARC程序介绍

1.1MARC程序结构特点

1.2MARC安装后的有关目录

1.3MARC程序手册的组成

1.4MARC的文件系统

1.5MARC软件运行

1.6运行MARC所需的Job文件的数据构成

2.弹塑性分析

2.1理论概要

2.2选项的使用方法

2.3应力分析的控制

2.3.1增量步迭代算法

2.3.2迭代求解控制（CONTROL选项）

3.热传导分析

3.1理论概要

3.2自动步长控制

3.3选项的使用方法

4.热应力分析

4.1理论概要

4.2自动步长控制

4.3选项的使用方法

5.实例分析

5.1弹塑性分析

5.2热传导分析

热传导分析的输入数据和结果

5.3热应力分析

热应力分析的输入数据和结果

附录

I有限单元法的基本步骤

n单元的使用方法及数值积分

m求解器

iv内存分配和空间需求

v重起动和后处理

VI用TYING实现节点自由度间的线性约束

VIISingularityRatio（奇异比）

1.MARC程序介绍

MARCAnalysisResearchCorporation（简称MARC）始创于

1967年，总部设在美国加州的PaloAlto，是全球第一家非线性有限

元软件公司。创始人是美国著名布朗大学应用力学系教授，有限元分

析的先驱公司在创立之初便独具慧眼，瞄准

PedroMarceloMARC

非线性分析这一未来分析发展的必然，致力于非线性有限元技术的研

究、非线性有限元软件的开发、销售和售后服务。对于学术研究机构,

MARC公司的一贯宗旨是提供高水准的CAE分析软件及其超强灵活

的二次开发环境，支持大学和研究机构完成前沿课题研究。对于广阔

的工业领域,MARC软件提供先进的虚拟产品加工过程和运行过程

的仿真功能，帮助市场决策者和工程设计人员进行产品优化和设计，

解决从简单到复杂的工程应用问题。经过三十余年的不懈努力，

MARC软件得到学术界和工业界的大力推崇和广泛应用，建立了它

在全球非线性有限元软件行业的领导者地位。

MARC公司的主要产品之一是通用的有限元分析软件

MARC/MENTA。包括前后处理界面MENTAT和求解器MARC。可

从下图了解MARC和MNETAT之间的数据传输关系：

MARC和MENTAT之间可分可合°MENTAT能够自动生成

MARC分析计算所需的模型数据文件（,dat\MARC分析后所

生成的结果文件即为后处理文件（上19或416）,可由MENTAT

读入后进行数据结果的图形显示。

本培训资料主要是从求解器分析的角度，介绍MARC的基本使用

方法。

1.1MARC程序结构特点

MARC是基于位移法的有限元程序，在非线性方面具有强大

的功能。程序按模块化编程，工作空间可根据计算机内存大小自

动进行调整。用户如果对精度要求较高，可选用双精度进行运算。

当单元数、节点数太多，内存不能满足需要时，程序能够自动利

用硬盘空间进行分析。在分析过程中，利用网格自适应和重划分

技术，能够变更单元的划分和节点数目。

MARC对于非线性问题采用增量解法，在各增量步内对非线

性代数方程组进行迭代以满足收敛判定条件。根据具体分析的问

题可采用不同的分析方法，如对于弹塑性分析和大位移分析可采

用切线刚度法，对于蠕变分析或热应力分析可采用初应变法。

单元刚度矩阵采用数值积分法生成。连续体单元及梁、板、

壳单元的面内区域采用高斯积分法，而梁、板、壳单元厚度方向

则采用任意奇数个点的Simpson积分法。应变-位移函数根据高

斯点来评价。程序计算、存贮单元所有积分点或单元中心点的应

力、应变、温度等。总刚矩阵，总质量矩阵等采用轮廓或稀疏存

贮法存贮。

输入数据可借助于Mentat界面生成，由用户填卡也可生成。

输入数据文件由四部分组成，其中分析功能的指示由参数选项组

确定，分析模型的内容由模型定义选项组确定，增量步数据由历

程定义选项组确定，图形输出由绘图定义选项组确定。

MARC程序拥有许多对用户开放的子程序即用户子程序,用

户能够根据各自需要用FORTRAN语言编制用户子程序,实现对

输入数据的修改、材料本构关系的定义、载荷条件、边界条件、

约束条件的变更，甚至扩展MARC程序的功能。

MARC程序拥有以下4个库

・单元库

•功能库

•分析库

・材料库

用MARC软件分析的每一个实际问题，都采用了这四种库的

一种之上的元素。用户根据各种具体的结构分析进行适当的选择,

下面简单介绍一下这四个库。

单元库

MARC的单元库提供了近150种单元,除了少部分单元

外，其它单元均可用于线性和非线性分析，分析中单元数和

单元类型可自由选择，不同类型单元可组合使用，如连接出

现不协调，可用MARC提供的多种标准连接约束来保证单元

间的一致性。四边形单元能够退化成三角形单元，六面体单

元能够退化成五面、四面体单元。

MARC单元库的主要单元包括：

•三维杆单元

•接触/摩擦单元

定方向接触

定距离接触

用拉格朗日乘子引入接触/摩擦约束条件

.弯管单元

考虑载面变形的弯管

•轴对称壳单元

直线单元

曲线单元

任意加载的3节点曲线单元

.梁单元

二维直线单元、曲线单元

三维直线单元、曲线单元

二维Timoshyenko曲梁单元

三维薄壁闭载面单元

三维薄壁开载面单元

壳梁组合、加固构件

考虑载面翘曲和扭转的薄壁开口梁单元

•平面应力单元

3节点、4节点、8节点等参单元

•平面应变单元

3节点、4节点、8节点等参单元

6节点、9节点半无限单元

•轴称实体元

3节点、4节点、8节点等参单元

考虑扭转的4节点单元

任意加载8节点单元

•广义平面应变单元

4(+2)节点、8(+2)节点等参单元

对于相对位移、相对转动采用附加自由度考虑弯曲的

轴对称单元

•考虑弯曲的轴对称单元

4节点、8节点单元

接触/摩擦单元

•三维薄膜元

4节点、8节点等参单元

•平板单元

3节点、4节点、8节点单元

•曲面壳单元

3节点、4节点、8节点薄壳单元

•三维壳单元

8节点薄壳/厚壳单元

4节点薄壳/厚壳单元

.三维实体单元

8节点、20节点等参单元

•不可压缩单元

5节点、8节点平面应变单元

一般平面应变单元

5节点、8节点轴对称单元

才虎扭转的轴对称单元

任意加载轴对称单元

采用Herrmann变分原理，将静水压力作为附加自由

度

•Rebar单元

8节点平面应变单元

一般平面应变单元

8节点三维单元

8节点轴对称单元

20节点三维单元

.三维剪切板单元

4节点三维单元

•热传导单元

2节点、3节点三维连接元

3节点、4节点、8节点轴对称单元

3节点、4节点、8节点轴对称单元

8节点。20节点三维固体单元

4节点、8节点壳元

2节点、3节点轴对称壳单元

•其它单元

8节点三维磁场单元

12节点磁场半无限单元

4节点电磁场平面单元

4节点电磁场轴对称单元

8节点电磁场三维单元

•减缩积分单元

高次等参单元减少积分点数，对于低次单元采用归约

积分，只有一个积分点

•常体积膨胀

适用低阶单元，用于模拟不可压缩行为，避免传统方

法中过于刚强硬、易于锁死的缺陷

•假定应变法

增强低阶单元模拟弯曲变形能力

.缆索单元

功能库

MARC功能库包含了对分析目标进行准确模拟、快速生

成输入数据、准确高效进行分析以及多种结果输出的众多功

能，具体如下：

•网格自动划分功能

二维单元自动划分

三维单元自动划分

壳坐标的生成

PIPELINE自动划分

增分法、分割法

•输入数据自动行生成功能

增分法

分割法

有规律的边界条件、载荷条件的生成

•几何约束条件

局部坐标系

壳单元局部坐标系

弹簧（线性、非线性）

Tying（节点自由度间的线性约束条件）

边界条件

•载荷条件

集中载荷、分布载荷、重力、离心力、向心力

热载、初始应力

•载荷增量的控制

机械载荷增量，指定位移增量自动调节、时间增量的自

动调节（动力分析、热传导分析、蠕变分析、粘弹性分

析）

•输入方法

自由格式

统一输入方法（关键词、数）

集名的定义

隐含值的利用

.带宽优化功能

Cuthill-Mckee算法

Sloan算法

从外部文件读内部节点和用户节点的对应表能够几种

优化方法组合使用

能够几种优化方法组合使用

.用户子程序

各种本构关系（弹性、塑性、粘弹性、蠕变等）

各向异性（弹性、塑性主轴的定义）

输入、输出（输入数据的生成、修正、输出数据的控制）

节点坐标值的生成、修改

单元划分

Tying条件式

绘图输出变量的定义

•外存的利用

必要时可使用外部文件

.刚度矩阵的部分再生成

适用于超大型问题

•重起动功能

•后处理功能

单元、节点坐标、各种应力、应变、位移、反力等

.输出控制功能

对于必须的输出数据进行详细控制

•adaptivemeshing功能

自动加密或减疏网格

.单元死活功能

分析过程中增加或减少单元

•Rezoning功能

网格退化的处理

分析结果的平滑化

•求解器

直接求解器

迭代EBE

迭代SPARSE

非对称求解器

且行求解器

分析库

MARC分析库包含许多分析类型,用户根据具体问题,

需作必要的选择，分析库具体内容如下：

•线性分析

.弹塑性分析

载荷按比例加一步达到初始屈服（SCALING）

自动增量

•蠕变分析

自动时间增量

•热应力分析

从后处理文件将热传导分析结果读入

自动增量

.粘弹性分析

自动增量

•大变形分析

更新拉格朗日法

屈曲特征值、屈曲模态

载荷一位移曲线

蠕变屈曲

牛顿-拉弗林迭代法

应变修正法

.有限鳖性应变分析

超弹性材料应变能函数

刚塑性分析

欧拉法

定常速度场

•有限应变弹性分析

超弹性材料的应变能函数

.断裂分析

J积分

Lorenzi求J积分方法

•裂纹扩展

裂纹承受拉伸载荷时材料软化，承受压缩载荷时裂纹闭

锁、压坏。

•动力分析

固有频率、固有模态因子

Lancos法、子空间法、模态迭加法、Houbolt法、中

心差分法、Newmark-Beta法

瑞利阻尼

数值阻尼

•频谱分析

适用于地震最大响应分析

•粘弹性体简谐振动分析

•流体-固体耦合分析

•轴对称结构承受任意载荷

傅立叶展开

•多次线性分析结果迭加

•热传导分析

稳态、瞬态热传导分析

对流换热边界条件（包括辐射）

热流、潜热

•热-电耦合分析

.润滑分析

•热-机械耦合分析

•自动塑性加工分析

弹塑性、刚塑性、应变速率相关本构关系、自动接触条

件

材料库

MARC材料库包含30多种材料本构模型,能够考虑材料

的线性和多种非线性材料特性的温度相关性、各向异性等。

具体如下：

•弹性

材料特性和温度相关

•塑性

VonMises屈服准则

Prandlt-Reuss流动准则

各向同性硬化准则

运动硬化准则

混合硬化准则

屈服应力、硬化模量和温度用关

T-T-T（时间、温度相关的相变）

•蠕变

等效应力、等效蠕变应变、温度及时间硬化蠕变模型

纯膨胀引起的蠕变模型

•ORNL推荐的模型

高温下金属循环塑性辅助本构关系

第10个循环的屈服应力作反向屈服应力

第100个循环的屈服应力作循环加载的屈服应力

塑性和蠕变相互作用

•和静水压力相关的屈服准则

岩土、冰等材料常用

线性Mohr-Coulomb屈服准则

抛物线型Mohr-Coulomb屈服准则

.亚弹性

任意的应力增量-应变增量关系式

•不可压缩材料

波松比等于0.5或接近0.5的材料用

Herrmann描述

Taylor各向异性本构关系

•超弹性

橡胶类材料

Mooney-Rivlin本构关系

Jamus-Green-Simpson本构关系

Neo-Hookean本构关系

•粘弹性

Maxwell应力松驰模型

Vogit蠕变模型

温度相关性

TRS（单纯的热粘弹性）材料

•大变形影响

第二类Piola-Kirchhoff应力

格林应变

•有限应变影响

其实应力

对数应变

.热膨胀系数

温度相关

•各向异性

弹性

塑性

热膨胀系数

•复合材料

破坏条件（最大应力、最大应变、Tsai-Wu、Hoffman）

的多种定义

•低抗拉材料

冰、带裂纹构件

温度相关性

.一般鳖性本构关系

屈服条件

动准则

1.2MARC安装后的有关目录

目录名内容

marck72/bin各个执行文件

marck72/lib各个程序的。切文件

marck72/commonMARCK72分析程序的公共块

marck72/toolsMARC程序的执行命令及各

个

程序执行模块生成的命令

marck72/demo例题输入文件及其用户子程序

marck72/user用户子程序的模板文件

1.3MARC程序手册的组成

MARC提供以下手册以供参考

A卷程序功能

第1章MARC系统

第2章程序初始化

第3章数据输入

第4章网格定义

第5章结构分析库

第6章非结构分析库

第7章材料库

第8章陲

第9章边界条件

第10章单元库

第U章非线性系统的求解过程

第12章结果输出

B卷单元库

B1介绍

B2单元分类

B3各种单元信息

C卷程序输入

C1程序概要

C2参数选项

C3模型定义选项

C4历程定义选项

C5REZONING

D卷用户子程序

D1介绍

D2定义载荷、边界条件和状态变量的子程序

D3定义材料各向异性和本构关系

D4定义粘塑性和本构关系

D5定义粘弹性

D6定义几何修正

D7定义输出结果

D8定义滑动轴承

D9定义用户的特殊输出要求

E卷例题集

E1介绍

E2线弹性分析

E3弹塑性和蠕变分析

E4大变形分析

E5热传导分析

E6动力问题

E7接触

E8高级课题

E9流体

E10灵敏度和优化

F卷理论文章

1.4MARC的文件系统

1.5MARC软件的运行

■在MENTAT界面中激活RUN>submitl

-以命令行的方式，按下列命令执行分析程序

(1)运行MARC

在UNIX操作系统下，采用〃run_marc〃命令进行运行,

〃run_marc〃具有下列许多选项,具体含义将在(2)中说明。

run_marc-jidjidname-ridridname-pidpidname

一sidsidname-userusername-prog

FORTRAN

PurposeFormat

UnitNumber

16PostoutputUnformatted

17PostnpulUnfamatted

18OptimizationFormaHed

19PostoutputFormatted

20PostnpulFormatted

21Dummy

22Flud^oiidinteraction-LanczosUnlormatted

23LanczosUnformatted

24TemperatureinputFormatted

25TemperatureinputUnformatted

31SubstructuresDataBaseUnlormatted

41PostoutputlordomaindecompositionUnformatted

42PostoutputlordomaindecoopositKXiFormatted

45DesignOptimlzabonFormatted

46DesignSensitivityorOptimizationUnkxmatted

49DefaultsFileFormatted

50RadiationView!actorFormatled

51PostFileLockFileFormatted

52DynamcControlFiloFormatted

97EXI1MessagesFormatted

98GlobalDefaultFileFonnaHed

99MatenalDataBaseFormatted

prognamebackbackvaluesave

savevalue-vfviewname-defdefname

ververifyvalue

"runmarc"运行时，屏幕响应如下:

run_marc-jide2xl-backno

programname:marck72

JobID:e2xl

Usersubroutinename

RestartfilejobID

SubstructureID

PostfilejobID

Savegeneratedmodule:no

Runjobinbackground:no

Aretheseparameterscorrect(yes,no,quit(no))?Yes

MARCJOBe2xlBEGINGEXECUTION

MEMORYINMAINPROGRAMSETTO=18000000

MARCPROGRAMEXECUTION

i

l

l

MARC的信息

MARCEXITNUMBER3004

在NT操作系统下，采用"run_marc〃命令进行运行,

run_marcw具有下列许多选项，具体含义将在（2）中说明°

run_marc-jidjidname-ridridname-pidpidname

-sidsidname-userusername-prog

progname-vfviewname-defdefname-

savesavevalue

"runmarc"运行时，屏幕响应如下：

run_marc-jide2xl

MARCK72Windows_NTVersion

programname:marck72

JobID:e2xl

Usersubroutinename:

RestartfilejobID:

SubstructureID:

PostfilejobID:

Savegeneratedmodule:no

ViewFactorID:

DefaultFileID:

MARCJOBe2xlBEGINGEXECUTION

MEMORYINMAINPROGRAMSETTO=

18000000

MARCPROGRAMEXECUTION

MARC的信息

MARCEXITNUMBER3004

(2)Mrun_marcM关键词

关键词变量内容

-prog*〃marc〃运行MARC

“mesh3d〃运行1MARC的mesh3d

簸

"plot"运行MARCPLOT

“Pipe"运行PIPE

(progname)用户生成的MARC执

行文件名

-jid(jidname)输入文件名

(jidname).dat必须存在

-rid(ridname)国亍MARC日瞳血皮件S

运行MESH3D,PIPE时

生成的网格文件名

-pid(pidname)运行MARC时热应力

计算的后处理文件名

运行PLOT时读入的后处

理文件名

-sid(sidname)子结构用文件名

-user(username)含有用户子程序的用

户文件<usemame>.f

生成可执行文件

<username>.marc

-back〃yes〃or〃y〃在后台执行

*〃no〃or〃n〃在前台执行

-save〃yos〃ory保留生成的执行文件

*〃no〃or不保留生成蝴，行文件

-vf(viewname)辐魄｛角藏懒敢件

-def(defname)用户缺省的输入数据

文件

-ver*"yes"or"y""run_marc〃检查命令

选项

〃no〃or〃run_marc〃不检直命

令选项

(★为隐含值)

(3)使用实例

1.run_marc-jide2xl

采用输入数据文件e2xl.dat进行分析

2.run_marc-jide2xl4-useru2xl4-savyu2x

14.f

输入数据文件为e2xl4.dat,执行文件为由用户文件

u2x14.f编译、连接形成的u2x14.marc且将执行文件

u2xl4.marc保留。

3.run_marc-jidex14a-progu2xl4

利用2已生成的执行文件u2xl4.marc进行分析

4.run_marc-jide3x2a-vno-bno

在前台运行且不检查输入命令的选项

5.run_marc-jide3x2b-ride3x2a利用4生成的重

起动文件e3x2a.t08继续分析e3x2b.dat

进行e5xlla的温度场分析

6.run_marc-jide5xlla

7.run_marc-jide5xllc-pide5xlla

利用6生成的后处理文件e5xlla.tl6,进行e5xllc

的热应力分析

1・6运行MARC所需的Job文件的数据构成

MARC程序的输入数据大致由以下几部分组成

Param«t»r

包括：

1.参数选项

2.模型定义选项

3.历程定义选项

4.绘图输出选项

如果要使用用户子程序，用户必须提供用FORTRAN语言编

写的用户子程序，且提供这些用户子程序所需的一些数据。

参数选项具有以下功能

1.确保分析所需的工作空间

2.注册使用的单元

3.指示使用的功能

4.对作业的空间进行分配（边界条件的最大个数、分布载荷

单元表中最大单元个数）

重起动时这部分数据不允许改动

这部分数据以‘END，为结束标记

模型定义选项具有以下功能

1.输入的单元节点编号及节点坐标

2.输入材料常数、几何形状参数

3.输入位移边界条件、载荷条件

4.输入其它进行模型分析时所需的数据及分析方法、收敛

判定方法、收敛判定条件等分析控制数据

除控制数据以外，其它数据在重起动时不允许改动

这部分数据以‘ENDOPTION'为结束标记

历程定义选项具有以下功能

1.载荷增量、位移增量、时间增量等增量数据的定义

2.输出控制、约束条件的变更数据输入

这部分数据以‘CONTINUE，为增量步结束标记,能够

填入多个历程定义数据选项，描述多步加载历程。

绘图输出选项具有以下功能

结构图（未变形1变形图、等应力线图、等应变线图绘制

数据的定义

各绘图输出选项以TNDPLOT结束。

输入数据顺序

分别按参数选项、模型定义选项、历程定义选项的顺序输

入，绘图输出选项根据要求的位置插入。

输入数据的规则及特点

1.以关键词的形式输入数据，根据分析所需，选择必要

的关键词输入°

2.输入数据有一定的顺序要求，应根据手册C卷要求的

顺序输入。

3.数据的输入基本按整数15、实数E10或E15的格式

输入。

4.除采用固定格式外，仍可用自由格式输入，此时数据

之间以逗号相隔。同一个输入文件能够固定格式

和自由格式混用，但每一个选项只能用一种格式。如

果某一行数据只输入一个整型数而又用自由格式输

入，则需在整数后加〃；'。

5.同一关键词能够使用多次，当数据重复定义时，程序

采用最后一次输入的数据。

6.载荷数据、边界条件数据处理成增量值，指定载荷状

态、位移状态按增量值输入。

项目表（List）输入方法

载荷条件、位移边界条件施加的节点和材料特性、几何形

状定义的单元，称为项目表。项目表的输入方法有新旧二种，

对于Kxx版均采用新的格式，下面仅介绍新格式的输入方法，

旧格式可参见MARC手册A卷、C卷。

1.将表内容全部直接列出

例：315241

2.用TO、BY形式将有关内容列出

例：5TO37BY8

3.用集名列出，集名所含项目在DEFINE选项中定义

例：MYSET

另外，可用AND.INTERSECT.EXCEPT将多个子集进行

组合

例：〃子集1AND子集2EXCEPT子集3〃，各

子集用之上3种方法定义

例：5TO37BY8ANDMYSETEXCEPT3

15241

输入、输出文件举例

下面通过一个简单例题对输入数据及输出结果文件作概要

说明。

例题：悬臂梁受端部载荷的问题

一端固定、长为10in、截面为1x1iM的悬臂梁,承

受端部800ps的集中力。材料的弹性模量为30000000

ps/iG,波松比0.3。屈服应力和塑性应变的关系见下图

有限元分析模型采用40个四节点四边形单元，单元类号

为11。梁的左端固支，右顶端施加800ps载荷。

输入数据

参数选项

必须有的选项

TITLE:标题

SIZING:单元总数、节点总数等

LARGEDISP:考虑几何非线性

DISTLOAD:载荷类型

ELEMENTS:单元类型

END:参数选项结束

模型定义选项

CONNECTIVITY定义单元节点编号

COORDINATES定义节点坐标

GEOMETRY输入单元截面积

ISOTROPIC定义材料常数

POINTLOAD定义节点集中载荷

FIXEDDISP定义边界条件

ENDOPTION模型定义选项结束

历程定义选项

CONTROL求解控制信息

AUTOLOAD载荷步长控制

POINTLOAD增量步载荷

CONTINUE历程选项结束

对于弹塑性问题的求解，程序要求输入历程定义选项进行增

量非线性分析。

输入数据文件内容如下:

titlejobl

S....MARCinputfileproducedbyMentat3.2

sizing10000004055110

elements3

processor111

$nolist

largedisp

allpoints

distloads001

setname3

end

solver

0000000

optimize9

connectivity

131672

232783

333894

4349105

53611127

63712138

73813149

839141510

9311161712

10312171813

11313181914

12314192015

16319242520

19323282924

20324293025

21326313227

22327323328

23328333429

24329343530

25331363732

26332373833

27333383934

28334394035

29336414237

30337424338

31338434439

32339444540

33341464742

34342474843

35343484944

36344495045

37346515247

38347525348

39348535449

40349545550

coordinates

355

10.00000-00.00000+00.00000+0

20.00000-02.50000-10.00000+0

30.00000-05.00000-10.00000+0

40,00000-07.50000-10.00000+0

50.00000-01.00000+00.00000+0

61.00000-00.00000+00.00000+0

71.00000-02.50000-10.00000+0

81.00000-05.00000-10.00000+0

91.00000-07.50000-10.00000+0

101.00000+01.00000+00.00000+0

112.000004-00.00000+00.00000+0

122.000004-02.50000-10.00000+0

132.000004-05.00000-10.00000+0

142.00000+07.50000-10.00000+0

152.000004-01.00000+00.00000+0

163.000004-00.00000+00.00000+0

173.000004-02.50000-10.00000+0

183.00000+05.00000-10.00000+0

193.000004-07.50000-10.00000+0

203.000004-01.00000+00.00000+0

214.00000+00.00000+00.00000+0

224.000004-02.50000-10.00000+0

234.00000+05.00000-10.00000+0

244.00000+07.50000-10.00000+0

254.00000+01.00000+00.00000+0

265.000004-00.00000+00.00000+0

275.00000+02.50000-10.00000+0

285.000004-05.00000-10.00000+0

295.000004-07.50000-10.00000+0

305.00000+01.00000+00.00000+0

316.000004-00.00000+00.00000+0

326.000004-02.50000-10.00000+0

336.00000+05.00000-10.00000+0

346.00000+07.50000-10.00000+0

356.000004-01.00000+00.00000+0

367.000004-00.00000+00.00000+0

377.000004-02.50000-10.00000+0

387.000004-05.00000-10.00000+0

397.00000+07.50000-10.00000+0

407.000004-01.00000+00.00000+0

418.000004-00.00000+00.00000+0

428.000004-02.50000-10.00000+0

438.000004-05.00000-10.00000+0

448.00000+07.50000-10.00000+0

458.000004-01.00000+00.00000+0

469.000004-00.00000+00.00000+0

479.000004-02.50000-10.00000+0

489.000004-05.00000-10.00000+0

499.00000+07.50000-10.00000+0

509.000004-01.00000+00.00000+0

511.000004-10.00000+00.00000+0

521.000004-12.50000-10.00000+0

531.00000>15,00000-10.00000+0

541.000004-17.50000-10.00000+0

551.00000+11.00000+00.00000+0

isotropic

Ivonmisesisotropic000Omateriall

3.00000+73.00000-17.33161-40.00000+02.00000+40.00000+00.00000+0

0.00000+0

1to40

workharddata

3010

2.0000000000+40.0000000000+0

2.5000000000+41.0900000000-1

3.0000000000+43.0500000000-1

geometry

1.00000+00.00000+01.00000+00.00000+00.00000+00.00000+0

1to40

fixeddisp

0.00000+0

1

1to5

0.00000+0

2

1

pointload

001

0.00000+00.00000+0

55

noprint

post

116171019200107

311

parameters

1.00000+01.00000+91.00000+21.00000+62.50000-15.00000-1

8.62500+02.00000+11.00000-41.0000-121.00000+01.00000-4

endoption

$................

$....startofloadcaseIcasel

titleIcasel

control

999991000010010

1.00000-10.00000+00.00000+00.00000+00.00000+00.00000+00.00000+0

0.00000+0

parameters

1.00000+01.00000+91.00000+21.00000+62.50000-15.00000-1

8.62500+02.00000+11.00000-41.0000-121.00000+01.00000-4

autoload

200

timestep

5.00000-2

pointload0

001

0.00000+0-8.00000+1

55

continue

$....endofloadcaseIcasel

$................

输出结果

运行结果有以下三种

•标准输出（程序输出，为.。仇文件）

.后处理输出（为由mentat可读的乜6或乜9文件）

.重起动输出（为.t08文件）

对于标准输出有以下内容

1.MARC公司的标记及程序版本

2.输入数据的回应及解释

对输入数据进行运行检查。如发现错误，便会在该行

数据前后指明。数据检查后输出对一些数据进行的解

释。

3.运行信息

输出分析的一些信息

a）整个模型上分布载荷及集中载荷的总和

b）奇异比

进行线性方程组求解时，输出的刚度矩阵性质指

标奇异比r能够指示浮点运算精度位数的降低程

度。MARC运行于多种不同的计算机系统，具有

不同的版本，但一般有以下规律：

1>r>10-4解的精度在容许范围

之内

10-4〉r>10-8解的精度位数降低

10-8>r>10-14矩阵非正定

c)CPU时间

MARC在运行中会输出下列关键时刻已用的CPU

时间

STARTOFINCREMENT开始增量

STARTOFASSEMBLY开始刚度集成

STARTOFMATRIXSOLUTION开始矩阵求解

ENDOFMATRIXSOLUTION结束矩阵求解

4.分析结果的输出

输出增量分析结束时的元素信息（应力、应变等1节

点信息（位移、等效节点力、反力等\

a）单元信息

各类单元高斯点应力、应变输出

•连续体单元

整体坐标系下的应力分量、主应力、静水压

力、等效应力

.壳单元

中面应力、沿厚度方向上各积分点的应力值

・梁单元

高斯点的轴向力、弯矩、扭矩

b）节点信息

•静力分析

位移增量、总位移、级节点力、反力、残余力

•动力分析

模态分析—固有频率

动力响应分析—总位移、速度、加速度、等

效节点力、反力、残余力

•热传导分析

节点温度

单元信息、节点信息的输出可用PRINTCHOICE选项

进行控制。

卜面通过悬臂梁受弯的分析，举例说明MARC运行的结果文

件.out的一般格式。

MM

+WW

MMMMMMMMMM

+WWWWWWWWWW

h/MMMMMMMMMMMMMM

MMM

+WWWWWWWWWWWV/WW

WWWW

MI7MMMMMMMMMMMMMMMM

MMMMMMMM

+WWWWWWWWWWWWWWWWW

WWWWWWWWW

MMMMMMMMMMMMMMMMMMMMM

MMMMMMMMMMMM

+WWVi/WWWWWWWWWWWWWWWWW

WWWWWWWWWWWWW

MMIVMMMMMMMMMMMMMMMMM

MMMMMMMMMMM

+WWWWWWWWWWWWWWWWV7WW

WWWWWWWWWWWW

MMMMMMMMMMMMMMMMM

MMMMMM

+WWWWWWWWWWWWWWWWW

wwwwww

MMIVMMMMMMMM

MMMM

+WWWWWWWWWWW

WWWW

MMMMMMM

+WWWWWWW

MMM

+WWW

MMMMMMM

+WWWWWWW

MMIVMMMMMMMM

MMMM

+WWWWWWWWWWW

wwww

MMMMMMMMMMMMMMMMM

MMMMMM

+WWWWWWWWWWWWWWWWW

wwwwww

MMIVMMMMMMMMMMMMMMMMM

MMMMMMMMMMM

+WWWWWWWWWWWWWWWWWWW

WWWWwwwwwwww

MMMMMMMMMMMMMMMMMMMMM

MMMMMMMMMMMM

+WWWWWWWWWWWWWWWWWV7WW

wwwwwwwwwwwww

MK/MMMMMMMMMMMMMMMM

MMMMMMMM

+WWWWWWWWWWWWWWWWW

wwwwwwwww

K/MMMMMMMMMMMMMM

MMM

+WWWWWWWWWWWWWW

WWWW

MMMMMMMMMM

+WWWWWWWWWW

MM

WW

MARC-PC<7-2-L<

(ciCOPYRIGHT1998MARCANALYSISRESEARCHCORPORATION,al

Irightsreserved

MARC-PCinputdata

page1

5101520253035404