短梁的約束以及荷載  約束左邊  右邊中部施加荷載

fixeddofs=1:2*(nely+1);

    alldofs=1:2*(nely+1)*(nelx+1);

    freedofs=setdiff(alldofs,fixeddofs);

    loaddof=2*(nely+1)*nelx+nely+2;

    F=sparse(loaddof,1,-1,2*(nely+1)*(nelx+1),1);

問題是如何施加均布荷載 以及多重荷載