非线性有限元作业任务老骆整编.pdf

非线性有限元作业任务老骆整编

**

1.轴对称问题的弹塑性分析

流程图：

开始

读入单元号，B矩阵位数，单刚位数，单元

节点号，单刚等各项参数

EN1存储单元节点号，局部坐标系转

换为全局坐标

Y

JD0

N

打印错误E=σ̅⁄ε̅调用子函数DEMATR求[D]

调用子函数BMATR求[B]

Y

R=1

N

[C]=[De][B]

[C]=[Dep][B]

[S]=∑[C]T[B]=∑[B]T[D][B]

切线刚度阵[EK]=[S][Q1]=∑[B]T[D][B]·JD·RN·H(I1)H(J1)

Y

返回各值

结束

非线性有限元作业任务老骆整编

**

解析解。厚壁筒受内压，采用Mises屈服准则

2

σ−σ=σ(1)

θγ√3s

经计算知，当t=12.5（12.91）时，材料处于弹塑性交界面。

弹性区为:t≤r≤b