2020A参数折磨人系列

clc,clear
format short g
% T15=173
% T6=198
% T7=230
% T89=257
% T1011=25
% v=78/60 %单位cm/s
T15=175
T6=195
T7=235
T89=255
T1011=25
v=70/60 %单位cm/s
dt=0.5 %每隔0.5 s焊接区域中心的温度
t=0:dt:435.5/v;
% ds=30.5
% di=5%interval
% darea=25 %炉前区域和炉后区域
T0=Temp(v,t,T15,T6,T7,T89,T1011)%各个时间点T0(i)表示第i个点在x=0的温度
% t=0对应T0下标1，t=i*0.5，对应下标i+1  t/0.5+1
len=length(T0);
f=[];
for i=1:len %下标i对应的是(i-1)*0.5*v这个点f=[f;(i-1)*0.5*v,T0(i)];
end
plot(f(:,1),f(:,2))
xlabel('x/cm','FontSize',10);
ylabel('u/度','FontSize',10,'Rotation',0);
title('炉内温度变化','FontSize',12); disp('----------------------------')
% 0.15mm 0.075
x=0.00075;
dx=0.00005;
m=x/dx;
m=15
n=length(t)-1 %n是区间段的个数，length是点的个数%已经在Temp函数(T0)中加上t=0
mm=m+1;
nn=n+1;
% % m=5
% % r=1;
k=1;
c=1;
rho=1;
h=1;
arpha=0.0000003;%k/(c*rho);
r=dt*arpha/(dx*dx)
mm
A=zeros(mm,mm);
A=(1+2*r)*eye(mm)+diag(-r*ones(1,mm-1),1)+diag(-r*ones(1,mm-1),-1);
A(1,1)=-1-dx*n/k;
A(mm,mm)=1;
A(mm,mm-1)=-1;
A(1,2)=1;
u=zeros(mm,nn); %每个时间对应一列u(:,1)=u(:,1)+25*ones(mm,1);% t=0时，电路板1~mm都是25°
for i=2:5%nnb=[]  ;
b=u(:,i-1);
b(1,1)=-dx*7*T0(i);%h/k
b(mm)=0;tmp=A\b;u(:,i)=u(:,i)+tmp;bAtmp
end
% u
% for i=2:length(t)
%     u(:,n+1)=u(:,n)+a^2*dt/dx^2*A*u(:,n);% a^2*(u关于x的二阶导)
%     u(1,n+1)=m1(n+1); %边界条件
%     u(end,n+1)=m2(n+1);
%
% end
% % u

Temp.m

function T=Temp(v,t,T15,T6,T7,T89,T1011)
%   t1=25/v  %炉前
%   t2=t1+(5*30.5+4*5)/v%1~5
%   t3=t2+(30.5+5)/v %小温区6
%   t4=t3+(30.5+5)/v  %小温区7
%   t5=t4+(2*30.5+2*5)/v%小温区8`9
% dis=[25,30.5,5,]
x1=25;%炉前区域
x2=x1+(30.5+5)*4+30.5;%小温区15及其之间4个间隙
x3=x2+5;%小温区5~小温区6之间间隙
x4=x3+30.5;%小温区6
x5=x4+5;%小温区6~小温区7之间间隙
x6=x5+30.5;%小温区7
x7=x6+5;%小温区7~小温区8之间间隙
x8=x7+30.5*2+5;%小温区89及其之间1个间隙
x9=x8+5;%小温区89~小温区1011之间间隙
x10=x9+30.5*2+5+25  %435.5 小温区1011、炉后区域之间间隙
% t=0:0.5:x10;T_air=25len=length(t);T=zeros(len-1,1);for i=1:len-1tt=i*0.5;d=tt*v;if(d<=x1)  k=(T15-T_air)/(25/v);T(i)=T_air+k*(tt-0);elseif(d>x1&&d<=x2) %小温区15及其之间间隙T(i)=T15;elseif(d>x2&&d<=x3)k=(T6-T15)/((x3-x2)/v);T(i)=T15+k*(tt-x2/v);elseif(d>x3&&d<=x4) T(i)=T6;elseif(d>x4&&d<=x5)%67之间的间隙k=(T7-T6)/((x5-x4)/v);T(i)=T6+k*(tt-x4/v);elseif(d>x5&&d<=x6)T(i)=T7;elseif(d>x6&&d<=x7)k=(T89-T7)/((x7-x6)/v);T(i)=T7+k*(tt-x6/v);elseif(d>x7&&d<=x8)T(i)=T89; elseif(d>x8&&d<x9)%小温区89~小温区1011之间间隙k=(T89-T1011)/((x9-x8)/v);T(i)=T89-k*(tt-x8/v);  elseif(d>x9&&d<=x10)T(i)=T1011;endendT=[25;T];
% T=?
% clc,clear
end

2020A参数折磨人系列相关推荐

华住数据库_华住内控人系列故事（四）技术领先篇——搭建大数据风险数据仓，实现自助取数...
我们不断的学习和引入好的"利器",以新的技术提升内控内审工作效率,工欲善其事必先利其器. 01.从"依赖,等待"到"靠自己",搭建大数据风险 ...
报表 labview_【LabVIEW懒人系列教程小白入门】1.3LabVIEW数据类型
上期作业解答: 如何对齐控件(动图) 第一种方法: 第二种方法: 第三种方法: 第四种方法: 大家尝试至少用第一种方法编写出该题. 作业讲解结束 labview中需要熟悉的控件类型有: 基本数据类型: ...
【LabVIEW懒人系列教程-小白入门】1.13LabVIEW程序结构之事件结构
上期1.2作业讲解: 下面是效果演示: 今日讲解程序框图中的事件结构运用方式, 事件结构必须要存在超时事件分支,否则程序会报错且无法执行顾名思义,事件结构就是满足某个事件执行对应事件的功能分支,事件 ...
【LabVIEW懒人系列教程-小白入门】1.15LabVIEW程序结构之跑马灯
上期1.14作业讲解: 本期教大家如何利用程序结构编写小程序<跑马灯> 程序要求,多数布尔灯依次排序,按照从左往右的顺序进行亮灭操作. 如图: 方法①/利用While循环组合顺序结构实现 ...
【LabVIEW懒人系列教程-小白入门】1.16LabVIEW程序结构之小试身手
上期1.15作业讲解: <分数序列求和> For循环组合移位寄存器实现计算: 本期带来两个习题,主要目的:运用前期所学知识点进行反复练习,提高对于函数运用的熟练度,了解Labview程序组 ...
折磨人小游戏之《一个都不能死》html5网页小游戏源代码
介绍: 今天给大家分享一款相当折磨人的小游戏,玩这个小游戏简直可以让人抓狂,在下也是玩了几局,最多也就活到20多秒就撑不下去了,脑子和手都不够用. 这款<一个都不能死>小游戏html源码可 ...
【LabVIEW懒人系列教程-小白入门】1.2LabVIEW前面板控件讲解
上次课程作业编写:d=a+b+c 前面板创建4个控件:数值型输入控件:a,b,c:数值型显示控件d ctrl+e 切换到程序框图界面,添加2个 "+"算法,然后连线,a,b,c,d ...
折磨人，中国人的最爱
喜欢在外人面前显示自己的气派大,朋友多,交际广,吃得开,有面子几乎是中国人几千年来与生俱来的天性. 其实中国人绝大多数并不会真正计较自己的朋友多一个或者少一个,关键在于不能让外人觉得你朋友少,&quo ...
笔记14 笨办法中折磨人的习题ex23.py，人类语言与计算机语言
笔记14 笨办法中折磨人的习题ex23.py,人类语言与计算机语言开始学习习题23了,但为了那个languages.text的下载,费了老大劲好像还是没有下载全.后来载下来了,那个代码文件录下来,却 ...

2020A参数折磨人系列

2020A参数折磨人系列相关推荐

最新文章

热门文章