abaqus单元节点选择与重构程序

clc; clear; nodfix= xlsread('d:\tdat.csv'); nodnum= xlsread('d:\node_num.csv'); [M N]=size(nodnum); nodnum1=zeros(M*N,1); [MM NN]=size(nodfix); k=0; for i=1:M for j=1:N nodnum1(k+j)=nodnum(i,j); end k=k+N; end nodn1=nodfix(:,1); len=M*N; nodloc=zeros(len,NN); for i= 1:len nodno=nodnum1(i); num=find(nodn1==nodno); nodloc(i,:)= nodfix(num,:); end prop = wblrnd(3.0e10,6,len,1); locprop= [nodloc prop]; pason=zeros(len,1); pason(:,:)=0.17; abaprop= [prop pason nodloc(:,2:end)]; save nodloc; % plot(te1(:,1),te1(:,2),'y-',te2(:,1),te2(:,2),'go') % % Get node number in the elements % elm= xlsread('d:\elm.csv'); % elm1=elm(:,2:end); % [M0 N0]= size(elm1); % arr_elm=zeros(M0, N0); % nodnum2=nodnum1(1:end-6); % for i=1:M0 % for j=1:N0 % nodno= elm1(i,j); % l= find(nodnum2==nodno); % arr_elm(i,j)=l; % end % end % % for obtain detT in the computation % E= 5.2785e+09; % v=0.17; % p=2400; % Le=0.06; % lamb= E*v/(1+v)/(1-2*v); % mu=E/2/(1+v); % cd=sqrt((lamb+2*mu)/p); % detT= Le/cd;