兰州交大MIDAS作业之拱桥空间分析计算报告

兰州交大MIDAS作业之拱桥空间分析计算报告
兰州交大MIDAS作业之拱桥空间分析计算报告

一、建模思路

1、定义材料和截面

按给定材料和截面数据定义相关材料和截面,其中虚设梁为指定车辆移动荷载而设。对虚设梁的材料和截面可输入任意的足够小的值以使其不致于对分析结果产生影响。

2.建模

先利用拱桥建模助手,输入矢高和跨径,生成拱肋。然后利用扩展线单元,选择“投影”并选中吊杆节点,选择相应材料、截面建立吊杆。全选复制生成另一拱肋及吊杆单元。扩展线单元生成横系梁及纵梁,桥面系斜撑的建立可先复制后镜像生成。最后建立拱上斜支撑和横向支撑及虚拟梁。(其中吊杆单元为桁架单元,其余均为梁单元)。建立完成后的模型如图1:

图1.全桥模型

3.边界条件

按照题目所给B1-B4约束条件对模型施加边界约束,支撑构件的两端对单元坐标系y、z轴为铰接条件,与主梁相连接的横系梁两端对单元坐标系y、z轴为铰接条件。吊杆与拱肋及主梁选择刚性连接。边界施加后如图2所示。

图2.边界条件

4. 车辆移动荷载和静力荷载

按要求输入人行道荷载和固定荷载,加载方式为梁单元均布加载。以虚拟梁为参照,进行车辆移动荷载加载。选择的公路工程技术标准(JTG B01-2003)的荷载,程序默认为公路-I 级荷载。移动荷载分析控制数据里需输入基频,因此需将自重转换为质量,并将恒载转为质量,求解出基频。经计算,基频为0.000442。程序运行结果如图3所示。(附:可能是程序破解狗有问题,程序得出的基频不正常,太偏小。)

图3.自振周期及频率

为求得吊杆初张力,先假设为18根吊杆初张力均为1KN。运行结果后再进行调试。初张力数值表如图4所示。

图4.初张力数据表

5.运行分析,查看结果

运行后,进行荷载组合。将自重、固定荷载及人行道荷载组合做为恒载状况。分别查看恒载和活载内力图及变形图。图形结果如图所示。

图5.恒载作用下内力图

图6.活载作用下内力图

图7.恒载作用下变形图

图8.活载作用下变形图

二、吊杆力计算及吊杆优化设计

1. 刚性支承连续梁法

所谓刚性支承连续法就是求一组恒载张拉力值, 使系梁与吊杆连接处节点在恒载和张拉力的作用下,在成桥状态下的位移为零,并且同时认为系梁内的弯矩为刚性支承连续梁弯矩。可见此处在确定吊杆的内力时, 主要以位移为控制目标(若为一次落架时同于零位移法) ,以保证最终成桥时的线型。

根据一次落架方式可以算出在恒载g 作用下而柔性吊杆初始张拉力为零时, 吊杆与主梁连接处节点的垂直位移,然后依次算出吊杆受单元力时对这些点的位移影响量ki, k 为节点编号, i 为吊杆编号, 于是在恒载和吊杆力的共同作用下, 以使各控制节点变位等于零为目标, 可以写出线性方程组:[ A] {X} + {△} = 0 ( 1)其中:{X} = { x1x2??xn }xi :表示第i根吊杆的张拉力矢量;{△} = {△1g△2g??△ng}△ig表示位移矢量。

2. 吊杆初张力求解

首先定义荷载组合,将固定荷载和各个吊杆初张力系数均设为1,如图9所示。

图9.荷载组合表

然后利用迈达斯的索力调整功能,以吊杆与主梁相交的各个节点的竖向位移Dz=0做为

控制目标。点击搜索按钮进行相关参数设置,系统自动进行求解。索力调整界面如图10所示。

图10.索力调整表

但是由于计算数据量庞大,电脑经常处于死机状态。故采取未知荷载系数法求解影响矩阵。相关数据见附图11。

图11.未知荷载系数法

生成的结果如图12所示,将其导出为excel文件。

图12.未知荷载系数法求解结果将结果整理制作表格如下:

3. 优化设计

(1)调索前后内力比较

将原模型复制供优化设计使用,将求解出来的初张力代替原来的假设张力。运行分析,查看位移变化是否满足要求。结果如下图:

图13.位移变化

由位移变化等值线图可知:主梁平面大部分节点位移均为0,中间两个单元位移稍大,为1.26e-4,近似为0。说明求解方法可行。

此时,查看恒载作用下内力图,如下图所示:

图14.恒载作用下内力图(确定吊杆初张力后)

图15.活载作用下内力图(确定吊杆初张力后)

图16.恒载和吊杆张力作用下内力图(确定吊杆初张力后)由以上图形数据可得知:在张拉力调整前,最大弯矩为1.9e+003KN*m;调整后仅在恒载作用下,弯矩降为1.22 e+003KN*m;吊杆张力与恒载共同作用下,弯矩进一步降为

4.74 e+002KN*m。这说明吊杆参与了结构受力,并减小了主梁相应的弯矩。

(2)吊杆截面尺寸优化

方法:做出吊杆内力图及内力表格,将表格内数据绘制图形查看吊杆内力趋势。按照标准抗拉强度1860MPa考虑,即将各个吊杆内力除以钢绞线强度值,即可得出各个吊杆优化后的面积。相关图表如下:

图17.吊杆内力图(确定吊杆初张力后)

图18.各个吊杆内力图(确定吊杆初张力后)

综合以上信息,可知除了40#单元和48#单元直径由原来的100mm变为22mm外,其余均变为20mm。

统计西安交大期末考试试题(含答案)

西安交大统计学考试试卷 一、单项选择题(每小题2 分,共20 分) 1.在企业统计中,下列统计标志中属于数量标志的是(C) A、文化程度 B、职业 C、月工资 D、行业 2.下列属于相对数的综合指标有(B ) A、国民收入 B、人均国民收入 C、国内生产净值 D、设备台数 3.有三个企业的年利润额分别是5000 万元、8000 万元和3900 万元,则这句话中有(B)个变量? A、0 个 B、两个 C、1 个 D、3 个 4.下列变量中属于连续型变量的是(A ) A、身高 B、产品件数 C、企业人数 D、产品品种 5.下列各项中,属于时点指标的有(A ) A、库存额 B、总收入 C、平均收入 D、人均收入 6.典型调查是(B )确定调查单位的 A、随机 B、主观 C、随意 D 盲目 7.总体标准差未知时总体均值的假设检验要用到(A ): A、Z 统计量 B、t 统计量 C、统计量 D、X 统计量 8.把样本总体中全部单位数的集合称为(A ) A、样本 B、小总体 C、样本容量 D、总体容量 9.概率的取值范围是p(D ) A、大于1 B、大于-1 C、小于1 D、在0 与1 之间 10.算术平均数的离差之和等于(A ) A、零 B、1 C、-1 D、2 二、多项选择题(每小题2 分,共10 分。每题全部答对才给分,否则不计分) 1.数据的计量尺度包括(ABCD ): A、定类尺度 B、定序尺度 C、定距尺度 D、定比尺度 E、测量尺度 2.下列属于连续型变量的有(BE ): A、工人人数 B、商品销售额 C、商品库存额 D、商品库存量 E、总产值 3.测量变量离中趋势的指标有(ABE ) A、极差 B、平均差 C、几何平均数 D、众数 E、标准差 4.在工业企业的设备调查中(BDE ) A、工业企业是调查对象 B、工业企业的所有设备是调查对象 C、每台设备是 填报单位D、每台设备是调查单位E、每个工业企业是填报单位 5.下列平均数中,容易受数列中极端值影响的平均数有(ABC ) A、算术平均数 B、调和平均数 C、几何平均数 D、中位数 E、众数 三、判断题(在正确答案后写“对”,在错误答案后写“错”。每小题1 分,共10 分) 1、“性别”是品质标志。(对)

迈达斯-截面特性值计算器

<图 1-(1)> 生成Plane 截面的过程 建立截面的轮廓 生成Plane 截面 利用网格进行计算

※注意事项 MIDAS/Civil和Gen数据库中提供的规则截面的抗扭刚度计算方法参见附录一。 对于MIDAS/Civil和Gen数据库中提供的规则截面,利用 MIDAS/Civil、Gen的截面特性计算功能计算截面特性值比SPC更好一些。 MIDAS/Civil和Gen数据库中提供的PSC截面,当用户输入的截面属于薄壁型截面时,应使用本截面特性值中的Line方式重新计算抗扭刚度,然后在截面特性值增减系数中对抗扭刚度进行调整。 对于Plane形式的截面,程序是通过有限元法来近似计算抗扭刚度的。在抗扭问题里使用的近似求解法有Ritz法(或者Galerkin法)、Trefftz法,所有的近似求解都与实际结果多少有点误差,其特征如下: J Ritz≤J Exact≤J Trefftz 像SPC一样利用有限元法近似地计算抗扭刚度时,通常使用Ritz法, 故其计算结果有可能比实际的抗扭刚度小。用户可通过加大网格划分密度方法来提高结果的精确度。 对于Line形式的截面, 如薄壁截面,线的厚度很薄时几乎可以准确地计算其抗扭刚度。但是如果是闭合截面(无开口截面),这种计算方式会导致其抗扭刚度的计算结果随着线厚度的增加而变小,所以对于不是薄壁截面的闭合截面应尽量避免使用Line的方式计算截面特性。 在SPC中对薄壁闭合截面,对闭合部分一定要使用model>closed loop>Register指定闭合。 SPC可以在一个窗口里任意的建立很多个截面,并分别进行分析,且可根据名称、位置、截面特性值等可以很方便地对截面进行搜索及排列。 <图2> 将DXF文件中的截面形状导入后,生成截面并进行排列

2013西安交大电路在线作业及答案

一、单选题(共20 道试题,共40 分。) V 1. 在RL串联的正弦电路中,如果电阻上电压为30伏,电感上电压为40伏,那么总电压为( B). A. 70伏 B. 50伏 C. 10伏 D. 50安 满分:2 分 2. 在电感电路中,当通入交流电时,电压与电流的关系是(B)。 A. 同相 B. 电压超前电流 C. 电流超前电压 D. 反相 满分:2 分 3. 网孔电流法中,以网孔电流作为电路变量时,需要按照(C)列出电路方程。 A. KCL B. 元件的VCR C. KVL D. 以上三个都需要 满分:2 分 4. 在动态电路中,不能发生跃变的是(B)。 A. 电压 B. 电容电压 C. 电容电压 D. 电流 满分:2 分 5. 一个元件的电压电流为关联参考方向,若电流I=-3A,元件产生的功率是12w,则电压=(C)V。 A. -4 B. -1/4 C. 4 D. 36 满分:2 分 6. 额定值为1W 10Ω的电阻器使用时通过电流的限额是(C)。 A. 0.1A B. 0.2A C. 0.316A D. 0.5A 满分:2 分 7. 现有额定值为220V、25W的电灯一只,准备串联一个小灯泡放在另一房间作为它的信号灯,应选择___A__规格的小灯泡? A. 6V,0.15A B. 6V,0.1A C. 1.5V,0.5A

D. 1.5V,0.3A 满分:2 分 8. 若RLC串联电路的谐振角频率为w0,则在角频率w>w0时电路呈现(B) A. 纯电阻性 B. 电感性 C. 电容性 D. 不能确定 满分:2 分 9.下列几项中,哪一对不是对偶元素的?(C) A. 电压-电流 B. 电阻-电导 C. 电容-电阻 D. 电容-电感 满分:2 分 10. 理想电感元件在某时刻其电压值为0,则该时刻电感上的储能(A) A. 不一定为0 B. 一定为0 C. 一定不为0 D. 无法确定 满分:2 分 11. 正弦稳态电路中, 一个无源单口网络的导纳Y=3 –j 4 Ω,电路的性质是(B)。 A. 纯阻性#感性 B. 纯容性 C. 容性 满分:2 分 12. R1和R2为两个串联电阻,已知R1=4R2,若R1上消耗的功率为1W,则R2上消耗的功率为(c)。 A. 5W B. 20W C. 0.25W D. 400W 满分:2 分 13. 在三相四线制电路中,通常所说的220伏和380伏指的是电压的(C )值 A. 平均 B. 最大 C. 有效 D. 瞬时 满分:2 分 14. 若C = 1F ,某时刻电容两端电压u为2V,则此时流过电容的电流=(D)。 A. 2A B. 0.5A C. -2A D. 不能确定 满分:2 分

西安交通大学计算方法B上机试题

1.计算以下和式:01421181 84858616n n S n n n n ∞ =?? =--- ?++++??∑ ,要求: (1)若保留11个有效数字,给出计算结果,并评价计算的算法; (2)若要保留30个有效数字,则又将如何进行计算。 (1)题目分析 该题是对无穷级数求和,因此在使用matlab 进行累加时需要一个累加的终止条件。这里令?? ? ??+-+-+-+= 681581482184161n n n n a n n ,则 ()()1.016 1 6855844864816114851384128698161 681581482184161148113811282984161111<< ? ??? ????? ??++++++???? ????? ??++++++=??? ????? ??+-+-+-+??? ????? ??+-+-+-+=+++n n n n n n n n n n n n n n n n a a n n n n n n 故近似取其误差为1+≈k a ε,并且有m -1m -111021 21 ?=?=≈+βεk a , (2)算法依据 使用matlab 编程时用digits 函数和vpa 函数来控制位数。 (3)Matlab 运行程序 %%保留11位有效数字 k1=11; s1=0;%用于存储这一步计算值 for n=0:50 a=(1/16^n)*(4/(8*n+1)-2/(8*n+4)-1/(8*n+5)-1/(8*n+6)); n1=n-1; if a<=0.5*10^(1-k1) break end end; for i=0:1:n1 t=(1/16^i)*(4/(8*i+1)-2/(8*i+4)-1/(8*i+5)-1/(8*i+6)); s1=s1+t; end s11=vpa(s1,k1); disp('保留11位有效数字的结果为:');disp(s11); disp('此时n 值为:');disp(n1); %%保留30位有效数字 clear all; k2=30;

现浇箱梁支架计算书 (midas计算稳定性)

---------------------------------------------------------------最新资料推荐------------------------------------------------------ 现浇箱梁支架计算书(midas计算稳定性) 温州龙港大桥改建工程满堂支架法现浇箱梁设计计算书计算:复核:审核:中铁上海工程局温州龙港大桥改建工程项目经理部2015年12月30日 1/ 24

温州龙港大桥改建工程现浇箱梁支架计算书目录 1 编制依据、原则及范围············· - 1 1.1 编制依据················· - 1 1.2 编制原则················· - 1 1.3 编制范围················· - 2 -2 设计构造··················· - 2 2.1 现浇连续箱梁设计构造··········· - 2 2.2 支架体系主要构造············· - 2 -3 满堂支架体系设计参数取值··········· - 8 3.1 荷载组合················· - 8 3.2 强度、刚度标准·············· - 9 3.3 材料力学参数··············· - 10 -4 计算····················· - 10 4.1 模板计算················· - 11 4.2 模板下上层方木计算············ - 11 4.3 顶托上纵向方木计算············ - 13 4.4 碗扣支架计算··············· - 14 4.5 地基承载力计算·············· - 18 -

西安交大数字图像处理第一次作业

数字图像处理第一次作业 姓名: 班级: 学号: 提交日期:2015年3月13日

摘要 本次报告首先简单阐述了BMP图像格式及其相关数据结构,随后主要完成了作业要求中关于图像处理与计算的各项任务。本次作业以Matlab 2014为平台,通过对lena.bmp,elain.bmp图像文件的编程处理,分别得到了lena.bmp图像的8到1级灰度逐级递减显示,lena.bmp图像的均值和方差,通过近邻、双线性和双三次插值法对lena.bmp进行4倍放大后得到的2048×2048尺寸图像,和对lena.bmp、elain.bmp图像分别进行水平偏移变换和旋转变换后的图像及其4倍插值放大图像。以上任务完成后均得到了预期的结果。

1.Bmp图像格式简介。 1.1 BMP格式概述 BMP(全称Bitmap)是Windows操作系统中的标准图像文件格式,可以分成两类:设备相关位图(DDB)和设备无关位图(DIB),使用非常广。它采用位映射存储格式,除了图像深度可选以外,不采用其他任何压缩,因此,BMP文件所占用的空间很大。BMP文件的图像深度可选lbit、4bit、8bit及24bit。BMP文件存储数据时,图像的扫描方式是按从左到右、从下到上的顺序。 1.2 BMP格式组成 典型的BMP图像文件由四部分组成: (1)位图头文件数据结构,它包含BMP图像文件的类型、显示内容等信息; (2)位图信息数据结构,它包含有BMP图像的宽、高、压缩方法,以及定义颜色等信息;(3)调色板,这个部分是可选的,有些位图需要调色板,有些位图,比如真彩色图(24位的BMP)就不需要调色板; (4)位图数据,这部分的内容根据BMP位图使用的位数不同而不同,在24位图中直接使用RGB,而其他的小于24位的使用调色板中颜色索引值。 1.3 BMP格式对应数据结构 BMP文件由文件头、位图信息头、颜色信息和图形数据四部分组成。 1.3.1 BMP文件头(14字节) BMP文件头数据结构含有BMP文件的类型、文件大小和位图起始位置等信息。 其结构定义如下: typedef struct tagBITMAPFILEHEADER { WORD bfType; // 位图文件的类型,必须为BM(1-2字节) DWORD bfSize; // 位图文件的大小,以字节为单位(3-6字节) WORD bfReserved1; // 位图文件保留字,必须为0(7-8字节) WORD bfReserved2; // 位图文件保留字,必须为0(9-10字节) DWORD bfOffBits; // 位图数据的起始位置,以相对于位图(11-14字节) // 文件头的偏移量表示,以字节为单位 } BITMAPFILEHEADER; 1.3.2 位图信息头(40字节) BMP位图信息头数据用于说明位图的尺寸等信息。 typedef struct tagBITMAPINFOHEADER{ DWORD biSize; // 本结构所占用字节数(15-18字节) LONG biWidth; // 位图的宽度,以像素为单位(19-22字节) LONG biHeight; // 位图的高度,以像素为单位(23-26字节) WORD biPlanes; // 目标设备的级别,必须为1(27-28字节) WORD biBitCount;// 每个像素所需的位数,必须是1(双色),(29-30字节) // 4(16色),8(256色)16(高彩色)或24(真彩色)之一 DWORD biCompression; // 位图压缩类型,必须是0(不压缩),(31-34字节) // 1(BI_RLE8压缩类型)或2(BI_RLE4压缩类型)之一 DWORD biSizeImage; // 位图的大小(其中包含了为了补齐行数是4的倍数而添加的空字节),以字节为单位(35-38字节) LONG biXPelsPerMeter; // 位图水平分辨率,每米像素数(39-42字节) LONG biYPelsPerMeter; // 位图垂直分辨率,每米像素数(43-46字节)

西安交大课程考试《政治经济学》在线作业

一、单选题(共26 道试题,共52 分。)V 1. 资本周围速度()、 A. 与周转时间成正比,与周转次数成反比 B. 与周转时间成反比,与周转次数成正比 C. 与周转时间成正比,与周转次数成正比 D. 与周转时间成反比,与周转次数成反比 标准答案:B 2. 把剩余价值看成是全部预付资本的产物,剩余价值就转化为()、 A. 成本价格 B. 利润 C. 生产价格 D. 利息 标准答案:B 3. 社会资本简单再生产的基本实现条件是()、 A. I(c+v+m)=I(c+△c)+II(c+△c) B. I(v+m)>IIc C. I(v+m)=IIc D. I(v+m)<IIc 标准答案:C 4. 划分固定资本和流动资本的根据是()、 A. 生产资本不同部分在价值形成中的作用不同 B. 生产资本不同部分在剩余价值生产中的作用不同 C. 生产资本不同部分的价值转移方式不同 D. 生产资本不同部分的价值周转方式不同 标准答案:D 5. 我国在发展对外经济关系中必须坚持的根本方针是()、 A. 独立自主自力更生 B. 自筹资金没有外债 C. 大量进口不怕逆差 D. 自给自足剩余出口 标准答案:A 6. 生产关系的基础是()、 A. 产品归谁所有 B. 人们在生产中的地位及相互关系 C. 生产力 D. 生产资料所有制 标准答案:D 7. 农业雇佣工人创造的剩余价值中,土地所有者获得的是()、 A. 全部剩余价值 B. 租金 C. 平均利润 D. 超额利润 标准答案:D 8. 社会主义制度的经济基础是()、 A. 股份制

B. 国家所有制 C. 多种经济形式 D. 生产资料公有制 标准答案:D 9. 资本主义经济危机的实质是()、 A. 生产绝对过剩 B. 生产相对不足 C. 生产相对过剩 D. 资本主义基本矛盾 标准答案:C 10. 国家垄断资本主义的实质是()、 A. 国家政权同私人垄断资本相结合 B. 国家直接掌握垄断资本 C. 国有资本和私人资本在企业内、外部的结合 D. 私人垄断资本利用国家机器来为其服务的手段 标准答案:D 11. G-W…W`-G`是()、 A. 货币资本循环公式 B. 生产资本循环公式 C. 商品资本循环公式 D. 流通资本循环公式 标准答案:A 12. 垄断资本主义阶段,资本主义对外经济关系的一个重要经济特征是()、 A. 商品输出 B. 原料输出 C. 资本输出 D. 劳动力输出 标准答案:C 13. 社会主义按劳分配的对象是()、 A. 生产资料 B. 个人消费品 C. 社会总产品 D. 生存资料 标准答案:B 14. 垄断资本主义阶段,价值规律发挥作用的形式是()、 A. 商品价格围绕价值上下波动 B. 商品价格采取垄断价格形式 C. 商品价格围绕生产价格上下波动 D. 价值规律已经不再起作用 标准答案:B 15. 职能资本家使用借贷资本所获得的平均利润可以分割为()、 A. 利息和企业利润 B. 产业利润和商业利润 C. 利息和超额利润

西安交通大学18年9月课程考试《土力学及地基基础》作业考核试题(100分)

【奥鹏】[西安交通大学]西安交通大学18年9月课程考试《土力学及地基基础》作业考核试题 试卷总分:100 得分:100 第1题,土的抗剪强度取决于土粒间() A、总应力 B、有效应力 C、孔隙水压力 D、黏聚力 正确答案:B 第2题,预制桩的最小间距为() A、2.5 D B、3 D C、3.5 D D、4 D 正确答案:B 第3题,在土的抗剪强度的影响因素中,最重要的是() A、剪切速率 B、应力历史 C、排水条件 D、应力状态 正确答案:C 第4题,绝对柔性基础在均匀受压时,基底反力分布图形简化为() A、矩形 B、抛物线形 C、钟形 D、马鞍形 正确答案:A 第5题,某黏性土的液性指数IL=0.5则该土的软硬状态为( ) A、硬塑 B、可塑 C、软塑 D、流塑 正确答案:B 第6题,无黏性土坡的稳定性应()

A、与坡高有关,与坡角无关 B、与坡角有关,与坡高无关 C、与坡高和坡角都无关 D、与坡高和坡角都有关 正确答案:B 第7题,计算自重应力时,地下水位以下的土层应采用( ) A、湿重度 B、饱和重度 C、有效重度 D、天然重度 正确答案:B 第8题,设地基的最终变形量为100mm,则当变形量为50mm时,地基的平均固结度为() A、30% B、50% C、80% D、100% 正确答案:B 第9题,在土的压缩性指标中() A、压缩系数a与压缩模量成正比 B、压缩系数a与压缩模量成反比 C、压缩系数越大,土的压缩性越低 D、压缩模量越小,土的压缩性越低 正确答案:B 第10题,利用静载荷试验确定地基承载力特征值时,同一土层参加统计的试验点数不应少于() A、3点 B、5点 C、8点 D、10点 正确答案:A 第11题,下列说法中,错误的是() A、压缩试验的排水条件为双面排水 B、压缩试验不允许土样产生侧向变形 C、在压缩试验中土样既有体积变形,也有剪切变形

西安交大数理统计作业(完整版)

第一章 1.1 X~N(μ,2 σ) 则X~N(μ, 2 n σ ),所以X-μ~N(0, 2 n σ ) P{X-μ <1}= P{ = 0.95 N(0,1),而(0.975) 1.96 Φ= 所以n最小要取[2 1.96x2σ]+1 1.2 (1)至800小时,没有一个元件失效 这个事件等价于P{ 123456 X X X X X X>800}的概率 由已知X服从指数分布,可求得P{ 123456 X X X X X X>800}=7.2 e-(2)至3000小时,所有六个元件都失效的概率 等价与P{ 123456 X X X X X X<3000}的概率 可求得P{ 123456 X X X X X X<3000}= 4.56 (1) e- - 1.5 2 1 () n i i X a = - ∑=2 1 [()()] n i i X X X a = -+- ∑ =22 111 ()2()()() n n n i i i i i X X X a X X X a === -+--+- ∑∑∑ 因为 1 () n i i X X = - ∑=0 所以2 1 () n i i X a = - ∑=22 11 ()() n n i i i X X X a == -+- ∑∑ =22 1 () n i nS X a = +- ∑ 所以当a=X时,2 1 () n i i X a = - ∑有最小值且等于2nS 1.6 (1)由 1 1n i i X X n= =∑

有等式的左边= 221 12n n i i i i X X n μμ==-+∑∑ 等式的右边= 22221122n n i i i i X X X nX nX nX n μμ==-++-+∑∑ = 22 2 2 211 22n n i i i i X nX nX nX X n μμ==-++-+∑∑ = 221 1 2n n i i i i X X n μμ==-+∑∑ 左边等于右边,结论得证。 (2) 等式的左边= 22 11 2n n i i i i X X X nX ==-+∑∑=221 n i i X nX =-∑ 等式的右边= 221 n i i X nX =-∑ 左边等于右边,结论得证。 1.7 (1)由11n n i i X X n ==∑ 及 22 1 1()n n i n i S X X n ==-∑ 有左边=1111111111()1111 n n n n n i i n i i i i X X X X X X n n n n ++++=====+=+++++∑∑∑ 111 ()111 n n n n n nX X X X X n n n ++= +=+-+++=右边 左边等于右边,结论得证。 (2)由 左边=12 21 11 1()1n n i n i S X X n +++==-+∑ 121111[()]11 n i n n n i X X X X n n ++==---++∑ 121111[()()]11 n i n n n i X X X X n n ++==---++∑ 12 2112 1121[()()()()]11(1) n i n i n n n n n i X X X X X X X X n n n +++==----+-+++∑

西安交通大学计算方法B大作业

计算方法上机报告 姓名: 学号: 班级:

目录 题目一------------------------------------------------------------------------------------------ - 4 - 1.1题目内容 ---------------------------------------------------------------------------- - 4 - 1.2算法思想 ---------------------------------------------------------------------------- - 4 - 1.3Matlab源程序----------------------------------------------------------------------- - 5 - 1.4计算结果及总结 ------------------------------------------------------------------- - 5 - 题目二------------------------------------------------------------------------------------------ - 7 - 2.1题目内容 ---------------------------------------------------------------------------- - 7 - 2.2算法思想 ---------------------------------------------------------------------------- - 7 - 2.3 Matlab源程序---------------------------------------------------------------------- - 8 - 2.4计算结果及总结 ------------------------------------------------------------------- - 9 - 题目三----------------------------------------------------------------------------------------- - 11 - 3.1题目内容 --------------------------------------------------------------------------- - 11 - 3.2算法思想 --------------------------------------------------------------------------- - 11 - 3.3Matlab源程序---------------------------------------------------------------------- - 13 - 3.4计算结果及总结 ------------------------------------------------------------------ - 14 - 题目四----------------------------------------------------------------------------------------- - 15 - 4.1题目内容 --------------------------------------------------------------------------- - 15 - 4.2算法思想 --------------------------------------------------------------------------- - 15 - 4.3Matlab源程序---------------------------------------------------------------------- - 15 - 4.4计算结果及总结 ------------------------------------------------------------------ - 16 - 题目五----------------------------------------------------------------------------------------- - 18 -

西安交大《微观经济学》习题三答案

《微观经济学》习题三答案 一、单项选择题(本大题共30小题,每小题2分,共60分) 1、经济学研究的基本问题是( D )。 A.如何在股市上赚钱 B.证明市场可以配置资源 C.选择最公平的收入分配方式 D.稀缺资源的配置与利用问题 2、一国生产可能线以内的一点表示( B )。 A.通货膨胀 B.失业或者资源没有被充分利用 C.该国技术水平提高 D.产品供求平衡 3、正常物品价格上升导致需求量减少的原因在于( C )。 A. 替代效应使需求量增加,收入效应使需求量减少 B. 替代效应使需求量增加,收入效应使需求量增加 C. 替代效应使需求量减少,收入效应使需求量减少 D. 替代效应使需求量减少,收入效应使需求量增加 4、已知一条线性需求曲线,a 点为AB 线段的中点,如图,则( B )。 A. b 点的需求价格弹性等于c 点的需求价格弹性 B. b 点的需求价格弹性大于1 C. b 点的需求价格弹性小于1 D. b 点的需求价格弹性等于1 5、若无差曲线上任何一点的斜率 dY/dX = - 1/ 2 ,这意味着消费者有更多的X 时,他愿意放弃( A )单位X 而获得一单位Y 。 A. 1/2 B.2 C.1 D.1.5 6、下列哪种情况不属消费者均衡的条件( D )。 A. MUx/Px=MUy/ Py = MUz/ Pz B.货币在每种用途上的边际效用相等 C.MUn=λPn D .各种商品的边际效用相等 7、在总产量、平均产量和边际产量的变化过程中,下列何者首先发生( A )。 A.边际产量下降 B.平均产量下降 C.总产量下降 D.B 和C 。 8、从原点出发与TC 曲线相切的直线的斜率是( D )。 P A b a c B Q

西交大计算方法上机报告

计算方法(B)实验报告 姓名: 学号: 学院: 专业:

实验一 三对角方程组Tx f =的求解 一、 实验目的 掌握三对角方程组Tx f =求解的方法。 二、 实验内容 求三对角方程组Tx f =的解,其中: 4 -1 -1 4 -1 -1 4 1 -1 4T ????????=?? ?? ???? , 3223f ?? ? ? ?= ? ? ??? 三、 算法组织 设系数矩阵为三对角矩阵 11222333111 b c a b c a b c a b c b n n n n T ---???????? =?????? ?????? 则方程组Tx f =称为三对角方程组。 设矩阵T 非奇异,T 可分解为T=LU ,其中L 为下三角矩阵,U 为单位上三角矩阵,记 1 1 212 313 1 1 1111 ,11n n n n n r l r l r L U l r l μμμμμ---???? ? ? ? ? ? ?== ? ? ? ? ? ? ? ? ? ?? ? ? ? 可先依次求出,L U 中的元素后,令Ux y =,先求解下三角方程组Ly f =得出 y ,再求解上三角方程组Ux y =。 追赶法的算法组织如下: 1.输入三对角矩阵T 和右端向量f ;

2.将Tx f =压缩为四个一维数组{}{}{}{}i i i i a b c d 、、、,{}{}{}i i i a b c 、、是T 的三对角线性方程组的三个对角,{}i d 是右端向量。将分解矩阵压缩为三个一维数组 {}{}{}i i i l r μ、、。 3.对T 做Crout 分解(也可以用Doolittle 分解)导出追赶法的计算步骤如下: 1111,b r c μ== for 2i n = 111, , ,i i i i i i i i i i i i i l a b a r r c y d l y μμ---==-==- end 4.回代求解x /n n n x y μ= for 11i n =- 1()/i i i i i x y c x μ+=- end 5. 停止,输出结果。 四、 MATLAB 程序 MATLAB 程序见附件1. 五、 结果及分析 实验结果为: (1.0000 1.0000 1.0000 1.0000)T x =

midas截面几何性质计算2

看大家对横向力分布系数计算疑惑颇多,特在这里做一期横向力分布系数计算教程(本教程讲的比较粗浅,适用于新手)。 总的来说,横向力分布系数计算归结为两大类(对于新手能够遇到的): 1、预制梁(板梁、T梁、箱梁) 这一类也可分为简支梁和简支转连续 2、现浇梁(主要是箱梁) 首先我们来讲一下现浇箱梁(上次lee_2007兄弟问了,所以先讲这个吧) 在计算之前,请大家先看一下截面 这是一个单箱三室跨径27+34+27米的连续梁,梁高1.55米,桥宽12.95米!! 支点采用计算方法为为偏压法(刚性横梁法) mi=P/n±P×e×ai/(∑ai x ai) 跨中采用计算方法为修正偏压法(大家注意两者的公式,只不过多了一个β) mi=P/n±P×e×ai×β/(∑ai x ai) β---抗扭修正系数β=1/(1+L^2×G×∑It/(12×E×∑ai^2 Ii) 其中:∑It---全截面抗扭惯距 Ii ---主梁抗弯惯距Ii=K Ii` K为抗弯刚度修正系数,见后 L---计算跨径 G---剪切模量G=0.4E 旧规范为0.43E P---外荷载之合力 e---P对桥轴线的偏心距 ai--主梁I至桥轴线的距离 在计算β值的时候,用到了上次课程https://www.360docs.net/doc/641845795.html,/thread-54712-1-1.html 我们讲到的计算截面几何性质中的抗弯惯矩和抗扭惯矩,可以采用midas计算抗弯和抗扭,也可以采用桥博计算抗弯, 或者采用简化截面计算界面的抗扭,下面就介绍一下这种大箱梁是如何简化截面的: 简化后箱梁高度按边肋中线处截面高度(1.55m)计算,悬臂比拟为等厚度板。 ①矩形部分(不计中肋): 计算公式:It1=4×b^2×h1^2/(2×h/t+b/t1+b/t2) 其中:t,t1,t2为各板厚度

西交计算方法A上机大作业

计算方法A 上机大作业 1. 共轭梯度法求解线性方程组 算法原理:由定理3.4.1可知系数矩阵A 是对称正定矩阵的线性方程组Ax=b 的解与求解二次函数1()2 T T f x x Ax b x =-极小点具有等价性,所以可以利用共轭梯度法求解1()2 T T f x x Ax b x = -的极小点来达到求解Ax=b 的目的。 共轭梯度法在形式上具有迭代法的特征,在给定初始值情况下,根据迭代公式: (1)()()k k k k x x d α+=+ 产生的迭代序列(1)(2)(3)x x x ,,,... 在无舍入误差假定下,最多经过n 次迭代,就可求得()f x 的最小值,也就是方程Ax=b 的解。 首先导出最佳步长k α的计算式。 假设迭代点()k x 和搜索方向()k d 已经给定,便可以通过()()()() k k f x d φαα=+的极小化 ()()min ()()k k f x d φαα=+ 来求得,根据多元复合函数的求导法则得: ()()()'()()k k T k f x d d φαα=?+ 令'()0φα=,得到: ()() ()()k T k k k T k r d d Ad α=,其中()()k k r b Ax =- 然后确定搜索方向()k d 。给定初始向量(0)x 后,由于负梯度方向是函数下降最快的方向,故第一次迭代取搜索方向(0) (0)(0)(0)()d r f x b Ax ==-?=-。令 (1)(0)00x x d α=+ 其中(0)(0)0(0)(0) T T r d d Ad α=。第二次迭代时,从(1) x 出发的搜索方向不再取(1)r ,而是选取(1) (1)(0)0d r d β=+,使得(1)d 与(0)d 是关于矩阵A 的共轭向量,由此可 求得参数0β:

大学计算机基础习题答案(西安交大)

习题参考答案 习题一 1.第一代计算机的主要部件是由(电子管和继电器)构成的。 2.未来全新的计算机技术主要指(光子计算机),(生物计算机)和(量子计算机)。3.按照Flynn分类法,计算机可以分为(单指令流单数据流),(单指令流多数据),(多指令流单数据流)和(多指令流多数据流)4种类型。 4.计算机系统主要由(硬件系统)和(软件系统)组成。 5.说明以下计算机中的部件是属于主机系统、软件系统、还是属于外部设备。 (1)CPU (主机系统) (2)内存条(主机系统) (3)网卡(主机系统) (4)键盘和鼠标(外设) (5)显示器(外设) (6)Windows操作系统(软件系统) 6.控制芯片组是主板的的核心部件,它由(北桥芯片)部分和(南桥芯片)部分组成。7.在计算机系统中设计Cache的主要目的是(提高存去速度)。 8.计算机各部件传输信息的公共通路称为总线,一次传输信息的位数称为总线的(宽度)。9.PCIE属于(系统)总线标准,而SATA则属于(硬盘接口或外设)标准。 10.在微机输入输出控制系统中,若控制的外部设备是发光二极管,最好选用的输入输出方法是(程序控制)方式;若控制的对象是高速设备,则应选则(DMA)控制方式。11.操作系统的基本功能包括(处理器管理或进程管理)、(文件管理)、(存储器管理)、(设备管理)和用户接口。 12.虚拟存储器由(主内存)和(磁盘)构成,由操作系统进行管理。13.CPU从外部设备输入数据需要通过(输入接口),向外设输出数据则需要通过(输出接口)。 14.简述CPU从外部设备输入数据和向外设输出数据的过程。 请参见教材第18页关于输入输出过程的描述。 15.普适计算的主要特点是(是一种无处不在的计算模式)。

迈达斯midas简支梁模型计算

第一讲 简支梁模型的计算 1.1 工程概况 20米跨径的简支梁,横截面如图1-1所示。 图1-1 横截面 1.2 迈达斯建模计算的一般步骤 后处理理处 前 第五步:定义荷载工况 第八步:查看结果 第七步:分析计算第六步:输入荷载 第四步:定义边界条件 第三步:定义材料和截面 第二步:建立单元第一步:建立结点 1.3 具体建模步骤 第01步:新建一个文件夹,命名为Model01,用于存储工程文件。这里,在桌面的“迈达斯”文件夹下新建了它,目录为C:\Documents and Settings\Administrator\桌面\迈达斯\模型01。 第02步:启动Midas Civil.exe ,程序界面如图1-2所示。

第03步:选择菜单“文件(F)->新项目(N)”新建一个工程,如图1-3所示。 图1-3 新建工程 第04步:选择菜单“文件(F)->保存(S)”,选择目录C:\Documents and Settings\Administrator\桌面\迈达斯\模型01,输入工程名“简支梁.mcb”。如图1-4所

示。 图1-4 保存工程 第05步:打开工程目录C:\Documents and Settings\Administrator\桌面\迈达斯\模型01,新建一个excel文件,命名为“结点坐标”。在excel里面输入结点的x,y,z 坐标值。如图1-5所示。 图1-5 结点数据 第06步:选择树形菜单表格按钮“表格->结构表格->节点”,将excel里面的数据拷贝到节点表格,并“ctrl+s”保存。如图1-6所示。

图1-6 建立节点 第07步:打开工程目录C:\Documents and Settings\Administrator\桌面\迈达斯\模型01,再新建一个excel文件,命名为“单元”。在excel里面输入单元结点号。如图 1-6所示。 图1-6 单元节点

西交大有限元原理及应用-大作业

有限元原理及工程应用 ——大作业 学院:机械工程学院 班级:硕4002班 小组成员:李追3114001089 陈草3114001080 2015.5.19

作业题目: 利用有限元方法对简支梁问题进行求解,梁的横截面为矩形,其约束情况如图1所示。 已知梁的几何尺寸和物理参数如下: (1)几何尺寸:长度40cm L =,截面尺寸2cm 02cm .b t ?=?; (2)物理参数:弹性模量70E =GPa ,泊松比0.3ν=,密度3 =2700kg/m ρ。 图1.梁及其横截面示意图 要求: (1) 至少划分五个节点(四个单元); (2) 给出单元节点信息; (3) 给出单元刚度矩阵和质量矩阵; (4) 给出总刚度矩阵和总质量矩阵; (5) 求出梁各界固有频率及振型(五阶); (6) 将所得结果与理论值进行对比,验证方法的可行性。

解:由有限元知识,根据Rayleigh-Ritz 法,解有限元分为四步:建立离散化、单元分析、形 成总体方程、解方程,具体步骤如下: (1)建立离散化 这里我们将矩形截面简支梁等分四等分,即分为六节点的五个杆单元,如图2所示: 每个单元尺寸40 cm=8cm 55 L l = =,这里只考虑杆在竖直平面的弯曲,每个节点只有y 方向位移和绕z 轴的旋转自由度。 (2)单元分析 构造一组Lagrange 插值基函数,在本节点值为1,其他节点值为0。从Rayleigh-Ritz 法可以看到,插值函数要p 次可微,最高阶导数出现在应变能表达式中;同样,我们可以这一原则适用于基函数的选择以及形状函数,否则我们将无法正确计算应变能当我们使用有限元逼近方法。 梁的弯曲问题,应变能计算公式: ()2 220,12L z v x t U EI dx x ???= ???? ? (1-1) 其中,E 为弹性模量,I z 为截面惯性矩。从公式可知,位移函数必须连续,并且二阶导数平方可积。 如图3,是一维杆单元模型,每个节点两个自由度,该单元含有四个自由度,即(,,,i zi j zj v v θθ)。本题中我们采用三次多项式插值函数: ()231234u x x x x αααα=+++ (1-2)

西安交通大学计算方法A实验报告

实验一 矩阵的分解 一、实验目的 掌握矩阵的分解原理和一般方法,学会利用矩阵分解直接求解线性方程组。 二、实验内容 求矩阵() 2020 =ij A α?的T LDL 分解与Cholesky 分解,其中 ,min(,),ij i i j i j i j α=?=? ≠? 。 三、问题分析 1. Cholesky 分解 Cholesky 分解是针对被分解矩阵为对称正定的情况给出的。 分解步骤如下: 11g =1111/y b g =,1111i i g g α= 2i n = ; DO 2j n = jj g = IF 0jj g < STOP ,JUMP TO (5) DO 1i j n =+ 1 1j ij ik kj k ij jj g g g g α-=??- ? ? ?=∑ ji ij g g = 1 1j i ik k k i jj b g y y g -=??- ? ? ?=∑ END DO END DO

2. T LDL 分解 T LDL 分解是针对Cholesky 分解中的开平方运算进行的改进。 分解步骤如下: 11i i r α=,1111/i i r r r =,11y b = 1i n = DO 2i n = DO j i n = 1 1i ij ij ik kj k r l r α-=??=- ??? ∑ /ji ij ii l r r = 1 1i i i ik k k y b l b -=??=- ??? ∑ END DO END DO 四、matlab 求解 分别写出T LDL 分解和Cholesky 分解的函数程序gaijinsqrt.m 和.cholesky m ,调用格 式如下: 1. [index,x,r]=gaijinsqrt(A,b) 参数说明: A 和b 分别是线性代数方程组Ax =b 的系数矩阵和右端向量;输出x 为解向量。 [index,x,g]=Cholesky(A,b) 参数说明: A 和b 分别是线性代数方程组Ax =b 的系数矩阵和右端向量;输出x 为解向量。 然后写出主程序2homework .m 如下: %生成矩阵A A=zeros(20,20); for i=1:20 for j=1:20 if i~=j if i>j A(i,j)=j; else A(i,j)=i; end

相关文档
最新文档