基础平面刚架MATLAB法度模范(国外英文资料).docVIP

基本平面刚架MATLAB程序(国外英文资料) The % plane just a MATLAB program 2003.9.16, 2007.2.28, 2009.10, 2009.10 2011.10, 2011, 2010 % * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * Percent variable declaration % NPOIN NVFIX NFPOIN NFPRES The number of points summed, the number of units, the number of constraints, the number of force nodes, the number of non-nodes % COORD LNODS YOUNG % structure node coordinates array, unit definition array, elastic modulus A % FPOIN FPRES FORCE FIXED The % nodal force array, the non-nodal force array, the total load vector, the constraint information array % HK DISP The total stiffness matrix, the node displacement vector % * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * Format short e % sets the output type Clear % clears the memory variable FP1 = fopen ( x.txt , rt) % open the initial data file % read in control data NELEM = fscanf (FP1, % d, 1); % unit number NPOIN = fscanf (FP1, % d, 1). % nodes NVFIX = fscanf (FP1, % d, 1) % number of constraints NFPOIN = fscanf (FP1, % d, 1). The number of nodes of the payload NFPRES = fscanf (FP1, % d, 1); % non-node load number YOUNG = fscanf (FP1, % f, 1); % elastic modulus % reads the structure information LNODS = fscanf (FP1, % f, [4, NELEM]) Unit definition: left, right node, area, moment of inertia (total NELEM group) COORD = fscanf (FP1, % f, [2, NPOIN]) % coordinates: x, y coordinates (total NPOIN group) FP1, % f, [4, NFPOIN] ] % node force (NFPOIN group) : the force node, the X direction force (to the right), % Y direction (up positive), M couple (counterclockwise) FPRES = fscanf (FP1, % f, [4, NFPRES]) % of the total % NFPRES group) : unit number, load type, load size, distance left % load type 1 - uniformly distributed load 2 - horizontal concentration 3 - longitudinal concentration FIXED = fscanf (FP1, % f, NVFIX) % constraint information: the corresponding displacement encoding (total NVFIX