1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45
| clear all, close all, clc,
Lx = 0.6; Nx = 60; hx = Lx / (Nx);%hx=hy
Ly = 0.3; Ny = 30; hy = Ly / (Ny);
H = 0.005000;
Vf = 0.500000;
Kxx = 1.0e-09;
Kxy = 0;
Kyy = 1.0e-010;
mu = 0.1804;
Number_nodes = (Nx+1) * (Ny+1);
Number_elem = (Nx ) * (Ny );
fid = fopen('wrong.dmp','w');
fprintf(fid,'Number of nodes : %5.0f \n',Number_nodes);
for j = 1:Ny+1;
for i = 1:Nx+1;
node_id = (j-1)*(Nx+1)+i;
x = (i-1)*hx;
y = (j-1)*hy;
z = 0;
NODE(node_id,:) = [node_id x y z];
fprintf(fid,'%6.0f %12.6f %12.6f %12.6f \n',NODE(node_id,:));
end
end
fprintf(fid,'Number of elements : %5.0f \n',Number_elem);
for j = 1:Ny;
for i = 1:Nx;
elem_id = (j-1)*(Nx )+i;
k1 = (j-1)*(Nx+1)+i;
k2 = k1+1;
k3 = k1+(Nx+1)+1;
k4 = k3-1;
ELEMENTS(elem_id,:) = [elem_id 4 k1 k2 k3 k4 H Vf Kxx Kxy Kyy];
fprintf(fid,'%6.0f %4.0f %6.0f %6.0f %6.0f %6.0f %32.6f %15.6f %15.4e %15.4e %15.4e \n',ELEMENTS(elem_id,:));
end
end
fprintf(fid,'Resin Viscosity model NEWTON \n');
fprintf(fid,'Viscosity : %15.6f \n',mu);
fprintf(fid,'# Done ************************** \n');
fclose(fid); |
Partager