adina提取内力

adina

中如何查看三维实体单元的截面轴力、弯矩与剪力？

1.

cut surface

方法

下面结合一个自由端作用集中荷载的三维悬臂梁实例，讲解如何计算某截面的轴

< br>力、弯矩与剪力。

1

）实例概况

截面尺寸为

0.10*0 .1

，

长度为

1

，

在自由端作用

2

个集

中力，

数值均为

1000

，

需要计算离自由端距离为

0.5

单位的 横截面上的轴力、

弯矩

与剪力，按照结构力学，该计算截面的轴 力为

0

，剪力为

2000

，弯矩为

1000.

下面

通 过

ADINA

程序验证上述数值的正确性。

2

）建模并求解

由于模型比较简单，

不详细讲解了 ，

需要说明的是，

坐标原点位移自由端截面最下

边。命令流如下，最终模型如下图：

*

DATABASE NEW SAVE=NO PROMPT=NO

FEPROGRAM ADINA

CONTROL FILEVERSION=V83

*

COORDINATES POINT SYSTEM=0

@CLEAR

1 0.00 0.00 0.00 0

2 0.00 0.100 0.00 0

3 0.00 0.100 0.100 0

4 0.00 0.00 0.100 0

@

*

SURFACE VERTEX NAME=1 P1=3 P2=4 P3=1 P4=2

*

VOLUME EXTRUDED NAME=1 SURFACE=1 DX=1.00,

DY=0.00 DZ=0.00 SYSTEM=0

PCOINCID=YES,

PTOLERAN=1.00E-05 NDIV=1 OPTION=VECTOR,

RATIO=1.00 PROGRESS=GEOMETRIC CBIAS=NO

*

FIXBOUNDARY SURFACES FIXITY=ALL

@CLEAR

6

'ALL'

@

*

LOAD FORCE NAME=1 MAGNITUD=1000.

FX=0.00,

FY=0.00 FZ=-1.00

*

APPLY-LOAD BODY=0

@CLEAR

1

'FORCE' 1

'POINT' 3 0 1 0.00 0 -1 0 0 0

'NO',

0.00 0.00 1 0

2

'FORCE' 1

'POINT' 4 0 1 0.00 0 -1 0 0 0

'NO',

0.00 0.00 1 0

@

*

MATERIAL ELASTIC NAME=1 E=2.00E+11

NU=0.300,

DENSITY=7800. ALPHA=0.00

MDESCRIP='NONE'

*

EGROUP THREEDSOLID NAME=1 DISPLACE=DEFAULT STRAINS=DEFAULT

MATERIAL=1,

RSINT=DEFAULT TINT=DEFAULT RESULTS=STRESSES DEGEN=NO

FORMULAT=0,

STRESSRE=GLOBAL INITIALS=NONE FRACTUR=NO CMASS=DEFAULT,