冲击响应谱斜率计算实例

disp(' ')

disp(' srs.m ver 2.0 July 3, 2006')

disp(' by Tom Irvine Email: tomirvine@aol.com') disp(' ')

disp(' This program calculates the shock response spectrum') disp(' of an acceleration time history, which is pre-loaded into Matlab.') disp(' The time history must have two columns: time(sec) & acceleration') disp(' ') %

clear t; clear y; clear yy; clear n; clear fn; clear a1; clear a2 clear b1; clear b2; clear jnum; clear THM; clear resp; clear x_pos; clear x_neg; %

iunit=input(' Enter acceleration unit: 1= G

桥梁计算书计算实例

前言

本设计是根据设计任务书的要求和《公路桥规》的规定,对高坎——上官伯段的高坎大桥进行方案比选和设计的。对该桥的设计，本着“安全、经济、美观、实用”的八字原则，本论文提出两种不同的桥型方案进行比较和选择：方案一为预应力混凝土简支梁桥，方案二为 拱桥。经由以上的八字原则以及设计施工等多方面考虑、比较确定预应力混凝土简支梁桥（锥形锚具）为推荐方案。

1.水文计算

1.1原始资料

1.1.1水文资料：

浑河发源于辽宁省新宾县的滚马苓，从东向西流过沈阳后，折向西南，至海城市三岔河与太子河相汇，而后汇入辽河。浑河干流长364公里，流域面积11085平方公里。本桥位上游45公里的大伙房水库，于1958年建成，该水库控制汇流面积5563平方公里，对沈阳地区的浑河洪峰流量起到很大的削减作用。根据水文部门的资料，建库前浑河的沈阳水文站百年一遇洪峰流量位11700立方米/秒，建库后百年一遇推算值为4780立方米/秒。浑河没年12月初开始结冰，次年3月开始化冻。汛期一般在7月初至9月上旬，河流无通航要求。桥为处河段属于平原区次稳定河段。

1.1.2设计流量

根据沈阳水文站资料，近50年的较大的洪峰流量如下： 大伙房水库建库前 1935年 5

潮流计算实例

潮流例题：根据给定的参数或工程具体要求（如图），收集和查阅资料；学习相关软件（软件自选：本设计选择Matlab进行设计）。

2.在给定的电力网络上画出等值电路图。 3.运用计算机进行潮流计算。 4.编写设计说明书。 一、设计原理

1． 牛顿-拉夫逊原理

牛顿迭代法是取x0 之后，在这个基础上，找到比x0 更接近的方程的跟，一步一步迭代，从而找到更接近方程根的近似跟。牛顿迭代法是求方程根的重要方法之一，其最大优点是在方程f(x) = 0 的单根附近具有平方收敛，而且该法还可以用来求方程的重根、复根。电力系统潮流计算，一般来说，各个母线所供负荷的功率是已知的，各个节点电压是未知的（平衡节点外）可以根据网络结构形成节点导纳矩阵，然后由节点导纳矩阵列写功率方程，由于功率方程里功率是已知的，电压的幅值和相角是未知的，这样潮流计算的问题就转化为求解非线性方程组的问题了。为了便于用迭代法解方程组，需要将上述功率方程改写成功率平衡方程，并对功率平衡方程求偏导，得出对应的雅可比矩阵，给未知节点赋电压初值，一般为额定电压，将初值带入功率平衡方程，得到功率不平衡量，这样由功率不平衡量、雅可比矩阵、节点电压不

1

潮流计算实例

平衡量（未知的）构成了误差

Al的声子谱计算

hhwj340

本文介绍计算Al的声子的计算的过程。使用的软件是ABINIT软

件包(version 7.4.2)。使用的赝势是模守恒赝势(TM型),交换关联能使用 GGA近似(PBE型)。计算声子时，应用的理论方法是线性响应理论。计算的具体参数见输入文件。

Step 1:对Al的结构进行驰豫,输入文件如下：

计算得到晶格常数为：7.6500879384(Bohr)

Step2:获得计算声子谱所有所需K点

将Step1中的晶格常数由7.6500879384代替，设置nstep=1，及nline=1,ngkpt12 12 12,简单运行之后即可以得到所有所需K点。K点可以在生成的out,或者DDB文件中截取，注意如果生成的K点多余50个，out文件中只列出前50个，所以最后在DDB文件中截取。

Step3:计算所有所需Ｋ点的DDB文件，输入文件如下：

qpt21 qpt22 qpt23 qpt24 qpt25 qpt26 qpt27 qpt28 qpt29 qpt30 qpt31 qpt32 qpt33 qpt34 qpt35 qpt36 qpt37 qpt38 qpt39 qpt40 qpt41 qpt42 qpt43 qpt44 qpt45

用VASP计算Ｈ原子的能量

氢原子的能量为 。在这一节中，我们用VASP计算H原子的能量。对于原子计算，我们可以采用如下的INCAR文件 PREC=ACCURATE

NELMDL = 5 make five delays till charge mixing ISMEAR = 0; SIGMA=0.05 use smearing method 采用如下的KPOINTS文件。由于增加K点的数目只能改进描述原子间的相互作用，而在单原子计算中并不需要。所以我们只需要一个K点。 Monkhorst Pack 0 Monkhorst Pack 1 1 1 0 0 0

采用如下的POSCAR文件 atom 1

15.00000 .00000 .00000 .00000 15.00000 .00000 .00000 .00000 15.00000 1 cart

0 0 0

采用标准的H的POTCAR 得到结果如下：

k-point 1 : 0.0000 0.0000 0.0

直线斜率的求法

直线的斜率是反映直线倾斜程度的特征量，在解决有关直线的方程问题中占据着重要的地位.下面例析直线斜率的几种常见求法，以期帮助同学们掌握斜率这一重要知识点.

一、 已知倾斜角定义求

例1 如图，菱形ABCD中， ADC 1200，分别求出BC、CD、AC、BD所在直线的斜率.

分析：准确的找出（或求出）所求直线的倾斜角是关键. 每一条直线都有唯一的倾斜角，直线与横坐标轴正半轴方向的夹角即为该直线的倾斜角.

解：因为在菱形ABCD中， ADC 1200，

所以， BAD 600， ABC 1200，

故kBC tan(1800

1200)=tan60 ；

因为CDPABPx轴，所以直线CD倾斜角为0，故kCD tan00 0; 00

又因为菱形的对角线是相应角的角平分线，

所以 BAC 300， DBA 600，

所以 DBx 1800 DBA 1200，

所以，kAC tan300

kBD tan1200 点评：由直线的倾斜角求斜率，必须正确利用直线的倾斜角与斜率的

000关系：k tan （其中 0,180 且 90）.要注意斜率k随着倾斜角

的变化而变化的趋势：当 00时，k 0；当00 900时，k为正且随着 的增大而增大；当 90

用VASP计算Ｈ原子的能量

氢原子的能量为 。在这一节中，我们用VASP计算H原子的能量。对于原子计算，我们可以采用如下的INCAR文件 PREC=ACCURATE

NELMDL = 5 make five delays till charge mixing ISMEAR = 0; SIGMA=0.05 use smearing method 采用如下的KPOINTS文件。由于增加K点的数目只能改进描述原子间的相互作用，而在单原子计算中并不需要。所以我们只需要一个K点。 Monkhorst Pack 0 Monkhorst Pack 1 1 1 0 0 0

采用如下的POSCAR文件 atom 1

15.00000 .00000 .00000 .00000 15.00000 .00000 .00000 .00000 15.00000 1 cart

0 0 0

采用标准的H的POTCAR 得到结果如下：

k-point 1 : 0.0000 0.0000 0.0

《化工设备设计基础》

课程设计计算说明书

学生姓名： 学 号： 所在学院： 专 业： 设计题目： 指导教师：

2006 年 月 日

目 录

一. 设计任务书…………………………………………………………2 二. 设计参数与结构简图………………………………………………4 三. 设备的总体设计及结构设计………………………………………5 四. 强度计算……………………………………………………………7 五. 设计小结…………………………………………………………..13 六. 参考文献…………………………………………………………..14

1

一、设计任务书

1、设计题目

根据《化工原理》课程设计工艺计算内容进行填料塔（板式塔）设计。各个同学按照自己的工艺参数确定设计题目

内蒙古工业大学本科毕业设计说明书

学校代码： 10128 学 号：200810107019

本科毕业设计说明书

（

题 目：基于MATLAB及VB的无限冲击响应

数字滤波器的设计

学生姓名：

学 院：机械学院 系 别：测控系

专 业：测控技术与仪器 班 级：测控08-2 指导教师： 高级实验师

二○一二年 六月

内蒙古工业大学本科毕业设计说明书

摘 要

随着信息时代和数字世界的到来，数字信号处理已成为当今一门极其重要的学科和技术领域。面对各种各样的原始信号, 若要过滤噪声信号，提取有用信号，这就需要使用滤波器。滤波器的作用主要是选择所需频带的信号内容而抑制不需要的其他频带的信号内容。无限冲击响应数字滤波器（IIR）属于经典数字滤波器的一种，使用MATLAB7.6设计数字滤波器,可以按照设计要求方便地调整设计参数,极大地减轻了设计的工作量,有利于滤波器设计。MATLAB7.6因其强大的数据处理功能被广泛应用于工程计算,其丰富的工具箱为工程计算提供了便利,利用MATLAB7.6信号处理工具箱可以快速有效地设计各种数字滤波器,设计简单方便。VisualBasic是由Microsoft公

判断脱丙烷塔进料热状况

在实际生产中，加入精馏塔的原料液可能有五种热状况：①温度低于泡点的冷液体；②泡点下的饱和液体；③温度介于泡点和露点之间的气液混合物；④露点下的饱和蒸汽；⑤温度高于露点的过热蒸汽。

（1）泡点法计算原料液的泡点温度

设T1=50℃时，由烃类的p-T-K列线图读得的p=2MPa，T1=50℃时各组分的

Ki值列于表3-10中，

表3-10 泡点温度计算表

Table 3-10 the table of bubble point temperature calculation

组成 C20 C3= C30 iC40 C4=-1 iC4= TC4=-2 nC40 C-C4= nC50 合计

Mol% 0.86 33.49 10.07 19.16 6.06 9.1 8.84 4.58 6.46 1.38 100

Ki(50℃) 2.5 1.05 0.92 0.46 0.46 0.5 0.46 0.33 0.2 0.13

） Kixi（50℃0.0215 0.351645 0.092644 0.088136 0.027876 0.0455 0.040664 0.015114 0.01292 0.001794

0.697793

Ki