2018A高温作业服设计 (显式 隐式差分代码)(12)

LEN1=int8(x(1)/dx(1)) 1; %记录每段接触点的位置，外界与第一段接触点的位置
LEN2=LEN1 x(2)/dx(2);%第一段与第二段接触点的位置
LEN3=LEN2 x(3)/dx(3);%第二段与第三段接触点的位置

LEN4=LEN3 x(4)/dx(4);%第三段与第四段接触点的位置
B=zeros(LEN4,1);
B(1,1)=2*dt*h(1)*Tout/(rho(1)*c(1)*dx(1));