ABAQUS梁计算例题部分含详细敕令流

ABAQUS计算指导0：应用梁单元计算简支梁的挠度

£1

1

0

0

0.1m

对于梁的分析可以使用梁单元、壳单元或是固体单元。Abaqus的梁单元需要设定线的方向，用选

中所需要的线后，输入该线梁截面的主轴1方向单位矢量(x,y,z),截面的主轴方向在截面Profile设定中

有规定。

注意：

因为ABAQUS软件没有UNDO功能，在建模过程中，应不时地将本题的CAE模型(阶段结果)

保存，以免丢失已完成的工作。

简支梁，三点弯曲，工字钢构件，结构钢材质，E=210GPa,〃=0.28,片7850kg/m3(在不计重力的

静力学分析中可以不要)。F=10kN,不计重力。计算中点挠度，两端转角。理论解：/=2.239X10-5m4,

w巾=2.769X10-3m,8边=2.077X104

文件与路径：

顶部下拉菜单File,SaveAsExpAbqOO。

-部件

1创建部件：Module,Part,CreatePart,

命名为Prat-l；3D,可变形模型，线，图形大约范围10(程序默认长度单位为m)。

2绘模型图：选用折线，从(0,0)-(2,0)-(4,0)绘出梁的轴线。

3退出：Done。

二性质

1创建截面几何形状：Module,Property,CreateProfile,

命名为Profile-1,选I型截面,按图输入数据，/=0.1,〃=0.2,仇=0.1,岳=0.1,«=0.01,1=0.01，口=0.01，

关闭。

2定义梁方向：Module,Property,AssignBeamOrientation,

选中两段线段，输入主轴1方向单位矢量(0,0,1)或(0,0,-1),关闭。

3定义截面力学性质：Module,Property,CreateSection,

命名为Section-1,梁，梁，截面几何形状选Profile-1,输入E=210e9(程序默认单位为N/m2,

GPa=l(PN/m2),

G=82.03e9,v=0.28,关闭。

4将截面的几何、力学性质附加到部件上：Modu1e,Property,AssignSection,

选中两段线段，将Section-1信息注入Part-1o

三组装

创建计算实体：Module,Assembly,顶部下拉菜单Instance,Create,

CreateInstance,以Prat-1为原形，用Independent方式生成实体。

四分析步

创建分析步：Module,Step,

CreateStep,命名为Step-1,静态Static,通用General。注释：无，时间：不变，非线性开关：

关。

五载荷

1施加位移边界条件：Modu1e,Load,CreateBoundaryCondition,

命名为BC-1,在分析步Siep-1中，性质：力学，针对位移和转角，Continue0选中梁左端，Don

e,约束〃I、“2、〃3、“RI、〃R2各自由度。

命名为BC-2,在分析步Step-1中，性质：力学，针对位移和转角，Continueo选中梁右端，Don

e,约束“2、"3、uRi>〃R2各自由度。

2创建载荷：Modu1e,Load,CreateLoad,

命名为Load-1,在分析步Step-1中，性质：力学，选择集中力ConcentratedForce,Continue

o选中梁中点，Done,施加产y(CF2)=・10000(程序默认单位为N)。

六网格

对实体Instance进行。

1撒种子：Modu1e,Mesh,顶部下拉菜单Seed,Instance,

G1obalSeeds,Approximategloba1size0.2全局种子大约间距0.2。

2划网格：Module,Mesh,顶部下拉菜单Mesh,Instance,yes。

七建立项目

1建立项目：Modu1e,Job,CreateJob,Instance,

命名为ExpAbqOO,选择完整分析，其余先不变，OK。

2计算：Modu1e,Job,JobManager,

Submit,可以用Monitor进行求解进程观察，算完。

八观察结果

1看图：

打开结果文件ExpAbqOO.odb,看位移彩图，学习修改图形参数。

2看数据：计算结果卬中=2.853XIO'm,0^=2.083X10-3。

学习显示节点号等信息，学习存储结果数据文件，学习查找节点、单元数据结果。

*Heading

**Jobname:ExpAbqOOintModelname:Model-1

*Preprint,echo=NO,mode1=NO,history=NO,contact=NO

**

**PARTS

*Part,name=Part-1

*EndPart

**

**ASSEMBLY

*Assembly,name=Assembly

*1nstance,name=Part-1-1,part=Part-1

*Node

1,0.0,0.0,0.0

21,4.0,0.0,0.0

*Ngen,,Nset=Partl

1,21,1,,,,,0.0,1.0,0.0

*Element,type=B31

1,1,2

*E1gen,Elset=Part1

1,20,1,1

**Section:Section-1Profile:Profile-1

*BeamGenera1Section,eIset=Par11,poisson=0.28,section=I

0.1,0.2,0.1,0.1,0.01,0.01,0.01

O.,O.,-1.

2.1e+11,8.203e+l0

*EndInstance

**

*Nset,nset=A,interna1,instance=Part-1-1

1,

*Nset,nset=B,interna1,instance=Part-1-1

21,

*Nset,nset=C,interna1,instance=Part-1-1

”,

*EndAssembly

**__________________________________________

**STEP:Step-1

*Step,name=Step-1

*Static

1.,1.,le-05,1.

**BOUNDARYCONDITIONS

**Name:BC-lType:Displacement/Rotation

*Boundary

A,1,5

**Name:BC-2Type:Displacement/Rotation

*Boundary

B,2,5

**LOADS

**Name:Load-1Type:Concentratedforce

*Cload

C,2,-10000.

**OUTPUTREQUESTS

*Restart,write,frequency=0

**FIELDOUTPUT:F-0utput-1

*Output,fieId,variable=PRESELECT

**HISTORYOUTPUT:H-Output-1

*Output,history,variable=PRESELECT

*EndStep

ABAQUS计算指导ExpAbqOI

ABAQUS计算指导1：应用梁单元计算平面结构

尸1廿

Bz,2

F

z,2

4mdp

D

勿R

门式框架，方管构件，结构钢材质，E=2lOGPa,"=0.28,p=7850kg/m3(在不计重力的静力学

分析中可以不要)。尸=1(X)()N,不计重力。计算最大弯矩，B点水平位移。理论解：/产8.617XI

0-7m4,h=2.887XIO^m4,

也=^T=°。15396m=肛=等=2500N.m

文件与路径

顶部下拉菜单File,SaveAsExpAbq010

一部件

1创建部件：Modu1e,Part,CreatePart,

命名为Prat/；3D,可变形模型，线，图形大约范围20(程序默认长度单位为m)。

2绘模型图：选用折线，从A(0,0)一8(0,4)—C(8,4)-0(8,0)绘出框架的轴线。

3退出：Done。

二性质

1创建截面几何形状：Module,Property,CreateProfile,

命名为Profile-1,选Box型截面，按图输入数据，Q=0.1,b=0.05,Thickness：uniform,Th

ickness：0.01,关闭。

2定义各段梁的方向：Modu1e,Property,AssignBeamOrientation,

选中左、右立柱，输入主轴1方向单位矢量(1,0,()),选中横梁，输入主轴1方向单位矢量(0,1,0),关闭。

3定义截面力学性质：Module,Property,CreateSection,

命名为Section-1,梁，梁，截面几何形状选Profile-1,输入E=210e9（程序默认单位为N/m2

GPa=109N/m2）,

G=82.03e9,v=0.28,关闭。

4将截面的几何、力学性质附加到部件上：Module,Property,AssignSection,

选中左右立柱和横梁，将Section-1信息注入Part-1o

三组装

ABAQUS计算指导ExpAbqOl

创建计算实体：Module,Assemb1y,顶部下拉菜单Instance,Create,

CreateInstance,以Prat-1为原形，用Independent方式生成实体。

四分析步

创建分析步：Module,Step,

CreateStep,命名为Step-1,静态Static,通用General。注释：无，时间：不变，非线性开关：

关。

五载荷

1施加位移边界条件：ModuIe,Load,CreateBoundaryCondition,

命名为BC-1,在分析步Step-1中，性质：力学，针对位移和转角。选中立柱两脚，约束全部自由度。

2创建载荷：Modu1e,Load,CreateLoad,

命名为Load-1,在分析步St叩-1中，性质：力学，选择集中力。选中立柱两顶点，施加F*=CF1=1

000(程序默认单位为N)。

六网格

对实体Instance进行。

1撒种子：Modu1e,Mesh,顶部下拉菜单Seed,Instance,

G1obalSeeds,Approximategloba1size0.2全局种子大约间距0.2。

2划网格：Module,Mesh,顶部下拉菜单Mesh,Instance,yes«

七建立项目

1建立项目：ModuIe,Job,CreateJob,Instance,

命名为ExpAbqOO,选择完整分析，其余先不变，OK。

2计算：Modu1e,Job,JobManager,

Submit,可以用Monitor进行求解进程观察，算完。

八观察结果

1看图：打开结果文件ExpAbqOl.odb,看位移彩图，学习修改图形参数。

2看数据：

显示节点号等信息。学习存储结果数据文件，学习查找节点、单元数据结果，“&=0.0157383。

3保存你的模型：将本题的CAE模型保存为ExpAbq()1.caeo

全令流

*Heading

**Jobname:ExpAbqO1inpModelname:Model-1

*Preprint,echo=NO,mode1=NO,history=NO,contact=NO

**

**PARTS

ABAQUS计算指导ExpAbqOI

*Part,name=Part-1

*EndPart

**

**ASSEMBLY

*Assembly,name=Assembly

*1nstance,name=Part-1-1,part=Part-1

*Node

1,0.0,().(),().()

21,0.0,4.0,0.0

61,8.0,4.0,0.0

81,8.0,0.0,0.0

*Ngen,,Nset=Partl

1,21,1,,,,,,,

21,61,1,,,,,,,

61,81,1,,,,,,,

*Element,type=B31

1,1,2

*EIgen,Elset=Part1

1,80,1,1

*EIset,elset=Part11,internal,generate

1,20,1

61,80,1

*EIset,elset=Part12,intemal,generate

21,60,1

**Region:(Section-l:V),(BeamOrientation:2)

**Section:Section-1Profile:Profile-1

*BeamGenera1Section,eIset=Partl1,poisson=0.28,section=BOX

0.1,0.05,0.01,0.01,0.01,0.01

2.1e+11,8.203e+l0

**Region:(Sectio(BeamOrientation:l)

*EIset,elset=__I2,interna1,generate

21,60,1

**Section:Section-1Profile:Profile-1

*BeamGenera1Section,eIset=Partl2,poisson=0.28,section=BOX

0.1,0.05,0.01,0.01,0.01,0.01

0.,l.,0.

2.1e+11,8.203e+l0

*EndInstance

*Nset,nset=AD,intemal,instance=Part-1-1

1,81

ABAQUS计算指导ExpAbqOI

*Nset,nset=BC,intemal,instance=Part-1-1

21,61

*EndAssembly

**STEP:Step-1

*Step,name=Step-1

*Static

1.,1.,Ie-05,1.

**BOUNDARYCONDITIONS

**Name:BC-lType:Displacement/Rotation

*Boundary

AD,1,6

**LOADS

**Name:Load-1Type:Concentratedforce

*Cload

BC,1,1000.

**OUTPUTREQUESTS

*Restart,write,frequency=0

**FIELDOUTPUT:F-Output-1

*0utput,fieId,variable=PRESELECT

**HISTORYOUTPUT:H-Output-1

*0utput,history,variable=PRESELECT

*EndStep

ABAQUS计算指导ExpAbqO2

ABAQUS计算指导2:应用梁单元计算平面结构

3(在不计重

力的静力学分析中可以不要)。F=1000N,此题要计入重力。计算水平梁中点下降位移。

文件与路径

顶部下拉菜单File,SaveAsExpAbq02«

-部件

创建部件，命名为Prat-1。

3D,可变形模型，线，图形大约范围2()(m)。

选用折线绘出整个图示屋架。

退出Done

二性质

1创建截面几何形状：Modu1e,Property,CreateProfi1e,

将截面(1)命名为Profile-1,选Box型截面，按图输入数据，关闭。直至完成截面(3)。

2定义各段梁的方向：

选中所有立杆，输入截面主惯性轴1方向单位矢量(1,0,0),选中横梁和斜杆，输入截面主轴1方向单位

矢量(0,1,0),关闭。还有好办法，请大家自己捉摸。

3定义截面力学性质：

将截面(1)Profile-1命名为Section-1,梁，梁，截面几何形状选Profile-1,输入E=210GPa,G

=82.03GPa,v=0.28,p=7850,关闭。直至完成截面(3)Section-3»

4将截面的几何、力学性质附加到部件上：

选中左右立柱和横梁，将各Section-l~3信息注入Part-1的各个杆件上，要对号入座。

5保存模型:

将本题的CAE模型保存为ExpAbq02.cae。

三组装

ABAQUS计算指导ExpAbqO2

创建计算实体，以Prat-1为原形，用Independent方式或Dependent生成实体。

四分析步

创建分析步，命名为Step-1,静态Static,通用General。

注释：无，时间：不变，非线性开关：关。

五载荷

1施加位移边界条件：

命名为BC-1,在分析步Step-1中，性质：力学，针对位移和转角。

选中立柱两脚，约束全部自由度。

2创建载荷：

命名为Load-1,在分析步Step-1中，性质：力学，选择集中力。

选中顶点，施加F产尸2=-1000(N)。

六网格

对部件Prat-1进行。

1撒种子：

针对部件，全局种子大约间距0.8。

2划网格：

针对部件，OK。

3保存你的模型：

将本题的CAE模型保存为ExpAbq02.cae.

七建立项目

命名为ExpAbqO2,选择完整分析，其余不变，OK。

八观察结果

1看图：

打开结果文件ExpAbq()2.odb,看位移图。

2看数据：

显示某节点结果数据，例如水平梁中点垂直位移。

3保存你的模型：

将本题的CAE模型保存为ExpAbq02.caeo

ABAQUS计算指导ExpAbqO3

ABAQUS计算指导3：应用梁单元计算空间结构

zm

oe

钢结构框架，圆管构件，材质E=210GPa,〃=0.28,0=7850kg/m3。尸=1000N,此题要计入重

力。计算最大下降位移。

文件与路径

顶部下拉菜单File,SaveAsExpAbq03»

一部件

创建部件，命名为Prat-1。

3D,可变形模型，线，图形大约范围20(m)。

选用折线绘出图示的后面一片刚架。

创建部件，命名为Prat-2。

3D,可变形模型，线，图形大约范围20(m)。

选用直线绘出长度为2m一线。

退出Done。

二组装

创建计算实体，以Prat-i为原形，用Dependent方式生成实体。向z轴间距2m拷贝出位于前面的另

一片刚架。

创建计算实体，以Prat-2为原形，用Dependent方式生成实体。转动至z轴方向，移动至适当位置；

拷贝另8个，全移动至适当位置。

融合全体，命名为Prat-3。

删除无用实体、部件。

三性质

1创建截面几何形状：

将截面(1)命名为Profile-1,选Pipe型截面，按图输入数据，关闭。

2定义各梁方向：

ABAQUS计算指导ExpAbqO3

选中所有Z方向杆，输入截面主惯性轴1方向单位矢量(1,0,0),选其余全部杆，输入截面主轴1方向单位

矢量(0,0,1),关闭。

3定义截面力学性质：

将截面Profile-1命名为Section-1,梁，梁，截面几何形状选Profile-1,输入E=2lOGPa,G=82.0

3GPa,v=0.28,p=7850,关闭。

4将截面的几何、力学性质附加到部件上：

选中所有梁，将Section-1信息注入Part-3的所有杆件上。

四分析步

创建分析步，命名为Step-1,静态Static,通用General。

注释：无，时间：不变，非线性开关：关。

五载荷

1施加位移边界条件：

命名为BC-1,在分析步St叩-1中，性质：力学，针对位移和转角。

选中结构左端4脚，约束全部自由度。

2创建载荷：

命名为Load-1,在分析步Step-1中，性质：力学，选择集中力。

选中结构右端上2点，施加尸尸尸2=■上00(N)。

六网格

对部件Prat-3进行。

1撒种子：

针对部件，全局种子大约间距0.3~0.5m。

2划网格：

针对部件，OK。

七建立项目

命名为ExpAbqO3,选择完整分析，其余不变，OK。

A观察结果

打开结果文件ExpAbq()3.odb,看位移图。

ABAQUS计算指导ExpAbqO4

ABAQUS计算指导4：应用壳单元计算开口薄壁截面梁的横力弯曲

ABAQUE的壳单元是在面上生成的，几何模型没有厚度，物理的厚度可以用定义壳单元的截面厚度

来实现。

3

用壳单元模拟槽形截面梁，梁的长度/=lm»材质E-210GPa,〃=0.28,p=7850kg/moF=100

0N,不计重力。计算最大Mises应力，最大绕z轴转角(附加扭转)。

文件与路径

顶部下拉菜单Fi】e,SaveAsExpAbq04。

一部件

创建部件，命名为Prat-1。

3D,可变形模型，壳，拓展，图形大约范围0.4(m)。

选用折线绘出图示的截面壳的中线，沿z方向拓展l(/=lm)。

-性质

1创建材料：

将材料命名为Material-1,输入£=210GPa,v=0.28,关闭。

2创建截面：

名为Section-1,壳，各向同性，材料名选Material-1,厚度0.002m,关闭。

3将截面的性质附加到部件上：

选中Prat-1,将Section-1信息注入Part-1。

三组装

创建计算实体，以Prat-1为原形，用Inependent方式生成实体。

四分析步

创建分析步，命名为Step-1,静态Static,通用General。

注释：无，时间：不变，非线性开关：打开非线性计算。

五载荷

1施加位移边界条件：

命名为BC-1,在分析步Step-1中，性质：力学，针对位移和转角。

选中结构后端3线，约束全部自由度。

ABAQUS计算指导ExpAbqO4

2创建载荷：

命名为Load-1,在分析步Step-1中，性质：力学，选择集中力。

选中槽钢的前端左上角点，施加尸产尸2=-100()(N)。

六网格

对部件Prat-1进行。

1撒种子：

针对部件，全局种子大约间距0.02(m)。

2划网格：

针对部件，OK。

七建立项目

命名为ExpAbq03,选择完整分析，其余不变，OK。

八观察结果

打开结果文件ExpAbq04.odb,看位移彩图。

在上述过程中不时地将本题的CAE模型保存为ExpAbq04.cae»

牛伞诵

"Heading

**Jobname:ExpAbq04inpModelname:Model-1

*Preprint,echo=NO,model=NO,history=NO,contact=NO

**

**PARTS

*Part,name=Part-1

*EndPart

**

**ASSEMBLY

*Assembly,name=Assembly

*Instance,name=Part-1-1,part=Part-1

*Node

1,0.1,0.08,0.0

6,0.0,0.08,0.0

14,0.0,-0.08,0.0

19,0.1,-0.08,0.0

1001,0.1,0.08,1.0

1006,0.0,0.08,1.0

1014,0.0,-0.08,1.0

1019,0.1,-0.08,1.0

*Ngen,,Nset=PartlA

ABAQUS计算指导ExpAbqO4

1,6,1,,,,,,,

6,14,1,,,,,,,

14,19,1,,,,,,,

*Ngen,,Nset=PartlB

1001,1006,1,,,,,,,

1006,1014,1,,,,,,,

1014,1019,1……，

Nset=Partl

PartiA,PartIB,50,20

*Element,type=S4R

1,1,2,22,21

*Elgen,Elset=Partl

1,18,1,1,50,20,20,1

**Section:Section-1

*ShellSection,elset=Part1,material=Material-1

0.002,5

*EndInstance

*Nset,nset=AD,internal,generate,instance=Part-l-l

1,19,1

*Nset,nset=BCtinternal,instance=Part-l-l

1(X)6

*EndAssembly

**

**MATERIALS

*Material,name=Material-l

.Density

7850,,

*Elastic

2.1e+U,0.28

**____________________________________________

**STEP:Step-1

*Step,name=Step-1,nlgeom=YES

*Static

1.,1.,Ie-05,1.

**BOUNDARYCONDITIONS

**Name:BC-1Type:Displacement/Rotation

“Boundary

AD,1,6

**LOADS

**Name:Load-1Type:Concentratedforce

*Cload

ABAQUS计算指导ExpAbqO4

BC,2,-1000.

**OUTPUTREQUESTS

"Restart,write,frequency=O

**FIELDOUTPUT:F-Output-1

*Output,field,variable=PRESELECT

**HISTORYOUTPUT:H-OutpuM

“Output,history,variable=PRESELECT

*EndStep

ABAQUS计算指导ExpAbqO5

ABAQUS计算指导5：应用壳单元计算轴对称问题

有内压的容器或弹性元件，其工作方式分为两类。第一类内压是通过调节而保持恒定的。对于恒压情况

,内压p不随变形改变。第二类是密闭容器。如空气弹簧、篮球等情况，施加外载荷后内压是随变形改

变的，采用常内压0无法模拟这一变化。另外，在计算时所称的“内压”实际上是内部压强与环境压

强的压差，这时造成变形的有效载荷。

用轴对称壳单元模拟镀青铜充压弹簧，材质E=35GPa,4=0.33。弹簧顶部承受外载荷g=0.5MPa,

内压p=0.1MPa(加载前初始内压与环境的压差)，不计重力。计算最大Mises应力，最大z轴位移。

如下是该题按常内压求解的算例，ExpAbq05,命令流为ExpAbq05l.inp,本题作为密闭容器的变

内压解法见命令流ExpAbq()52.inp,因为变形很大，两种情况的求解结果大不一样。

文件与路径

顶部下拉菜单File,SaveAsExpAbq05»

一部件

创建部件，命名为Prat-1。

轴对称，可变形模型，线，图形大约范围l(m)。

选用线和圆弧绘出图示的截面壳的中线。造1层，拷贝5层，修正后完成Part-1。

二性质

1创建材料：

将材料命名为Material-1,输入E=35GPa,v=0.33,关闭。

2创建截面：

名为Section-1,壳，各向同性，材料名选Material-1,厚度0.005m,关闭。

3将截面的性质附加到部件上：

选中Prat-1,将Section-1信息注入Part-1。

三组装

ABAQUS计算指导ExpAbqO5

创建计算实体，以Prat-1为原形，用Independent方式生成实体。

四分析步

创建分析步，命名为Step-1,静态Static,通用General。

注释：无，时间：不变，非线性开关：关。

五载荷

1施加位移边界条件：

命名为BC-1,在分析步Step中，性质：力学，针对位移和转角。

选中轴上两线端点0、C,约束自由度X。

命名为BC-2,在分析步Step-1中，性质：力学，针对位移和转角。

选中底线。4,约束自由度y。

2创建载荷：

命名为Load-1,在分析步Step-1中，性质：力学，选择压强，选中顶线，选择外侧，值=0.5e6。

命名为Load-2,在分析步Step-1中，性质：力学，选择压强，选中各线，选择内侧，值=0.1e6。

六网格

对部件Prat-1进行。

1撒种子：

针对部件，全局种子大约间距().()05m。

2划网格：

针对部件，OK。

七建立项目

命名为ExpAbqO5,选择完整分析，其余不变，OK。

八观察结果

打开结果文件ExpAbq05.odb,看位移彩图。

在上述过程中不时地将本题的CAE模型保存为ExpAbq05.cae»

翁令瀛

*Heading

**Jobname:ExpAbqO51Modelname:Model-1

*Preprint,echo=NO,mode1=NO,history=NO,contact=NO

**

**PARTS

*Part,name=Part-1

*EndPart

**

**ASSEMBLY

*Assembly,name=Assembly

ABAQUS计算指导ExpAbqO5

*1nstance,name=Part-1-1,part=Part-1

*Node

1,0.0,0.0

4,*1，,0.一2,，().0

61,0.2,0.04

81,0.1,().04

101,0.1,0.08

121,0.2,0.08

141,0.2,0.12

161,0.1,0.12

181,0.1,0.16

201,0.2,0.16

221,0.2,0.2

241,0.1,0.2

261,0.1,0.24

281,0.2,0.24

3()1,0.2,0.28

321,0.1,0.28

341,0.1,0.32

361,0.2,0.32

381,0.2,0.36

421,0.0,0.36

*Ngen,,Nset=Partl

1,41,1♦，，，，，，

*Ngen,Line=C,Nset=Partl

41,61,1,,0.20,0.02,0.00,0.00,0.00,1.00

*Ngen,,Nset=Partl

61,81,1,,,,,,,

*Ngen,Line=C,Nset=Partl

81,101,1,,0.10,0.06,0.00,0.00,0.00,-1.00

*Ngen,,Nset=Partl

101,121,1,,,,,,,

*Ngen,Line=C,Nset=Partl

121,141,1,,0.20,0.10,0.00,0.00,0.00,1.00

*Ngen,,Nset=Partl

141,161,1,,,,,,,

*Ngen,Line=C,Nset=Par11

161,181,1,,0.10,0.14,0.00,0.00,0.00,-1.00

*Ngen,,Nset=Partl

181,201,1,,,,,,,

*Ngen,Line=C,Nset=Partl

ABAQUS计算指导ExpAbqO5

201,221,1,,0.20,0.18,0.00,0.00,0.00,1.00

*Ngen,,Nset=Partl

221,241,1,,,,,,,

*Ngen,Line=C,Nset=Partl

241,261,1,,0.10,0.22,0.00,0.00,0.00,-1.00

*Ngen,,Nset=Partl

261,281,1,,,,,,,

*Ngen,Line=C,Nset=Partl

281,30L1,,0.20,0.26,0.00,0.00,0.00,1.00

*Ngen,,Nset=Partl

301,321,1,,,,,,,

*Ngen,Line=C,Nset=Partl

321,341,1,,0.10,0.30,0.00,0.00,0.00,-1.00

*Ngen,,Nset=Partl

341,361,1,,,,,,,

*Ngen,Line=C,Nset=Partl

361,381,1,,0.20,0.34,0.00,0.00,0.00,1.00

*Ngen,,Nset=Partl

381,421,1,,,,,,,

*Element,type=SAX1

1,1,2

*ELGEN,ELSET=Part1

1,420,1,1

**Section:Section-1

*She11Section,elset=Part1,material=Material-1

0.005,5

*EndInstance

*Nset,nset=DUX,internaI,instance=Part-1-1

1,421

*Nset,nset=DUY,interna1,instance=Part-1-1

1,41,1

*EIset,elset=SurfQE,internal,instance=Part-1-1,generate

380,420,1

*Surface,type=ELEMENT,name=SurfQS,interna1

SurfQE,SNEG

*EIset,elset=SurfPE,internal,instance=Part-1-1,generate

1,420,1

*Surface,type=ELEMENT,name=SurfPS,interna1

SurfPE,SPOS

*EndAssembly

ABAQUS计算指导ExpAbqO5

**

**MATERIALS

*MateriaI,name=Material-1

*Elastic

3.5e+10,0.33

**__________________________________________

**STEP:Step-1

*Step,name=Step-1

*Static

1.,1.,Ie-05,1.

**BOUNDARYCONDITIONS

**Name:BC-lType:Displacement/Rotation

*Boundary

DUX,1,1

**Name:BC-2Type:Displacement/Rotation

*Boundary

DUY,2,2

**LOADS

**Name:Load-1Type:Pressure

*Dsioad

SurfQS,P,500000.

**Name:Load-2Type:Pressure

*Dsioad

SurfPS,R100000.

**

**OUTPUTREQUESTS

*Restart,write,frequency=0

**FIELDOUTPUT:F-0utput-1

*0utput,fieId,variable=PRESELECT

**HISTORYOUTPUT:H-Output-1

*0utput,history,variable=PRESELECT

*EndStep

*Heading

**Jobname:ExpAbqO52Modelname:Model-1

*Preprint,echo=NO,mode1=NO,history=NO,contact=NO

**

*PHYSICALCONSTANTS,ABSOLUTEZERO=-273.15

**

**PARTS

ABAQUS计算指导ExpAbqO5

*Part,name=Part-1

*EndPart

**

**ASSEMBLY

*Assembly,name=Assembly

*1nstance,name=Part-1-1,part=Part-1

*Node

1,0.0,0.0

41,0.2,0.0

61,0.2,0.04

81,0.1,0.04

101,0.1,0.08

121,0.2,0.08

141,0.2,0.12

161,0.1,0.12

181,0.1,0.16

201,0.2,0.16

221,0.2,0.2

241,0.1,0.2

261,0.1,0.24

281,0.2,0.24

301,0.2,0.28

321,0.1,0.28

341,0.1,0.32

361,0.2,0.32

381,0.2,0.36

421,0.0,0.36

*Ngen,,Nset=Partl

1,41,1，，，，，，，

*Ngen,Line=C,Nset=Partl

41,61,1,,0.20,0.02,0.00,0.00,0.00,1.00

*Ngen,,Nset=Partl

61,81,1，，，，，，，

*Ngen,Line=C,Nset=Par11

81J01,1,,0.10,0.06,0.00,0.00,0.00,-1.00

*Ngen,,Nset=Partl

101,121,1，，，，，，，

*Ngen,Line=C,Nset=Partl

121,141,1,,0.20,0.10,0.00,0.00,0.00,1.00

*Ngen,,Nset=Partl

141,161,1，，，，，，，

ABAQUS计算指导ExpAbqO5

*Ngen,Line=C,Nset=Partl

161,181,1,,0.10,0.14,0.00,0.00,0.00,-1.00

*Ngen,,Nset=Partl

181,201,1,,,,,,,

*Ngen,Line=C,Nset=Partl

201,221,1,,0.20,0.18,0.00,0.00,0.00,1.00

*Ngen,,Nset=Partl

221,241,1,,,,,,,

*Ngen,Line=C,Nset=Partl

241,261,1,,0.10,0.22,0.00,0.00,0.00,-1.00

*Ngen,,Nset=Partl

261,281,1,,,,,,,

*Ngen,Line=C,Nset=Par11

281,30h1,,0.20,0.26,0.00,0.00,0.00,1.00

*Ngen,,Nset=Partl

301,321,1,,,,,,,

*Ngen,Line=C,Nset=Partl

321,341,1,,0.10,0.30,0.00,0.00,0.00,-1.00

*Ngen,,Nset=Partl

341,36hl,,,,,,,

*Ngen,Line=C,Nset=Partl

361,381,1,,0.20,0.34,0.00,0.00,0.00,1.00

*Ngen,,Nset=Partl

381,421,1,,,,,,,

*EIement,type=SAX1

1,1,2

*ELGEN,ELSET=Part1

1,420,1,1

*Node

1000,0.00,0.18,0.00

*Nset,nset=CAVITY,intemal

1000,

**Section:Section-1

*She11Section,elset=Part1,material=Material-1

0.005,5

**FLUIDELEMENTS

*ELEMENT,TYPE=FAX2,ELSET=FLUID

501,1,2

*ELGEN,ELSET=FLUID

501,420,1,1

*FLUIDPROPERTY,REFNODE=CAVITY,TYPE=HYDRAULIC,

ABAQUS计算指导ExpAbqO5

ELSET=FLUID,AMBIENT=100000.00

*FLUIDDENSITY,PRESSURE=100000.00,TEMPERATURE=27.0

1.0432,

*EndInstance

**

*Nset,nset=CAVITY,internal,instance=Part-1-1

1000,

*Nset,nset=DUX,interna1,instance=Part-1-1