有限单元法作业提交全部打包

蔡永昌老师班物体，根据对称型取1/4作为分析模型，如图(b)所示。划分为4个三角形单元，6个结点。3如果受到如图(c)1.1Inputdata.mfunction%%输入数据 %输出数据 % %globalgMatgNdtgEltgEInfo%Material E,Poissonratio, gMat=[10.00.0%Node gNdt=[02;01;11；00;10;20];%Element node1,node2,gElt=[123;245;253;35 forie=1:ieltgEInfo(ie)=1;%Boundary type=1--fixeddisp;type=2--point *, gBct=[10-1.00.00.00.00.021.00.00.00.00.00.01；21.00.00.00.00.01；41.01.00.00.00.01；50.01.00.00.00.00.01；61.00.00.00.00.0

nodal cement:x, 0000 3[-0.0879-由[Bm,A]=Tri3BMtr(xe)、De=DMtrElas(Mat)得到单元1的单元代入{ε}=[B]{a^e}和{σ}=[D]{ε},1的应变和应力。令xe=[02;01;11],=00010010001001AMat=\[1,0\];De=DMtrElas(Mat)De=000000\{ε\}=\[-0.0879;-\{σ\}=\[-0.0879;-1.将垂直于1-6边上的均布荷载转化为等效结点荷载。在In-putdata.m中定义：gBct=\[1-0.5-0.50000 3-1-10002；6-0.5-0.500002;110000021000001;41100001;5010001;6010000E=1000v=0.25。function%处理施加上的荷载边界条件(局部坐标系下%输入数据 %输出数据 %globalgNdtgBctglobalgFA%%node1=node2gBct(ib+1,1);%gBct1、22个结点Cor1(1:2,1)=Cor2(1:2,1gNdt(node2,:)相应的结点坐len=sqrt(dDx*dDx+dDy*dDy);%val10=gBct(ib,4:5);%结点1的局部坐标下的荷载 %结点2的局部坐标下val10=val20trf'*val20'转换到整体坐标系下%val2=(val20/3.0+val10/6.0)*len计算等效forgFA(irowgFA(irow)+val1(ind)组集第gFA(irowgFA(irow)+val2(ind)组集第functionFe=eLoopBoundaryLoad(ie)globalgMatgNdtgEltgEInfogBctfunctiongnodenum=findnode(gNdt)fori=1:iefor%forn=1:tfunction%%输入数据 %输出数据 globalgBctglobalgKAgFAiBct=size(gBct,1);forib=1:iBctiNode=iType=ifiType==2 forind=1:2gFA(irow)=gFA(irow)+gBct(ib,ind+1);if forind=1:2if(gBct(ib,ind+1)==1) functionF=Triload(ie,Fe)globalgEltgNdtfori=1:EN globalgMatgNdtgEltgEInfogBctglobalgKAgFAglobal %InputDatatry;%11结点作用分布荷载%InputData_Canttry;%138结点作用分