西安交通大学——温度场数值模拟(matlab)

温度场模拟matlab代码：

clear,clc,clf

L1=8;L2=8;N=9;M=9;% 边长为8cm的正方形划分为8*8的格子

T0=500;Tw=100; % 初始和稳态温度

a=0.05; % 导温系数

tmax=600;dt=0.2; % 时间限10min和时间步长0.2s

dx=L1/(M-1);dy=L2/(N-1);

M1=a*dt/(dx^2);M2=a*dt/(dy^2);

T=T0*ones(M,N);

T1=T0*ones(M,N);

t=0;l=0;k=0;

Tc=zeros(1,600);% 中心点温度，每一秒采集一个点

for i=1:9

for j=1:9

if(i==1|i==9|j==1|j==9)

T(i,j)=Tw;% 边界点温度为100℃

else

T(i,j)=T0;

end

end

end

if(2*M1+2*M2<=1) % 判断是否满足稳定性条件

while(t

t=t+dt;

k=k+1;

for i=2:8

for j=2:8

T1(i,j)=M1*(T(i-1,j)+T(i+1,j))+M2*(T(i,j-1)+T(i,j+1))+(1-2*M1-2*M2)*T(i,j);

end

end

for i=2:8

for j=2:8

T(i,j)=T1(i,j);

end

end

if(k==5)

l=l+1;

Tc(l)=T(5,5);

k=0;

end

end

i=1:9;j=1:9;

[x,y]=meshgrid(i); figure(1);

subplot(1,2,1);

mesh(x,y,T(i,j))% 画出10min 后的温度场 axis tight;

xlabel('x','FontSize',14);ylabel('y','FontSize',14);zlabel('T/℃','FontSize',14) title('1min 后二维温度场模拟图','FontSize',18) subplot(1,2,2);

[C,H]=contour(x,y,T(i,j)); clabel(C,H);axis square;

xlabel('x','FontSize',14);ylabel('y','FontSize',14); title('1min 后模拟等温线图','FontSize',18) figure(2); xx=1:600;

plot(xx,Tc,'k-','linewidth',2)

xlabel('时间/s','FontSize',14);ylabel('温度/℃','FontSize',14);title('中心点的冷却曲线','FontSize',18)

else disp('Error!') % 如果不满足稳定性条件，显示“Error ！” end

实验结果：

时间/s

温度/℃

中心点的冷却曲线

x

1min

后二维温度场模拟图

T /℃

x

y

1min 后模拟等温线图

x

5min 后二维温度场模拟图

T /℃

x

y

5min 后模拟等温线图

x

10min

后二维温度场模拟图

T /℃

x

y

10min 后模拟等温线图

x

10min 后二维温度场模拟图（不满足稳定性条件）

y

T /℃

21

时间/s

温度/℃

中心点的冷却曲线（不满足稳定性条件）

现代信号处理Matlab仿真——例611

例6.11 利用卡尔曼滤波估计一个未知常数 题目： 设已知一个未知常数x 的噪声观测集合，已知噪声v(n)的均值为零， 方差为 ，v(n)与x 不相关，试用卡尔曼滤波估计该常数 题目分析： 回忆Kalman 递推估计公式 由于已知x 为一常数，即不随时间n 变化，因此可以得到： 状态方程： x(n)=x(n-1) 观测方程： y(n)=x(n)+v(n) 得到A(n)=1,C(n)=1, , 将A(n)=1,代入迭代公式 得到：P(n|n-1)=P(n-1|n-1) 用P(n-1)来表示P(n|n-1)和P(n-1|n-1),这是卡尔曼增益表达式变为 从而 2v σ1??(|1)(1)(1|1)(|1)(1)(1|1)(1)()()(|1)()[()(|1)()()]???(|)(|1)()[()()(|1)](|)[()()](|1)H w H H v x n n A n x n n P n n A n P n n A n Q n K n P n n C n C n P n n C n Q n x n n x n n K n y n C n x n n P n n I K n C n P n n --=----=----+=--+=-+--=--2()v v Q n σ=()0w Q n =(|1)(1)(1|1)(1)()H w P n n A n P n n A n Q n -=----+21 ()(|1)[(|1)]v K n P n n P n n σ-=--+22(1)()[1()](1)(1)v v P n P n K n P n P n σσ-=--=-+

维导热物体温度场的数值模拟

传热大作业 二维导热物体温度场的数值模拟（等温边界条件） 姓名： 班级： 学号：

墙角稳态导热数值模拟（等温条件） 一、物理问题 有一个用砖砌成的长方形截面的冷空气空道，其截面尺寸如下图所示，假设在垂直于纸面方向上冷空气及砖墙的温度变化很小，可以近似地予以忽略。在下列两种情况下试计算： （1）砖墙横截面上的温度分布； （2）垂直于纸面方向的每米长度上通过砖墙的导热量。外矩形长为，宽为；内矩形长为，宽为。 第一种情况：内外壁分别均匀地维持在0℃及30℃； 第二种情况：内外表面均为第三类边界条件，且已知： 外壁：30℃，h1=10W/m2·℃, 内壁：10℃，h2= 4 W/m2·℃ 砖墙的导热系数λ= W/m·℃ 由于对称性，仅研究1/4部分即可。 二、数学描写 对于二维稳态导热问题，描写物体温度分布的微分方程为拉普拉斯方程

02222=??+??y t x t 这是描写实验情景的控制方程。 三、方程离散 用一系列与坐标轴平行的网格线把求解区域划分成许多子区域，以网格线的交点作为确定温度值的空间位置，即节点。每一个节点都可以看成是以它为中心的一个小区域的代表。由于对称性，仅研究1/4部分即可。依照实验时得点划分网格： 建立节点物理量的代数方程 对于内部节点，由?x=?y ，有 )(411,1,,1,1,-+-++++=n m n m n m n m n m t t t t t 由于本实验为恒壁温，不涉及对流，故内角点，边界点代数方程与该式相同。

设立迭代初场，求解代数方程组。图中，除边界上各节点温度为已知且不变外，其余各节点均需建立类似3中的离散方程，构成一个封闭的代数方程组。以C t 000 为场的初始温度，代入方程组迭代，直至相邻两次内外传热值之差小于，认为已达到迭代收敛。 四、编程及结果 1) 源程序 #include <> #include <> int main() { int k=0,n=0; double t[16][12]={0},s[16][12]={0}; double epsilon=; double lambda=,error=0; double daore_in=0,daore_out=0,daore=0; FILE *fp; fp=fopen("data3","w"); for (int i=0;i<=15;i++) for (int j=0;j<=11;j++) { if ((i==0) || (j==0)) s[i][j]=30; if (i==5) if (j>=5 && j<=11) s[i][j]=0; if (j==5) if (i>=5 && i<=15) s[i][j]=0; } for (int i=0;i<=15;i++)

Matlab仿真实例-卫星轨迹

卫星轨迹 一．问题提出 设卫星在空中运行的运动方程为： 其中是k 重力系数(k=401408km3/s)。卫星轨道采用极坐标表示，通过仿真，研究发射速度对卫星轨道的影响。实验将作出卫星在地球表面（r=6400KM ，θ=0）分别以v=8KM/s,v=10KM/s,v=12KM/s 发射时，卫星绕地球运行的轨迹。 二．问题分析 1．卫星运动方程一个二阶微分方程组，应用Matlab 的常微分方程求解命令ode45求解时，首先需要将二阶微分方程组转换成一阶微分方程组。若设，则有： 2．建立极坐标如上图所示，初值分别为：卫星径向初始位置，即地球半径：y(1,1)=6400；卫星初始角度位置：y(2,1)=0；卫星初始径向线速度：y(3,1)=0；卫星初始周向角速度：y(4,1)=v/6400。 3．将上述一阶微分方程及其初值带入常微分方程求解命令ode45求解，可得到一定时间间隔的卫星的径向坐标值y(1)向量；周向角度坐标值y(2)向量；径向线速度y(3)向量；周向角速度y(4)向量。 4．通过以上步骤所求得的是极坐标下的解，若需要在直角坐标系下绘制卫星的运动轨迹，还需要进行坐标变换，将径向坐标值y(1)向量；周向角度坐标值y(2)向量通过以下方程转换为直角坐标下的横纵坐标值X,Y 。 5．卫星发射速度速度的不同将导致卫星的运动轨迹不同，实验将绘制卫星分别以v=8KM/s ，v=10KM/s ，v=12KM/s 的初速度发射的运动轨迹。 三．Matlab 程序及注释 1．主程序 v=input('请输入卫星发射速度单位Km/s ：\nv=');%卫星发射速度输入。 axis([-264007000-1000042400]);%定制图形输出坐标范围。 %为了直观表达卫星轨迹，以下语句将绘制三维地球。 [x1,y1,z1]=sphere(15);%绘制单位球。 x1=x1*6400;y1=y1*6400;???????-=+-=dt d dt dr r dt d dt d r r k dt r d θ θθ2)(2 22222θ==)2(,)1(y r y ?????????????**-=**+*-===)1(/)4()3(2)4()4()4()1()1()1()3()4()2() 3()1(y y y dt dy y y y y y k dt dy y dt dy y dt dy ???*=*=)] 2(sin[)1(Y )]2(cos[)1(X y y y y

matlab绘制温度场

通过在室内的某些位置布置适当的节点，采集回来室内的温湿度以及空气质量等实际参数。首先对室内空间建模，用一个无限细化的三维矩阵来模拟出室内的温度分布情况，针对采集回来的数据，采用插值法和适当次数的拟合函数的拟合，得出三维矩阵的实际值的分布，最后结合matlab软件绘制出计算出的温度场的三维图像。 一．数据的采集与处理 因为影响人的舒适感的温度层只是室内的某一高度范围内的温度，而温度传感器虽然是布置在一个平面内，但是采用插值法和拟合函数法是可以大致再现出影响人的舒适感的温度层的温度变化的。同时，在构建出的三维模型中，用第三维表示传感器层面的温度。 在传感器层面，传感器分布矩阵如下： X=【7.5 36.5 65.5】（模型内单位为cm） Y=【5.5 32.5 59.5】 Z=【z1 z2 z3； z4 z5 z6； z7 z8 z9；】（传感器采集到的实时参数） 采用meshgrid（xi，yi，zi，…）产生网格矩阵； 首先按照人的最小温度分辨值，将室内的分布矩阵按照同样的比例细化，均分，使取值点在坐标一定程度上也是接近于连续变化的，从而才能最大程度上使处理数据得来的分布值按最小分辨值连续变化！ 根据人体散热量计算公式：C=hc（tb-Ta） 其中hc为对流交换系数； 结合Gagge教授提出的TSENS热感觉指标可以计算出不同环境下人的对环境温度变化时人体温度感知分辨率，作为插值法的一个参考量，能使绘制出的温度场更加的符合人体的温度变化模式。 例如按照10cm的均差产生网格矩阵（实际上人对温度的分辨率是远远10cm大于这个值的，但是那样产生的网格矩阵也是异常庞大的，例如以0.5cm为例，那么就可以获得116*108=12528个元素，为方便说明现已10cm为例）： [xi yi]=meshgrid(7.5:10:65.5,5.5:10:59.5) xi = 7.5000 17.5000 27.5000 37.5000 47.5000 57.5000 7.5000 17.5000 27.5000 37.5000 47.5000 57.5000 7.5000 17.5000 27.5000 37.5000 47.5000 57.5000 7.5000 17.5000 27.5000 37.5000 47.5000 57.5000 7.5000 17.5000 27.5000 37.5000 47.5000 57.5000 7.5000 17.5000 27.5000 37.5000 47.5000 57.5000

西安交通大学——温度场数值模拟(matlab)

温度场模拟matlab代码： clear,clc,clf L1=8;L2=8;N=9;M=9;% 边长为8cm的正方形划分为8*8的格子 T0=500;Tw=100; % 初始和稳态温度 a=0.05; % 导温系数 tmax=600;dt=0.2; % 时间限10min和时间步长0.2s dx=L1/(M-1);dy=L2/(N-1); M1=a*dt/(dx^2);M2=a*dt/(dy^2); T=T0*ones(M,N); T1=T0*ones(M,N); t=0;l=0;k=0; Tc=zeros(1,600);% 中心点温度，每一秒采集一个点 for i=1:9 for j=1:9 if(i==1|i==9|j==1|j==9) T(i,j)=Tw;% 边界点温度为100℃ else T(i,j)=T0; end end end if(2*M1+2*M2<=1) % 判断是否满足稳定性条件 while(t

end i=1:9;j=1:9; [x,y]=meshgrid(i); figure(1); subplot(1,2,1); mesh(x,y,T(i,j))% 画出10min 后的温度场 axis tight; xlabel('x','FontSize',14);ylabel('y','FontSize',14);zlabel('T/℃','FontSize',14) title('1min 后二维温度场模拟图','FontSize',18) subplot(1,2,2); [C,H]=contour(x,y,T(i,j)); clabel(C,H);axis square; xlabel('x','FontSize',14);ylabel('y','FontSize',14); title('1min 后模拟等温线图','FontSize',18) figure(2); xx=1:600; plot(xx,Tc,'k-','linewidth',2) xlabel('时间/s','FontSize',14);ylabel('温度/℃','FontSize',14);title('中心点的冷却曲线','FontSize',18) else disp('Error!') % 如果不满足稳定性条件，显示“Error ！” end 实验结果： 时间/s 温度/℃ 中心点的冷却曲线

MATLAB实现通信系统仿真实例

补充内容：模拟调制系统的MATLAB 仿真 1.抽样定理 为了用实验的手段对连续信号分析，需要先对信号进行抽样（时间上的离散化），把连续数据转变为离散数据分析。抽样（时间离散化）是模拟信号数字化的第一步。 Nyquist 抽样定律：要无失真地恢复出抽样前的信号，要求抽样频率要大于等于两倍基带信号带宽。 抽样定理建立了模拟信号和离散信号之间的关系，在Matlab 中对模拟信号的实验仿真都是通过先抽样，转变成离散信号，然后用该离散信号近似替代原来的模拟信号进行分析的。 【例1】用图形表示DSB 调制波形)4cos()2cos(t t y ππ= 及其包络线。 clf %%计算抽样时间间隔 fh=1;%%调制信号带宽(Hz) fs=100*fh;%%一般选取的抽样频率要远大于基带信号频率，即抽样时间间隔要尽可能短。 ts=1/fs; %%根据抽样时间间隔进行抽样，并计算出信号和包络 t=(0:ts:pi/2)';%抽样时间间隔要足够小，要满足抽样定理。 envelop=cos(2*pi*t);%%DSB 信号包络 y=cos(2*pi*t).*cos(4*pi*t);%已调信号 %画出已调信号包络线 plot(t,envelop,'r:','LineWidth',3); hold on plot(t,-envelop,'r:','LineWidth',3); %画出已调信号波形 plot(t,y,'b','LineWidth',3); axis([0,pi/2,-1,1])% hold off% xlabel('t'); %写出图例 【例2】用图形表示DSB 调制波形)6cos()2cos(t t y ππ= 及其包络线。 clf %%计算抽样时间间隔 fh=1;%%调制信号带宽(Hz) fs=100*fh;%抽样时间间隔要足够小，要满足抽样定理。 ts=1/fs; %%根据抽样时间间隔进行抽样

二维导热物体温度场的数值模拟

金属凝固过程计算机模拟题目：二维导热物体温度场的数值模拟 Solidworks十字接头的传热分析 作者：张杰 学号：S2******* 学院：北京有色金属研究总院 专业：材料科学与工程 成绩： 2015 年12 月

二维导热物体温度场的数值模拟 图1 二维均质物体的网格划分 用有限差分法模拟二维导热物体的温度场，首先将二维物体划分为如图1所示的网格，x ?与y ?可以是不变的常量,即等步长,也可以是变量(即在区域内的不同处是不同的),即变步长?如果区域内各点处的温度梯度相差很大,则在温度变化剧烈处,网格布得密些,在温度变化不剧烈处,网格布得疏些?至于网格多少,步长取多少为宜,要根据计算精度与计算工作量等因素而定? 在有限的区域内,将二维不稳定导热方程式应用于节点 ,)i j （可写成: ,2222 ,i j P P p i j T T T C x y ρλτ?????=+ ?????? ,1 , ,()i j P P P i j i j T T T οτττ+-???= +? ????? () , 1 , , 1 ,22 2()i j P P P P i j i j i j T T T T x x x ο+--+??? =+? ????? () , ,1 , ,122 2()i j P P P P i j i j i j T T T T y y y ο+--+???=+? ?????τ?、x ?、y ? 当τ?、x ?、y ?较小时，忽略()οτ?、2()x ο?、2 ()y ο?项。当x y ?=?时， 即x 、y 方向网格划分步长相等?最后得到节点 ,)i j （的差分方程: ()1 , ,0 1 , 1 , ,1 ,1 ,4P P P P P P P i j i j i j i j i j i j i j T T F T T T T T ++-+-=++++- 式中:() 02 p F C x λτ ρ?= ??

基于生死单元的激光熔覆温度场数值模拟

基于生死单元的激光熔覆温度场数值模拟 基于生死单元的激光熔覆温度场数值模拟 摘要:计算了不同激光功率条件下粉末颗粒到达基底前的温升，并以粉末颗粒到达基底前的温度为初始条件。用生死单元法研究了单通道和多通道激光熔覆温度场。利用熔池的大小和形态，验证了模型的可靠性。结果表明，粉末颗粒的温升与激光功率呈线性关系。单个包层的温度变化是锯齿状的。温升过程近似为直线，温降曲线近似为双曲线。在多通道熔覆过程中，温度场呈微椭圆形。节点上的热循环经过一个逐渐增加的峰值。峰值温度最终趋于稳定。0系列 激光熔覆根据送粉工艺不同可分为两种类型，即粉末预置法和同步送粉法。本发明具有易于自动控制、激光能量吸收率高、无内部气孔的优点。特别是对于覆层金属陶瓷，覆层的抗裂性可以显著提高，并且硬质陶瓷相可以均匀地分布在覆层中。有广阔的应用空间。国内学者利用ANSYS [1-4对激光熔覆过程的温度场和应力场进行了大量的研究工作。目前，利用ANSYS模拟激光熔覆温度场的研究没有考虑激光束与粉末的相互作用。事实上，激光束首先作用于粉末。除了损失的能量，部分激光束被包覆粉末吸收。另一部分通过粉末被基质吸收。除了直接吸收激光束能量，基质还吸收从粉末转移到基质的能量。因此，有必要在仿真前弄清激光能量的分布，使所建立的模型更接近实际，仿真结果更有说服力。本文将粉末在到达基体前吸收能量后的温升作为初始温度场加载到基体上。同时，利用有限元分析软件

ANSYS中的生死单元技术模拟了熔覆单元的生长过程。高斯体热源加载基体吸收的能量，模拟送粉激光熔覆的温度场分布。在此基础上，模拟了多道次激光熔覆的温度场，研究了多道次激光熔覆的温度场。当屏蔽激光时， 1粉末到达基体前的温度为 粉末。它还吸收部分激光能量，从而提高其温度。事实上，粒子直接吸收激光辐射能量并发射辐射能量，而不考虑等离子体的影响(能量密度低于105W /cm2)。在空气中，粉末颗粒也因空气对流而耗散能量，并且颗粒也相互加热。这些能量在总能量中的比例非常小。目前，关于粉体颗粒温升的模型很少。此外，有必要在模型[5]中建立假设条件。为了便于计算，模型中假设: (1)气体-粉末射流中粉末颗粒的体积分数很低，并且受到激光反射、折射、颗粒离子间相互加热和束屏蔽等的影响。可以忽略。(2)粉末颗粒是半径为rP的球体。由于粉末颗粒足够小，它们被认为是能量计算中的一个点。颗粒的导热性是无限的，即粉末颗粒的温度被认为是均匀的，并且在光接收表面和背光表面之间没有差异。(3)粉末颗粒仅吸收光接收表面上的能量，但是外部辐射发生在整个球体的表面上。(4)粉末不吸收来自基质的光反射。基于上述假设，粉末颗粒的温升可以根据颗粒的能量方程来计算。这个方程是一个非线性方程。利用Matlab软件，采用迭代法求解方程。当激光功率P=2 kW时，方程的解在1500 ~ 1600k范围内，因此初始值被设置为t = 1500k，并且通过迭代发现方程的一个实根是t = 1570k。改变激光功率，获得了当

MATLAB 画等温线

测量到不同坐标点的高度值,如何用matlab画三维图 附上部分数据： A=[-210.6627 -33391.1192 5.0273 -221.3052 -33387.7415 4.5969 -210.9391 -33393.0068 5.5647 -221.8901 -33390.7396 5.0077 -211.384 -33394.7093 5.6505 -222.6117 -33392.778 5.0554 -212.7074 -33397.5459 5.7381 -225.8973 -33397.5869 5.5587]; 解：代码在matlab2009a版以上均可运行。 A=[-210.6627 -33391.1192 5.0273 -221.3052 -33387.7415 4.5969 -210.9391 -33393.0068 5.5647 -221.8901 -33390.7396 5.0077 -211.384 -33394.7093 5.6505 -222.6117 -33392.778 5.0554 -212.7074 -33397.5459 5.7381 -225.8973 -33397.5869 5.5587]; xData = A(:,1); yData = A(:,2); zData = A(:,3); fitresult = fit( [xData, yData], zData, 'linearinterp'); figure( 'Name', '三维图' ); plot( fitresult, [xData, yData], zData ); xlabel( 'x' ); ylabel( 'y' ); zlabel( 'z' ); grid on view( -53, 50 );

激光焊接温度场数值模拟讲解

第２４卷第２期 ２ＯＯ 焊接学报 ｖ０１．２４Ａｐｒｉｌ Ｎｏ．２２０Ｏ３ ３年４月ＴＲＡＮＳＡＣ’ｎ０ＮＳ０ＦＴＨＥＣＨＩＮＡ碍砸ＬＤＩＮＧＩＮＳｎＴＵｒｎＯＮ 激光焊接温度场数值模拟 薛忠明，顾 兰， 张彦华 （北京航空航天大学机械工程及自动化学院。北京１０００８３） 摘要：深入分析了激光焊接小孔传热模型的特点，在此基础上选取合适的热源形式，研究了移动线热源和高斯分布热源作用下，准稳态与瞬态激光焊接温度场。利用ＭＡＴ－ＬＡＢ软件及ＡＮｓＹＳ有限元分析程序对激光焊接温度场分别进行了计算及模拟，并且将两种分析结果进行了比较。最后还将有限元的模拟值与实测值进行了对比分析，进一步验证了小孔模型与高斯热源在激光焊接温度场模拟中的适用性。关键词：激光焊接；温度场；有限元；ＡＮｓＹｓ 中围分类号：１嘶６ Ｏ 文献标识码：Ａ文章编号：０２５３—３６０ｘ（２００３）０１—７９—０４薛忠明 序言 实测值进行了对比分析，验证了小孔模型与高斯热源在激光焊接温度场模拟中的适用性（板厚≤４ｍｍ）。 激光焊接是利用高能量密度的激光束作为热源的一种高效精密的焊接方法。激光焊接具有高能量密度、可聚焦、深穿透、高效率、高精度、适应性强等优点，广泛应用于航空航天、汽车、微电子、轻工业、医疗及核工业等要求高精度和高质量的焊接领域。 １ 激光焊接中的小孔传热模型 当激光功率密度达到１０６Ｗ，／ｃｍ２时，激光能量 由于激光焊接是一巾陕速而不均匀的热循环过 程，焊缝附近出现很大的温度梯度，因此在焊后的结构中也会出现不同程度的残余应力和变形，这些都成为影响焊接结构质量和使用性能的重要因素。准确地认

识焊接热过程，对焊接结构力学分析、显微组织分析以及最终的焊接质量控制具有重要意义。 ２０世纪７０年代以来，国外很多学者对激光焊接机理进行了深入的研究，提出了蒸汽小孔模型。考虑熔池形状以及熔池中金属的流动和热流分布，考虑电子密度、离子化程度、等离子体对入射激光的吸收系数和激光焊接工艺参数对熔深的影响，建立了不同的能量吸收模型”。。这些研究偏向于应用物理和量子力学的研究领域，在实际工程分析中存在一定的局限性。在国内，有关激光焊接机理以及激光焊接温度场与力学场的数值模拟方面的研究正在引起重视。 作者深入分析了激光焊接小孔传热模型，在此基础上选取合适的热源形式，研究了移动线热源和高斯分布热源作用下，准稳态与瞬态激光焊接温度场。利用ＭＡｌｌＡＢ软件及ＡＮＳＹｓ有限元程序对激光焊接温度场分别进行了计算及模拟，并且将两种分析结果进行了比较。最后还将有限元的模拟值与 收稿日期：２００２—０７—１２ 向工件输入的速率远大于传导、对流、辐射散热的速率，材料表面产生汽化而形成小孔，激光能量是通过小孔而进行转换和传递的。 激光焊接中熔池与小孔的几何特征如图ｌ所示。焊件表面被加热、熔化、蒸发，在蒸汽压力的作用下形成小孔，当小孔产生的蒸汽压力与熔池中液体金属的静应力达到平衡时，小孔是稳定存在的‘“。 固１Ｈｇ．１ 激光焊接熔池与小孔几何特征囤 Ｇ岫ｅｔｒ萱ｃｆ嘲ｕｒ嚣０ｆｍｏｌｔｅｎ ａｎｄｋｅｙｈ０Ｉｅｉｎ ｐｏｏｌ Ｊ∞ｅｒｗｅⅫ咂ｇ 激光焊接中，小孔与工件作相对运动，运动过程 中的动量扩散和热量扩散的相对程度由佩克莱特准 万方数据 ８０

MATLAB电路仿真实例

题14.14 图(a)所示电路，已知 V )2cos(15S t u =二端口网络阻抗参数矩阵 Ω?? ????=46j 6j 10Z 求ab 端戴维南等效电路并计算电压o u 。 u -+o u 图题14.14 (一)手动求解： 将网络N 用T 型电路等效，如图(b)所示 S U +-o U 等效阻抗 Ω=-+-?+ -=4.6j615j6j6)15(6j 6j 4i Z 开路电压 V 2j302 15j6j6105j6OC =?∠?+-+=U V 1482 18.3j46.42j3j4j4Z j4OC o ?∠=+?=?+=U U i

所以 )1482cos(18.3o ?+=t u V (二)Matlab 仿真： ⒈分析：本次仿真需输入各阻抗Zl 、Z1、Z2、Z3、Z4以及激励源Us 的参数值，仿真结果需输出开路电压Uoc 、等效阻抗Zi 以及电感两端电压U0的幅值和相位信息，并绘制Uoc ，U0的值随时间变化的波形曲线。其中各元件与原图的对应关系如下图所示： ⒉编辑M 文件的源程序如下： clear %清空自定义变量 z1=4-6j;z2=6j;z3=10-6j;z4=5;us=15*exp(j*0);zl=4j;%输入各元件参数 zi=z1+(z2*(z3+z4)/(z2+z3+z4));%等效阻抗zi 的计算表达式uoc=us*z2/(z2+z3+z4);%开路电压uoc 的计算表达式u0=zl/(zi+zl)*uoc;%电感两端电压uo 的计算表达式disp('The magnitude of zi is'); %在屏幕上显示“The magnitude of zi is ”disp(abs(zi)) %显示等效阻抗zi 的模disp('The phase of zi is'); %在屏幕上显示“The phase of zi is ”disp(angle(zi)*180/pi)%显示等效阻抗zi 的辐角 disp('The magnitude of uoc is'); %在屏幕上显示“The magnitude of uoc is ” disp(abs(uoc))%显示开路电压uoc 的模

matlab绘制温度场复习过程

m a t l a b绘制温度场

通过在室内的某些位置布置适当的节点，采集回来室内的温湿度以及空气质量等实际参数。首先对室内空间建模，用一个无限细化的三维矩阵来模拟出室内的温度分布情况，针对采集回来的数据，采用插值法和适当次数的拟合函数的拟合，得出三维矩阵的实际值的分布，最后结合matlab软件绘制出计算出的温度场的三维图像。 一．数据的采集与处理 因为影响人的舒适感的温度层只是室内的某一高度范围内的温度，而温度传感器虽然是布置在一个平面内，但是采用插值法和拟合函数法是可以大致再现出影响人的舒适感的温度层的温度变化的。同时，在构建出的三维模型中，用第三维表示传感器层面的温度。 在传感器层面，传感器分布矩阵如下： X=【7.5 36.5 65.5】（模型内单位为cm） Y=【5.5 32.5 59.5】 Z=【z1 z2 z3； z4 z5 z6； z7 z8 z9；】（传感器采集到的实时参数） 采用meshgrid（xi，yi，zi，…）产生网格矩阵； 首先按照人的最小温度分辨值，将室内的分布矩阵按照同样的比例细化，均分，使取值点在坐标一定程度上也是接近于连续变化的，从而才能最大程度上使处理数据得来的分布值按最小分辨值连续变化！ 根据人体散热量计算公式：C=hc（tb-Ta） 其中hc为对流交换系数；

结合Gagge教授提出的TSENS热感觉指标可以计算出不同环境下人的对环境温度变化时人体温度感知分辨率，作为插值法的一个参考量，能使绘制出的温度场更加的符合人体的温度变化模式。 例如按照10cm的均差产生网格矩阵（实际上人对温度的分辨率是远远10cm 大于这个值的，但是那样产生的网格矩阵也是异常庞大的，例如以0.5cm为例，那么就可以获得116*108=12528个元素，为方便说明现已10cm为例）： [xi yi]=meshgrid(7.5:10:65.5,5.5:10:59.5) xi = 7.5000 17.5000 27.5000 37.5000 47.5000 57.5000 7.5000 17.5000 27.5000 37.5000 47.5000 57.5000 7.5000 17.5000 27.5000 37.5000 47.5000 57.5000 7.5000 17.5000 27.5000 37.5000 47.5000 57.5000 7.5000 17.5000 27.5000 37.5000 47.5000 57.5000 7.5000 17.5000 27.5000 37.5000 47.5000 57.5000 yi = 5.5000 5.5000 5.5000 5.5000 5.5000 5.5000 15.5000 15.5000 15.5000 15.5000 15.5000 15.5000 25.5000 25.5000 25.5000 25.5000 25.5000 25.5000 35.5000 35.5000 35.5000 35.5000 35.5000 35.5000 45.5000 45.5000 45.5000 45.5000 45.5000 45.5000 55.5000 55.5000 55.5000 55.5000 55.5000 55.5000 产生网格矩阵之后，就可以在测得的实时数据的基础上，通过相关的温度场的专业的估算函数，以及相关的数值处理函数来估计整个分布面（有最小的分辨率）上的温度了。即在这些函数的基础之上，对参数进行一些必要的处理。

路基温度场数值模拟及变化规律研究

Open Journal of Transportation Technologies 交通技术, 2018, 7(1), 11-19 Published Online January 2018 in Hans. https://www.wendangku.net/doc/2a7777303.html,/journal/ojtt https://https://www.wendangku.net/doc/2a7777303.html,/10.12677/ojtt.2018.71002 Research on the Numerical Simulation and Change Rules of the Subgrade Temperature Field Lei Xu, Yunliang Li, Lun Ji, Yiqiu Tan School of Transportation Science and Engineering, Harbin Institute of Technology, Harbin Heilongjiang Received: Jan. 1st, 2018; accepted: Jan. 15th, 2018; published: Jan. 22nd, 2018 Abstract Numerical analysis model of the subgrade temperature field was established. Based on the pro-gram ANSYS, the distribution and the time-varying properties of the subgrade temperature field were analyzed. Results show that for the vertical subgrade temperature field, the temperature in-creases with the depth increasing, and the subgrade frost depth at the coldest time in January is about 2.0 m. The ambient temperature affects the horizontal subgrade temperature field within a scope of about 2.0 m, which is the same as the subgrade frost depth. Temperature change trend of various positions in the subgrade during a year basically agrees with the ambient temperature change trend. The temperature of the pavement surface is basically the same as the ambient tem-perature, while the ambient temperature affects a little on the subgrade temperature, and soil in the depths of the subgrade keeps permafrost, or seasonal frozen. Keywords Subgrade, The Temperature Field, Numerical Analysis, ANSYS 路基温度场数值模拟及变化规律研究 徐垒，李云良，纪伦，谭忆秋 哈尔滨工业大学，交通科学与工程学院，黑龙江哈尔滨 收稿日期：2018年1月1日；录用日期：2018年1月15日；发布日期：2018年1月22日 摘要 建立了路基温度场的数值分析模型，基于ANSYS软件分析了路基温度场的分布规律及时变特性。研究表

二维导热物体温度场的数值模拟教程文件

二维导热物体温度场的数值模拟

金属凝固过程计算机模拟题目：二维导热物体温度场的数值模拟 Solidworks十字接头的传热分析 作者：张杰 学号： S2******* 学院：北京有色金属研究总院 专业：材料科学与工程 成绩： 2015 年 12 月

二维导热物体温度场的数值模拟 图1 二维均质物体的网格划分 用有限差分法模拟二维导热物体的温度场，首先将二维物体划分为如图1所示的网格，x ?与y ?可以是不变的常量,即等步长,也可以是变量(即在区域内的不同处是不同的),即变步长?如果区域内各点处的温度梯度相差很大,则在温度变化剧烈处,网格布得密些,在温度变化不剧烈处,网格布得疏些?至于网格多少,步长取多少为宜,要根据计算精度与计算工作量等因素而定? 在有限的区域内,将二维不稳定导热方程式应用于节点 ,)i j （ 可写成: ,2222 ,i j P P p i j T T T C x y ρλτ?????=+ ?????? ,1 , ,()i j P P P i j i j T T T οτττ+-???= +? ????? () , 1 , , 1 ,22 2()i j P P P P i j i j i j T T T T x x x ο+--+???=+? ????? () , ,1 , ,122 2()i j P P P P i j i j i j T T T T y y y ο+--+???=+? ?????τ?、x ?、y ? 当τ?、x ?、y ?较小时，忽略()οτ?、2()x ο?、2 ()y ο?项。当x y ?=?时，即x 、y 方向网格划分步长相等?最后得到节点 ,)i j （的差分方程: ()1 , ,0 1 , 1 , ,1 ,1 ,4P P P P P P P i j i j i j i j i j i j i j T T F T T T T T ++-+-=++++- 式中: () 02 p F C x λτ ρ?= ??

matlab-SIMULINK仿真实例资料

二并联杆数控螺旋面钻头尖刃磨机的机构仿真 一、仿真原理一、实训题目：全自动洗衣机控制系统 实训目的及要求： 1、掌握欧姆龙PLC的指令，具有独立分析和设计程序的能力 2、掌握PLC梯形图的基本设计方法 3、培养分析和解决实际工程问题的能力 4、培养程序设计及调试的能力 5、熟悉传输带控制系统的原理及要求 实训设备:： 1、OMRON PLC及模拟实验装置1台 2、安装CX-P编程软件的PC机1台 3、PC机PLC通讯的RS232电缆线1根 实训内容： 1、分析工艺过程，明确控制要求 （1）按下启动按扭及水位选择开关，相应的显示灯亮，开始进水直到高（中、低）水位，关水。 （2）2秒后开始洗涤。 （3）洗涤时，正转30秒停2秒；然后反转30秒停2秒。 （4）循环5次，总共320秒，然后开始排水。排水后脱水30秒。 图1 全自动洗衣机控制 2、统计I/O点数并选择PLC型号 输入：系统启动按钮一个，系统停止按钮一个，高、中、低水位控制开关三个，高、中、低液位传感器三个，以及排水液位传感器一个。

输出：进出水显示灯一盏，高、中、低水位显示灯各一盏，电机正、反转显示灯各一盏，排水、脱水显示灯灯各一盏。 PLC的型号：输入一共有9个，考虑到留有15％～20％的余量即9×（1＋15％）＝10.35，取整数10，所以共需10个输入点。输出共有8个，8×（1＋15％）＝9.2，取整数9，所以共需9个输出点。可以选OMRON公司的CPM1A/CPM2A 型PLC就能满足此例的要求。 3、I/O分配 表1 全自动洗衣机控制I/O分配表 输入输出 地址名称地址名称 00000 启动系统按钮01000 排水显示灯 00001 高水位选择按钮01001 脱水显示灯 00002 中水位选择按钮01002 进、出水显示灯 00003 低水位选择按钮01003 高水位显示灯 00004 排水液位传感器01004 中水位显示灯 00005 停止系统按钮01005 低水位显示灯 00006 高水位液位传感器01006 电机正转显示灯 00007 中水位液位传感器01007 电机反转显示灯 00008 低水位液位传感器 4、PLC控制程序设计及分析 实现功能：当按下按钮00000,中间继电器20000得电并自锁,按下停止按钮00005，中间继电器20000掉电。中间继电器20000为系统总启动。 实现功能：当按下按钮00001,中间继电器20001得电并自锁；当中间继电器20002、20003、20004、20007任意一个为ON，或按下停止按钮00005，或01000、01001为ON时，中间继电器20001掉电。

实验2 Matlab绘图操作

实验2 Matlab 绘图操作 实验目的： 1、 掌握绘制二维图形的常用函数； 2、 掌握绘制三维图形的常用函数； 3、 掌握绘制图形的辅助操作。 实验内容： 1. 设sin .cos x y x x ?? =+ ??+?? 23051，在x=0~2π区间取101点，绘制函数的曲线。 2. 已知： y x =21，cos()y x =22，y y y =?312，完成下列操作： （1） 在同一坐标系下用不同的颜色和线性绘制三条曲线； （2） 以子图形式绘制三条曲线； （3） 分别用条形图、阶梯图、杆图和填充图绘制三条曲线。 3. 已知：ln(x x e y x x ?+≤??=??+>??2 0102 ，在x -≤≤55区间绘制函数曲线。 4. 绘制极坐标曲线sin()a b n ρθ=+，并分析参数a 、b 、n 对曲线形状的影响。 5．在xy 平面内选择区域[][],,-?-8888， 绘制函数z =的三种三维曲面图。 6. 用plot 函数绘制下面分段函数的曲线。 ,(),,x x f x x x x ?+>? ==??+

（1）.y x =-205 （2）sin()cos ,sin()sin x t t t y t t π=?≤≤? =?303 详细实验内容： 1.设sin .cos x y x x ?? =+ ??+?? 23051，在x=0~2π区间取101点，绘制函数的曲线。 >> x=(0:2*pi/100:2*pi); >> y=(0.5+3*sin(x)/(1+x.^2))*cos(x); >> plot(x,y) 2.已知： y x =2 1，cos()y x =22，y y y =?312，完成下列操作： （1）在同一坐标系下用不同的颜色和线性绘制三条曲线； >> x= linspace(0, 2*pi, 101); >> y1=x.*x; >> y2=cos(2x); >> y3=y1.*y2; plot(x,y1,'r:',x,y2,'b',x,y3, 'ko')

振动流化床内温度场和流场的数值模拟

硕士学位论文开题报告及论文工作计划书 课题名称振动流化床内温度场和流场的数值模拟 学号1000614 姓名张 专业机械设计及理论 学院机械工程与自动化 导师张 选题时间2011 年09 月01 日 东北大学研究生院 2011年10月20日

填表说明 1、本表一、二、三、四、五项在导师指导下如实填写。 2、学生在通过开题后一周内将该材料交到所在学院、研究所。 3、学生入学后第三学期应完成论文开题报告，按有关规定，没有完成开题报告的学生不能申请论文答辩。

一、立论依据 课题来源、选题依据和背景情况、课题研究目的、理论意义和实际应用价值 （一）课题来源和背景情况 通常人们把热物理方法去湿的过程称为“干燥”，其特征是采用加热、陷湿、减压或其他能量传递的方式使物料中的湿分产生挥发、冷凝、升华等相交过程与物体分离以达到去湿目的。由于干燥的产品便于加工、运输、贮存和使用，干燥在国民经济的各部门有着广泛的应用。化学工业的产品如肥料、染料、无机盐到医药工业、粮食、食品、饲料的生产过程均离不开干燥，产品经过干燥以后具有良好的扩散性、均性。正确地完成干燥过程有利于保证和改进产品的质量，同时对提高生产效率，促进国民经济的发展有十分重要的作用。现代干燥技术在国民生产中应用的程度与一个国家的综合国力和国民生活质量的水平密切相关，从某种意义上说，它标志着这个国家国民经济和社会文明的发展程度。现在中国的经济处于飞速的发展期，各行业也处于发展的转型期，干燥设备制造业也亟需进行改进和优化，以便适应大环境的发展。 传统的工业生产普遍采用的干燥技术主要有：厢式干燥、隧道干燥、转筒干燥、转鼓干燥、带式干燥、盘式连续干燥、卧式桨叶式干燥、流化床干燥、改型流化床干燥、喷动床干燥、喷雾干燥、气流干燥、真空冷冻干燥、太阳能干燥、微波和高频干燥、红外热辐射干燥等。此外在各个行业，例如谷物、水果和蔬菜、木材、茶叶、乳品、中药材等行业也有适合自身特点的专有干燥技术。这些传统的干燥技术发展历史较长、成熟可取，在我国以及世界已经得到广泛的应用。近些年来，国际上涌现出一批新型的干燥技术，作为代表的有：脉冲燃烧干燥、对撞流干燥、冲击穿透干燥、声波场干燥、超临界流体干操、过热蒸汽干操、接触吸附干燥等等。这些新技术相对于传统干燥技术在机理上有一定的突破，但在工业化应用方面仍有待于完善。 流化床干燥是现代干燥技术的一种，是60年代发展起来的一种干燥技术，目前在化工、轻工、医药、食品以及建材等方面都得到了广泛的应用。由于干燥过程中固体颗粒悬浮在干 燥介质中，因而流体与固体接触面较大，热容量系数可达8000～25000 3 /() KJ m h C? ??（按干 燥器总体积计算），又由于物料剧烈搅动，大大减小了气膜阻力，因而热效率较高，可达60%～80%（干燥结合水时为30%～40%）。流化床干燥装置密封性能好，传动机械又不接触物料，因而不会有杂质混入，这对要求纯度高的制药工业来说也是十分重要。