有限单元法作业非线性+选读.docVIP

几何非线性大作 和弧长法程序设计 系（所）：建筑工程系 学 号：1432055 姓 名：焦 联 洪 培养层次：专业硕士 指导老师：吴 明 儿 2015年6月19日 几何非线性大作业 Newton-Raphson法） 用荷载增量法（Newton-Raphson法）编写几何非线性程序： （1）用平面梁单元，可分析平面杆系 （2）算例：悬臂端作用弯矩。悬臂梁最终变形形成周长为悬臂梁长度的圆。 1.1 Newton-Raphson算法基本思想 图1.1 Newton-Raphson算法基本思想 1.2 悬臂梁参数 ：L=2m, D=0.03m, A=7.069E-4m2, I=3.976E-08m4 ,E=2.0E11N/m2 图 悬臂梁单元信息悬臂梁分成10个单元，如图1.2所示 2.1 MATLAB输入信息 材料信息 单元信息 约束信息（0为约束，1为放松） 荷载信息（FX,FY,M） 节点信息 2.2 求解过程 梁弯成圆形：理论弯矩M=EIY=24981.944N.m ,直径为0.642m 运用ABAQUS和MATLAB进行求解对比： 图1.3 加载图 图1.4 ABAQUS变形图 图1.5 MATLAB变形曲线 ABAQUS和MATLAB变形对比，最终在理论荷载作用下都弯成了一个圆，其直径为0.64716m，与理论值相对比值为（0.64716-0.642）/0.642=4.非常接近。 2.3 加载点荷载位移曲线 图1.5 加载点Y方向的荷载位移曲线 加载点的最大竖向位移分别为1.4525m和1.45246m，相对比值（1.4525-1.45246）/1.45246=2.75395E-05。完全相同，说明计算结果很好。 几何非线性大作业 1) 用平面梁单元，可分析平面杆系结构 2) 算例 (1)受集中荷载的拱：考察拱的矢跨比、荷载位置对荷载位移曲线的影响。 (2)其他有复杂平衡路径的结构 3) 将结果与相关文献进行对比 1.1 弧长法基本思想 图2.1 弧长法基本思想 1.2 拱基本参数 L=100m, A=0.32m2 ,I=1m4 ,E=1.0e7N/m2，F=-5000N，拱曲线 y=5×sin(3.1415926*x/L)分成25个单元，如图2所示 图 拱单元信息 2.1 MATLAB输入信息 材料信息 单元信息 约束信息（0为约束，1为放松） 荷载信息（FX,FY,M） 节点信息 2.2 运用ANSYS和MATLAB进行求解对比（两端铰接） ANSYS中模型 图2.3 ANSYS模型 图2.4 MATLAB和ANSYS变形图 2.3 加载点荷载位移曲线 图2.5 加载点荷载位移曲线 ANSYS求得的极限承载力3042.53,对应位移3.00142 MATLAB求得的极限承载力3043.8, 对应位移3.0768 相对误差分别为0.0417%,2.45%，模拟效果比较好。 2 不同矢跨比（1/20，3/40，1/10，3/20）的荷载位移曲线 计算 图2.6 荷载位移曲线 图2.7 荷载位移曲线 表1 各矢跨比下拱结构的极限荷载 参数 矢高 极值点 F(N) 位移（）N) 位移（） 5mm 3043.8 3.0768 1765.2 7.0816 7.5mm 7623.3 4.0335 -595.82 11.21 10mm 14974 5.4026 -6408.1 14.886 20mm 39791 9.4831 -63049 30.513 从表中初步得出：随着增加仍然呈现跳跃失稳的形式承载能力有大幅度的在的力出反向，相当于有一个拉力在阻止拱结构发生跳跃失稳矢跨比越大，拱越不容易发生跳跃失稳。拱的矢跨比一定范围后，拱将发生复杂的不同于的失稳形式。 2）MATLAB与ANSYS计算结果对比 图2.8 ANSYS和MATLAB对比荷载位移曲线 表2 各矢跨比下拱结构的极限荷载 参数 矢高 F(N)MAT 位移（）N)ANA 位移（））） 5mm 3043.8 3.0768 3042.53 3.00142 0.04 2.45 7.5mm 7623.3 4.0335 7624.91 3.96303 -0.02 1.75 10mm 14974 5.4026 14974.3 5.3157 0.00 1.61