jueves, 29 de marzo de 2012

Propiedades Geométricas de Secciones

Programa para  la Casio fx-4500P, que calcula el área, centroides, inercias, etc. de un polígono cualquiera expresada por sus coordenadas.

Mcl:N"NºVERTICES":Defm N*2
I=0:Lbl 0:I=I+1:I<=N=>{V}:X[I]=V"X":{W}:X[N+1+I]=W"Y":Goto 0 ∆
X[I]=X[1]:X[N+1+I]=X[N+2]
A=0:B=0:C=0:D=0:E=0:F=0:I=0:Lbl 1:I=I+1:I<=N=>A=A+0.5*(X[I+1]+X[I]*(X[N+I+1]-X[N+I+2])
B=B+(X[I+1]-X[I])*((X[N+I+2]^2+X[N+I+1]^2+X[N+I+1]*X[N+I+2])/6
C=C+(X[N+I+1]-X[N+I+2])*((X[I+1]^2+X[I]^2+X[I]*X[I+1])/6)
D=D+(X[I+1]-X[I])*((X[N+I+1]^3+X[N+I+2]^3+X[N+I+1]*X[N+I+2]^2+X[N+I+1]^2*X[N+I+2])/12)
E=E+X[N+I+1]-X[N+I+2])*((X[I]^3+X[I+1]^3+X[I]*X[I+1]^2+X[I]^2*X[I+1])/12)
F=F+(X[I+1]-X[I])*((X[I+1]*(3*X[N+I+2]^2+2*X[N+I+2]*X[N+I+1]+X[N+I+1]^2)+X[I]*(X[N+I+2]^2+2*X[N+I+2]*X[N+I+1]+3*X[N+I+1]^2))/24
Goto 1 ∆
A:"AREA="▲G=C/A:"XG="▲H=B/A:"YG="▲D:"IX="▲E:"IY="▲"PIXY="▲F-A*G*H:"PIG="▲
J=D-A*H^2:"IXG="▲K=E-A*G^2:"IYG="▲
√(J/A):"RXG="▲:√(K/A):"RYG="▲
Defm 0

DEFINICIÓN DE VARIABLES:
IG,YG = CENTRO GEOMÉTRICO DEL POLÍGONO
IX,IY = MOMENTO DE INERCIA CON RESPECTO A LOS EJES X,Y
PIXY = PRODUCTO DE INERCIA CON RESPECTO A LOS EJES X,Y
IXG,IYG = MOMENTO DE INERCIA CON RESPECTO A LOS EJES X,Y CENTROIDALES
PIG = PRODUCTO DE INERCIA CON RESPECTO A LOS EJES X,Y CENTROIDALES
RXG,RYG = RADIO DE GIRO CON RESPECTO A LOS EJES X,Y CENTROIDALES