几何非线性大作业荷载增量法和弧长法程序设计汇编完整.pdfVIP

几何非线性大作业荷载增量法

和弧长法程序设计汇编

(可以直接使用，可编辑优秀版资料，欢迎下载）

几何非线性大作业荷载增量法

和弧长法程序设计

系（所）：建筑工程系

学号：

姓名：焦联洪

培养层次：专业硕士

指导老师：吴明儿

2021年6月19

日

一、几何非线性大作业（Newton-Raphson法）

用荷载增量法（Newton-Raphson法）编写几何非线性程

序：

（1）用平面梁单元，可分析平面杆系

（2）算例：悬臂端作用弯矩。悬臂梁最终变形形成周长

为悬臂梁长度的圆。

1.1Newton-Raphson算法基本思想

图1.1Newton-Raphson算法基本思想

1.2悬臂梁参数

基本参数：L=2m,D=0.03m,A=7.069E-4m2,

I=3.976E-08m4,E=2.0E11N/m2

图1.2悬臂梁单元信息

将悬臂梁分成10个单元，如图1.2所示

2.1MATLAB输入信息

材料信息单元信

息

约束信息（0为约束，1为放松）荷

载信息（FX,FY,M）

节点信息

2.2求解过程

运用ABAQUS和MATLAB进行求解对比：

图1.3加载图

图1.4ABAQUS变形图

图1.5MATLAB变形曲线

ABAQUS和MATLAB变形对比，最终在理论荷载作用

下都弯成了一个圆，其直径为0.64716m，与理论值相对比

值为：（0.64716-0.642）/0.642=0.00804.非常接近。

2.3加载点荷载位移曲线

图1.5加载点Y方向的荷载位移曲线

加载点的最大竖向位移，相对比值（1.4525-1.45246）

/1.45246=2.75395E-05。完全相同，说明MATLAB的计算

结果很好。

二、几何非线性大作业（弧长法）

用弧长法编写几何非线性程序，分析荷载位移全过程曲线：

1)用平面梁单元，可分析平面杆系结构

2)算例

(1)受集中荷载的拱：考察拱的矢跨比、荷载位置对荷载

位移曲线的影响。

(2)其他有复杂平衡路径的结构

3)将结果与相关文献进行对比

1.1弧长法基本思想

图2.1弧长法基本思想

1.2拱基本参数

242

拱参数：L=100m,