16 lines
820 B
Mathematica
16 lines
820 B
Mathematica
|
function ke=ele_stiff_matrix(element,node,ele_id,A,D,t)
|
|||
|
|
i=ele_id;
|
|||
|
|
a1=node(element(i,2),1)*node(element(i,3),2)-node(element(i,3),1)*node(element(i,2),2);
|
|||
|
|
a2=node(element(i,3),1)*node(element(i,1),2)-node(element(i,1),1)*node(element(i,3),2);
|
|||
|
|
a3=node(element(i,1),1)*node(element(i,2),2)-node(element(i,2),1)*node(element(i,1),2);
|
|||
|
|
b1=node(element(i,2),2)-node(element(i,3),2);
|
|||
|
|
b2=node(element(i,3),2)-node(element(i,1),2);
|
|||
|
|
b3=node(element(i,1),2)-node(element(i,2),2);
|
|||
|
|
c1=-node(element(i,2),1)+node(element(i,3),1);
|
|||
|
|
c2=-node(element(i,3),1)+node(element(i,1),1);
|
|||
|
|
c3=-node(element(i,1),1)+node(element(i,2),1);
|
|||
|
|
|
|||
|
|
B(:,:,i)=1/(2*A)*[b1,0,b2,0 ,b3,0;0,c1,0,c2,0,c3;c1,b1,c2,b2,c3,b3];
|
|||
|
|
%<EFBFBD><EFBFBD>ԪB<EFBFBD><EFBFBD><EFBFBD><EFBFBD>
|
|||
|
|
|
|||
|
|
ke=B(:,:,i)'*D*B(:,:,i)*t*A; %<EFBFBD><EFBFBD>Ԫ<EFBFBD>նȾ<EFBFBD><EFBFBD><EFBFBD>
|