matlab有限元编程实例梁,平面梁单元MATLAB有限元程序.txt

平面梁单元MATLAB有限元程序.txt

function beam2nclear allclose allclearclose BEAM2 -disp;disp PROGRAM BEAM2 ;disp Beam Bending Analysis ;disp T.R.Chandrupatla and A.D.Belegundu ;disp;Data;Bandwidth;Stiffness;ModifyForBC;BandSolver;ReactionCalc;Output; function Data -function Data;global NN NE NM NDIM NEN NDNglobal ND NL NCH NPR NMPC NBWglobal X NOC F AREA MAT SMI Sglobal PM NU U MPC BT REACTglobal CNSTglobal TITLE FILE1 FILE2global LINP LOUTglobal NQdispblanks1;FILE1 Data File Name ,s;LINP FILE1,r;FILE2 Output Data File Name ,s;LOUT FILE2,w;DUMMY fgetsLINP;TITLE fgetsLINP;DUMMY fgetsLINP;TMP str2numfgetsLINP;NN, NE, NM, NDIM, NEN, NDN dealTMP1,TMP2,TMP3,TMP4,TMP5,TMP6;NQ NDN * NN;DUMMY fgetsLINP;TMP str2numfgetsLINP;ND, NL, NMPC dealTMP1,TMP2,TMP3;NPR1; E- Coordinates -DUMMY fgetsLINP;for I1NNTMP str2numfgetsLINP;N, XN,dealTMP1,TMP21NDIM;end- Connectivity -DUMMY fgetsLINP;for I1NETMP str2numfgetsLINP;N,NOCN,, MATN, SMIN .dealTMP1,TMP21NEN, TMP2NEN, TMP3NEN;end- Specified Displacements -DUMMY fgetsLINP;for I1NDTMP str2numfgetsLINP;NUI,,UI, dealTMP1, TMP2;end- Component Loads -DUMMY fgetsLINP;F zerosNQ,1;for I1NLTMP str2numfgetsLINP;N,FNdealTMP1,TMP2;end- Material Properties -DUMMY fgetsLINP;for I1NMTMP str2numfgetsLINP;N, PMN, dealTMP1, TMP2NPR1;end- Multi-point Constraints B1*QiB2*QjB0if NMPC 0DUMMY fgetsLINP;for I1NMPCTMP str2numfgetsLINP;BTI,1, MPCI,1, BTI,2, MPCI,2, BTI,3 .dealTMP1,TMP2,TMP3,TMP4,TMP5;endendfcloseLINP; function Bandwidth -function Bandwidth;global NN NE NM NDIM NEN NDNglobal ND NL NCH NPR NMPC NBWglobal X NOC F AREA MAT SMI Sglobal PM NU U MPC BT STRESS REACTglobal CNSTglobal TITLE FILE1 FILE2global LINP LOUT- Bandwidth uation -NBW 0;for N1NENABS NDN*absNOCN, 1 - NOCN, 2 1;if NBW 0SNR, NC SNR, NC SEI, J;endendendendendend function ModifyForBC -function ModifyForBC;global NN NE NM NDIM NEN NDNglobal ND NL NCH NPR NMPC NBWglobal X NOC F AREA MAT SMI Sglobal PM NU U MPC BT REACTglobal CNSTglobal NQ- Decide Penalty Parameter CNST -CNST 0;for I 1NQif CNST I2; IR I2; endIC absI2 - I1 1;SIR, IC SIR, IC CNST * BTI, 1 * BTI, 2;FI1 FI1 CNST * BTI, 1 * BTI, 3;FI2 FI2 CNST * BTI, 2 * BTI, 3;end function BandSolver -function BandSolver;global NN NE NM NDIM NEN NDNglobal ND NL NCH NPR NMPC NBWglobal X NOC F AREA MAT SMI Sglobal PM NU U MPC BT REACTglobal CNSTglobal NQ- Equation Solving using Band Solver -dispSolving using Band Solverbansol.m;F bansolNQ,NBW,S,F; function ReactionCalc -function ReactionCalc;global NN NE NM NDIM NEN NDNglobal ND NL NCH NPR NMPC NBWglobal X NOC F AREA MAT SMI Sglobal PM NU U MPC BT REACTglobal CNSTfor I 1NDN NUI;REACTI CNST * UI - FN;end function Output -function Output;global NN NE NM NDIM NEN NDNglobal ND NL NCH NPR NMPC NBWglobal X NOC F AREA MAT SMI Sglobal PM NU U MPC BT REACTglobal CNSTglobal TITLE FILE1 FILE2global LINP LOUT dispsprintfOutput for Data from file sn,FILE1;fprintfLOUT,Output for Data from file sn,FILE1;dispTITLE;fprintfLOUT,sn,TITLE;disp Node X-Displ Rotation;fprintfLOUT, Node X-Displ Rotationn;I1NN; print a matrixdispsprintf 4d 15.4E 15.4En,I,F2*I-1,F2*I;fprintfLOUT, 4d 15.4E 15.4En,I,F2*I-1,F2*I;- Reaction Calculation -dispsprintf DOF Reaction;fprintfLOUT, DOF Reactionn;for I 1NDN NUI;R CNST * UI - FN;dispsprintf 4d 15.4E,N,REACTI;fprintfLOUT, 4d 15.4En,N,REACTI;endfcloseLOUT;dispsprintfThe Results are available in the text file s, FILE2;