function edelta=elementdelta(ie,delta,element)
% 計算單元的變形列向量
i=element(ie,1);
j=element(ie,2);
edelta(1,1)=delta((i-1)*3+1);
edelta(2,1)=delta((i-1)*3+2);
edelta(3,1)=delta((i-1)*3+3);
edelta(4,1)=delta((j-1)*3+1);
edelta(5,1)=delta((j-1)*3+2);
edelta(6,1)=delta((j-1)*3+3);
return