OpenFOAM计算实例

目录

目录 (1)

1．文献综述 (3)

2．模拟计算方案与过程 (3)

2.1几何模型与数学模型 (3)

2.2计算区域网格划分 (3)

2.3网格文件的导入与转换 (4)

2.3.1计算文件设置 (5)

2.3.2初边值条件设置 (6)

2.3.3物理模型参数设置 (6)

2.4模拟计算 (7)

3．计算结果与分析 (7)

3.1 速度分析 (7)

3.2 两相区域与相界面分析 (10)

3.3 压力分布分析 (12)

4．体会与启示 (15)

5．参考文献 (15)

6．计算文件附录 (15)

水-水蒸气两相相界面的数值模拟

摘要：在OpenFoam平台上对类似散列靶的铅铋合金流动情况的水-水蒸气两相流动相界面进行了数值模拟，根据不同边界条件设置情况做了三组计算，并从计算时间，稳定状态和瞬态三个方面对其做了对比分析。结果表明三组计算在达到稳态后的压力和速度分布以及相界面是基本一致的，区别在于达到稳定所需时间和达到稳态之前的瞬时状态不同。

1．文献综述

2．模拟计算方案与过程

2.1几何模型与数学模型

2.2计算区域网格划分

计算区域为变截面圆管形，即为轴对称三维区域，所以为了节省计算资源和计算时间，在不影响计算结果的情况下将其进行简化，只计算整个圆周中的5度旋转角的楔形区域。因此计网格划分中有楔形网格生成。

OpenFOAM软件有自带的前处理模块blockMesh来生成计算网格，但是blockMesh模块在linux平台下可视化程度不够高，需要通过命令操作。另一方面，windows平台下的前处理软件发展已经十分成熟，且其极高的可视化程度使得windows平台下的各前处理软件已经发展的比较成熟，使用比较方便。综合考虑方便与效率后决定在windows平台下运行Gambit软件进行网格划分，再将网格模型导入OpenFoam软件进行求解计算。

在Gambit软件中划分模型网格，具体过程为：先在Gambit软件中绘制二维平面集合体，如图3所示，再将图3所示的平面集合体分区域划分网格，对于中间部分由于实际过程中运动比较复杂，需要对网格进行加密，所以在绘制网格的时候此部分的网格绘制比较细密。

（a）整体网格

（b）入口段网格

（c）中间段网格

图3 模型网格

Gambit软件中绘制平面网格多种组合，具体情况如表1所示。

表1 Gambit软件中网格划分组合

其中，“√”表示可以采用的组合方法，为了获得较高质量的网格，选择的组合为Quad+Map，所划分的网格全为平行四边形，共划分网格数为22648。

对于由Gambit划分的网格模型，需要在Gambit软件中指定边界类型，随后再在求解器中进行具体边界条件设置。表2所示为网格模型边界类型设置情况。

在Gambit软件中将网格划分完毕并设置边界类型之后，导出网格文件。以其中一个网格数为22648的网格为例，将其命名为xt-ads.2.3w。

2.3网格文件的导入与转换

Gambit软件绘制的网格文件默认情况下是以fluent软件为目标求解器的，所以讲其绘制的网格文件导入OpenFoam软件时需要进行网格转换操作。值得注意的一点是，OpenFoam 软件只对三维网格进行求解计算，对于导入的二维网格，其将在第三维方向上指定网格数为一个，也即指定原始二维网格为在第三维方向为一层的三维网格。在读取二维Fluent网格

时，OpenFoam将在第三维方向上生成两个与原始二维平面平行且几何形状完全一致的两个平面，将其分别命名为frontAndBackPlanes_Pos和frontAndBackPlanes_Neg，而将原始的二维平面类型指定为empty类型，在求解计算中该平面不参与计算[5]。

Gambit和Fluent软件中默认的长度单位为毫米，而OpenFoam软件中默认长度单位为米，因此在将Gambit绘制的网格模型导入到OpenFoam软件中时要注意单位的转换。单位的转换通过scale命令来实现。

由于计算求解的问题不是简单的二维或者三维问题，而是旋转轴对称问题，这就需要将三维的单层网格转换为楔形网格。对于三维单层网格向楔形网格的转换，典型安装模式下的OpenFoam软件时不能实现的，需要新编译一个MakeAxialMesh模块来实现。MakeAxialMesh 模块的功能即为将平面二维网格转换为轴对称网格。

在进行二维平面网格向轴对称网格的抓换过程中，会在楔形体的轴附近产生多需的边，这就需要用collapseEdges命令来清除这些多余的边，collapseEdges命令的成功使用需要在system文件夹下放置如下三个文件：collapseDict，collapseEdges.c和meshQualityDict，文件具体内容见附件。

Gambit软件绘制的网格模型导入到OpenFoam软件并转换为符合求解条件的网格的操作步骤如下：

(1)建立文件夹2D.2.3w，将网格文件xt-ads.2.3w.msh放置于此文件夹下，并在此文件

夹下新建立文件夹system，system文件夹下需要有控制文件ControlDict。

(2)在Ubuntu终端进入到步骤（1）中建立的2D.2.3w文件夹，依次输入以下命令：

①fluentMeshToFoam xt-ads.2.3w.msh -scale 0.001

②makeAxialMesh -axis axis123 -wedge frontAndBackPlanes

③collapseEdges

其中，在分别执行上诉三条命令后会在文件夹2D.2.3w内生分别成新的文件夹constant、2e-5和4e-5。在执行完三条命令之后需要将这就需要将constant和2e-5文件夹删除，并将文件夹4e-5更名为constant。

(3)经过第（2）步后通过命令checkMesh检查网格，如果检查结果为Mesh OK，则表

明已经成功地将Gambit软件划分的二维平面网格转换为了OpenFoam软件能够使

用的轴对称网格。

2.3.1计算文件设置

OpenFoam软件的计算求解需要具备符合具体求解要求的文件结构，对于一般的模拟计算，其求解计算的算例文件夹下需有三个文件夹：0，constant和system。其中0文件夹为初始场与边界场文件夹，即包含0时刻时的各边界条件，如压力和速度等；constant文件夹包含一些固定不变的参数，如网格信息，物理模型及其参数，物理边界定义等；system文件夹为系统文件夹，包含求解计算的运行文件，离散格式和代数方程求解器等文件。图4

为计算的文件结构。

图4

算例文件结构

文件的详细设置情况参见附件。

2.3.2 初边值条件设置 本次计算的0文件夹中包含alpha.water ，epsilon ，k ，nut ，p_rgh 和U 文件，分别为水的体积分数文件，湍动和湍动能文件，初始压力和初始速度文件。Epsilon 和k 的值需要通过如下公式求解计算而得：

k =2

3

(u inlet T i )2

ε=

C μ34k 3

2l

此处的l 取值为散列靶纵向尺寸长度，即取l =0.7m 。初边值条件设置情况如表3所示。

表3 初边值条件

2.3.3 物理模型参数设置

物理模型参数的设置主要包括重力加速度和物性参数设置。由于模型网格模型重力加速

度方向为y轴负方向，所以重力加速度值应设置为（0，-9.81，0）m/s2；物性参数的设置中将水的气蚀压强设置为2370Pa；表面张力设置为0.07kg/s2。计算求解过程中使用k- epsilon 湍流模型。

2.3.4离散格式与求解器设置

对控制方程的离散，需要根据实际问题分别对时间项，对流项和压力项以及耗散项等分别设置离散格式。本文所所求解的问题为两相流动问题，相应离散格式选择入下表4所示。

Euler离散格式对时间为一阶隐式差分，Gauss linear离散格式对梯度项的离散精度为2阶，Gauss linear corrected离散格式对拉普拉斯项离散精度为2阶。对于散度项（耗散项），需要对控制方程的具体项设置不同的高斯离散格式。

对控制方程的离散格式的求解，需要对不同的参数方程选择不同的求解器。本次模拟计算需要求解压力场，速度场以及组分百分数分布场。但是由于流动时湍流问题，需要将速度场的求解与湍流动能场耦合求解，因此对压力场采用线性求解器（Linear solver）进行求解，对速度，湍流动能耦合场和液态水百分数分布场选择光滑求解器（smooth solvers）求解。此外，PCG和PBiCG线性求解器分别用于求解对称和非对称问题，本次求解的问题为轴对称问题，所以选择PCG线性求解器。

2.4模拟计算

主要模拟计算了网格数为22648的模型，共做了三组有效模拟，设置情况如表5所示。计算中需要对比计算区域初始速度场对计算结果的影响以及不同压力出口边界条件设置得到的结果。

表5 三组算例的设置情况

除去表5所示差异外，三组计算文件的其他条件设置完全一致。

3．计算结果与分析

在求解计算过程中，三组算例的自动时间步长是基本一致的，稳定计算时时间步长在3x10-5秒附近波动。

由库劳特数（courant）的定义

?t|U|

Co=

库劳特数和网格间距一定时，时间步长与运动速度成反比。

模拟计算结果分别以速度、压力和组分百分数来进行评价。图5至图7分别为三组计算的速度、压力和组分百分数分布情况随时间变化的情况。

3.1 速度分析

图5所示为三组计算的速度分布随时间的变化情况。

（a）case1速度分布

（b）case2速度分布

（c）case3速度分布图5 不同时刻速度分布图

从图5可以看出，液态水气化后其上部蒸汽速度较大，且远大于主流区液态水的运动速度。蒸汽运动速度在靠近顶部压力入口处为最大，越往下越小。而中部靠右侧的蒸汽运动速度则和水的运动速度差别不大，各处的值也差别不大。对比分析可以发现，在0.5秒之前，三个case的蒸汽运动速度变化都较大，但在0.5秒以后就都趋于稳定状态。Case1和case2的运动速度在0.3秒左右就达到稳定状态，而case3的运动速度在0.4秒左右才达到稳定，相对于case1和case2滞后了0.1秒。

3.2 两相区域与相界面分析

图6所示为三组计算的两相分布随时间的变化情况。

（a）case1两相分布

（b）case2两相分布

（c）case3 两相分布

图6 两相分布随时间变化情况

从图6可以看出，三个算例在最终达到稳定后的两相分布是基本一致的，两个气态水蒸气区域分别分布在压力入口附近和中部两侧，且两个相界面的相对位置也都基本一致。但三者在达到稳定状态之前的瞬态过程是不一致的，且达到稳定状态所需的时间也是不一致的。

对比case1和case3可以发现，二者差异主要在于1.0秒之前的两相分布情况。case1中两相分布在1.5秒左右达到稳定状态，中部只有靠近外壁面处有气化核心，气化速率较case3的快，而case3在中部不仅在靠近外壁面处出现了气化核心，而且中部靠近轴心处也出现了气化核心，随着时间增加，靠轴心的气化核心逐渐消失，而靠近外壁面的气化核心逐渐长大并逐渐与case1的气化核心趋于一致；对比case1和case2可以发现，case1相对于case2能够较快地达到稳定状态，case1达到稳定所需时间为1.5秒，而case2则需要2.5秒。Case1中部气化核心相界面只经过一次波动即达到稳定，而case3经过多次逐渐衰减的振动才达到稳定。

对比图6和文献[3]的计算结果发现，到达稳态后图6中重心相界面与轴线相交处的角度较小，显得比较尖细，而文献[3]中这一角度较大。这一现象可能是水的表面张力的设置差异造成的

3.3 压力分布分析

图7所示为三组计算的压力分布随时间的变化情况。

（a）case1压力分布

（b）Case2压力分布

（c）case3压力分布

图7 压力分布随时间变化情况

从图7可知，流动区域中速度入口区域附近的压力一直是最大的，case1与case3的压力在大约1.5秒时刻达到稳定，而case2的压力在2.5秒左右才达到经稳定，这和图6中相界面的稳定过程是类似的。

将三组模拟计算结果与文献[1]的结果进行对比，发现三组算例在达到稳定后的两相分布情况与文献提供的计算结果都是基本一致的，而case1和case3的差别仅仅在于达到稳定之前的瞬态分布，内部初始速度场的设置并没有影响到稳态时的计算结果与达到稳定所需时间，而case1和case3达到稳定所需时间比case2的短的多，而稳定状态则基本一致。所以出于

经济与高效的考虑，应该采用case1或case3的设置来进行计算，即将压力出口设置为总压。4．体会与启示

5．参考文献

[1] A.G. Class, D. Angeli, A. Batta, etc. XT-ADS Windowless spallation target thermohydraulic design &

experimental setup, Journal of Nuclear Materials, 415, 378-384, 2011.

[2] A. Batta, A.G. Class, Numerical investigations on Geometrical Designs of the Windowless XT ADS

spallation target, Proceedings of ICAPP 2007, Nice, France, May 13-18, 2007.

[3] A. Batta, A.G. Class, Free surface modeling and simulation of the water experiment for the XT ADS

spallation target, Proceedings of ICAPP’ 08, Anaheim, CA USA, June 8-12, 2008.

[4] F. Roelofs, N.B. Siccama, H. Jeanmart, K. van Tichelen, M. Dierckx, P. Schuurmans, Application of a

controlled swirl in the XT-ADS spallation target, Proceedings of ICAPP’ 08, Anaheim, CA USA, June 8-12, 2008.

[5] OpenFoam User Guide.

6．计算文件附录

潮流计算(matlab)实例计算

潮流例题：根据给定的参数或工程具体要求（如图），收集和查阅资料；学习相关软件（软件自选：本设计选择Matlab进行设计）。 2.在给定的电力网络上画出等值电路图。 3.运用计算机进行潮流计算。 4.编写设计说明书。 一、设计原理 1．牛顿-拉夫逊原理 牛顿迭代法是取x0 之后，在这个基础上，找到比x0 更接近的方程的跟，一步一步迭代，从而找到更接近方程根的近似跟。牛顿迭代法是求方程根的重要方法之一，其最大优点是在方程f(x) = 0 的单根附近具有平方收敛，而且该法还可以用来求方程的重根、复根。电力系统潮流计算，一般来说，各个母线所供负荷的功率是已知的，各个节点电压是未知的（平衡节点外）可以根据网络结构形成节点导纳矩阵，然后由节点导纳矩阵列写功率方程，由于功率方程里功率是已知的，电压的幅值和相角是未知的，这样潮流计算的问题就转化为求解非线性方程组的问题了。为了便于用迭代法解方程组，需要将上述功率方程改写成功率平衡方程，并对功率平衡方程求偏导，得出对应的雅可比矩阵，给未知节点赋电压初值，一般为额定电压，将初值带入功率平衡方程，得到功率不平衡量，这样由功率不平衡量、雅可比矩阵、节点电压不平衡量（未知的）构成了误差方程，解误差方程，得到节点电压不平衡量，节点电压加上节点电压不平衡量构成新的节点电压初值，将新的初值带入原来的功率平衡方程，并重新形成雅可比矩阵，然后计算新

的电压不平衡量，这样不断迭代，不断修正，一般迭代三到五次就能收敛。 牛顿—拉夫逊迭代法的一般步骤： （1）形成各节点导纳矩阵Y。 （2）设个节点电压的初始值U和相角初始值e 还有迭代次数初值为0。 （3）计算各个节点的功率不平衡量。 （4）根据收敛条件判断是否满足，若不满足则向下进行。 （5）计算雅可比矩阵中的各元素。 （6）修正方程式个节点电压 （7）利用新值自第（3）步开始进入下一次迭代，直至达到精度退出循环。 （8）计算平衡节点输出功率和各线路功率 2．网络节点的优化 １）静态地按最少出线支路数编号 这种方法由称为静态优化法。在编号以前。首先统计电力网络个节点的出线支路数，然后，按出线支路数有少到多的节点顺序编号。当由n 个节点的出线支路相同时，则可以按任意次序对这n 个节点进行编号。这种编号方法的根据是导纳矩阵中，出线支路数最少的节点所对应的行中非零元素也２）动态地按增加出线支路数最少编号在上述的方法中，各节点的出线支路数是按原始网络统计出来的，在编号过程中认为固定不变的，事实上，在节点消去过程中，每消去一个节点以后，与该节点相连的各节点的出线支路数将发生变化(增加，减少或保持不变)。因此，如果每消去一个节点后，立即修正尚未编号节点的出线支路数，然后选其中支路数最少的一个节点进行编号，就可以预期得到更好的效果，动态按最少出线支路数编号方法的特点就是按出线最少原则编号时考虑了消去过程中各节点出线支路数目的变动情况。 3．MATLAB编程应用 Matlab 是“Matrix Laboratory”的缩写，主要包括：一般数值分析，矩阵运算、数字信号处理、建模、系统控制、优化和图形显示等应用程序。由于使用Matlab 编程运算与人进行科学计算的思路和表达方式完全一致，所以不像学习高级语言那样难于掌握，而且编程效率和计算效率极高，还可在计算机上直接输出结果和精美的图形拷贝，所以它的确为一高效的科研助手。 二、设计内容 1.设计流程图

《分布式计算、云计算与大数据》习题参考解答

第1章分布式计算概述 一、选择题 1，CD 2，ABC 3，ABCD 4，ACD 二、简答题 1，参考1.1.1和节 2，参考1.1.2节 3，分布式计算的核心技术是进程间通信，参考1.3.2节 4，单播和组播 5，超时和多线程 三、实验题 1.进程A在进程B发送receive前发起send操作 进程A进程B 发出非阻塞send操 作，进程A继续运行 发出阻塞receive操 作，进程B被阻塞进程B在进程A发起send前发出receive操作

发出非阻塞send 操作，进程A 继续运行 发出阻塞receive 操作，进程B 被阻塞 收到进程A 发送的数据，进程B 被唤醒 2. 进程A 在进程B 发送receive 前发起send 操作 进程A 进程B 发出阻塞send 操作， 进程A 被阻塞 发出阻塞receive 操作，进程B 被阻塞 进程B 在进程A 发起send 前发出receive 操作

发出阻塞send操作，进程A被阻塞 发出阻塞receive操作，进程B 被阻塞 收到进程A发送的数据，进程B 被唤醒 收到进程B返回的数 据，进程A被唤醒 3.1).在提供阻塞send操作和阻塞receive操作的通信系统中在提供非阻塞send操作和阻塞receive操作的通信系统中2)．P1，P2，P3进程间通信的顺序状态图 m1 m1 m2 m2 第2章分布式计算范型概述 1.消息传递，客户-服务器，P2P，分布式对象，网络服务，移动代理等 2.分布式应用最广泛最流行的范型是客户-服务器范型，参考节

3.分布式应用最基本的范型是消息传递模型，参考节 4.参考节，P2P应用有很多，例如Napster，迅雷，PPS网络电视等 5.参考节 6.参考节 7.略 8.消息传递模式是最基本的分布式计算范型，适用于大多数应用；客户-服务器范型是最 流行的分布式计算范型，应用最为广泛；P2P范型又称为对等结构范型，使得网络以最有效率的方式运行，适用于各参与者地位平等的网络；分布式对象范型，是抽象化的远程调用，适用于复杂的分布式计算应用等。 9.略 10.中间件又称为代理，中间件为参与对象提供内容抽象，隐藏对象引用，起到中介作用。 11.略 第3章 Socket编程与客户服务器应用开发 一、填空题 1.数据包socket，流式socket 2.无连接方式，面向连接方式 3.数据层，业务层，应用层 4.迭代服务器和并发服务器 5.有状态服务器和无状态服务器 二、简答题 1.API：Application Programming Interface,应用程序编程接口，是一些预先定义 的函数，目的是提供应用程序与开发人员基于某软件或硬件得以访问一组例程的能 力，而又无需访问源码，或理解内部工作机制的细节 Socket API：套接字应用程序编程接口，适用于进程间通信的套接字应用程序编程 接口

计算机潮流计算程序代码

G(5,5)=6.250;B(5,5)=-18.750;G(5,2)=-5.00;B(5,2)=15.000; G(5,3)=-1.250;B(5,3)=3.750;G(5,4)=0;B(5,4)=0; G(5,1)=0;B(5,1)=0; G(2,5)=-5.00;B(2,5)=15.000;G(2,2)=10.834;B(2,2)=-32.500;G(2,3)=-1.677;B(2,3)=5.00; G(2,4)=-1.667;B(2,4)=5.00;G(2,1)=-2.500;B(2,1)=7.500; G(3,5)=-1.250;B(3,5)=3.750;G(3,2)=-1.667;B(3,2)=5.00;G(3,3)=12.917;B(3,3)=-38.750; G(3,4)=-10.000;B(3,4)=30.000;G(3,1)=0;B(3,1)=0; G(4,5)=0;B(4,5)=0;G(4,2)=-1.667;B(4,2)=5.000;G(4,3)=-10.000;B(4,3)=30.000; G(4,4)=12.917;B(4,4)=-38.750;G(4,1)=-1.250;B(4,1)=3.750; G(1,5)=0;B(1,5)=0;G(1,2)=-2.500;B(1,2)=7.500;G(1,3)=0;B(1,3)=0; G(1,4)=-1.250;B(1,4)=3.750;G(1,1)=3.750;B(1,1)=-11.250; Y=G+j*B delt(1)=0; delt(2)=0; delt(3)=0; delt(4)=0; u(1)=1.0; u(2)=1.0; u(3)=1.0; u(4)=1.0; ps(2)=0.2;qs(2)=0.20;ps(3)=-0.45;qs(3)=-0.15; ps(4)=-0.4;qs(4)=-0.05; ps(1)=-0.6;qs(1)=-0.1; k=1;precision=1; N1=4; while precision>0.0001 delt(5)=0;u(5)=1.06; for m=1:N1 for n=1:N1+1 pt(n)=u(m)*u(n)*(G(m,n)*cos(delt(m)-delt(n))+B(m,n)*sin(delt(m)-delt(n))); qt(n)=u(m)*u(n)*(G(m,n)*sin(delt(m)-delt(n))-B(m,n)*cos(delt(m)-delt(n))); end pi(m)=sum(pt);qi(m)=sum(qt); dp(m)=ps(m)-pi(m); dq(m)=qs(m)-qi(m); end for m=1:N1 for n=1:N1 if m==n H(m,m)=-qi(m)-u(m)^2*B(m,m); N(m,m)=pi(m)+u(m)^2*G(m,m); J(m,m)=pi(m)-u(m)^2*G(m,m); L(m,m)=qi(m)-u(m)^2*B(m,m); JJ(2*m-1,2*m-1)=H(m,m); JJ(2*m-1,2*m)=N(m,m); JJ(2*m,2*m-1)=J(m,m); JJ(2*m,2*m)=L(m,m); else H(m,n)=u(m)*u(n)*(G(m,n)*sin(delt(m)-delt(n))-B(m,n)*cos(delt(m)-delt(n))); J(m,n)=-u(m)*u(n)*(G(m,n)*cos(delt(m)-delt(n))+B(m,n)*sin(delt(m)-delt(n))); N(m,n)=-J(m,n);L(m,n)=H(m,n); JJ(2*m-1,2*n-1)=H(m,n);JJ(2*m-1,2*n)=N(m,n); JJ(2*m,2*n-1)=J(m,n); JJ(2*m,2*n)=L(m,n); end

造价案例分析 第六讲

第六讲工程价款结算与竣工决算 内容提要 1．现行建筑安装工程价款结算方法； 2．竣工决算； 重点难点 1．现行建筑安装工程价款结算方法； 2．工程预付款及计算； 3，工程价款调整方法； 4．竣工决算； 5．新增资产的构成及其价值的确定 内容讲解 第六章工程款结算与竣工决算 本章考核知识点有： (1)工程价款借款与。支付。 (2)资金使用计划编制与投资偏差分析。 (3)竣工决算内容与编制方法。 (4)新增资产的分类及其价值的确定。 本章在2003年考试中占1题，约占20分。2004年的考核重点有： (1)掌握工程预付款及起扣点的计算。熟悉预付备料款的限额计算及扣回方法。例如，备料款的扣回可以从未施工工程尚需的主材及构件的价值相当于备料款数额时起扣，从每次结算价款中，按材料比重抵扣工程价款，竣工前全部扣清。其次，也可以在承包商完成金额累计达到合同总价的10％后，由承包商开始向发包方还款，发包方从每次应付给承包商的金额中扣回工程预付款，发包方至少应在合同规定的完工日期前三个月将工程预付款的总计金额按逐次分摊的方法扣回。 (2)掌握工程进度款的支付方法，以按月结算为例，施工企业可以在旬末或月中向建设单位提出预支工程款账单，预支一句或半月的工程款，月终再提出工程款结算帐单和已完工程月报表，收取当月工程款。 (3)熟悉工程造价指数调整法、调值公式法等工程价款价差调整方法o (4)掌握新增固定资产价值的确定方法。其中，其他费用的分摊应按单项工程以一定比例分摊。例如，对建设单位管理费由建筑工程、安装工程、需安装设备价值总额按比例分摊，土地征用费、勘察设计费等只按建筑工程造价比例分摊。 (5)熟悉固定资产、流动资产、无形资产、递延资产的含义及常见代表性费用。 第一节工程价款的结算 一、工程价款的主要计算方式 工程价款的结算，合同中可根据不同的情况采取不同的形式。教材介绍了四种方式。 (1)按月结算。这是一种以分部分项工程为对象，实行旬末或月中预交、月终结算、竣工后清算的结算办法。按月结算是案例分析中重点要掌握的结算办法。

电力系统分析潮流计算例题

电力系统的潮流计算 西安交通大学自动化学院 2012.10 3．1 电网结构如图3—11所示，其额定电压为10KV 。已知各节点的负荷功率及参数： MVA j S )2.03.0(2 +=， MVA j S )3.05.0(3+=， MVA j S )15.02.0(4+= Ω+=)4.22.1(12j Z ，Ω+=)0.20.1(23j Z ，Ω+=)0.35.1(24j Z 试求电压和功率分布。 解：（1）先假设各节点电压均为额定电压，求线路始端功率。 0068.00034.0)21(103.05.0)(2 2223232232323j j jX R V Q P S N +=++=++=?0019.00009.0)35.1(10 15.02.0)(2 2 224242242424j j jX R V Q P S N +=++=++=?

则： 3068.05034.023323j S S S +=?+= 1519.02009.024424j S S S +=?+= 6587.00043.122423' 12 j S S S S +=++= 又 0346 .00173.0)4.22.1(106587.00043.1)(2 2 212122'12'1212j j jX R V Q P S N +=++=++=? 故： 6933.00216.112'1212 j S S S +=?+= （2） 再用已知的线路始端电压kV V 5.101 =及上述求得的线路始端功率 12 S ，求出线 路 各 点 电 压 。

kV V X Q R P V 2752.05 .104.26933.02.10216.1)(11212121212=?+?=+=? kV V V V 2248.101212=?-≈ kV V V V kV V X Q R P V 1508.100740.0) (24242 2424242424=?-≈?=+=? kV V V V kV V X Q R P V 1156.101092.0) (23232 2323232323=?-≈?=+=? （3）根据上述求得的线路各点电压，重新计算各线路的功率损耗和线路始端功率。 0066.00033.0)21(12.103.05.02 2 223j j S +=++=? 0018.00009.0)35.1(15 .1015.02.02 2 224j j S +=++=? 故 3066.05033.023323j S S S +=?+= 1518.02009.024424j S S S +=?+= 则 6584.00042.122423' 12 j S S S S +=++= 又 0331.00166.0)4.22.1(22 .106584.00042.12 2 212j j S +=++=? 从而可得线路始端功率 6915.00208.112 j S +=

潮流计算matlab程序

clear; %各节点参数：节点编号，类型，电压幅值，电压相位，注入有功，注入无功%类型：1＝PQ节点，2＝PV节点，3＝平衡节点 %本程序中将最后一个节点设为平衡节点 R_1=[1 1 1.0 0 0.2 0.2j; 2 1 1.0 0 -0.45 -0.15j; 3 1 1.0 0 -0.45 -0.05j; 4 1 1.0 0 -0.6 -0.1j; 5 3 1.0 0 0 0]; %支路号首端节点末端节点支路导纳 R_2=[1 5 2 1.25-3.75j; 2 2 3 10.00-30.00j; 3 3 4 1.25-3.75j; 4 1 4 2.50-7.50j; 5 1 5 5.00-15.00j; 6 1 2 1.667-5.00j]; n=5;L=6;%需要改变的到此为止 i=0;j=0;a=0;precision=1;k=0; Y=zeros(n,n);u=zeros(1,n);delt=zeros(1,n);P=zeros(1,n);Q=zeros(1,n); G=[];B=[];PP=[];uu=[];U=[];dp=[];dq=[]; for a=1:L i=R_2(a,2); j=R_2(a,3); Y(i,j)=-R_2(a,4); Y(j,i)=Y(i,j); end for a=1:n for b=1:n if a~=b Y(a,a)=Y(a,a)+Y(a,b); end end end for i=1:n for j=1:n if i==j Y(i,j)=-Y(i,j); end end end Y %形成导纳矩阵 for i=1:n for j=1:n G(i,j)=real(Y(i,j));

2002年造价工程师《工程造价案例分析》试题

2002年造价工程师执业资格考试 工程造价案例分析考试试题 试题一：（25分） 某建设项目有关资料如下： 1.项目计算期10年，其中建设期2年。项目第3年投产，第5年开始达到100％设计生产能力。 2.项目固定资产投资9000万元（不含建设期贷款利息和固定资产投资方向调节税），预计8500万元形成固定资产，500万元形成无形资产。固定资产年折旧费为673万元，固定资产余值在项目运营期末收回，固定资产投资方向调节税税率为0。 3.无形资产在运营期8年中，均匀摊入成本。 4.流动资金为1000万元，在项目计算期末收回。 5.项目的设计生产能力为年产量1.1万吨，预计每吨销售价为6000元，年销售税金及附加按销售收入的5％计取，所得税率为33％。 6.项目的资金投入、收益、成本等基础数据，见表1.1。 表1.1建设项目资金投入、收益及成本表单位：万元 7.还款方式：在项目运营期间（即从第3年到第10年）按等额本金法偿还，流动资金贷款每年付息。长期贷款利率为6.22％（按年付息），流动资金贷款利率为3％。 8.经营成本的80％作为固定成本。 问题： 1.计算无形资产摊销费。 2.编制借款还本付息表，把计算结果填入答题纸表1.2中（表中数字按四舍五入取整，表1.3、表1.4同）。 3.编制总成本费用估算表，把计算结果填入答题纸表1.3中。 4.编制项目损益表，把计算结果填入答题纸表1.4中。盈余公积金提取比例为10％。

5.计算第7年的产量盈亏平衡点（保留两位小数）和单价盈亏平衡点（取整），分析项目盈利能力和抗风险能力。 试题二：（25分） 某房地产开发公司对某一地块有两种开发方案。 A方案：一次性开发多层住宅45000平方米建筑面积，需投入总成本费用（包括前期开发成本、施工建造成本和销售成本，下同）9000万元，开发时间（包括建造、销售时间，下同）为18个月。 B方案：将该地块分为东、西两区分两期开发。一期在东区先开发高层住宅36000平方米建筑面积，需投入总成本费用8100万元，开发时间为15个月。二期开发时，如果一期销路好，且预计二期销售率可达100％（售价和销量同一期），则在西区继续投入总成本费用8100万元开发高层住宅36000平方米建筑面积；如果一期销路差，或暂停开发，或在西区改为开发多层住宅22000平方米建筑面积，需投入总成本费用4600万元，开发时间为15个月。两方案销路好和销路差时的售价和销量情况汇总于表2.1。 根据经验，多层住宅销路好的概率为0.7，高层住宅销路好的概率为0.6。暂停开发每季损失10万元。季利率为2％。

调值公式法在公路工程价格调整中的应用

调值公式法在公路工程价格调整中的应用 摘要：通过工程实例，介绍了公路工程价格调整的必要性，以及调值公式法在价格调整中的应用，并且介绍了在确定固定因子和调价成本因素的权重系数过程中应考虑的因素，用实际数据演示了调值公式法的计算过程，为公路从业人员掌握调值公式法提供帮助。 主题词：公路价格调整调值公式法必要性应用 采用调价公式法对合同进行调整是国际性招标项目的惯例。根据FIDIC通用合同条款第70条第一款的规定，在合同履行期间，随着劳务、材料或影响工程施工成本的任何其他事项的价格涨落而引起施工成本增减，根据合同专用条件规定的调价公式给予价格调整计算，在合同价格基础上增加或减少相应的金额。对于国内公路招标项目，交通部《公路工程国内招标文件范本》（交公路发[2003]94号）合同条款 70.1规定：凡是合同预期工期在24个月以上者，在合同执行期内由于人工和材料的价格涨落因素，应对合同价格进行调整，调价时，应按公式计算，每年进行一次调整。国家发改委等七部委发布的工程建设项目施工招标投标办法》（[2003]第30号令）第30条规定施工招标项目工期超过12个月的，招标文件中可以规定工程造价指数体系，价格调整因素和调整方法。 一、公路工程价格调整的必要性 对于合同工期较长（一般规定为超过12个月）的公路工程，采用价格调整可以由业主和承包商公平合理地分担价格意外风险，减少承包商由于材料价格上涨引起的风险和损失，有利于承包商为保证工程质量和工程进度而进行必要的投入。并且使业主能获得比较真实可靠的报价，以及在工程决算时能在一个合理的价格水平上承受工程费用，确保双方能顺利履行合同。笔者所参建的新疆奎屯—赛里木湖高等级公路工程为世界银行贷款项目，合同工期为2002年7月10日至2005年10月20日，时间跨度达四年。在几年中，受国内政策、伊拉克战争和国际原油价格以及供需多方面的影响，公路建设项目主要材料价格变化明显，涨幅较大，如钢筋：2002年平均价格（Ⅰ级钢筋占25%，Ⅱ级钢筋占75%）为2601元/吨，2004年平均价格为3860元/吨，涨幅达48.4%；交通石油沥青：2002年平均价格为1644元/吨，2004年平均价格为1935元/吨，较2002年上涨17.7%；0#柴油：2002年平均价格为3.00元/公斤，2004年平均价格为3.71元/公斤，较2002年上涨23.67%，尤其是钢材、汽柴油和沥青处于持续上涨态势。对于如此巨大的材料价格涨幅，所有施工单位均难以承受，亏损严重，公司总部垫资维持，因此造成2004年上半年奎赛公路项目整体进度滞后。针对此种情况，新疆自治区交通厅、高管局和奎赛项目执行办公室积极采取措施，对固定因子以及参与价格调整的因素和价格指数进行了优化和补充，并与世行进行了谈判、确认，使其更符合实际情况，从而扭转了不利局面，使工程顺利竣工通车。 二、公路工程价格调整的方法 价格调整的计算方法主要有：①、实际价格结算法②、按主材计算价差③、主材按抽料计算价差，其他材料按系数计算材差④、竣工调价系数法⑤、调值公式法。由于前三种计算方法存在工程实际材料消耗数量统计不准确，各合同段实际配合比以及材料用量计算有差异等弊端，因此目前极少应用，而第四种方法又缺乏工程进行过程的动态调整，因此调值公式法是目前应用最广泛的，国际和国内招标中推荐的调价方法。 三、调值公式法的计算方法 1、首先确定参与价格调整的因素（或指标）。 选择对工程投资、工程成本影响较大且投入数量较多，并且预计在合同工期内有可能发生较大价格变动的主要材料和成本费用因素作为调价指标。如新疆维吾尔自治区国道218

基于MATLAB的潮流计算源程序代码(优.选)

%*************************电力系统直角坐标系下的牛顿拉夫逊法潮流计算********** clear clc load E:\data\IEEE014_Node.txt Node=IEEE014_Node; weishu=size(Node); nnum=weishu(1,1); %节点总数 load E:\data\IEEE014_Branch.txt branch=IEEE014_Branch; bwei=size(branch); bnum=bwei(1,1); %支路总数 Y=(zeros(nnum)); Sj=100; %********************************节点导纳矩阵******************************* for m=1:bnum; s=branch(m,1); %首节点 e=branch(m,2); %末节点 R=branch(m,3); %支路电阻 X=branch(m,4); %支路电抗 B=branch(m,5); %支路对地电纳 k=branch(m,6); if k==0 %无变压器支路情形 Y(s,e)=-1/(R+j*X); %互导纳 Y(e,s)=Y(s,e); end if k~=0 %有变压器支路情形 Y(s,e)=-(1/((R+j*X)*k)); Y(e,s)=Y(s,e); Y(s,s)=-(1-k)/((R+j*X)*k^2); Y(e,e)=-(k-1)/((R+j*X)*k); %对地导纳 end Y(s,s)=Y(s,s)-j*B/2; Y(e,e)=Y(e,e)-j*B/2; %自导纳的计算情形 end for t=1:nnum; Y(t,t)=-sum(Y(t,:))+Node(t,12)+j*Node(t,13); %求支路自导纳 end G=real(Y); %电导 B=imag(Y); %电纳 %******************节点分类************************************* * pq=0; pv=0; blancenode=0; pqnode=zeros(1,nnum); pvnode=zeros(1,nnum); for m=1:nnum; if Node(m,2)==3 blancenode=m; %平衡节点编号 else if Node(m,2)==0 pq=pq+1; pqnode(1,pq)=m; %PQ 节点编号 else if Node(m,2)==2 pv=pv+1; pvnode(1,pv)=m; %PV 节点编号 end end end end %*****************************设置电压初值********************************** Uoriginal=zeros(1,nnum); %对各节点电压矩阵初始化 for n=1:nnum Uoriginal(1,n)=Node(n,9); %对各点电压赋初值 if Node(n,9)==0;

调值公式实例

掌握建筑工程竣工结算 竣工调值公式法： 小案例: 1.某项目2008年5月1日签署合同价1000万,合同约定固定系数0.2。2009年5月1日结算时发现,钢材涨价13%,水泥涨价16%,其中钢材和水泥分别占调值部分的比例分别为25%和30%。求结算价。 答： )8.0%45100 100 8.0%301001168.0%251001132.0(1000??+??+??+ =P 2.某项目 2008年5月1日签署合同价1000万,合同约定固定系数0.2。2009年5月 1日结算时发现,钢材涨价13%,水泥涨价16%，人工涨价22％,其中钢材、水泥和人工分别占调值部分的比例分别为25%、30%和35％。求结算价。 答： ) 8.0%10100100 8.0%351001228.0%301001168.0%251001132.0(1000??+??+??+??+=P 例:某项目有效合同价为2000万元，第一个分项工程土石方费用为460万元，第二个分项工程钢筋混凝土为1540万元，各指标费用占土石方、钢筋混凝土两个分项工程比例如下表，不调值的费用占工程价款的15%。 各分项工程费用比例、有关物价指标及指数表

若2001年9月完成的工程价款占合同价的10%，试计算该月的调价额为多少？ [解 ]各项参加调值的费用占工程价值的比例为： 人工费权重系数： 机械费权重系数： 钢筋权重系数： 水泥权重系数： 取上一个月的物价指数 2001年9月的调价额为： P197【案例2A320054—2】 练习题 案例一 背景：某项工程业主与承包商签订了工程承包合同，合同中含有两个子项工程，估算工程量甲项为2300m 3，乙项为3200m 3，经协商甲项单价为180元/m 3，乙项为160元/m 3。承包合同规定： 1.开工前业主应向承包商支付工程合同价20%的预付款； 2.业主自第一个月起，从承包商的工程款中，按5%的比例扣留质量保修金； 3.当子项工程累计实际工程量超过估算工程量10%时，可进行调价，调价系数为0.9；

电力系统分析潮流计算课程序设计及其MATLAB程序设计

电力系统分析潮流计算程序设计报告题目：13节点配电网潮流计算 学院电气工程学院 专业班级 学生姓名 学号 班内序号 指导教师房大中 提交日期 2015年05月04日

目录 一、程序设计目的 (2) 二、程序设计要求 (4) 三、13节点配网潮流计算 (4) 3.1主要流程................................................................................................ 错误！未定义书签。 3.1.1第一步的前推公式如下（1-1）-（1-5）： .................................. 错误！未定义书签。 3.1.2第二步的回代公式如下（1-6）—（1-9）： ................................ 错误！未定义书签。 3.2配网前推后代潮流计算的原理 (7) 3.3配网前推后代潮流计算迭代过程 (8) 3.3计算原理 (9) 四、计算框图流程 (10) 五、确定前推回代支路次序....................................................................................... 错误！未定义书签。 六、前推回代计算输入文件 (11) 主程序： (11) 输入文件清单： (12) 计算结果： (13) 数据分析： (13) 七、配电网潮流计算的要点 (14) 八、自我总结 (14) 九、参考文献 (15) 附录一 MATLAB的简介 (15)

《分布式计算、云计算与大大数据》习题参考解答

《分布式计算、云计算与大数据》习题解答参考第1章分布式计算概述 一、选择题 1，CD 2，ABC 3，ABCD 4，ACD 二、简答题 1，参考1.1.1和1.1.2节 2，参考1.1.2节 3，分布式计算的核心技术是进程间通信，参考1.3.2节 4，单播和组播 5，超时和多线程 三、实验题 1.进程A在进程B发送receive前发起send操作 进程A进程B 发出非阻塞send操 作，进程A继续运行 发出阻塞receive操 作，进程B被阻塞 进程B在进程A发起send前发出receive操作

发出非阻塞send 操作，进程A 继续运行 发出阻塞receive 操作，进程B 被阻塞 收到进程A 发送的数据，进程B 被唤醒 2. 进程A 在进程B 发送receive 前发起send 操作 进程A 进程B 发出阻塞send 操作， 进程A 被阻塞 发出阻塞receive 操作，进程B 被阻塞 进程B 在进程A 发起send 前发出receive 操作

发出阻塞send操作，进程A被阻塞 发出阻塞receive操作，进程B 被阻塞 收到进程A发送的数据，进程B 被唤醒 收到进程B返回的数 据，进程A被唤醒 3.1).在提供阻塞send操作和阻塞receive操作的通信系统中 receive operation send operation t=1 在提供非阻塞send操作和阻塞receive操作的通信系统中

t=1 receive operation send operation 2)．P1，P2，P3进程间通信的顺序状态图 m1 m1 m2 m2 第2章分布式计算型概述 1.消息传递，客户-服务器，P2P，分布式对象，网络服务，移动代理等 2.分布式应用最广泛最流行的型是客户-服务器型，参考2.2节 3.分布式应用最基本的型是消息传递模型，参考2.1节 4.参考2.3节，P2P应用有很多，例如Napster，迅雷，PPS网络电视等 5.参考2.4节 6.参考2.7节 7.略 8.消息传递模式是最基本的分布式计算型，适用于大多数应用；客户-服务器型是最流行 的分布式计算型，应用最为广泛；P2P型又称为对等结构型，使得网络以最有效率的方

潮流计算程序

最大负荷时 %本程序的功能是用牛顿——拉夫逊法进行潮流计算 % B1矩阵：1、支路首端号；2、末端号；3、支路阻抗；4、支路对地电纳 % 5、支路的变比；6、支路首端处于K侧为1，1侧为0 % B2矩阵：1、该节点发电机功率；2、该节点负荷功率；3、节点电压初始值% 4、PV节点电压V的给定值；5、节点所接的无功补偿设备的容量 % 6、节点分类标号 clear; n=12;%input('请输入节点数：n='); nl=14;%input('请输入支路数：nl='); isb=1;%input('请输入平衡母线节点号：isb='); pr=0.00001;%input('请输入误差精度：pr='); B1=[1 2 3.3+11.48i 0 1 0; 1 4 3.87+15.375i 0 1 0; 1 6 5.8+18.8i 0 1 0; 2 3 2.39+34.95i 0 1.02 1; 4 5 1.59+23.3i 0 1.02 1; 2 4 4.67+18.34i 0 1 0; 6 7 2.39+24.35i 0 1.02 1; 8 9 2.39+24.35i 0 1.02 1; 10 11 2.39+24.35i 0 1.02 1; 8 10 1.7+11.2i 0 1 0; 4 12 4.2+13.9i 0 1 0; 6 12 4.2+14.64i 0 1 0; 8 12 2.2+13.9i 0 1 0; 10 12 2.4+15.2i 0 1 0]; %input('请输入由支路参数形成的矩阵： B1='); B2=[0 0 231 231 0 1; 0 0 220 0 0 2; 0 100+75i 220 0 0 2; 0 0 220 0 0 2; 0 100+75i 220 0 0 2; 0 0 220 0 0 2; 0 150+112.5i 220 0 0 2; 0 0 220 0 0 2; 0 150+112.5i 220 0 0 2; 0 0 220 0 0 2; 0 150+112.5i 220 0 0 2; 450 0 231 231 0 3];%input('请输入各节点参数形成的矩阵： B2='); Y=zeros(n);e=zeros(1,n);f=zeros(1,n);V=zeros(1,n);sida=zeros(1,n);S1= zeros(nl); %－－－－－－－修改部分－－－－－－－－－－－－ ym=1;

2020年一级造价工程师《建设工程造价案例分析(土木建筑工程、安装工程)》 章节题库(工程结算与决算

2020年一级造价工程师《建设工程造价案例分析（土木建筑工程、安装工程）》章节题库 第六章工程结算与决算下 案例十二 某住宅楼工程合同价款为1000万元，工程价款按季进行动态结算，结算时人工单价、机械台班单价不予调整，主要材料采用调值公式法给承包商以调价补偿。合同报价日期为2010年3月份，开工日期为4月1日，计划年底完工。各季度主要材料价格指数如表6-19所示。 表6-19各季度主要材料价格指数 各季度的产值情况如表6-20所示。 表6-20各季度产值表 并且合同中规定了以下内容： （1）开工前业主预付承包商20%工程款，预付款在第三、四季度结算时平均扣回；

（2）业主从第二季度起，按承包商每季度工程款的3%扣留保留金，保修期满无质量问题后退还承包商。 实际施工中发生以下情况： （1）4月份由于公网停电使得工期拖延2天，承包商损失2万元；在基础开挖过程中，因工程受特大暴雨袭击，导致承包商费用损失3万元，工期延误5天。 （2）5月份为赶回由于特大暴雨袭击导致的工期延误，工程师发布加速施工命令，双方协商加速施工费用通过工期奖罚加以补偿，在本季度结算，提前完工一天，奖励5000元；结果承包商加速施工使工程按原定计划完成了第二季度任务。 （3）7月份承包商租赁的某种施工机械因运输问题推迟运到现场，造成承包商窝工损失5万元，工期延误4天。 （4）10月份工程师要求对某隐蔽工程再做剥露检验，经查不合格，导致承包商延误工期1天，损失5000元。 【问题】 1．承包商对实际施工中发生的各事项可提出工期索赔几天？费用索赔多少？ 答：各事项索赔天数及费用索赔如下： （1）4月份因公网停电是业主责任，可索赔2天，费用索赔2万；特大暴雨属客观原因，工期可索赔5天，费用不可索赔。 （2）加速施工的工期、费用不可索赔。 （3）机械问题是承包商责任，工期、费用不可索赔。 （4）隐蔽工程质量问题是承包商责任，工期、费用不可索赔。

matlab潮流计算

附录1 使用牛顿拉夫逊法进行潮流计算的Matlab程序代码 % 牛拉法计算潮流程序 %----------------------------------------------------------------------- % B1矩阵：1、支路首端号；2、末端号；3、支路阻抗；4、支路对地电纳 % 5、支路的变比；6、支路首端处于K侧为1，1侧为0 % B2矩阵：1、该节点发电机功率；2、该节点负荷功率；3、节点电压初始值% 4、PV节点电压V的给定值；5、节点所接的无功补偿设备的容量% 6、节点分类标号：1为平衡节点（应为1号节点）；2为PQ节点；3为PV节点； %------------------------------------------------------------------------ clear all; format long; n=input('请输入节点数:nodes='); nl=input('请输入支路数:lines='); isb=input('请输入平衡母线节点号:balance=');

pr=input('请输入误差精度:precision='); B1=input('请输入由各支路参数形成的矩阵:B1='); B2=input('请输入各节点参数形成的矩阵:B2='); Y=zeros(n);e=zeros(1,n);f=zeros(1,n);V=zeros(1,n);sida=zeros(1,n);S1=zeros(nl); %------------------------------------------------------------------ for i=1:nl %支路数 if B1(i,6)==0 %左节点处于1侧 p=B1(i,1);q=B1(i,2); else %左节点处于K侧 p=B1(i,2);q=B1(i,1); end Y(p,q)=Y(p,q)-1./(B1(i,3)*B1(i,5)); %非对角元 Y(q,p)=Y(p,q); %非对角元 Y(q,q)=Y(q,q)+1./(B1(i,3)*B1(i,5)^2)+B1(i,4); %对角元K侧 Y(p,p)=Y(p,p)+1./B1(i,3)+B1(i,4); %对角元1侧 end %求导纳矩阵

第三章简单电力系统的潮流计算汇总

第一章 简单电力系统的分析和计算 一、 基本要求 掌握电力线路中的电压降落和功率损耗的计算、变压器中的电压降落和功率损耗的计 算；掌握辐射形网络的潮流分布计算；掌握简单环形网络的潮流分布计算；了解电力网络的简化。 二、 重点内容 1、电力线路中的电压降落和功率损耗 图3-1中，设线路末端电压为2U 、末端功率为222~jQ P S +=，则 （1）计算电力线路中的功率损耗 ① 线路末端导纳支路的功率损耗： 222 2* 222~U B j U Y S Y -=?? ? ??=? ……………（3-1） 则阻抗支路末端的功率为： 222~~~Y S S S ?+=' ② 线路阻抗支路中的功率损耗： ()jX R U Q P Z I S Z +'+'==?2 2 22222 ~ ……（3-2） 则阻抗支路始端的功率为： Z S S S ~ ~~21?+'=' ③ 线路始端导纳支路的功率损耗： 2121* 122~U B j U Y S Y -=?? ? ??=? …………（3-3） 则线路始端的功率为： 111~ ~~Y S S S ?+'= ~~~图3-3 变压器的电压和功率 ~2 ? U （2）计算电力线路中的电压降落 选取2U 为参考向量，如图3-2。线路始端电压 U j U U U δ+?+=2 1 其中 2 2 2U X Q R P U '+'= ? ； 222U R Q X P U '-'=δ ……………（3-4） 则线路始端电压的大小： ()()2 221U U U U δ+?+= ………………（3-5） 一般可采用近似计算： 2 2 2221U X Q R P U U U U '+'+ =?+≈ ………………（3-6）

潮流计算设计

毕业设计（论文）题目配电网潮流计算与程序设计 完成日期年月日

目录 1绪论 (4) 1.1潮流计算的目的和意义 (4) 1.2潮流计算的发展 (4) 1.3本文的主要工作 (6) 2 配电网络模型 (6) 2.1线路模型 (6) 2.2变压器的模型 (9) 2.3 负荷模型 (13) 3基于前推回代法的潮流计算 (14) 4基于matlab的程序设计 (15) 4.1 MATLAB简介 (15) 4.2 程序设计 (16) 4.3 算例 (20) 5 总结 (24)

配电网潮流计算与程序设计 摘要：在电力系统的正常运行中，随着用电负荷的变化和系统运行方式的改变，网络中的损耗也将发生变化。要严格保证所有的用户在任何时刻都有额定的电压是不可能的，因此系统运行中个节点出现电压的偏移是不可避免的。为了保证电力系统的稳定运行，要进行潮流调节。 Abstract:In the normal operation of the power system, with the change of power load and the change of system operation mode, the network losses will also change. To strictly ensure that all users at any time have a rated voltage is impossible, so the node appears in the operation of the system voltage deviation is inevitable. In order to guarantee the stable operation of power system, must carry on the trend of the regulation. 关键词: 电网潮流计算前推后代法MATLAB Keywords:Power flow calculation Push the offspring before MATLAB 前言 电力是以电能作为动力的能源。发明于19世纪70 年代，电力的发明和应用掀起了第二次工业化高潮。成为人类历史18世纪以来，世界发生的三次科技革命之一，从此科技改变了人们的生活。既是是当今的互联网时代我们仍然对电力有着持续增长的需求，因为我们发明了电脑、家电等更多使用电力的产品。不可否认新技术的不断出现使得电力成为人们的必需品。 20世纪出现的大规模电力系统是人类工程科学史上最重要的成就之一，是由发电、输电、变电、配电和用电等环节组成的电力生产与消费系统。它的功能是将自然界的一次能源通过发电动力装置转化成电能，再经输电、变电和配电将电能供应到各用户。为实现这一功能，电力系统在各个环节和不同层次还具有相应的信息与控制系统，对电能的生产过程进行测量、调节、控制、保护、通信和调度，以保证用户获得安全、经济、优质的电能。 电力工业发展初期，电能是直接在用户附近的发电站（或称发电厂）中生产的，各发电站孤立运行。随着工农业生产和城市的发展，电能的需要量迅速增加，而热能资源和水能资源丰富的地区又往往远离用电比较集中的城市和工矿区，为了解决这个矛盾，就需要在动力资源丰富的地区建立大型发电站，然后将电能远距离输送给电力用户。同时，为了提高供电的可靠性以及资源利用的综合经济性，又把许多分散的各种形式的发电站，通过送电线路和变电所联系起来。这种由发电机、升压和降压变电所，送电线路以及用电设备有机连接起来的整体，即称为电力系统。 随着电力系统及在线应用的发展，计算机网络已经形成，为电力系统的潮流计算提供了物质基础。电力系统潮流计算是电力系统分析计算中最基本的内容，也是电力系统运行及设计中必不可少的工具。根据系统给定的运行条件、网络接线及元件参数，通过潮流计算可以确定各母线电压的幅值及相角、各元件中流过