西南交大工程计算实践作业
西南交通大学结构分析计算机程序与应用课程设计题目参考
请大家写成报告的形式。
要求：
1、主要的计算步骤截图（白底，黑图、黑字）
2、计算结果
3、计算结果分析 报告提交时间另行通知。
习题题目
1、图示交叉梁系的各梁与周边斜交
45°，梁端在周边为铰支，在结
点1、2、3、4、5 处各作用有垂直于结构平面的集中荷载20kN ，设各梁28/101.2m kN E ⨯=,45104m I -⨯=，241050m A -⨯=，。
绘出各梁最终弯矩图。
2、试用桁架与刚架分别建立下图的计算模型，P=10Kn, 比较计算结果。
3. 用结构分析有限元程序计算图2.40所示连续梁。
设72310kN/m E =⨯，截面形状为矩形，尺寸25cm ×50cm ，热膨胀系数51001-⨯=α.。
上侧温度升高25°，下侧温度升高10°。
要求分别计算出温度变化、均布荷载和集中力单独作用时的内力，并计算三者同时作用的组合内力。
4 计算错误！未找到引用源。
示下承式连续梁拱组合结构内力图。
拱轴线为抛物线，其方程为()x l x l
f
Z -=24，杆件截面0.6m 0.6m ⨯,弹性模量822.010kN/m E =⨯。
12m
30m
12m
6m
50kN
70kN
80kN
4m
5、空间刚架结构。
计算三层框架结构，层高3m ，楼板厚15cm ；柱尺寸为0.5⨯0.5m ，梁截面为0.2⨯0.4m ，计算楼面荷载为3kN/m 2时的内力。
723.010kN/m E =⨯.。
西南交通大学数学实验作业全面版
实验课题一基础编程第一大题：编程完成下列计算1. 当x = 3, x =2π 时，求1sin()xy x e =+ 的值。
%第一大题 %1x=[3,2*pi];y1=sin(x)+exp(x) %{ y1 =1517/75 31594/59 %}2. 用冒号法作等差数列x = 2,4,6,8,10求对应的函数22y x =+%2x=2:2:10;y2=x.^2+sqrt(2*x) %{ y2 =6 3841/204 7143/181 68 3761/36 %}3. 已知：22,35,a b c e π===计算：31sin cos();532tan()cot .3a y b c a y b ⎛⎫=+⨯ ⎪⎝⎭⎛⎫= ⎪⎝⎭%3a=2*pi,b=35,c=exp(2); y31=sin(a/5)+cos(b)*c y32=tan(b)*cot(a/3) %{ y31 =-4060/709y32 =-1019/3725 %}4. 将数据格式转换成有理格式后，清屏后重新输出a ,b ,c ,y 31,y 32（提示：参数选项或format rational ，清屏clc ） %4format rationalclc5.查看工作空间已有变量及信息。
(提示：打开变量信息窗口或whos)%5whos%{Name Size Bytes Class AttributesA 3x3 72 doubleA1 3x3 72 doubleA2 1x1 8 doubleA3 3x3 72 doubleS 21x2 336 doubleX 1x21 168 doubleY 1x21 168 doublea 1x1 8 doublea1 1x1 8 doublea11 1x1 8 doublea2 1x1 8 doublea21 1x1 8 doublea3 1x1 8 doublea31 1x1 8 doubleb 1x1 8 doublec 1x1 8 doubles 1x1 8 doublex 1x2 16 doubley1 1x2 16 doubley2 1x5 40 doubley31 1x1 8 doubley32 1x1 8 doubley71 1x1 8 doubley72 1x1 8 double%}6.a1=-6.28 a2=7.46 a3=5.37将a1,a2,a3分别向零取整后赋给a11,a21,a31。
西南交通大学数值分析上机实习
目录解题： (1)题目一： (1)1.1计算结果 (1)1.2结果分析 (1)题目二： (2)2.1计算结果 (2)2.2结果分析 (3)题目三： (4)3.1计算结果 (4)3.2结果分析 (5)总结 (5)附录 (6)Matlab程序： (6)题目一： (6)第一问Newton法： (6)第二问Newton法： (6)第一问Steffensen加速法： (7)第二问Steffensen加速法： (7)题目二 (8)1、Jacobi迭代法 (8)2、Causs-Seidel迭代法 (8)题目三： (9)题目一：分别用牛顿法，及基于牛顿算法下的Steffensen 加速法(1)求ln(x +sin x )=0的根。
初值x0分别取0.1, 1,1.5, 2, 4进行计算。
(2)求sin x =0的根。
初值x0分别取1，1.4，1.6, 1.8，3进行计算。
分析其中遇到的现象与问题。
1.1计算结果求ln(x +sin x )=0的根，可变行为求解x-sinx-1=0的根。
1.2结果分析从结果对比我们可发现牛顿—Steffensen 加速法比牛顿法要收敛的快，牛顿法对于初值的选取特别重要，比如第（1）问中的初值为4的情况，100次内没有迭代出来收敛解，而用Steffensen 加速法，7次迭代可得；在第（2）问中的初值为1.6的情况，收敛解得31.4159，分析其原因应该是x x f cos )('=，x0=1.62π≈,0)('≈x f ；迭代式在迭代过程中会出现分母趋近于0，程序自动停止迭代的情况，此时得到的x 往往非常大，而在第一问中我们如果转化为用x+sinx=1,则可以收敛到结果。
用雅格比法与高斯－赛德尔迭代法解下列方程组Ax=b，研究其收敛性，上机验证理论分析是否正确，比较它们的收敛速度，观察右端项对迭代收敛有无影响。
(1)A行分别为A1=[6,2,-1],A2=[1,4,-2],A3=[-3,1,4]；b1=[-3,2,4]T,b2=[100,-200,345]T,(2) A行分别为A1=[1,0,8,0.8],A2=[0.8,1,0.8],A3=[0.8,0.8,1]；b1=[3,2,1]T,b2=[5,0,-10]T,(3)A行分别为A1=[1,3],A2=[-7,1]；b=[4,6]T2.1计算结果初值均为0矩阵带入(1)A行分别为A1=[6,2,-1],A2=[1,4,-2],A3=[-3,1,4]；b1=[-3,2,4]T,b2=[100,-200,345]T2) A行分别为A1=[1,0,8,0.8],A2=[0.8,1,0.8],A3=[0.8,0.8,1]；b1=[3,2,1]T,b2=[5,0,-10]TT2.2结果分析ρ小于1,故方程组雅可比迭代收第一小题的经计算谱半径为5427B(=).0敛。
西南交通大学限修课数学实验题目及答案三
西南交通大学限修课数学实验题目及答案三实验课题三向量与曲线绘图第一个主要问题：向量的创建和操作用元素输入法创建向量x11=(2c58c171-83259)x11=[2-58-171-83259]用冒号输入法创建向量x12=(246810121416182022)x12=2:2:22用等分取值法创建向量x13，其初值为0，终值为2π，共20个元素.X13=linspace（0,2*PI，20）使用随机输入法x14x14=rand（1,8）创建8维行向量用随机输入法创建6维整数列向量x15x15=fix(rand(6,1)*100)取向量X11的绝对值大于3的元素，形成向量x16 x16=X11（abs（X11）>3）求空间两点间距离m1(5,?4,?9)、m2(8,?6,?3)d=norm([863]-[549])向量的线性运算：x18=4 x X11+7x12 x18=4+X11+7*x12做向量的数量积x19=x11x12.x19=点（x11，x12）分别取x11与x12的前三个元素做向量的叉积赋给x10.x10=cross(x11([1:3]),x12([1:3]))第二大题：曲线绘图：构造坐标向量并绘制单词“天”的图形（首先给出组成单词的数据点的坐标）axis([0,6,0,6])x=[1155115533]；y=[1551133551]；直线（x，y）绘制向量y=[4553235678]的图形。
y=[4553235678];plot(y)数据数组x23=（0.10.110.12…10），函数Y23=30/x23，绘制函数曲线。
x23=0.1:0.01:10；y23=30x23；绘图（x23，y23）数据数组x24为区间[-5,5]上等分的30个点列，绘出函数y24=5x24cos(x24)的曲线图。
x24=linspace(-5,5,30);y24=5*x24.*cos(x24);绘图（x24，y24）数据数组x25是从[-2?，2?]中，取50个点,在同一块图形窗口绘出蓝色、数据点o、实线线型的y25=sin(x25)和红色、数据点*、虚线线型的z25=cos(x25)。
西南交大《工程数学Ⅰ》1-4次离线作业
工程数学Ⅰ第1次离线作业三、主观题(共15道小题)29.求5元排列52143的逆序数。
解答：在排列52143中，排在5之后，并小于5的数有4个；排在2之后，并小于2的数有1个；排在1之后，并小于1的数有0个；排在4之后，并小于4的数有1个。
所以30.计算行列式解答：容易发现D的特点是：每列（行）元素之和都等于6，那么，把二、三、四行同时加到第一行，并提出第一行的公因子6，便得到由于上式右端行列式第一行的元素都等于1，那么让二、三、四行都减去第一行得31.求行列式中元素a和b的代数余子式。
解答：行列式展开方法==32.计算行列式解答：容易发现D的特点是：每列元素之和都等于6，那么，把二、三、四行同时加到第一行，并提出第一行的公因子6，便得到由于上式右端行列式第一行的元素都等于1，那么让二、三、四列都减去第一列，第一行就出现了三个零元素，即33.设，求解答：34.，求解答：35.求矩阵X使之满足解答：36.解矩阵方程，其中解答：首先计算出，所以A是可逆矩阵。
对矩阵（A，B）作初等行变换所以所以秩（A）= 4。
37.解答：38.求向量组解答：设39.求解非齐次线性方程组解答：对增广矩阵施行初等行变换化成简单阶梯形矩阵40.设解答：若41.设，求A的特征值和特征向量。
解答：42.求一个正交矩阵P，将对称矩阵化为对角矩阵。
解答：43.已知二次型，问：满足什么条件时，二次型 f 是正定的；满足什么条件时，二次型 f 是负定的。
解答：二次型 f 的矩阵为计算 A 的各阶主子式得工程数学Ⅰ第2次离线作业三、主观题(共14道小题)30.判断（1）；（2）是否是五阶行列式 D5 中的项。
解答：（1）是；（2）不是；31.设求的根。
解答：行列式特点是：每行元素之和都等于 a+b+c+x，那么，把二、三、四列同时加到第一列，并提出第一列的公因子a+b+c+x，便得到二、三、四列-a依次减去第一列的-a、-b、-c倍得32.计算四阶行列式解答：D的第一行元素的代数余子式依次为由行列式的定义计算得33.用克莱姆法则解方程组解答：34.解答：35.解答：36.用初等行变换把矩阵化为阶梯形矩阵和简单阶梯形矩阵。
西南交通大学新学期《工程数学I》在线作业一
西南交《工程数学I》在线作业一
如果矩阵A满足A^2=A，则( )
A:A=0
B:A=E
C:A=0或A=E
D:A不可逆或A-E不可逆
参考选项：D
A、B均为n阶方阵，则必有
A:det(A)det(B)=det(B)det(A)
B:det(A+B)=det(A)+det(B)
C:(A+B)的转置=A+B
D:(AB)的转置＝A的转置乘以B的转置
参考选项：A
设A为n阶方阵，r(A)<n，下列关于齐次线性方程组Ax=0的叙述正确的是
（）
A:Ax=0只有零解
B:Ax=0的基础解系含r(A)个解向量
C:Ax=0的基础解系含n-r(A)个解向量
D:Ax=0没有解
参考选项：C
n阶行列式的展开式中共有（）项
A:n
B:n^2
C:n!
D:n(n+1)/2
参考选项：C
设3阶实对称矩阵A的特征值分别为2，0，-3，则（）
A:|A|≠0
B:A负定
C:A正定
D:|A|=0
参考选项：D
设A，B均为n阶方阵,则等式(A+B)(A-B) = A2-B2成立的充分必要条件是( )．A:A=E
B:B=O
C:A=B
D:AB=BA
1。
西南交大土木实习报告
一、实习目的与要求本次实习是西南交通大学土木工程专业教学计划中的一项重要实践教学环节，旨在将所学的理论知识与实际工程实践相结合，提高学生的实际操作能力和工程管理能力。
通过实习，我们期望达到以下目的和要求：1. 了解土木工程的基本施工流程，熟悉施工现场的各项工作；2. 掌握土木工程相关材料的性能和施工技术；3. 培养团队协作精神和工程实践能力；4. 提高对土木工程行业现状和发展趋势的认识。
二、实习时间与地点本次实习时间为2021年6月1日至2021年6月30日，实习地点为成都市某大型住宅小区建设项目。
三、实习内容1. 施工现场考察在实习期间，我们深入施工现场，对住宅小区的建设项目进行了全面考察。
主要内容包括：（1）了解施工现场的组织架构和工程管理流程；（2）熟悉施工现场的安全措施和环境保护措施；（3）观察施工现场的施工进度和施工质量；（4）了解施工现场的施工设备和材料。
2. 施工技术学习在实习期间，我们学习了以下施工技术：（1）地基基础施工技术；（2）主体结构施工技术；（3）装饰装修施工技术；（4）园林景观施工技术。
3. 工程管理实践在实习期间，我们参与了以下工程管理工作：（1）施工进度管理：了解施工进度计划，掌握施工进度控制方法；（2）施工质量管理：了解施工质量控制标准，掌握施工质量检验方法；（3）施工安全管理：了解施工现场安全规定，掌握施工现场安全措施；（4）施工成本管理：了解施工成本构成，掌握施工成本控制方法。
四、实习心得体会1. 实践出真知。
通过本次实习，我深刻认识到理论知识与实际工程实践相结合的重要性。
在实际工程中，我们需要灵活运用所学知识，解决实际问题。
2. 团队协作精神。
在实习过程中，我们与团队成员密切配合，共同完成了各项工作任务。
这使我明白了团队协作在工程实践中的重要性。
3. 安全意识。
在施工现场，安全至关重要。
通过本次实习，我更加重视安全意识，遵守施工现场的各项规定。
4. 职业素养。
西南交通大学计算机程序设计基础C++实验九
实验 9 实验报告教学班级：_制药一班______ 学生学号：xxxxxx 学生姓名：刘祥龙实验日期：__6.1_______ 实验机房：____x7407____ 助教老师：__________一、实验目的掌握对一组数据进行插入、删除、修改、以及查找的常用算法。
二、实验任务1.建立一个一维数组a，完成以下任务，输入输出格式参见下图。
（1）利用rand()函数产生10个10~100之间的整数，存入数组a中并输出；（2）使用冒泡排序法，将数组a中的数据按降序排列，并输出；（3）使用折半查找法在数组a中查找是否有键盘输入的指定数据，并输出查找结论（没找到或在第X个位置找到）。
2.建立一个一维数组a，完成以下任务，输入输出格式参见下图。
（1）自动生成斐波那契数列前10项的值，存入数组a，并输出。
（斐波那契数列前两项的值为1，从第三项开始，每一项的值是前两项值的和）；（2）将键盘输入的一个数据插入数组a，使得数组a中的元素仍然保持有序；（3）删除数组a中的重复数据，并输出删除后的结果；（4）求出数组a所有元素的平均值，并输出；（5）删除数组a中所有低于平均值的元素，并输出删除后的结果。
三、实验结果（源程序+ 注释+ 运行界面截图）（要求：（1）提供源程序；（2）源程序包含必要的注释，包括：在程序头部给出本程序功能的注释、所定义变量的用途注释、程序块的功能注释、子函数的功能注释等；（3）运行结果的窗口截图，程序运行结果中输出的第一行内容必须为：程序作者的学号和姓名。
）1.#include<iostream>;using namespace std;int main(){cout<<"xxxxxxxxx"<<endl;int a[10],i,j,t,top=0,mid,bot=9,x;cout<<"10个10到100间的随机整数为："<<endl;for(i=0;i<10;i++){a[i]=10+rand()%(100-10+1);cout<<a[i]<<" ";}cout<<endl;//使用冒泡法，将数组a按降序排列，并输出for(i=0;i<9;i++){for(j=0;j<9;j++)if(a[j]<a[j+1]){t=a[j];a[j]=a[j+1];a[j+1]=t;}}cout<<"按从大到小排列的结果为："<<endl;for(i=0;i<10;i++)cout<<a[i]<<" ";cout<<endl;//使用折半查找完成指定数据的查找，并输出是否查找到的信息cout<<"请输入要查找的数：";cin>>x;while(top<=bot){mid=(top+bot)/2;if(a[mid]==x){break;}elseif(a[mid]<x)bot=mid-1;elseif(a[mid]>x)top=mid+1;}if(top<=bot)cout<<x<<"在第"<<mid+1<<"个位置。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
工程计算实践第一次作业
班级：专隧1班
姓名：XXXXX
学号：xxxxxxxxxx
2017年4月
1.工程背景
选取一相似实际工程，该隧道为单洞双线隧道。
区域为中低山区，地形起伏较大，山坡自然坡度约30~60度，相对高差300~860m左右，主要岩层地层岩性较单一，地表零星分布有第四系坡残积层，主要为花岗岩、花岗闪长岩、全~弱风化，全风化层厚度约为20~30m，强风化层厚3~5m，下伏弱风化层。
埋深约为30m。
掌子面围岩主要为强风化花岗岩、花岗闪长岩。
掌子面揭露呈黄褐色，岩体破碎，节理裂隙较发育，裂缝密集且部分张开；毛开挖面能够自稳，偶有掉块现象。
初期支护喷射35cm厚，C25混凝土，钢架间距60cm。
每三棍在拱部1400范围内施做一环4m长小导管，间距30cm。
二衬厚度60cm, C30混凝土。
2.计算模型建立
根据隧道支护结构尺寸信息，采用ANSYS建模，仅考虑二次衬砌受力，二次衬砌承担50%的荷载。
2.1计算参数
二衬C30混凝土以及V级围岩的参数见表1-1。
根据《铁路隧道设计规范》，可以计算出深埋隧道围岩的垂直和水平均布荷载，如表1-2所示。
对于竖向和水平的均布荷载其等效节点力分别近似取节点两相邻单元水平或垂直投影长度的一半衬砌计算宽度这一面积范围内的分布荷载总和。
结构自重直接通过ANSYS 添加密度施加。
2. 2模拟步骤
2.2.1建立模型
打开ANSYS，进入前处理器，输入如下命令可初步建立模型。
图1 模型示意图
2.2.2为模型赋予属性
模型材料选择为梁单元(beam3 )，截面宽度lm，高度0.4m，弹性模量32Gpa泊松比0.2，密度2500。
同时选择弹簧单元（combin14）设置其相关物理和几何参数，如图2所示
2.2.3网格划分
分别对用梁单元（BEAM3）弹簧单元（combin14）对其进行网格划分。
用杆单元来模拟弹簧单元。
2.2.4施加约束及荷载
施加节点约束，分别对弹簧单元施加X,Y方向的约束，约束值设置为0，同时在
隧道底部约束住隧道水平方向的位移，约束值设置为0。
:如图3所示
图3 约束设置图2.2.5施加节点荷载
节点荷载计算方法
节点荷载计算示意图如下
图4 单元ij受力图
衬砌结构承受主动垂直线性荷载和水平线性荷载，计算时将其转换为节点荷载。
假设单元ij受力如图1-10所示荷载作用，则节点i和节点j的等效节点荷载可分别由如下的荷载阵列求出。
命令流如下
*SET,_Q1, 99.404e3
*SET,_Q2, 99.404e3
*SET,_E1, 24.851e3
*SET,_E2, 24.851e3
*SET,_E3, 24.851e3
*SET,_E4, 24.851e3
*GET,_E_MAX,ELEM,,NUM,MAX
*GET,_E_MIN,ELEM,,NUM,MIN
*GET,_E_COUNT,ELEM,,COUNT
*GET,_N_MAX,NODE,,NUM,MAX
*GET,_N_MIN,NODE,,NUM,MIN
*GET,_N_COUNT,NODE,,COUNT
*GET,_LOC_X_MAX,NODE,,MXLOC,X
*GET,_LOC_X_MIN,NODE,,MNLOC,X
*GET,_LOC_Y_MAX,NODE,,MXLOC,Y
*GET,_LOC_Y_MIN,NODE,,MNLOC,Y
*SET,_WIDTH,_LOC_X_MAX-_LOC_X_MIN
*SET,_HEIGHT,_LOC_Y_MAX-_LOC_Y_MIN
*SET,_CENT_X,(_LOC_X_MAX+_LOC_X_MIN)/2
*SET,_CENT_Y,(_LOC_Y_MAX+_LOC_Y_MIN)/2
*SET,_I_ND,_N_MIN
*SET,_I_EL,_E_MIN
*DIM,_NF_INFO,ARRAY,_N_COUNT,4
*DIM,_ND_INDEX,ARRAY,_N_MAX
*DO,I,1,_N_COUNT,1
*IF,NX(_I_ND),EQ,_LOC_X_MAX,THEN
_CENT_Y=NY(_I_ND)
*ENDIF
*IF,NY(_I_ND),EQ,_LOC_Y_MAX,THEN
_CENT_X=NX(_I_ND)
*ENDIF
_NF_INFO(I,1)=_I_ND
_ND_INDEX(_I_ND)=I
_I_ND=NDNEXT(_I_ND)
*ENDDO
*DO,I,1,_E_COUNT,1
_NI=NELEM(_I_EL,1)
_NJ=NELEM(_I_EL,2)
_LOC_XI=NX(_NI)-_LOC_X_MIN
_LOC_YI=NY(_NI)-_LOC_Y_MIN
_LOC_XJ=NX(_NJ)-_LOC_X_MIN
_LOC_YJ=NY(_NJ)-_LOC_Y_MIN
*GET,_E_CENT_X,ELEM,_I_EL,CENT,X
*GET,_E_CENT_Y,ELEM,_I_EL,CENT,Y
*IF,_E_CENT_X,GT,_CENT_X,THEN
*IF,_E_CENT_Y,GT,_CENT_Y,THEN
_E_I=-_E3-(_E4-_E3)*(1-_LOC_YI/_HEIGHT)
_Q_I=-_Q1-(_Q2-_Q1)*_LOC_XI/_WIDTH
_E_J=-_E3-(_E4-_E3)*(1-_LOC_YJ/_HEIGHT)
_Q_J=-_Q1-(_Q2-_Q1)*_LOC_XJ/_WIDTH *ELSE
_E_I=-_E3-(_E4-_E3)*(1-_LOC_YI/_HEIGHT)
_Q_I=0
_E_J=-_E3-(_E4-_E3)*(1-_LOC_YJ/_HEIGHT)
_Q_J=0
*ENDIF
*ELSE
*IF,_E_CENT_Y,GT,_CENT_Y,THEN
_E_I= _E1+(_E2-_E1)*(1-_LOC_YI/_HEIGHT)
_Q_I=-_Q1-(_Q2-_Q1)*_LOC_XI/_WIDTH
_E_J= _E1+(_E2-_E1)*(1-_LOC_YJ/_HEIGHT)
_Q_J=-_Q1-(_Q2-_Q1)*_LOC_XJ/_WIDTH
*ELSE
_E_I=_E1+(_E2-_E1)*(1-_LOC_YI/_HEIGHT)
_Q_I=0
_E_J=_E1+(_E2-_E1)*(1-_LOC_YJ/_HEIGHT)
_Q_J=0
*ENDIF
*ENDIF
_FX_I=0.05*(3*_E_J+7*_E_I)*ABS(_LOC_YJ-_LOC_YI)
_FY_I=0.05*(7*_Q_I+3*_Q_J)*ABS(_LOC_XJ-_LOC_XI)
_FX_J=0.05*(7*_E_J+3*_E_I)*ABS(_LOC_YJ-_LOC_YI)
_FY_J=0.05*(3*_Q_I+7*_Q_J)*ABS(_LOC_XJ-_LOC_XI)
_NF_INFO(_ND_INDEX(_NI),2)=_NF_INFO(_ND_INDEX(_NI),2)+_FX_I _NF_INFO(_ND_INDEX(_NI),3)=_NF_INFO(_ND_INDEX(_NI),3)+_FY_I _NF_INFO(_ND_INDEX(_NJ),2)=_NF_INFO(_ND_INDEX(_NJ),2)+_FX_J _NF_INFO(_ND_INDEX(_NJ),3)=_NF_INFO(_ND_INDEX(_NJ),3)+_FY_J ! MZ
_I_EL=ELNEXT(_I_EL)
*ENDDO
*DO,I,1,_N_COUNT,1
F,_NF_INFO(I,1),FX,_NF_INFO(I,2)
F,_NF_INFO(I,1),FY,_NF_INFO(I,3)
! MZ
*ENDDO
! FREE THE MEMORY
*DEL,_ND_INDEX,,NOPR
*DEL,_NF_INFO,,NOPR
3.求解
施加重力9.8，同时指定分析类型为静态分析，然后求解，求解完成，查看，节
点轴力，弯矩，剪力图如下图，5，6，7所示
轴力图(单位:N)
弯矩图(单位:N/m)
剪力图(单位:N}m)
其中最大压力在303号单元处，其值为0.49828E+07，最大拉力在296处，其值为0.12476E+07，最大负弯矩在296号单元处，其值为0.43405E+06，最大正弯矩在283号单元处，其值为0.79314E+06，最大剪力发生在296号单元处，其值为0.60868E+06。
计算各断面最不利位置的偏心距，去年权系数，按破损截断法验算衬砌截面强度，有：K=N极限/N，式中，N极限—截面的极限承载力（轴力）N—截面的实际轴力K—强度安全系数对于截面极限承载力N极限的计算，依轴力偏心距NMe的大小有两种情况：
对混凝土矩形截面构建，当e≤0.20h时，承载能力由抗压强度控制：
N极限=φαRαbh
当e0>0.20h时，抗拉强度也控制承载能力。
内力最大四个特征点验算：。