53练习题三角形单元塑性程序.docVIP

习题： 程序及解答 主程序 子程序 2.1 输入数据 function InputData global gMat gNdt gElt gEInfo gBct %材料属性 % E Poisson ration gravity ,yield stress, *, gMat [座机电话号码 0.3 0 17.32e3 0]; %结点坐标 % gNdt [0 0 40 0 1.111 0 2.222 0 3.333 0 4.444 0 5.556 0 6.667 0 7.778 0 8.889 0 10 0 …… （数量较多，不一一罗列） 370 3 4 90 3 370 90 1 3]; % 单元材料号 ielt size gElt,1 ; %单元总数 for ie 1:ielt gEInfo ie 1; end %边界条件 % type 1--fixed disp; 1--结点固定；0--结点自由） % type 2--point load % node x y *, *, *, *, type gBct [ 83 1 1 0 0 0 0 1; 84 1 1 0 0 0 0 1; 85 1 1 0 0 0 0 1; 86 1 1 0 0 0 0 1; 87 1 1 0 0 0 0 1; 88 1 1 0 0 0 0 1; 89 1 1 0 0 0 0 1; 90 1 1 0 0 0 0 1; 47 1 1 0 0 0 0 1; 42 0 -11000 0 0 0 0 2]; 2.2 初始化变量值 function AssignInitValues global gMat gNdt gElt gBct gEInfo global gKA gFA global gEData gNData iNdt size gNdt,1 ;%总结点数 iElt size gElt,1 ;%总单元数 gKA zeros 2*iNdt,2*iNdt ; gFA zeros 2*iNdt,1 ; for ie 1:iElt iMatNum gEInfo ie ; vdMat gMat iMatNum,: ; gEData ie .strs_tot zeros 1,3 ; gEData ie .strs_inc zeros 1,3 ; gEData ie .strs_yds vdMat 4 ;%初始屈服应力 gEData ie .strs_yield vdMat 4 ;%屈服应力 gEData ie .strs_equ 0.0;%等效应力 gEData ie .strn_equ 0.0;%等效塑性应变 end for in 1:iNdt gNData in .disp_tot zeros 1,2 ;%结点in的总位移 gNData in .disp_inc zeros 1,2 ;%结点in的增量位移 end 2.3 增量法求解塑性方程 function IncreMethod global gNdt global gKA gFA global gEData gNData %处理荷载边界条件（计算总荷载） LoopBoundaryLoad; % idivid 50;%荷载分段数 Fdiv gFA/idivid; %初始化向量 nfre_max size gFA,1 ; Fcur zeros nfre_max,1 ; %开始增量计算 for idiv 1:idivid %在屏幕上显示载荷分段数 sDisplay sprintf Current load divide number:%d,idiv ; disp sDisplay ; %当前增量步总荷载 Fcur Fcur+Fdiv; iNdt size gNdt,1 ; gKA zeros 2*iNdt,2*iNdt ; %计算并组集当前步的整体刚度矩阵 LoopCalcGlobalStif; %荷载列阵初始化为0 gFA zeros nfre_max,1 ; %计算前一步应力的等效荷载（初始应力部分） LoopElePrevForce; %处理位移边界条件 LoopBoundaryDisp; %用于本步计算的未平衡荷载 Ferr Fcur+gFA; %solve the equation Disp gKA\Ferr; % nlen size gNdt,1 ; for in 1:nle