哈工大测试大作业——传感器综合运用——题目四
传感器综合运用
1、 设计题目
如图所示工件，在生产线的30°滑道上自上而下滑落，要求在滑动过程中检测工件厚度，并且计数。图中4mm尺寸公差带为10。
图1.测量工件
2、 厚度检测传感器的选择
电容传感器是把被测的机械量，如位移、压力等转换为电容量变化的传感器。它的敏感部分就是具有可变参数的电容器。其最常用的形式是由两个平行电极组成、极间以空气为介质的电容器（见图）。若忽略边缘效应，平板电容器的电容为εA/δ，式中ε为极间介质的介电常数，A为两电极互相覆盖的有效面积,δ为两电极之间的距离。δ、A、ε 三个参数中任一个的变化都将引起电容量变化，并可用于测量。因此电容式传感器可分为极距变化型、面积变化型、介质变化型三类。极距变化型一般用来测量微小的线位移或由于力、压力、振动等引起的极距变化。面积变化型一般用于测量角位移或较大的线位移。介质变化型常用于物位测量和各种介质的温度、密度、湿度的测定。与电阻式或电感式传感器相比，电容传感器具有四大优点:
(l)分辨力高，常用于精密测量;
(2)动态响应速度快，可以直接用于某些生产线上的动态测量;
(3)从信号源取得的能量少，有利于发挥其测量精度;
(4)机械结构简单，易于实现非接触式测量。因此电容传感器在精密测量中占有重要的地位。
此外，电容器传感器还具有结构简单,价格便宜，灵敏度高，零磁滞，真空兼容，过载能力强，动态响应特性好和对高温、辐射、强振等恶劣条件的适应性强等优点。
因此，在本题中选择电容传感器作为厚度检测传感器。
3、 电容传感器的检测原理
电容式传感器可分为面积变化型、极距变化型、介质变化型三类，下面将分述其检测原理。
1、面积变化型电容传感器
这一类传感器输出特性是线性的，灵敏度是常数。这一类传感器多用于检测直线位移、角位移、尺寸等参量。测量装置如图2所示。
图2.变面积式电容传感器
其电容量计算公式为：
式中 L－外圆筒与内圆柱覆盖部分的长度
外圆筒内半径与内圆柱外半径
灵敏度
2、极距变化型电容传感器
极距变化型电容传感器一般用来测量微小的线位移或由于力、压力、振动等引起的极距变化。测量装置如图3所示。
图3.变极距式电容传感器示意图及其特性曲线
其电容量计算公式为：
式中 S－极板相对覆盖面积；
d －极板间的距离；
ε0－真空介电常数，ε0 =8.85×10-12F.m-1
εr －极板间介质的相对介电常数；
ε－介电常数，ε = ε0 εr 。
灵敏度
非线性误差：
3、介质变化型电容传感器
这种传感器由一对互相平行的金属极板A、B组成，其相对面积为S，极板间隙为，两极板间放置有厚度的金属待测工件。测量装置如图4所示。
图4.变介电常数式电容传感器
其电容量计算公式为：
式中 —总的输出电容量(F)；
—以待测工件为介质的电容量(F)；
—以空气为介质的电容量(F);
又
，
所以得到，
式中:S—极板相对面积;
—空气介质常数，;
—待测工件介质介电常数;
—两极板间距；
—非金属薄膜介质的厚度。
灵敏度，可以看出x/(a-x)越大， r越小，
灵敏度越大，非线性越小。
当待测工件进入电容两极板后，改变了范围内的介电常数，从而改变了电容量。这样电容量C与待测工件厚度形成了确定的关系，这就是变介电常数式电容传感器测厚的基本原理。
4、 电容传感器的设计
电容传感器为两圆柱形测量头，以其两个端面分别作为电容的一极，两者之间形成电容C。若把它看作平行板电容器，则电容器两极板之间的电力线明显地要产生边缘效应。要使C 与间距之间的关系按无穷大平行板电容器来处理，传感器结构上需采取较为特殊的设计。即在测量头外面加一测头保护层，保护屏蔽层与测量头之间电绝缘。传感器采用4Cr13 制作，因4Cr13 的线膨胀系数较低。绝缘层是聚四氟乙烯材料，它的绝缘强度非常高。保护层接电源地。接地后此保护层有利于克服测量头与周围导体的寄生电容的影响，防止外来信号的干扰。与传感器连接的电缆线采用低噪声的同轴屏蔽电缆。由于连接测量头的电缆、屏蔽层与测量头之间形成的电容很大，将把待测电容全部淹没，无法将有用信号拾取出来，而电容传感器的电容很小，又是高阻抗元件，所以测量头屏蔽和传输电缆电容的影响问题是电容传感器的技术关键，因此对传感器需采取屏蔽、在测量线路的前置级需采取驱动电缆技术。5、 测量电路
电容式传感器的常用测量电路有变压器电桥电路，差动脉冲调制电路，双T桥电路，运算放大器是电路，此外电容式传感器的测量电路还有普通交流电桥电路和紧耦合电感臂电桥电路等。由于运算放大器的放大倍数K非常大, 而且输入阻抗Z i很高。运算放大器的这一特点可以使其作为电容式传感器的比较理想的测量电路。此处选用该种测量电路。测量电路图如图5所示。
图5.电容式传感器测量电路
由运算放大器工作原理可得
式中： Cx为电容式传感器是输出信号电压
C0为固定电容是交流电源电压
如果传感器是一只平板电容,则C x=εS/d ,
代入上式, 有。
此式说明运算放大器的输出电压与极板间距离d呈线性关系。运算放大器电路解决了单个变极板间距离式电容传感器的非线性问题，但要求Zi及K足够大。为保证仪器精度, 还要求电源电压的幅值和固定电容C 值稳定。
6、 测量方案简图
采取以上措施以后，再经过传感器线性校正以后，其测量误差可以达到0.001‰，即测量
4mm的工件最大测量误差在4~5左右，满足测量要求。测量方案简图如下：
图6.测量方案简图
七、计数器的设计
1、设计方案
光电计数器采用光电传感器构成的光电门实现对通过光电门的物体进行计数，是一种非接触式计数，在部分场合有着其无比的优越性，从
而使其广泛应用于工业生产、实时监测、自动化控制等领域。本设计为实现光电计数器的功能，采用模数结合的电路，以红外对射光电传感器为传感器件。电路主要分为信号采集电路、两位十进制计数电路、 数码显示电路，分别实现对通过光电门的物体感应、计数、显示。
为实现以上要求，本设计方案采用红外发射管采用直流供电，接收对管判断是否有物体通过光电门，并且当物体通过光电门时输出一个高电平，触发后面的加法计数器，使其加一，为简单起见，计数器为一组BCD码输出，输出由BCD-七段数码管译码器译码，输入至数码管显示。另此设计还有脉冲电路为红外发光二极管供电，使其发射红外脉冲，从而可以提高瞬时功率而使平均功率满足正常工作要求，从而可以加大光电门的宽度，以通过较大的物体。脉冲发射电路中采用的振荡器采用555 构成的多谐荡器，其结构简单，便于调试。
2、 系统框图
本设计采系统结构如下图所示：
图7.系统框图
由图可见，系统可分为信号采集、计数、显示功能块。
3、单元电路设计
（1）信号采集电路设计
该电路主要由脉冲（方波）发生电路，光电转换电路，信号滤波比较电路组成。 脉冲发生电路为555构成的多谐振荡器，振荡频率大约为30KHz，555定时器是一种模拟和数字功能相结合的中规模集成器件，它内部包括两个电压比较器，三个等值串联电阻，一个RS触发器，一个放电管T及功率输出级，它提供两个基准电压，比较器的参考电压由三只5KΩ的电阻器构成分压。A1和A2的输出端控制RS触发器状态和放电管开关状态。当输入信号超过时，触发器复位，555的输出端3脚输出低电平，同时放电，开关管导通；当输入信号自2脚输入并低时，555的3脚输出高电平，同时放电，开关管截止。RD是复位端，当其为 0时，555输出低电平。平时该端开路或接VCC。Vc是控制电压端（5脚），平时输出2VCC作为比较器A1的参考电平，当35脚外接一个输入电压，即改变了比较器的参考电平，从而实现对输出的另一种控制，在不接外加电压时，通常接一个0.01μf的电容器到地，起滤波作用，以消除外来的干扰，以确保参考电平的稳定。T为放电管，当T导通时，将给接于脚7的
电容器提供低阻放电电路。电路如下图：
图8.信号采集电路
图中增加了两个导引二极管，使充放电的电流流过指定的电阻。（2）光电转换电路设计
光电转换电路采用红外发光二极管和光敏三极管作为光电转换传感器，采用典型电路，如下图：
图9.光电转换电路
从 555 多谐振荡器输出的信号控制小功率三极管的导通与截止，从而控制发光二极管发光与否。接收电路由光敏三极管上接电源并加一下拉电阻实现。发光二极管的正常工作压降大约是1.2V，电流大约是
20mA，从而可以估算出限流电阻大约是150Ω，但是由于由脉冲信号控制，工作时间大约占20%，因此瞬时功率可以提高5倍，即限流电阻减小80%，大约是30Ω。光敏三极管在无光照射下，电阻可达输兆欧姆，有足够的关照时，电阻降到几千欧姆，因此其分压电阻可以取几十到几百千欧，具体需根据实际情况调整，这里采用了典型值33KΩ。图中节点3输出含有直流分量的信号，其交流部分大致同步于前面脉冲输出信号。（3）滤波电路设计
光电转换电路输出信号在有、无物件遮挡住光线时，都输出含有直流交流分量的信号。因此通过滤波电路，可以将两种信号转换为不同电压的直流信号，再通过一定阈值的比较器，可以将两种信号转换为数字电路中的高低电平，从而控制后续电路。电路如下图所示：
图10.滤波电路设计
电路中，R5和C3构成滤波电路，为了减小纹波，一般电阻、电容取得较大。但电容和电阻取得过大又会使的在两种输入信号切换时，电路反应迟钝，这里选用了一般值。比较电路采用LM258构成，3脚输入参考电压，即是阈值电压，2脚输入待处理信号。阈值的选取需根据实际情况设定，是电路能够准确的区分两种不同的输入信号电压，因此采用一电位器从电源分压获得，其中电容C4为纹波滤除电容。 综合上述三个单元电路，信号采集电路可以实现抗干扰，防止背景光和瓶子抖动产生计数误差， 电路情况为当有物体通过时，遮挡住红外光线，电路输出高电平，平时输出 低电平。
（4）计数电路
计数电路主要采用计数器统计信号采集电路输出的脉冲个数，实现对物件计数的功能。为了使电路简单化，在此选用74LS190为此电路的
计数器。74LS190为可预置十进制同步加减计数器，其预置是异步的，当置入控制端 LD 为低电平时，不管时钟端CP状态如何，其计数是同步的，靠输出端即预置成与数据输入端相一致的状态。CP同时加在四个触发器上而实现的。当计数控制端CT为低电平时，在CP上升沿作用下同时变化，从而消除了异步计数器中出现的计数尖峰。当计数方式控制U/D 为低电平时进行加计数；当U/D为高电平时进行减计数。只有在CP为高电平时CT和U/D才可以跳变。其有超前进位功能，当计数上溢或下溢时，进位/借位输出端CO/BO输出一个宽度等于脉冲周期的高电平脉冲；行波时钟输出端RC输出一个宽度等于CP低电平部分的低电平脉冲。利用RC端可级联成N位同步计数器。当采用并行时钟控制时，则将RC接到后一级CT；当采用并行CT控制时，则将RC接到后一级CP。
该电路是实现将计数电路的计数值以数字方式显现出来。计数电路输出两组BCD码，为了使电路简单，应选用BCD码——七段数码管译码驱动器。在此选用了常用的 74LS48 共阴数码管驱动器，选取共阴八段数码管，但小数点位不用。74LS48为内部有上拉电阻的BCD—七段数码管译码器/驱动器，输出端为高电平有效，可用于驱动共阴数码管。
哈工大2008-2009年春季学期传热学试题A
哈工大08/09学年春季学期 一、名词解释（20分） 1、导热系数 2、热边界层 3、辐射强度 4、灰体 二、分析论述与回答问题（30分） 1、写出傅里叶导热定律表达式，并说明式中各量和符号的物理意义。 2、简述在对流传热研究中，引入边界层理论的意义。 3、写出努谢尔数Nu与毕渥数Bi的表达式并比较异同。 4、太阳能集热器采用选择性表面涂层，它对太阳辐射的吸收效率为0.9，它本身的 发射率为0.3，这一现象是否违背基尔霍夫定律？为什么？ 5、厚度等于δ的常物性无限大平板，初始温度均为t0，过程开始后，左侧有一定热 流密度q w的热源加热，右侧与低温流体t f相接触（t0>t f），表面传热系数h等于常数，所有物性参数已知，写出改导热问题的数学描写。
三、如图所示的二维稳态导热物体，其导热系数λ为常数，边界面与环境发生对流换热， 环境温度为t F ，边界面对流换热表面传热系数为h ，网格划分如下图所示，试建立数值求解节点温度t 4，t 5，t 6的离散方程。 四、在一个特殊应用中，空气流过一个热的表面，其边界层温度分布可近似为 s s T-T =1-exp(Pr )T -T u y v ∞∞-，其y 是离开表面的垂直距离，普朗特数Pr =0.7u a ∞=是一个无量纲的流体物性。如果来流温度T 400K ∞=，表面温度s T 300K =，且-15000m u v ∞=，求表面热流密度是多少？（10分） （空气导热系数，330K 时，0.0263W/(m K)λ=?；400K 时，0.0339W/(m K)λ=?）
λ=?导热系数的热绝缘层，六、在太空中飞行的宇宙飞船，表面贴有厚0.15m、0.045W/(m K) 而外表面的黑度为0.04，设飞船内空气温度为20℃，空气与内壁间的对流换热系数为6W/(m.K)， 试求飞船外表面的温度。（假设宇宙空间温度为0K；忽略飞船壁面的导热热阻）（15分）
哈工大概率论参考答案习题
习 题 一 1．写出下列随机试验的样本空间及下列事件中的样本点： （1）掷一颗骰子，记录出现的点数. A =‘出现奇数点’； （2）将一颗骰子掷两次，记录出现点数. A =‘两次点数之和为10’，B =‘第一次的点数，比第二次的点数大2’； （3）一个口袋中有5只外形完全相同的球，编号分别为1,2,3,4,5；从中同时取出3只球，观察其结果，A =‘球的最小号码为1’； （4）将,a b 两个球，随机地放入到甲、乙、丙三个盒子中去，观察放球情况，A =‘甲盒中至少有一球’； （5）记录在一段时间内，通过某桥的汽车流量，A =‘通过汽车不足5台’，B =‘通过的汽车不少于3台’。 解 （1）123456{,,,,,}S e e e e e e =其中i e =‘出现i 点’1,2,,6i =， 135{,,}A e e e =。 （2）{(1,1),(1,2),(1,3),(1,4),(1,5),(1,6)S = (2,1),(2,2),(2,3),(2,4),(2,5),(2,6) (3,1),(3,2),(3,3),(3,4),(3,5),(3,6) (4,1),(4,2),(4,3),(4,4),(4,5),(4,6) (5,1),(5,2),(5,3),(5,4),(5,5),(5,6) (6,1),(6,2),(6,3),(6,4),(6,5),(6,6)}； {(4,6),(5,5),(6,4)}A =； {(3,1),(4,2),(5,3),(6,4)}B =。 （3）{(1,2,3),(2,3,4),(3,4,5),(1,3,4),(1,4,5),(1,2,4),(1,2,5)S = (2,3,5),(2,4,5),(1 {(1,2,3),(1,2,4),(1,2,5),(1,3,4),(1,3,5),(1,4,5)}A = （4）{(,,),(,,),(,,),(,,),(,,),(,,),S ab ab ab a b a b b a =--------- (,,),(,,,),(,,)}b a a b b a ---，其中‘-’表示空盒； {(,,),(,,),(,,),(,,),(,,)}A ab a b a b b a b a =------。 （5）{0,1,2,},{0,1,2,3,4},{3,4,}S A B ===。 2．设,,A B C 是随机试验E 的三个事件，试用,,A B C 表示下列事件： （1）仅A 发生； （2）,,A B C 中至少有两个发生；
传热学MATLAB温度分布大作业完整版
传热学大作业（第四章） 姓名：张宝琪学号：03110608 一、题目及要求 1.各节点的离散化的代数方程 2.源程序 3.不同初值时的收敛快慢 4.上下边界的热流量（λ=1W/(m℃）） 5.计算结果的等温线图 6.计算小结 题目：已知条件如下图所示： 二、方程及程序 (1)各温度节点的代数方程 ta=(300+b+e)/4 ; tb=(200+a+c+f)/4; tc=(200+b+d+g)/4; td=(2*c+200+h)/4 te=(100+a+f+i)/4; tf=(b+e+g+j)/4; tg=(c+f+h+k)/4 ; th=(2*g+d+l)/4
ti=(100+e+m+j)/4; tj=(f+i+k+n)/4; tk=(g+j+l+o)/4; tl=(2*k+h+q)/4 tm=(2*i+300+n)/24; tn=(2*j+m+p+200)/24; to=(2*k+p+n+200)/24; tp=(l+o+100)/12 (2)源程序 【G-S迭代程序】 【方法一】 函数文件为： function [y,n]=gauseidel(A,b,x0,eps) D=diag(diag(A)); L=-tril(A,-1); U=-triu(A,1); G=(D-L)\U; f=(D-L)\b; y=G*x0+f; n=1; while norm(y-x0)>=eps x0=y; y=G*x0+f; n=n+1; end 命令文件为： A=[4,-1,0,0,-1,0,0,0,0,0,0,0,0,0,0,0; -1,4,-1,0,0,-1,0,0,0,0,0,0,0,0,0,0; 0,-1,4,-1,0,0,-1,0,0,0,0,0,0,0,0,0;
哈工大概率论与数理统计课后习题答案 一
·1· 习 题 一 1．写出下列随机试验的样本空间及下列事件中的样本点： （1）掷一颗骰子，记录出现的点数. A =‘出现奇数点’； （2）将一颗骰子掷两次，记录出现点数. A =‘两次点数之和为10’，B =‘第一次的点数，比第二次的点数大2’； （3）一个口袋中有5只外形完全相同的球，编号分别为1,2,3,4,5；从中同时取出3只球，观察其结果，A =‘球的最小号码为1’； （4）将,a b 两个球，随机地放入到甲、乙、丙三个盒子中去，观察放球情况，A =‘甲盒中至少有一球’； （5）记录在一段时间内，通过某桥的汽车流量，A =‘通过汽车不足5台’，B =‘通过的汽车不少于3台’。 解 （1）123456{,,,,,}S e e e e e e =其中i e =‘出现i 点’1,2,,6i = ， 135{,,}A e e e =。 （2）{(1,1),(1,2),(1,3),(1,4),(1,5),(1,6)S = (2,1),(2,2),(2,3),(2,4),(2,5),(2,6) (3,1),(3,2),(3,3),(3,4),(3,5),(3,6) (4,1),(4,2),(4,3),(4,4),(4,5),(4,6) (5,1),(5,2),(5,3),(5,4),(5,5),(5,6) (6,1),(6,2),(6,3),(6,4),(6,5),(6,6)}； {(4,6),(5,5),(6,4)}A =； {(3,1),(4,2),(5,3),(6,4)}B =。 （3）{(1,2,3),(2,3,4),(3,4,5),(1,3,4),(1,4,5),(1,2,4),(1,2,5)S = (2,3,5),(2,4,5),(1,3,5)} {(1,2,3),(1,2,4),(1,2,5),(1,3,4),(1,3,5),(1,4,5)}A = （4）{(,,),(,,),(,,),(,,),(,,),(,,),S ab ab ab a b a b b a =--------- (,,),(,,,),(,,)}b a a b b a ---，其中‘-’表示空盒； {(,,),(,,),(,,),(,,),(,,)}A ab a b a b b a b a =------。 （5）{0,1,2,},{0,1,2,3,4},{3,4,}S A B === 。 2．设,,A B C 是随机试验E 的三个事件，试用,,A B C 表示下列事件： （1）仅A 发生； （2）,,A B C 中至少有两个发生；
流体力学大作业
《计算流体力学》课程大作业 作业内容：3-4人为小组完成数值模拟，在第8次课上每组进行成果展示，并在课程结束后每组上交一份纸质版报告。 数值模拟实现形式：自编程或者使用任意的开源、商业模型。 成果展示要求：口头讲述和幻灯片结合的方式，每组限时10分钟（8分钟讲述，2分钟提问和讨论）。 报告要求：按照期刊论文的思路和格式进行撰写（包括但不限于如下内容：摘要、绪论\引言、数值模型简介、数值结果分析\讨论、结论、参考文献）。 （以下题目二选一） 题目一：固定单方柱扰流问题 根据文章《Interactions of tandem square cylinders at low Reynolds numbers》中的实验进行数值模拟，完成但不局限于如下工作： （1）根据Fig. 2 中的雷诺数和方柱排列形式，进行相同雷诺数不同间距比情况下的方柱绕流数值模拟，并做出流线图和Fig.2中的结果对比。 （2）根据Fig. 3 中的雷诺数和方柱排列形式，进行相同雷诺数后柱不同转角情况下的方柱绕流数值模拟，并做出流线图和Fig.3中的结果对比。 （3）根据Fig. 12, 13 中的雷诺数和方柱间距比的设置进行数值模拟，作出频率、斯特劳哈尔数、阻力系数随雷诺数变化的折线并与图中对应的折线画在同一坐标系下比较。 （中共有4条折线，对应4种不同的方柱排列形式下的物理参数随雷诺数变化的规律，仅需选取单柱模型和其中一种双柱模型进行数值模拟，共计16个工况）。 题目二：溃坝问题 根据文章《Experimental investigation of dynamic pressure loads during dam break》中的实验进行数值模拟，完成但不局限于如下工作： （1）分别完成二维、三维的溃坝的数值建模，讨论二维、三维模型的区别。 （2）分别将二维、三维溃坝的数值模拟结果和Fig. 7,10中各时刻的自由面形态进行对比，并分别观测溃坝前端水舌的位置随时间的变化，其结果和Fig. 12 种的各试验结果放在同一坐标系下进行对比。 （3）根据实验设置数值观测点，分别观测与实验测点相对应的数值观测点上的水体高度、压力随时间的变化曲线，并和Fig.16, 18,21,30,31,32,33,35中的实验结果进行对比。
传热学MATLAB温度分布大作业完整版
东南大学能源与环境学院 课程作业报告 作业名称：传热学大作业——利用matlab程序解决热传导问题 院系：能源与环境学院 专业：建筑环境与设备工程 学号： 姓名： 2014年11月9日
一、题目及要求 1.原始题目及要求 2.各节点的离散化的代数方程 3.源程序 4.不同初值时的收敛快慢 5.上下边界的热流量（λ=1W/(m℃）） 6.计算结果的等温线图 7.计算小结 题目：已知条件如下图所示： 二、各节点的离散化的代数方程 各温度节点的代数方程 ta=(300+b+e)/4 ; tb=(200+a+c+f)/4; tc=(200+b+d+g)/4; td=(2*c+200+h)/4 te=(100+a+f+i)/4; tf=(b+e+g+j)/4; tg=(c+f+h+k)/4 ; th=(2*g+d+l)/4 ti=(100+e+m+j)/4; tj=(f+i+k+n)/4; tk=(g+j+l+o)/4; tl=(2*k+h+q)/4
tm=(2*i+300+n)/24; tn=(2*j+m+p+200)/24; to=(2*k+p+n+200)/24; tp=(l+o+100)/12 三、源程序 【G-S迭代程序】 【方法一】 函数文件为： function [y,n]=gauseidel(A,b,x0,eps) D=diag(diag(A)); L=-tril(A,-1); U=-triu(A,1); G=(D-L)\U; f=(D-L)\b; y=G*x0+f; n=1; while norm(y-x0)>=eps x0=y; y=G*x0+f; n=n+1; end 命令文件为： A=[4,-1,0,0,-1,0,0,0,0,0,0,0,0,0,0,0; -1,4,-1,0,0,-1,0,0,0,0,0,0,0,0,0,0; 0,-1,4,-1,0,0,-1,0,0,0,0,0,0,0,0,0;
传热学
镁合金激光-TIG复合热源焊接热源模型 学院：材料学院 专业：材料加工工程 学号： 姓名： 指导教师： 江苏科技大学 2015年4 月11 日
镁合金激光—TIG复合焊接热源模型与热过程 1 前言 镁合金被称为“21世纪绿色工程材料”。镁合金是目前被国内外重新认识并积极开发的一种轻量化材料,具有低密度、高比强度、阻尼减震性好、易机械加工以及良好的可回收性等优点。高效合理的镁合金焊接方法将大大推动镁合金的发展与应用。激光--电弧复合热源焊接具有高速、高效、接头质量优异等特点,目前正在被国内外的研究者日益关注。对这一过程的焊接数值模拟研究有助于更深层次地理解过程的物理机制,从而实现指导焊接工艺、控制焊接质量的目的。目前，YAG激光--TIG复合热源焊接AE31B镁合金已经被证明是一种可行而且高质量的焊接工艺[1], 迫切需要数值模拟工作对这一过程进行指导,并通过数值模拟更深层次的理解复合热源焊接这一过程。 但目前复合热源的数值模拟工作开展的却非常有限。其中一个主要原因是复合热源焊接热源模型一直解决得。首先，高能束激光焊接的热源模型虽然经过线热源、面热源、柱状热源乃至双椭球体热源的变迁,始终没有得到很好的解决; 其次激光、电弧两热源之间存在着一定的物理机制, 需要考虑热源之间的能量影响关系。 在复合热源焊接工艺研究的基础上,结合镁合金材料特点,建立了基于旋转高斯体热源与高斯面热源相结合的复合热源模型:高能束激光热源由旋转高斯体热源描述;TIG电弧则由高斯面热源描述。热源模型的建立充分考虑了过程的物理特点与热源间的能量增强效应。 1.1激光--电弧复合热源焊接概况 激光--TIG电弧复合热源焊接的特点是YAG激光、TIG电弧这两种不同物理性质与能量传输机制的热源同时作用于焊接区。这种方法克服了单独采用激光和单独采用TIG电弧焊接的缺点,并且两种热源相互藕合获得了更大能量形式。其原理如图1.1。其在实践中的优点却是非常明显的:速度快,桥接能力强,焊接变形小,焊接过程稳定,焊接质量和效率高等[2-4]。
运筹学大作业 哈工大
课程名称：对偶单纯形法 一、教学目标 在对偶单纯形法的学习过程中，理解和掌握对偶问题；综合运用线性规划和对偶原理知识对对偶单纯形法与单纯形法进行对比分析，了解单纯形法和对偶单纯形法的相同点和不同点，总结出各自的适用范围；掌握对偶单纯形法的求解过程；并能运用对偶单纯形法独立解决一些运筹学问题。 二、教学内容 1) 对偶单纯形法的思想来源(5min) 2) 对偶单纯形法原理(5min) 3) 总结对偶单纯形法的优点及适用情况(5min) 4) 对偶单纯形法的求解过程(10min) 5) 对偶单纯形法例题(15min) 6) 对比分析单纯形法和对偶单纯形法(10min) 三、教学进程： 1）讲述对偶单纯形法思想的来源： 1954年美国数学家C.莱姆基提出对偶单纯形法（Dual Simplex Method ）。单纯形法是从原始问题的一个可行解通过迭代转到另一个可行解，直到检验数满足最优性条件为止。对偶单纯形法则是从满足对偶可行性条件出发通过迭代逐步搜索原始问题的最优解。在迭代过程中始终保持基解的对偶可行性，而使不可行性逐步消失。因此在保持对偶可行性的前提下，一当基解成为可行解时，便也就是最优解。 2）讲述对偶单纯形法的原理 A.对偶问题的基本性质 依照书第58页，我们先介绍一下对偶问题的六个基本性质： 性质一：弱对偶性 性质二：最优性。如果 x j （j=1...n)原问题的可行解，y j 是其对偶问题可 行解，且有 ∑=n j j j x c 1 =∑=m i i i y b 1 ，则x j 是原问题的最优解，y j 是其对偶问题的最
优解。 性质三：无界性。如果原问题（对偶问题）具有无界解，则其对偶问题（原问题）无可行解。 性质四：强对偶性。如果原问题有最优解，则其对偶问题也一定有最优解。 性质五：互补松弛型。在线性规划问题的最优解中，如果对应某一约束条件的对偶变量值为零，则该约束条件取严格等式；反之如果约束条件取严格不等式，则其对应的对偶变量一定为零。 性质六：线性规划的原问题及其对偶问题之间存在一对互补的基解，其中原问题的松弛变量对应对偶问题的变量，对偶问题的剩余变量对应原问题的变量；这些互相对应的变量如果在一个问题的解中是基变量，则在另一问题的解中是非基变量；将这对互补的基解分别代入原问题和对偶问题的目标函数有z=w. B.对偶单纯形法（参考书p64页） 设某标准形式的线性规划问题，对偶单纯形表中必须有c j -z j ≤0（j=1...n),但b i (i=1...m)的值不一定为正，当对i=1...m ，都有b i ≥0时，表中原问题和对偶问题均为最优解，否则通过变换一个基变量，找出原问题的一个目标函数值较小的相邻的基解。 3）为什么要引入对偶单纯形法 从理论上说原始单纯形法可以解决一切线性规划问题，然而实际问题中，由于考虑问题的角度不同，变量设置的不同，便产生了原问题及其对偶问题，对偶问题是原问题从另外一个角度考虑的结果。用对偶单纯形法求解线性规划问题时，当约束条件为“≥”时，不必引入人工变量，使计算简化。 例如，有一线性规划问题： min ω =12 y 1 +16y 2 +15 y 3 约束条件 ?? ?? ???≥=≥+≥+0)3,2,1(3522 423121 i y y y y y i
流体力学 大作业
一.选择题 1.牛顿内摩擦定律适用于（）。 A．任何流体B．牛顿流体C．非牛顿流体 2.液体不具有的性质是（）。 A．易流动性B．压缩性C．抗拉性D．粘滞性 3连续介质假定认为流体（）连续。 A．在宏观上B．在微观上C．分子间D．原子间 4.在国际单位制中流体力学基本量纲不包括（）。 A．时间B．质量C．长度D．力． 5.在静水中取一六面体，作用在该六面体上的力有（） A．切向力、正压力B．正压力C．正压力、重力D．正压力、切向力、重力 6．下述哪些力属于质量力( ) A．惯性力B．粘性力C．弹性力D．表面张力E．重力 7．某点存在真空时，（）（） A．该点的绝对压强为正值B．该点的相对压强为正值c.该点的绝对压强为负值D．该点的相对压强为负值 8.流体静压强的（）。 A．方向与受压面有关B．大小与受压面积有关B．大小与受压面方位无关 9.流体静压强的全微分式为（）。 A．B．C． 10.压强单位为时，采用了哪种表示法（）。 A．应力单位B．大气压倍数C．液柱高度 11.密封容器内液面压强小于大气压强，其任一点的测压管液面（）。A．高于容器内液面B．低于容器内液面C．等于容器内液面 12.流体运动的连续性方程是根据（）原理导出的。 A.动量守恒 B. 质量守恒 C.能量守恒 D. 力的平衡 13. 流线和迹线重合的条件为（）。
A.恒定流 B.非恒定流 C.非恒定均匀流 14.总流伯努利方程适用于（）。 A.恒定流 B.非恒定流 C.可压缩流体 15. 总水头线与测压管水头线的基本规律是：（）、（） A.总水头线总是沿程下降的。 B.总水头线总是在测压管水头线的上方。 C.测压管水头线沿程可升可降。 D.测压管水头线总是沿程下降的。 16 管道中液体的雷诺数与（）无关。 A. 温度 B. 管径 C. 流速 D. 管长 17.. 某圆管直径d=30mm，其中液体平均流速为20cm/s。液体粘滞系数为0.0114cm3/s，则此管中液体流态为（）。 A. 层流 B. 层流向紊流过渡 C.紊流 18.等直径圆管中紊流的过流断面流速分布是（）A呈抛物线分布B. 呈对数线分布 C.呈椭圆曲线分布 D. 呈双曲线分布 19.等直径圆管中的层流，其过流断面平均流速是圆管中最大流速的（） A 1.0倍B.1/3倍C. 1/4倍D. 1/2倍 20.圆管中的层流的沿程损失与管中平均流速的（）成正比. A. 一次方 B. 二次方 C. 三次方 D. 四次方 21..圆管的水力半径是( ) A. d/2 B. d/3 C. d/4 D. d/5. 22谢才公式中谢才系数的单位是（）A. 无量纲B. C. D. . 23. 判断层流和紊流的临界雷诺数是（） A.上临界雷诺数 B.下临界雷诺数 C.上下临界雷诺数代数平均 D.上下临界雷诺数几何平均 24.. 对于管道无压流，当充满度分别为（）时，其流量和速度分别达到最大。 A. 0.5, 0.5 B. 0.95, 0.81 C. 0.81, 081 D. 1.0, 1.0 25.对于a, b, c三种水面线，下列哪些说法是错误（）（） A．所有a、c型曲线都是壅水曲线，即，水深沿程增大。B．所有b型曲线都是壅水曲线，即，水深沿程增大。C．所有a、c型曲线都是降水曲线，即，水深沿程减小。C．所有b型曲线都是降水曲线，即，水深沿程减
哈工大传热学作业答案
一维非稳态导热计算 4-15、一直径为1cm,长4cm 的钢制圆柱形肋片，初始温度为25℃，其后，肋基温度突然升高到200℃，同时温度为25℃的气流横向掠过该肋片，肋端及两侧的表面传热系数均为 100。试将该肋片等分成两段（见附图），并用有 限差分法显式格式计算从开始加热时刻起相邻4个时刻上的温度分布（以稳定性条件所允许的时间间隔计算依据）。已知＝43W/(m.K)，。（提示：节点4的离散方程可按端面的对流散热与从节点3到节点4的导热相平衡这一条件列出）。 解：三个节点的离散方程为： 节点2： 节点3： 节点4： 。 以上三式可化简为： 稳定性要求，即 。 ，代入得： ， 如取此值为计算步长，则： ，。 于是以上三式化成为： )./(2 K m W λs m a /10333.12 5 -?=()()12223212222/2444k k k k k k k f t t t t t t d d d d x h t t c x x x πππλλπρτ+????????---++?-=?? ? ? ? ???????????? ()()12224323333/2444k k k k k k k f t t t t t t d d d d x h t t c x x x πππλλπρτ+????????---++?-=?? ? ? ? ???????????? () 22344/244k k k f t t d d h t t x ππλ????-=- ? ?????? 12132222 43421k k f a a h a h t t t t t x x cd x cd τττττρρ+????????????? =+++-- ? ? ? ????????????13243222 43421k k f a a h a h t t t t t x x cd x cd τττττρρ+????????????? =+++-- ? ? ? ??????????? ?()4322k k f xh t t xht λλ+?=+?2 3410a h x cd ττ ρ??- -≥?2341/a h x cd τρ???≤+ ????5 54332.25810 1.33310c a λρ-===??5253 1.33310410011/8.898770.020.013 2.258100.0999750.0124s τ-??????≤+== ???+??5221.333108.898770.29660.02a x τ-???==?5441008.898770.110332.258100.01h cd τρ???==??1132 20.29660.29660.1103k k f t t t t +?++=12430.29660.296620.1103k k k f t t t t ++?+=34 0.97730.0227k k f t t t +=
哈工大2015年概率统计试题及答案
2015年哈工大概率统计试题 一、填空题（每小题3分，共5小题，满分15分） 1．设()()0.7P A P B +=，且,A B 只发生一个的概率为0.5，则,A B 都发生的概率为 ________________ ． 2．设随机变量X 的概率密度为???<≥=0 ,00e )(-x x x f x X ，，则随机变量X Y e =的概率密度为 ()Y f y = ______________ _ _ ． 3．设随机变量, X Y 的相关系数为0.5，220,2EX EY EX EY ====，则 2()E X Y +=． 4．生产一个零件所需时间2(,)X N μσ ，观察25个零件的生产时间得 5.5x =秒，样本 标准差 1.73s =秒，则μ的置信度为0.95的置信区间为________________ __． 5．设随机变量, X Y 相互独立，且均服从区间[]0,3上的均匀分布，则 {max(,)1}P X Y ≤=______ ． 注：可选用的部分数值：0.050.0250.025(24) 1.7109, (24) 2.0639, (25) 2.0595,t t t === .95.0645.1975.096.1=Φ=Φ）（，）（ 二、选择题（每小题3分，共5小题，满分15分） 1．设()01,P B <<(|)(|)1P A B P A B +=，则 （A ）,A B 互不相容.（B ）,A B 互为对立事件. （C ）,A B 相互独立.（D ）,A B 不独立.【】 2．下列函数可作为随机变量的分布函数的是 （A ）()2 1,1F x x x =-∞<<+∞+．（B ）, 0() 1 0, 0 x x F x x x ?≥? =+??
运筹学课程课件
中原工 Zhongyuan University of Technology Course Syllabus December 5 Note: Each student should print out this course syllabus and bring it to each class session. COURSE TITLE: Principles of Marketing: A Global Perspective INSTRUCTOR: William Teng, Ph.D., CFA CONTACT: wyteng@https://www.wendangku.net/doc/fb14859350.html, FACULTY BIO: Dr. William Teng received his Ph.D. in Economics and Finance from the University of Memphis in Tennessee. Dr. Teng also holds the designations of Chartered Financial Analyst (CFA). Dr. Teng has more than fifteen years of teaching experience at both the undergraduate and the graduate levels. He has published research papers in the International Journal of Service Science, Information Technology Journal, and Economics System. TEXTBOOK: Keegan, Warren, Global Marketing, 8th edition, Prentice-Hall, 4 LEARNING OUTCOMES: Upon successful course completion, students should be able to: 1. Identify the principles of marketing and explain the impact these principles have on the global economic, social/cultural, legal/political, and regulatory environment. 2. Identify regional economic markets and explain how to qualify and quantify potential opportunities using research, segmentation, and targeting techniques. 3. Explain how marketing ‘mix’ decisions – product, price, physical distribution, promotion – impact a global marketing strategy. 4. Explain the strengths and weaknesses of a company’s global marketing plan.
哈工大工程流体力学(二)试题
1.沿程阻力， 2.时间平均压强， 3.水力短管，5.翼弦 6.点汇， 7.旋涡强度， 8.速度势函数， 9.水力粗糙管，10.紊流 1.局部阻力， 2.时间平均流速， 3.水力长管，，5.翼弦 1．6.点源，7.涡线，8.流函数，9.水力光滑管，10.层流 2．水击现象、边界层 3．入口起始段、攻角、空气动力翼弦 1.简述边界层的特点 2.何谓述叶栅理论中的正问题和反问题 二、简答题（10分） 1. 在机翼理论中，如何利用保角变换法解决机翼绕流问题的 2.试推求有压管路产生水击时压强最大升高值的计算公式， 并说明减小水击的措施。（10分） 二、简答题 1．试分析流体流经弯管时局部阻力产生的具体原因是什么？（8分） 2．结合流体对圆柱体的有环量绕流，分析升力是如何产生的？（7分） 3．简述粘性流体绕物体流动时压差阻力产生的原因。 4．简述水击现象的物理过程，并说明减少水击现象的措施。 5．简述曲面边界层的分离现象 三、推求边界层的动量积分关系式（15分） 四、推求边界层的微分方程（普朗特边界层方程）
四、试推导说明圆柱外伸管嘴出流流量大于同直径薄壁小孔口的出流流量（10分） 三．推导理想流体平面有势流动中偶极流的速度势函数和流函数。（15分） 说明速度势函数的存在条件，并证明速度势函数的特性 说明流函数的存在条件，并证明流函数的特性 四．流体在长为l 的水平放置的等直径圆管中作定常流动，若已知沿程损失因数为λ，管壁切应力为τ，断面平均流速为V ； 试证明：28 V λ τρ= 。 （15分） 试推导二元旋涡的速度和压强分布 试证明旋涡理论中的斯托克斯定理 试证明速度环量保持不变的汤姆逊定理 三、推导、证明题 1.试推导圆管层流流动的速度分布规律，并求： （1）断面平均流速 （2）动能修正因数 (15分) 五、用突然扩大使管道的平均流速从1V 减到2V ，如图所示，如果 cm d 51=及1V 一定，试求使测压管液柱差h 成为最大值的2V 及2d 为若 干？并求m ax h 是多少？（10分）
概率论课程期末论文大作业
《概率论与数理统计》论文题目：正态分布及其应用 学院：航天学院 专业：空间科学与技术 姓名：黄海京 学号：1131850108
正态分布及其应用 摘要：正态分布（normal distribution），是一个在数学、物理及工程等领域都非常重要的概率分布，在统计学的许多方面有着重大的影响力。正态分布有极其广泛的实际背景, 例如测量误差, 人的生理特征尺寸如身高、体重等 ,正常情况下生产的产品尺寸:直径、长度、重量高度,炮弹的弹落点的分布等, 都服从或近似服从正态分布，以及确定医学参考值范围，药品规格，用量等。可以说,正态分布是自然界和社会现象中最为常见的一种分布, 一个变量如果受到大量微小的、独立的随机因素的影响, 那么这个变量一般是一个正态随机变量。 关键词：正态分布， 一、正态分布的由来 正态分布（normal distribution）又名高斯分布（Gaussian distribution）。正态分布概念是由德国的数学家和天文学家Moivre于1733年受次提出的，但由于德国数学家Gauss率先将其应用于天文学家研究，故正态分布又叫高斯分布，高斯这项工作对后世的影响极大，他使正态分布同时有了“高斯分布”的名称，后世之所以多将最小二乘法的发明权归之于他，也是出于这一工作。 正态分布是一个在数学、物理及工程等领域都非常重要的概率分布，在统计学的许多方面有着重大的影响力。若随机变量X服从一个数学期望为μ、标准方差为σ2的高斯分布，记为：则其概率密度函数为正态分布的期望值μ决定了其位置，其标准差σ决定了分布的幅度。因其曲线呈钟形，因此人们又经常称之为钟形曲线。我们通常所说的标准正态分布是μ= 0,σ= 1的正态分布。 二、正态分布的特性 1. 正太分布的曲线特征 正态曲线呈钟型，两头低，中间高，左右对称，曲线与横轴间的面积总等于1。 （1）集中性：正态曲线的高峰位于正中央，即均数所在的位置。 （2）对称性：正态曲线以均数为中心，左右对称，曲线两端永远不与横轴相交。 （3）均匀变动性：正态曲线由均数所在处开始，分别向左右两侧逐渐均匀下降。
哈工大工程流体力学样本
《工程流体力学》综合复习资料 一、判断题 1、根据牛顿内摩擦定律, 当流体流动时, 流体内部内摩擦力大小与该处的流 速大小成正比。 2、一个接触液体的平面壁上形心处的水静压强正好等于整个受压壁面上所有 各点水静压强的平均值。 3、流体流动时, 只有当流速大小发生改变的情况下才有动量的变化。 4、在相同条件下, 管嘴出流流量系数大于孔口出流流量系数。 5、稳定( 定常) 流一定是缓变流动。 6、水击产生的根本原因是液体具有粘性。 7、长管是指运算过程中流速水头不能略去的流动管路。 8、所谓水力光滑管是指内壁面粗糙度很小的管道。 9、外径为D, 内径为d的环形过流有效断面, 其水力半径为 4d D- 。 10、凡是满管流流动, 任何断面上的压强均大于大气的压强。 二、填空题 1、某输水安装的文丘利管流量计, 当其汞-水压差计上读数cm h4 = ?, 经过的流量为s L/ 2, 分析当汞水压差计读数cm h9 = ?, 经过流量为L/s。 2、运动粘度与动力粘度的关系是v=u/p , 其国际单位是厘斯(mm2/s) 。 3、因次分析的基本原理是: 因次和谐的原理 ; 具体计算方法分为两 种。 4、断面平均流速V与实际流速u的区别是。 5、实际流体总流的伯诺利方程表示式为 , 其适用条件是。 6、泵的扬程H是指扬程, m。 7、稳定流的动量方程表示式为。
8、计算水头损失的公式为与。 9、牛顿内摩擦定律的表示式τ=μγ , 其适用范围是是指在温度不变 的条件下, 随着流速梯度的变化, μ值始终保持一常数。 10、压力中心是指作用在物体上的空气动力合力的作用点。 三、简答题 1、稳定流动与不稳定流动。---流体在管道内或在窑炉系统中流动时, 如果任 一截面上的流动状况(流速、压强、重度、成分等)都不随时间而改变, 这种流动就称为稳定流动; 反之, 流动各量随着时间而改变, 就称为不稳定流动。实际上流体(如气体, 重油等)在管道内或窑炉系统中流动时, 只要波动不太大, 都能够视为稳定流动。 2、 产生流动阻力的原因。---直管阻力: 流体流经直管段时, 由于克服流体的粘滞性及与管内壁间的磨擦所产生的阻力。有粘管壁, 其壁面的流动速度降为0. 局部阻力: 流体流经异形管或管件时, 由于流动发生骤然变化引起涡流所产生的能量损失。 3、串联管路的水力特性。---串联管路无中途分流和合流时, 流量相等, 阻力 叠加。串联管路总水头损失等于串联各管段的水头损失之和, 后一管段的流量等于前一管段流量减去前管段末端泄出的流量。 4、如何区分水力光滑管和水力粗糙管, 两者是否固定不变? ---在紊流中存在 层流底层, 当层流底层厚度δl＞5Δ时, 粗糙高度几乎全被层流底层淹没, 管壁对紊流区流体的影响很小, 这与流体在完全光滑的管道中流动类似, 这种情况的管子叫做水力光滑管。当层流底层厚度δl＜0.3Δ时, 管壁上几乎所有的凸峰都暴露在紊流中, 紊流去的流体质点与凸峰相互碰撞, 阻力增加, 此时的管子叫做水利粗糙管。 5、静压强的两个特性。---1.静压强的方向是垂直受压面, 并指向受压面。2. 任一点静压强的大小和受压面方向无关, 或者说任一点各方向的静压强均相等。
论传热学与本专业的联系
《传热学》课程大作业 题目：____ 论传热学与飞行器制造工程的关系___________ 姓名：____ _ __________ 学号：_____ ________ 授课教师：_______ ___ _ ___ 哈尔滨工业大学 2015年5月4日
论传热学与飞行器制造工程的关系 摘要：本文主要介绍传热学的研究内容，研究简史，飞行器制造工程专业的研究内容与方向。其中重点介绍传热学对飞行器制造工程专业相关研究的影响。 关键词：传热学；飞行器制造工程 一传热学研究简史及研究内容 传热学是研究物体内部或物体与物体之间由温度差引起热量传递过程的学科。飞行器及其推进系统的发展提出了大量的传热学问题。传热的基本方式有导热、对流传热和辐射传热。传热不仅是常见的自然现象，而且广泛存在于工程技术领域。提高锅炉的蒸汽产量，防止燃气轮机燃烧室过热、减小内燃机气缸和曲轴的热应力、确定换热器的传热面积和控制热加工时零件的变形等，都是典型的传热问题。 传热学作为学科形成于19世纪。在热对流方面，英国科学家牛顿于1701 年在估算烧红铁棒的温度时，提出了被后人称为牛顿冷却定律的数学表达式，不过它并没有揭示出对流换热的机理。 对流换热的真正发展是19世纪末叶以后的事情。1904年德国物理学家普朗特的边界层理论和1915年努塞尔的因次分析，为从理论和实验上正确理解和定量研究对流换热奠定了基础。1929年，施密特指出了传质与传热的类同之处。 二飞行器制造工程研究内容 飞行器制造工程专业以一般机械制造工程为基础，广泛吸收各种先进技术和科学理论的成果，针对飞行器的特点研究各种制造方法的机理和应用，探求制造过程的规律，合理利用资源，经济而高效率地制造先进优质飞行器的一门技术科学。 本专业主要以机械设计制造为基础，涉及机械工程、电机工程、电子技术、计算机技术、材料科学、管理工程、控制工程和系统工程等许多科学技术领域。各种新结构、新元件、新材料、新工艺、新方法的应用，正在加速整个飞行器制造工程的发展。设计-结构-材料-工艺技术的最佳配合将是飞行器制造工程中的一个新趋向。 三传热学对飞行器制造工程的影响 本专业涉及到各类制造加工技术，普通加工中切削加工占大多数，加工工件
哈工大概率论与数理统计课后习题答案四
习 题 四 1．一个袋子中装有四个球，它们上面分别标有数字1,2,2,3，今从袋中任取一球后不放回，再从袋中任取一球，以,X Y 分别表示第一次，第二次取出的球上的标号，求(,)X Y 的 分布列. 解 (,)X Y 的分布列为 其中 (1,1)(1)(1|1)0P X Y P X P Y X ======= (1,2)(1)(2|1)P X Y P X P Y X ====== 121436 =?= 余者类推。 2．将一枚硬币连掷三次，以 X 表示在三次中出现正面的次数，以Y 表示三次中出现正 面次数与出现反面次数之差的绝对值，试写出(,)X Y 的分布列及边缘分布列。 解 一枚硬币连掷三次相当于三重贝努里试验，故1 ~(3,).2 X B 331 ()(),0,1,2,32 k P X k C k ===，于是(,)X Y 的分布列和边缘分布为 其中 (0,1)(0)(1|0)0P X Y P X P Y X =======，
13 313(1,1)(1)(1|1)()128 P X Y P X P Y X C =======?=， 余者类推。 3．设(,)X Y 的概率密度为 1 (6),02,24, (,)80,.x y x y f x y ?--<<<
哈尔滨工程大学传热学大作业数值计算matlab程序内容
传热学作业数值计算
数值计算matlab程序内容：>> tw1=10; % 赋初值tw2=20; c=1.5; p2=20; p1=c*p2; L2=40; L1=c*L2; deltaX=L2/p2; a=p2+1; b=p1+1; ti=ones(a,b)*5; m1=ones(a,b); m1(a,2:b-1)=zeros(1,b-2); m1(2:a,1)=zeros(a-1,1); m1(2:a,b)=zeros(a-1,1); m1(1,:)=ones(1,b)*2; k=0; max1=1.0; tn=ti; while(max1>1e-6) max1=0; k=k+1; for i=1:1:a for j=1:1:b
m=m1(i,j); n=ti(i,j); switch m case 0 tn(i,j)=tw1; case 1 tn(i,j)=0.25*(tn(i,j+1)+tn(i,j-1)+tn(i+1,j)+tn(i-1,j)); case 2 tn(i,j)=tw1+tw2*sin(pi*(j-1)/(b-1)); end er=abs(tn(i,j)-n); if er>max1 max1=er; end end end ti=tn; end k ti max1 t2=ones(a,b); %求解析温度场 for i=a:-1:1 for j=1:1:b y=deltaX*(a-i); x=deltaX*(j-1); t2(i,j)=tw1+tw2*sin(pi*x/L1)*(sinh(pi*y/L1))/(sinh(pi*L2/L1)); end end t2 迭代次数k =706 数值解温度场ti