当前位置:文档之家› ABAQUS经验总结

ABAQUS经验总结

ABAQUS经验总结
ABAQUS经验总结

例如要取缸套的主推力侧上一列节点的径向位移(缸套轴线平行于Z轴),然后按照Z坐标画出径向位移曲线。总的方法是先选上要分析的节点,然后按空间位置顺序记录节点编号,根据编号生成一个Path,再以这个path为横坐标,以

要分析的量(如径向位移)为纵坐标画X-Y图。

1. 选出要分析区域节点:这列节点共有几十个,一个一个选太慢,打开要分析的odb文件,选择Display Group中的【Replace Selected】按钮,在目标类型中选择【Nodes】,然后结合视图方向、框选方法(矩形、圆形、多边形框)和选择方式(individually或者by angle)选出需要的节点集合,例子中的一列节点只有Z坐标不同,因此把视图方向调整到Z轴垂直于屏幕,这时一列点在屏幕上变成了一个点,在这一个点的位置上框选就可以得到一列点了。在individually 选择方式下,【shift+框选】是将选到的加入到已有选择中,【Ctrl+框选】是将选到的从已有选择中去除。例如要在一个曲面上选一列节点,可以先用【by angle】方式将整个曲面选上,然后调整到合适的视图方向上切换到【individually】

方式,用【Ctrl+框选】去掉多余的节点。

2. 按顺序记录节点编号:Path对顺序很敏感,节点号顺序的调整会改变最后plot的数据点顺序。要按照Z坐标大小顺序画出径向位移曲线,就必须按顺序记录节点编号。上一步选好节点按中键确定后,屏幕上是空白,因为ABAQUS 不会显示单独的节点,打开显示节点编号开关(在【Common Plot Options->Labels】中),这样屏幕上会显示出要选的那些节点的编号,按顺序记录下来。号码之间用逗号分隔,冒号表示连续和间隔,例如:1.2(1号和2号节点),1:10

(表示1,2,…,10),1:9:2(表示1,3,5,…,9)。

3. 生成目标表格和曲线,两种方法:a)用上面排好的节点序列生成一个Path,然后基于这个path生成一个X-Y图

(Create X-Y data -> Path),在【X-Y Data Manager】中双击生成的图线可以得到相应的二维表格;b)菜单【Report->Field Output】,这种方法将二维表格数据写入一个文件,本身不能出图,但是可以对计算结果求和,因此可以用来求接触力、压力的合力(但要注意这里的求和是数值相加,不是矢量求和,因此只适用于各个力方向基本一致的

情况下的合力估计)。

有时候历史变量输出的太多会出现

“The number of history output requests (22466) in this analysis step has exceeded the maximum value of 10000 specified by the Abaqus environment variable

'max_history_requests.'”错误,解决办法是:

在abaqus_v6.env中加一句“max_history_requests=0”即可。

abaqus_v6.env文件在C:\ABAQUS\6.7-1\site\下。

比如两个接触面之间有0.1mm的缝隙,在压力作用下发生接触,只要将Contact Controls下的稳定因子设为1即可正常计算(Interaction模块,主菜单Interaction->Contact Controls->Create->Stabilization->Automatic

Stabilization->factor设为1);

对于缝隙不均匀的情况,比如一端宽0.1mm,一端宽0.2mm,减小网格尺寸和稳

定因子对计算结果几乎没有影响;

[04/08] 关于ABAQUS历史输出变量

在帮助文档Analysis User’s Manual中4.2.2节Nodal variables部分定义了节点输出(用*Node Output关键字表示)可以包括的变量,从中可以看到应力(S)是不能作为节点计算结果输出的,只能作为单元积分点计算结果输出。原因是节点应力是根据单元积分点应力计算结果外插得到的,因此节点应力只能在场变量(而且是节点所在单元有输出)的情况下才能输出。如果在abaqus/cae中强行定义节点应力输出,INP文件中不会有相应的关键字,如果在INP文件中强行定义,求解时处理INP文件阶段会出错退出。

由于应力历史变量只能在单元积分点上输出,而在abaqus/cae中不能生成基于element型的set,除非直接使用网格模型做Part,但这样的模型往往不好操作,变通的方法是写由abaqus/cae生成INP文件,再手工修改这个INP文件,在其中添加一个element型的set,

然后让其输出应力历史变量。具体步骤:

1.记住要输出的单元编号,例如为27,在INP文件的*Assembly部分定义一个

element型的set:

*Elset, elset=Set-Name, instance=Part-Name

27,

2.在*End Step前加上:

*Output, history, frequency=10

*Element Output, elset=Set-Elem-forMeshPart

MISES,

[03/29] ABAQUS子程序VDLOAD学习笔记

[折叠]

使用的INP文件(beam.inp):

*Heading

load was added in cae

** Job name: Job-CAE Model name: Explicit-NoMove

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

**

** PARTS

**

*Part, name=beam

*Node

1, 0., 0.

2, 50., 0.

3, 100., 0.

4, 150., 0.

5, 200., 0.

6, 250., 0.

7, 300., 0.

8, 350., 0.

9, 400., 0.

10, 450., 0.

11, 500., 0.

*Element, type=B21

1, 1, 2

2, 2, 3

3, 3, 4

4, 4, 5

5, 5, 6

6, 6, 7

7, 7, 8

8, 8, 9

9, 9, 10

10, 10, 11

*Elset, elset=_PickedSet2, internal, generate

1, 10, 1

*Elset, elset=_PickedSet3, internal, generate

1, 10, 1

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

*Beam Section, elset=_PickedSet2, material=Material-1, poisson = 0.3, temperature=GRADIENTS, section=RECT

10., 5.

0.,0.,-1.

*End Part

**

**

** ASSEMBLY

**

*Assembly, name=Assembly

**

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

*End Instance

**

*Elset, elset=forLoad, instance=beam-1

6, 7

*Nset, nset=_PickedSet12, internal, instance=beam-1

1,

*Nset, nset=_PickedSet13, internal, instance=beam-1

11,

*Surface, type=ELEMENT, name=forLoadSPos, internal

forLoad, SPOS

*End Assembly

**

** MATERIALS

**

*Material, name=Material-1

*Damping, alpha=5.

*Density

7.8e-09,

*Elastic

210000., 0.3

**

** BOUNDARY CONDITIONS

**

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

*Boundary

_PickedSet12, 1, 1

_PickedSet12, 2, 2

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

*Boundary

_PickedSet13, 2, 2

** ----------------------------------------------------------------

**

** STEP: Step-1

**

*Step, name=Step-1

*Dynamic, Explicit

, 0.5

*Bulk Viscosity

0.06, 1.2

**

** LOADS

**

** Name: Load-1 Type: Concentrated force

*Dsload

forLoadSPos, PNU, 1.

**

** OUTPUT REQUESTS

**

*Restart, write, number interval=1, time marks=NO

**

** FIELD OUTPUT: F-Output-1

**

*Output, field, time interval=0.005

*Node Output

U,

*Element Output, directions=YES

S,

**

** HISTORY OUTPUT: H-Output-1

**

*Output, history, variable=PRESELECT

*End Step

使用的for文件(Gas.for):

subroutine vdload (

1 nblock, ndim, stepTime, totalTime,

1 amplitude, curCoords, velocity, dirCos, jltyp, sname,value )

C

include 'vaba_param.inc'

C

dimension curCoords(nblock,ndim), velocity(nblock,ndim),

1 dirCos(nblock,ndim,ndim), value(nblock)

character*80 sname

C

open(unit=17,file='e:\aa.txt',status='old',position='append')

write(17,*)'stepTime= ',stepTime

write(17,*)'ndim= ',ndim

write(17,*)'curCoords= ',curCoords

write(17,*)'jltyp= ',jltyp

value = 1

write(17,*)'sname= ',sname

write(17,*)'value= ',value

write(17,*)' '

close(17)

return

end

命令行:abaqus job=beam user=gas interactive

分析:

ABAQUS/Explicit每个增量步(Increment)都会调用用户子程序,这里就是gas.for,因此上面打开/关闭文件的动作会进行很多次,为避免后面打开文件将原有的内容删掉,就必须显式地指明打开的是一个已存在的文件(status='old'),写文件的位置在原文件末尾(position='append')。nblock表示载荷作用区域内对象的个数,当载荷类型为Pressure,作用在一些单元的面上时,nblock就是单元的数量。

stepTime表示当前时间,就是每个增量步发生时候的时间,时间起点为当前载

荷步,如果要总时间,用totalTime。

当模型为二维时,ndim值为2,三维时ndim值为3。

curCoords(nblock, ndim)是一个二维数组。curCoords(ki, kj)表示当前增量步时间下,序号为ki的点的位置坐标的序号为kj的坐标分量。在二维梁模型中,“点”表示有载荷作用的单元的中点, ki取值范围为1~nblock,kj取值范围

为1~ndim。

jltyp表示载荷的类型,与INP文件中的定义相一致,具体参见ABAQUS文档。sname表示当载荷类型为surface-based load(jltyp=0)时,载荷作用面的名称,与INP文件中的定义相一致,具体参见ABAQUS文档。

[03/24] ABAQUS用户子程序安装过程

软件版本

Microsoft.Visual.Studio.2005.Professional.Edition.DVD

https://www.doczj.com/doc/c36155360.html,piler.Pro.v9.1.032

ABAQUS6.7

安装顺序

原则是保证VC2005在Intel Fortran之前装,因为这样装Intel Fortran与VC2005集成的组件才有效,可以按VC2005-Intel Fortran-ABAQUS或者ABAQUS - VC2005-Intel Fortran顺序装,VC2005只装VC++即可。

安装过程

分为三个步骤:先装软件,再修改环境变量,最后是验证安装是否成功。

1)机器上已经安装了abaqus6.7。安装VC2005时选择custom方式,只装VC++,其他的诸如SQL Server之类的都不需要装,也不需要更改

VC2005默认的安装路径;安装Intel Fortran时选择所有项目;

2)更改环境变量PATH、INCLUDE和LIB,方法是在各自原有的内容后增

加:

PATH:C:\Program Files\Microsoft Visual Studio 8\Common7\IDE;

C:\Program Files\Microsoft Visual Studio 8\VC\bin;

C:\Program Files\Microsoft Visual Studio 8\Common7\Tools;

C:\Program Files\Microsoft Visual Studio 8\SDK\v2.0\bin;

C:\Program Files\Microsoft Visual Studio 8\VC\VCPackages;

C:\Program Files\Microsoft Visual Studio 8\VC\PlatformSDK\bin

LIB:C:\Program Files\Microsoft Visual Studio 8\VC\lib;C:\Program Files\Microsoft Visual Studio 8\SDK\v2.0\Lib;C:\Program Files\Microsoft

Visual Studio 8\VC\PlatformSDK\Lib

INCULDE:C:\Program Files\Microsoft Visual Studio

8\VC\include;C:\Program Files\Microsoft Visual Studio

8\VC\PlatformSDK\Include;

注:LINK_F90这个变量不需要增加;上面的都假设VC2005装在了C盘默认目录

下。

3)验证:

a)用记事本生成下面两个文件:

bbb.for文件

SUBROUTINE DLOAD(F,KSTEP,KINC,TIME,NOEL,NPT,LAYER,KSPT,COORDS,

1 JLTYP,SNAME)

C

INCLUDE 'ABA_PARAM.INC'

C

DIMENSION TIME(2), COORDS (3)

CHARACTER*80 SNAME

WRITE(6,*)'INCREMENT NUMBER',KINC

WRITE(6,*)'STEP TIME',TIME(1)

WRITE(6,*)'ELEMENT NUMBER',NOEL

WRITE(6,*)'INTEGRATION POINT NUMBER',NPT

WRITE(6,*)'COORDS:X',COORDS(1)

F=1.0*COORDS (1)

WRITE(6,*)'FORCE',F

WRITE(6,*)' '

RETURN

END

beam.inp文件

*Heading

** Job name: Job-1 Model name: beam

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

**

** PARTS

**

*Part, name=PART-1

*End Part

**

** ASSEMBLY

**

*Assembly, name=Assembly

**

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

*Node

1, 0., 0.

2, 20., 0.

3, 40., 0.

4, 60., 0.

5, 80., 0.

6, 100., 0.

*Element, type=B31

1, 1, 2

2, 2, 3

3, 3, 4

4, 4, 5

5, 5, 6

*Elset, elset=BEAM, generate

1, 5, 1

** Region: (Section-1-BEAM:BEAM), (Beam Orientation:BEAM)

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

*Beam Section, elset=BEAM, material=STEEL, temperature=GRADIENTS, section=RECT

0.2, 5.

0.,0.,-1.

*End Instance

*Nset, nset=ENDS, instance=PART-1-1

1, 6

*Nset, nset=_M4, internal, instance=PART-1-1

6,

*Nset, nset=_M5, internal, instance=PART-1-1

1,

*End Assembly

**

** MATERIALS

**

*Material, name=STEEL

*Elastic

210000., 0.3

*ELSET,ELSET=PY

BEAM

**

** BOUNDARY CONDITIONS

**

** Name: Disp-BC-1 Type: Symmetry/Antisymmetry/Encastre

*Boundary

_M4, ENCASTRE

** ----------------------------------------------------------------

**

** STEP: Step-1

**

*Step, name=Step-1

*Static

**

** LOADS

**

** Name: CFORCE-1 Type: Concentrated force

*DLOAD

PY,PYNU

**

** OUTPUT REQUESTS

**

**

** FIELD OUTPUT: F-Output-1

**

*Output, field, variable=PRESELECT

**

** FIELD OUTPUT: F-Output-2

**

*Output, field

*Element Output

SF,

**

** HISTORY OUTPUT: H-Output-1

**

*Output, history

*Node Output, nset=ENDS

CF1, CF2, CF3, CM1, CM2, CM3, RF1, RF2

RF3, RM1, RM2, RM3, U1, U2, U3, UR1

UR2, UR3

*El Print, freq=999999

*Node Print, freq=999999

*End Step

b)在Abaqus Command提示符后输入:

abaqus job=beam user=bbb interactive

运行完毕并没有提示出错后,当前目录下应该出现文件名为beam,扩展名不同的几个文件,打开beam.dat文件,可以找到bbb.for文件中WRITE语句指定的输出内容,详见笔记“ABAQUS用户子程序学习小结”

[01/19] 用ABAQUS进行多体分析的步骤

[折叠]

第一阶段:总体分析

1.画草图:上面需要表明的信息包括(参考[石亦平2006]图8-19):

a)各个部件(Part)的名称;

b)各个部件上每一个参考点的名称和位置坐标,一般就是与其他部

件连接的点,部件包括柔体、刚体和显示刚性体(本身为柔体,通过定义一个参考点,并施加显示体约束生成);

c)各个连接单元的名称;

d)每个连接单元要用到的两个参考点的名称;

e)每个连接单元要用到的局部坐标系的名称、位置坐标和方向(便

于后面输入坐标直接建立局部坐标系);

2.自由度分析:保证系统中不会出现过约束或者欠约束,步骤参考[石

亦平2006]p261“连接单元类型和边界条件”;

3.建立每个部件:位置由第1.b)步给出;

4.定义材料属性:显示刚性体改为柔体,可以不定义材料属性;

5.装配:一般情况下各个part的位置都是装配后的位置,如果需要移

动,可以用【Instance-Translate】来移动;

6.划网格;

第二阶段:刚体分析

目的:检验系统中是否存在过约束和欠约束;

原则:刚体加显示约束,未来柔体(在第三阶段要转化为柔体的Part)加刚体约

束。

7.定义各个参考点:按照第1.b)步生成各个参考点,如果模型比较复杂,

最好为每个参考点建立一个Set,这样后面建立约束和连接单元都会

很方便;

8.建立约束:

a)建立柔体与参考点(如A)之间的显示体约束,以生成各个显示

刚性体;

b)如果体上除A外还有其他参考点(如B),各点之间要建立刚体

约束,B的Region Type选“Tie(nodes)”;

9.定义连接单元的基准坐标系:坐标原点无所谓,关键是坐标轴方向;

10.定义连接单元:ABAQUS中建立运动副的方法是:有Part A和Part B,

质心参考点分别为ac和bc,在运动副处二者的参考点分别为arp和

brp,首先建立(ac和arp)、(bc和brp)的刚体约束,然后建立arp

和brp的连接单元。对于未来柔体,不要建立arp和brp,而要使用体

上已经存在的点(Vertex),这样第三阶段直接删除刚体约束即可进

行柔体分析。

a)定义连接属性:建立各个Connector Section;

b)定义连接对(Wire):就是建立连接各对参考点的Wire:

c)建立连接单元(Connector Assignment):在这里配对连接属性

和wire,同时指定第一个参考点使用的坐标系;

11.确定分析步类型和数量:如果不关心分析目标的动态响应(例如进行

强度分析),则使用静态分析步,但几何非线性要设为On;

12.设置历史变量输出:domain设为代表某个运动副的Set,输出变量设

为Connector中的CRF和CU;

13.定义边界条件(和载荷):按照第2步中的分析定义边界条件;

14.提交作业求解;

15.后处理:

a)用【Query-Probe Values】,输出场变量设为U,可以看到各个

参考点的位移量(参考点被作为装配体(Assembly)的节点(node),

【Probe Values-Key-in label-Part instance:Assembly】可以依次输

入各个参考点求其坐标,配合【Step/Frame】,可以看到各个迭

代步中参考点的位移),如果仅进行刚体分析,没有施加载荷,

也没有过约束,各个部件的node的位移值都不会显示;

b)【Result-History Output】可以看到第12步中设置的输出变量结

果,如CU是运动副上两个点的相对位移,CRF是两个点的反作用

力等等;

第三阶段:柔体分析

16.例如要分析Part A和Part B上的应力分布,则删除Part A和Part B上

的刚体约束,给Part A和Part B建立各自的Set,在场变量输出中增加

Set-Part-A和Set-Part-B的相关变量(例如应力、位移等),重新提交

分析。

[2007/10/29] ABAQUS用户子程序学习小结

1.1 FORTRAN语言不区分大小写,文件扩展名为for,正文从第7列开始,第1~5列是标号区,第六列是续行标志区。

1.2 FORTRAN语言中的“I-N规则”:I、J、K、L、M、N开头的为整型变量,其他开头为实型变量;

2 DIMENSION COORDS(3)表示声明一个含3个元素的数组,下标分别为1、2、3,访问形式为COORDS(n),n为1

~3;

3 子程序(*.for)文件中如何输出调试信息:WRITE(6,*)'COORDS(1)',COORDS(1),在*.dat文件中可看到输出,如果

希望WRITE输出到msg文件中,则写为WRITE(7,*)'COORDS...;

4 用户子程序DLOAD中COORDS数组的含义:COORDS(1)也是一个数组,存贮单元集合中所有单元积分点的X坐

标,COORDS(2)存贮Y坐标,相应INP文件中的写法为:

*DLOAD

PY,PYNU

其中PY为单元集合名称,定义方法为:

*Elset, elset=BEAM, generate

1, 5, 1

...

*ELSET,ELSET=PY

BEAM

5 DLOAD中F的定义方法:

F只有定义在单元积分点上才有效,例如:

F=1.0*COORDS (1)

附一个简单实例:

beam.inp文件:

*Heading

** Job name: Job-1 Model name: beam

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

**

** PARTS

**

*Part, name=PART-1

*End Part

**

** ASSEMBLY

**

*Assembly, name=Assembly

**

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

*Node

1, 0., 0.

2, 20., 0.

3, 40., 0.

4, 60., 0.

5, 80., 0.

6, 100., 0.

*Element, type=B31

1, 1, 2

2, 2, 3

3, 3, 4

4, 4, 5

5, 5, 6

*Elset, elset=BEAM, generate

1, 5, 1

** Region: (Section-1-BEAM:BEAM), (Beam Orientation:BEAM)

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

*Beam Section, elset=BEAM, material=STEEL, temperature=GRADIENTS, section=RECT

0.2, 5.

0.,0.,-1.

*End Instance

*Nset, nset=ENDS, instance=PART-1-1

1, 6

*Nset, nset=_M4, internal, instance=PART-1-1

6,

*Nset, nset=_M5, internal, instance=PART-1-1

1,

*End Assembly

**

** MATERIALS

**

*Material, name=STEEL

*Elastic

210000., 0.3

*ELSET,ELSET=PY

BEAM

**

** BOUNDARY CONDITIONS

**

** Name: Disp-BC-1 Type: Symmetry/Antisymmetry/Encastre

*Boundary

_M4, ENCASTRE

** ----------------------------------------------------------------

**

** STEP: Step-1

**

*Step, name=Step-1

*Static

**

** LOADS

**

** Name: CFORCE-1 Type: Concentrated force

*DLOAD

PY,PYNU

**

** OUTPUT REQUESTS

**

**

** FIELD OUTPUT: F-Output-1

**

*Output, field, variable=PRESELECT

**

** FIELD OUTPUT: F-Output-2

**

*Output, field

*Element Output

SF,

**

** HISTORY OUTPUT: H-Output-1

**

*Output, history

*Node Output, nset=ENDS

CF1, CF2, CF3, CM1, CM2, CM3, RF1, RF2

RF3, RM1, RM2, RM3, U1, U2, U3, UR1

UR2, UR3

*El Print, freq=999999

*Node Print, freq=999999

*End Step

bbb.for文件

SUBROUTINE DLOAD(F,KSTEP,KINC,TIME,NOEL,NPT,LAYER,KSPT,COORDS,

1 JLTYP,SNAME)

C

INCLUDE 'ABA_PARAM.INC'

C

DIMENSION TIME(2), COORDS (3)

CHARACTER*80 SNAME

WRITE(6,*)'COORDS(3)',COORDS(3)

F=1.0*COORDS (1)

RETURN

END

运行方法:

在Abaqus Command提示符后输入:

abaqus job=beam user=bbb interactive

阅读(199)

分享(0)

评论(2)

分类:CAE

发表于16:40

[2007/10/22] 推荐一本有限元方面的好书

阅读(85)

分享(0)

评论(0)

分类:CAE

发表于14:33

上页

14/16

下页

[2007/10/16] Abaqus/cae生成inp文件的结构

[折叠]

*Part, name=Part1

*Node

*Element

*Nset

*Elset

*Solid Section(定义截面属性)

*End Part

(重复以上*Part节,定义下一个Part)

*Assembly name=Assm1

*Instance, name=A

(Assembly中的各个Instance在此定义)

*Nset, nset=Set1

*Elset, elset=Set1

(重复这2行,定义下一个Set,似乎只用Nset就可以定义)

*Elset, elset=_Surf-1_S1

1, 100, 1 (数据行,定义这个Elset)

*Surface, type=ELEMENT, name=Surf-1

_Surf-1_S1, S1

(重复这4行,定义下一个Surface,第4行是用“元素集合+面标识号”的方法定义Surface)

*End Assembly

*Material (定义材料特性)

*Surface Interaction (定义接触属性)

*Boundary

Set-1, ENCASTRE (定义边界条件)

*Contact Pair, interaction=IntProp-1

Surf-1, Surf-2 (定义接触对)

*Step, name=Step1

......

*End Step

(重复,定义下一个Step)

[2007/10/16] ABAQUS中如何定义非均布时变载荷

[折叠]

Analytical Field、Amplitude与Load结合,可以在一个面上定义任意形式的非均布时变载荷;

1. 菜单操作:Load或Interaction模块:Tools->Analytical Field;

2. 坐标系的定义

a) 在【Edit Expression Field】对话框中的【Local system】定义要使用的坐标系;

b) 表达式中的坐标分量都以此坐标系原点为基准;

3. 在Abaqus Document(6.6/6.7)【ABAQUS/CAE User's Manual->Part VI: Using Toolsets->Chapter39

The Analytical Field toolset】中对Analytical Field的使用进行了详细的说明;

4. 构建表达式使用Python语法,可使用的运算符都包括在【Edit Expression Field】对话框的【Operators】

文本框内;

5. 在面上施加非均布时变压力的步骤是:

a) 在Load模块中用Analytical Field定义好表达式,默认取名为AnalyticalField-1;

b) 在【Tools->Amplitude】中创建一个随时间变化的函数曲线(详见Chapter 38 The Amplitude toolset),

默认名称为Amp-1;

c) Create Load,在【Distribution】中选择AnalyticalField-1,在【Magnitude】中输入一个数,在【Amplitude】中选择Amp-1,即可完成对该压力的定义;【Distribution】中的默认选项Uniform表示均布压力,【Amplitude】中的

Ramp表示按时间线性加载;

6. 定义好的压力不是以函数形式写入INP文件中,而是将面上单元中心点的位置坐标代入函数表达式,算

出该点上的压力值,再以*Dload关键字写入INP文件中;

7. 圆柱坐标(坐标分量为R、Th、Z,其中Th是角度)和球坐标(坐标分量为R、Th、P,其中Th和P

都是角度)中,Th和P的范围都是(-pi,pi);

8. 表达式中大小写表示不同含义;

9. 【View->Assembly Display Options->Attribute->Symbol】中可以控制载荷图标的尺寸和密度,对于一些

简单的分布形式,如按单一坐标分量线性分布(3*y之类),可以直接在图形显示区看到载荷的分布特点;

abaqus压杆屈曲分析

a b a q u s压杆屈曲分析 Revised by Petrel at 2021

压杆屈曲分析1.问题描述 在钢结构中,受压杆件一般在其达到极限承载力前就会丧失稳定性,所以失稳是钢结构最为突出的问题。压杆整体失稳形式可以是弯曲、扭转和弯扭。钢构件在轴心压力作用下,弯曲失稳是常见的失稳形式。影响轴心受压构件整体稳定性的主要因素为纵向残余应力、初始弯曲、荷载初偏心及端部约束条件等。实际的轴心受压构件往往会存在上述的一种或多种缺陷,导致构件的稳定承载力降低。 本文利用abaqus对一定截面不同长细比下的H型钢构件进行屈曲分析,通过考虑材料非线性、几何非线性并引入初弯曲,得出构件发生弯曲失稳的极限荷载。通过比较不同长细比下的弯曲失稳的临界荷载得出构件荷载位移曲线,并与《规范》中的构件曲线相比较。钢构件的截面尺寸如图1-1所示。 构件的材料特性:,, 图1-1 2.长细比计算 通过计算截面几何特性,截面绕y轴的回转半径为,长细比取值及杆件长度见表1: 表1 50 60 80 100 120 150 180 (m) 1.92 2.30 3.07 3.84 4.60 5.76 6.90 3.模型分析

ABAQUS非线性屈曲分析的方法有riks法,generalstatics法(加阻尼),或者动力法。非线性屈曲分析采用riks算法实现,可以考虑材料非线性、几何非线性已及初始缺陷的影响。其中,初始缺陷可以通过屈曲模态、振型以及一般节点位移来描述。 利用abaqus进行屈曲分析,一般有两步,首先是特征值屈曲分析,此分析为线性屈曲分析,是在小变形的情况进行的,也即上面提到过的模态,目的是得出临界荷载(一般取一阶模态的eigenvalue乘以所设定的load)。其次,就是后屈曲分析,此步一般定义为非线性,原因在于是在大变形情况进行的,一般采用位移控制加修正的弧长法,可以定义材料非线性,以及几何非线性,加上初始缺陷,所以也称为非线性屈曲分析。此步分析,为了得到极限值,需要得出荷载位移曲线的下降段。缺陷较小的结构初始位移变形较小,在极值点突变,而初始缺陷较大的结构,载荷位移曲线较平滑。 4.建模计算过程 建模计算过程以长细比为50的构件为例,其余构件建模计算过程与之类似。 4.1buckle分析 1在buckle分析中创建part模块,创建的模型为三位可变形壳体单元,截面参数见图1-1,构件长度1.92。如图4-1示 图4-1 2定义材料特性及截面属性并将其赋予单元。材料定义为弹塑性,泊松比0.3,屈服强度,弹性模量;腹板和翼缘板为壳单元,厚度分别为0.008和0,01。材料定义见图4-2

基于ABAQUS的钢管轴心受压非线性屈曲分析

一.问题描述 在钢结构中,受压杆件一般在其达到极限承载力前就会丧失稳定性,所以失稳是钢结构最为突出的问题。压杆整体失稳形式可以是弯曲、扭转和弯扭。钢构件在轴心压力作用下,弯曲失稳是常见的失稳形式。而影响轴心受压构件整体稳定性的主要因素为纵向残余应力、初始弯曲、荷载初偏心及端部约束条件等。实际的轴心受压构件往往会存在上述的一种或多种缺陷,导致构件的稳定承载力降低。 本文主要针对任意轴对称的圆形钢管截面,利用ABAQUS有限元非线性分析软件,对其在轴心受压情况下进行特征值屈曲分析和静态及动态的非线性屈曲分析(考虑材料弹塑性和初始缺陷的影响)。通过考虑材料非线性、几何非线性并引入初弯曲,得出构件发生弯曲失稳的极限荷载,并且由弯曲失稳的临界荷载得出的构件荷载位移曲线。同时再进行非线性分析时,需要施加初始扰动,以帮助非线性分析时失稳,可以通过特征值屈曲分析得到的初始弯曲模态来定义初始缺陷;最后由可以将特征值屈曲分析得到的临界荷载作为非线性屈曲分析时所施加荷载的参考。 二.结构模型 用ABAQUS中的壳单元建立轴心受压模型,采用SI国际单位制(m)。 1.构件的材料特性: E= 2.0×1011N m2,μ=0.3, f y=2.35×

108N m2,ρ=7800kg m3,钢管半径:60mm,厚度:3mm,长度:2.5m。 2.钢管的截面尺寸及钢管受到的约束和荷载施加的模型图如图2-1及图2-2所示。 图2-1 图2-2 三.建模步骤(Buckle分析) (1)创建部件 在创建part模块中命名构件的名字为gang guan,创建的模型为三维可变形壳体单元,如图3-1所示。截面参数见图2-1,构件长度2.5m。 图3-1

采用ABAQUS进行屈曲后屈曲和破坏分析

| w w w .3d s .c o m | ? D a s s a u l t S y s t èm e s | Buckling, Postbuckling, and Collapse Analysis with Abaqus | w w w .3d s .c o m | ? D a s s a u l t S y s t èm e s | Day 1 ?Lecture 1Basic Concepts and Overview ?Workshop 1Buckling and Postbuckling Analyses of a Crane Structure ?Lecture 2 Finite Element Formulation ?Lecture 3Finite Element Implementation in Abaqus ?Lecture 4Eigenvalue Buckling Analysis ?Workshop 2Eigenvalue Buckling of a Ring Subjected to External Pressure ?Workshop 3 Elastic Buckling of Ring-Supported Cylindrical Shell under Hydrostatic Pressure

| w w w .3d s .c o m | ? D a s s a u l t S y s t èm e s | Buckling, Postbuckling, and Collapse Analysis with Abaqus Day 2 ?Lecture 5 Regular and Damped Static Solution Procedures for Postbuckling Analyses ?Workshop 4Nonlinear Buckling of Ring-Supported Cylindrical Shell under Hydrostatic Pressure ?Workshop 5Static Buckling Analysis of a Circular Arch ?Lecture 6Modified Riks Static Solution Procedure for Postbuckling Analyses ?Workshop 5Static Buckling Analysis of a Circular Arch (continued)?Lecture 7Dynamic Analysis Solution Procedures for Postbuckling Analyses ?Workshop 5Static Buckling Analysis of a Circular Arch (continued)?Workshop 6Tube Crush Dynamic Analysis ?Lecture 8Putting It All Together… ?Workshop 7Capstone Workshop: Lee’s Frame Buckling Problem ?Workshop 8 Buckling and Postbuckling Analyses of a Stiffened Panel | w w w .3d s .c o m | ? D a s s a u l t S y s t èm e s | Legal Notices The Abaqus Software described in this documentation is available only under license from Dassault Systèmes and its subsidiary and may be used or reproduced only in accordance with the terms of such license. This documentation and the software described in this documentation are subject to change without prior notice. Dassault Systèmes and its subsidiaries shall not be responsible for the consequences of any errors or omissions that may appear in this documentation. No part of this documentation may be reproduced or distributed in any form without prior written permission of Dassault Systèmes or its subsidiary.? Dassault Systèmes, 2011. Printed in the United States of America Abaqus, the 3DS logo, SIMULIA and CATIA are trademarks or registered trademarks of Dassault Systèmes or its subsidiaries in the US and/or other countries. Other company, product, and service names may be trademarks or service marks of their respective owners. For additional information concerning trademarks, copyrights, and licenses, see the Legal Notices in the Abaqus 6.11 Release Notes and the notices at: https://www.doczj.com/doc/c36155360.html,/products/products_legal.html.

本人学习abaqus五年的经验总结 让你比做例子快十倍

第二章 ABAQUS 基本使用方法 [2](pp15)快捷键:Ctrl+Alt+左键来缩放模型;Ctrl+Alt+中键来平移模型;Ctrl+Alt+右键来旋转模型。 ②(pp16)ABAQUS/CAE 不会自动保存模型数据,用户应当每隔一段时间自己保存模型以避免意外 丢失。 [3](pp17)平面应力问题的截面属性类型是Solid(实心体)而不是Shell(壳)。ABAQUS/CAE 推荐的建模方法是把整个数值模型(如材料、边界条件、载荷等)都直接定义在几 何模型上。 载荷类型Pressure 的含义是单位面积上的力,正值表示压力,负值表示拉力。 [4](pp22)对于应力集中问题,使用二次单元可以提高应力结果的精度。 [5](pp23)Dismiss 和Cancel 按钮的作用都是关闭当前对话框,其区别在于:前者出现在包含只读数 据的对话框中;后者出现在允许作出修改的对话框中,点击Cancel 按钮可关闭对话框,而不保存 所修改的内容。 [6](pp26)每个模型中只能有一个装配件,它是由一个或多个实体组成的,所谓的“实体”(instance) 是部件(part)在装配件中的一种映射,一个部件可以对应多个实体。材料和截面属性定义在部件 上,相互作用(interaction)、边界条件、载荷等定义在实体上,网格可以定义在部件上或实体上, 对求解过程和输出结果的控制参数定义在整个模型上。 [7](pp26) ABAQUS/CAE 中的部件有两种:几何部件(native part)和网格部件(orphan mesh part)。 创建几何部件有两种方法:(1)使用Part 功能模块中的拉伸、旋转、扫掠、倒角和放样等特征来直 接创建几何部件。(2)导入已有的CAD 模型文件,方法是:点击主菜单 File→Import→Part。网 格部件不包含特征,只包含节点、单元、 面、集合的信息。创建网格部件有三种方法:(1)导入 ODB 文件中的网格。(2)导入INP 文件中的网格。(3)把几何部件转化为网格部件,方法是:进 入Mesh 功能模块,点击主菜单Mesh→Create Mesh Part。 [8](pp31)初始分析步只有一个,名称是initial,它不能被编辑、重命名、替换、复制或删除。在初 始分析步之后,需要创建一个或多个后续分析步,主要有两大类:(1)通用分析步(general analysis step)可以用于线性或非线性分析。常用的通用分析步包含以下类型:—Static, General: ABAQUS/Standard 静力分析 —Dynamics, Implicit: ABAQUS/Standard 隐式动力分析 —Dynamics, Explicit: ABAQUS/ Explicit 显式动态分析

abaqus压杆屈曲分析78112

压杆屈曲分析 1.问题描述 在钢结构中,受压杆件一般在其达到极限承载力前就会丧失稳定性,所以失稳是钢结构最为突出的问题。压杆整体失稳形式可以是弯曲、扭转和弯扭。钢构件在轴心压力作用下,弯曲失稳是常见的失稳形式。影响轴心受压构件整体稳定性的主要因素为纵向残余应力、初始弯曲、荷载初偏心及端部约束条件等。实际 2 压杆截面尺寸(单位:m) 图1-1 2.长细比计算 通过计算截面几何特性,截面绕y轴的回转半径为i y=0.0384m ,长细比取

值及杆件长度见表1: 表1 3.模型分析 ABAQUS非线性屈曲分析的方法有riks法,general statics法(加阻尼),或者动力法。非线性屈曲分析采用riks算法实现,可以考虑材料非线性、几何非线性已及初始缺陷的影响。其中,初始缺陷可以通过屈曲模态、振型以及一般节点位移来描述。 利用abaqus进行屈曲分析,一般有两步,首先是特征值屈曲分析,此分析为线性屈曲分析,是在小变形的情况进行的,也即上面提到过的模态,目的是得出临界荷载(一般取一阶模态的eigenvalue乘以所设定的load)。其次,就是后屈曲分析,此步一般定义为非线性,原因在于是在大变形情况进行的,一般采用位移控制加修正的弧长法,可以定义材料非线性,以及几何非线性,加上初始缺陷,所以也称为非线性屈曲分析。此步分析,为了得到极限值,需要得出荷载位移曲线的下降段。缺陷较小的结构初始位移变形较小,在极值点突变,而初始缺陷较大的结构,载荷位移曲线较平滑。 4.建模计算过程 建模计算过程以长细比为50的构件为例,其余构件建模计算过程与之类似。 4.1 buckle分析 1 在buckle分析中创建part模块,创建的模型为三位可变形壳体单元,截面参数见图1-1,构件长度1.92。如图4-1示

abaqus压杆屈曲分析63758

压杆屈曲分析 1.问题描述 在钢结构中,受压杆件一般在其达到极限承载力前就会丧失稳定性,所以失稳是钢结构最为突出的问题。压杆整体失稳形式可以是弯曲、扭转和弯扭。钢构件在轴心压力作用下,弯曲失稳是常见的失稳形式。影响轴心受压构件整体稳定性的主要因素为纵向残余应力、初始弯曲、荷载初偏心及端部约束条件等。实际的轴心受压构件往往会存在上述的一种或多种缺陷,导致构件的稳定承载力降低。 本文利用abaqus 对一定截面不同长细比下的H 型钢构件进行屈曲分析,通过考虑材料非线性、几何非线性并引入初弯曲,得出构件发生弯曲失稳的极限荷载。通过比较不同长细比下的弯曲失稳的临界荷载得出构件荷载位移曲线,并与《规范》中的构件曲线相比较。钢构件的截面尺寸如图1-1所示。 构件的材料特性: E =2.0×1011 N m 2? ,μ=0.3 , f y =3.45×108N m 2? 压杆截面尺寸(单位:m)

图1-1 2.长细比计算 通过计算截面几何特性,截面绕y轴的回转半径为i y=0.0384m ,长细比取值及杆件长度见表1: 表1 λ50 60 80 100 120 150 180 ι(m) 1.92 2.30 3.07 3.84 4.60 5.76 6.90 3.模型分析 ABAQUS非线性屈曲分析的方法有riks法,general statics法(加阻尼),或者动力法。非线性屈曲分析采用riks算法实现,可以考虑材料非线性、几何非线性已及初始缺陷的影响。其中,初始缺陷可以通过屈曲模态、振型以及一般节点位移来描述。 利用abaqus进行屈曲分析,一般有两步,首先是特征值屈曲分析,此分析为线性屈曲分析,是在小变形的情况进行的,也即上面提到过的模态,目的是得出临界荷载(一般取一阶模态的eigenvalue乘以所设定的load)。其次,就是后屈曲分析,此步一般定义为非线性,原因在于是在大变形情况进行的,一般采用位移控制加修正的弧长法,可以定义材料非线性,以及几何非线性,加上初始缺陷,所以也称为非线性屈曲分析。此步分析,为了得到极限值,需要得出荷载位移曲线的下降段。缺陷较小的结构初始位移变形较小,在极值点突变,而初始缺陷较大的结构,载荷位移曲线较平滑。 4.建模计算过程

基于ABAQUS复合材料薄壁圆筒的屈曲分析

基于ABAQUS复合材料薄壁圆筒的屈曲分析 由于玻璃钢复合材料的薄壁圆筒结构具有强度高、重量轻、刚度大、耐腐蚀,电绝缘及透微波等优点,目前已广泛应用于航空航天和民用领域中。工程中广泛使用的这些薄壁圆筒,当它们受压缩、剪切、弯曲和扭转等荷载作用时,最常见的失效模式为屈曲。因此,为了保证结构的安全,需要进行屈曲分析。 对结构进行屈曲分析,涉及到较复杂的弹(塑)性理论和数学计算,要通过求解高阶偏微分方程组,才能求解失稳临界荷载,而且只有少数简单结构才能求得精确的解析解。因此,只能采用能量法、数值方法和有限元方法等近似的分析方法进行分析。近20年来,随着计算机和有限元方法的迅猛发展,形成了许多的实用分析程序,提高了对复杂结构进行屈曲分析的能力和设计水平。ABAQUS 就是其中的杰出代表。 1.屈曲有限元理论 有限元方法中,对结构的屈曲失稳问题的分析方法主要有两类:一类是通过特征值分析计算屈曲载荷,另一类是利用结合Newton—Raphson迭代的弧长法来确定加载方向,追踪失稳路径的几何非线性分析方法,能有效分析高度非线性屈曲和后屈曲问题。 1.1线性屈曲 假设结构受到的外载荷模式为P0。,幅值大小为λ,结构内力为Q,则静力平衡方程应为 λP0=λQ 进一步考察结构在(λ+△λ)P0载荷作用下的平衡方程,得到 {[K E]+[K S(S+λ△S)]+[K G(u?+λu?)]}△u?=△λP0由于结构达到保持稳定的临界载荷时有△λ,代入上式得 {[K E]+λ[K S△σ]+K G(△u?)}△u?=0 该方程对应的特征值问题为 det{[K E]+λ[K S△σ]+K G(△u?)}=0 如果忽略几何刚度增量的影响,屈曲分析的方程又可进一步简化为 det{[K E]+λ[K S△σ]}=0 该方程即为求解线性屈曲的特征值方程。λ为屈曲失稳载荷因子,(△u?)为结构失稳形态的特征向量。

ABAQUS非线性屈曲分析步骤

ABAQUS6.7非线性屈曲分析步骤 riks法,或者general statics法(加阻尼),或者动力法 一共三种方法, 【问】在aba中能实现非线性屈曲分析吗?在step中选定line- perturbation下的各项,其Nlgeom都为Off,是不是意味着是进行不了啊? 【答】 line-perturbation应该是特征值屈曲分析,只能是线性的,要想进行非线性屈曲分析要引入初始缺陷 ABAQUS中非线性屈曲分析采用riks算法实现,可以考虑材料非线性、几何非线性已经初始缺陷的影响。其中,初始缺陷可以通过屈曲模态、振型以及一般节点位移来描述。 no.1:利用abaqus进行屈曲分析,一般有两步,首先是特征值屈曲分析,此分析为线性屈曲分析,是在小变形的情况进行的,也即上面提到过的模态,目的是得出临界荷载(一般取一阶模态的eigenvalue乘以所设定的load),且需要在inp 文件中,作如下修改 *node file,global=yes *End Step 此修改目的在于:在下一步后屈曲分析所需要的初始缺陷的节点输出为.fil文件。no.2:其次,就是所谓的后屈曲分析,此步一般定义为非线性,原因在于是在大变形情况进行的,一般采用位移控制加修正的弧长法,可以定义材料非线性,以及几何非线性,加上初始确定,所以也称为非线性屈曲分析。此步分析,为了得到极限值,需要得出荷载位移曲线的下降段,除了采用位移控制以及弧长法设定外,需在所得到的inp文件中,嵌入no.1中的.fil节点数据。修改如下: *IMPERFECTION(缺陷), FILE=results_file(此文件名为.fil), STEP=step(特征

ANSYS与ABAQUS稳定性分析比较

ANSYS与ABAQUS稳定性分析比较(转载-来自结构工程师崔家春的个人空间)其实,这些东西很简单,大多数朋友应该都比较了解。但是作为整个稳定性分析的一部分,觉得还是整理一下吧,也算是对后来者又抛了一块砖。 算例描述: 为了能体现出一般性,我故意找了一个比较大的结构。这是一个单层网壳结构,最大尺寸在90m左右,杆件长度在1.13m-3.63m之间,截面形式为箱型截面;构件布置见下图。荷载任意挑选一个标准组合(具体是哪个不记得,只是验证软件单元特征,没有关系)。 在ANSYS软件中分别采用BEAM44、BEAM188和BEAM189进行计算。分析结果见下文。 ANSYS BEAM44分析结果 E1E2E3E4E5 N1 6.10 6.367.117.438.03 N2 6.08 6.347.087.407.99 N3 6.08 6.347.087.407.98 N4 6.08 6.347.087.397.98 备注:表格中N1、N2分别代表每根构件采用1、2个单元;E1、E2代表第1、2阶屈曲荷载因子; ANSYS BEAM188分析结果 E1E2E3E4E5 N1 6.817.098.158.619.35 N2 6.25 6.527.347.698.34 N3 6.15 6.427.197.538.14 N4 6.12 6.387.147.478.07 N5 6.10 6.367.127.448.04 N6 6.09 6.357.107.438.02 N7 6.09 6.357.107.428.01 N8 6.08 6.347.097.418.00 N9 6.08 6.347.097.417.99 N10 6.08 6.347.097.407.99 ANSYS BEAM189分析结果 E1E2E3E4E5 N1 6.10 6.367.127.448.05 N2 6.07 6.337.087.407.98 N3 6.07 6.337.087.397.98 N4 6.07 6.337.087.397.97

abaqus压杆屈曲分析报告

压杆屈曲非线性分析 专业:结构工程 姓名:刘耀荣 学号:2110150113

压杆屈曲分析 1.问题描述 在钢结构中,受压杆件一般在其达到极限承载力前就会丧失稳定性,所以失稳是钢结构最为突出的问题。压杆整体失稳形式可以是弯曲、扭转和弯扭。钢构件在轴心压力作用下,弯曲失稳是常见的失稳形式。影响轴心受压构件整体稳定性的主要因素为纵向残余应力、初始弯曲、荷载初偏心及端部约束条件等。实际的轴心受压构件往往会存在上述的一种或多种缺陷,导致构件的稳定承载力降低。 本文利用abaqus 对一定截面不同长细比下的H 型钢构件进行屈曲分析,通过考虑材料非线性、几何非线性并引入初弯曲,得出构件发生弯曲失稳的极限荷载。通过比较不同长细比下的弯曲失稳的临界荷载得出构件荷载位移曲线,并与《规范》中的构件曲线相比较。钢构件的截面尺寸如图1-1所示。 构件的材料特性: E =2.0×1011 N m 2? ,μ=0.3 , f y =3.45×108N m 2? 图1-1 压杆截面尺寸(单位:m)

2.长细比计算 通过计算截面几何特性,截面绕y轴的回转半径为i y=0.0384m ,长细比取值及杆件长度见表1: 表1 3.模型分析 ABAQUS非线性屈曲分析的方法有riks法,general statics法(加阻尼),或者动力法。非线性屈曲分析采用riks算法实现,可以考虑材料非线性、几何非线性已及初始缺陷的影响。其中,初始缺陷可以通过屈曲模态、振型以及一般节点位移来描述。 利用abaqus进行屈曲分析,一般有两步,首先是特征值屈曲分析,此分析为线性屈曲分析,是在小变形的情况进行的,也即上面提到过的模态,目的是得出临界荷载(一般取一阶模态的eigenvalue乘以所设定的load)。其次,就是后屈曲分析,此步一般定义为非线性,原因在于是在大变形情况进行的,一般采用位移控制加修正的弧长法,可以定义材料非线性,以及几何非线性,加上初始缺陷,所以也称为非线性屈曲分析。此步分析,为了得到极限值,需要得出荷载位移曲线的下降段。缺陷较小的结构初始位移变形较小,在极值点突变,而初始缺陷较大的结构,载荷位移曲线较平滑。 4.建模计算过程 建模计算过程以长细比为50的构件为例,其余构件建模计算过程与之类似。 4.1 buckle分析 1 在buckle分析中创建part模块,创建的模型为三位可变形壳体单元,截面参数见图1-1,构件长度1.92。如图4-1示

ABAQUS非线性屈曲分析步骤

ABAQUS 6.7非线性屈曲分析步骤 riks法,或者general statics法(加阻尼),或者动力法 一共三种方法, 【问】在aba中能实现非线性屈曲分析吗?在step中选定line- perturbation 下的各项,其Nlgeom都为Off,是不是意味着是进行不了啊? 【答】 line-perturbation应该是特征值屈曲分析,只能是线性的,要想进行非线性屈曲分析要引入初始缺陷 ABAQUS中非线性屈曲分析采用riks算法实现,可以考虑材料非线性、几何非线性已经初始缺陷的影响。其中,初始缺陷可以通过屈曲模态、振型以及一般节点位移来描述。 no.1:利用abaqus进行屈曲分析,一般有两步,首先是特征值屈曲分析,此分析为线性屈曲分析,是在小变形的情况进行的,也即上面提到过的模态,目的是得出临界荷载(一般取一阶模态的eigenvalue乘以所设定的load),且需要在inp文件中,作如下修改 *node file,global=yes *End Step 此修改目的在于: 在下一步后屈曲分析所需要的初始缺陷的节点输出为.fil文件。 no.2:其次,就是所谓的后屈曲分析,此步一般定义为非线性,原因在于是在大变形情况进行的,一般采用位移控制加修正的弧长法,可以定义材料非线性,以及几何非线性,加上初始确定,所以也称为非线性屈曲分析。此步分析,为了得到极限值,需要得出荷载位移曲线的下降段,除了采用位移控制以

及弧长法设定外,需在所得到的inp文件中,嵌入no.1中的.fil节点数据。修改如下: *IMPERFECTION(缺陷), FILE=results_file(此文件名为.fil), STEP=step(特征值分析步名), 1(模态),2e-3(模态的比例因子,此值一般取杆件的1%,壳体厚度1%)此修改一般加在boundary之后step之前。 Re: 新手请教非线性屈曲中如何加初始扰动? 6.2.4 Unstable collapse and postbuckling analysis Rik法用于跳越失稳问题的研究,也可以用于分支屈曲的后屈曲研究。分支屈曲的后屈曲分析不能直接在分支屈曲后面研究,而是要给一个初始缺陷,使力学响应呈连续状态(非线性) 7.6.1 Introducing a geometric imperfection into a model 定义初始缺陷 Abaqus用三种手段定义初始缺陷,根据分支屈曲模型取一个线性组合,根据静力分析结果,直接指定。除非初始缺陷已经知道,一般采用第一种方法。 第一步,特征屈曲分析,Write the eigenmodes in the default global system to theresults file as nodal data (“Output to the data and results files,”Section 4.1.2).第二步,将这些特征屈曲模态添加到perfect几何体中,作为初始缺陷,where is the mode shape and is the associated scale factor. 一般来说系数w,在第一阶模态最大,而且w一般取结构几何参数的倍数,如壳的厚度的 0.1倍,等。 第三步,用rik法进行分析。

相关主题
相关文档 最新文档