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