============================================================= ============================================================= ======== DICCIONARIO DE VARIABLES DEL PROGRAMA MEF1D ====== ============================================================= ============================================================= Interpretacion de Variables: ____________________________ Elementos...................... Elementos Finitos VPs (Variables Primarias)...... Desplazamientos Nodales VSs (Variables Secundarias).... Tensiones TI (Terminos Independientes)... Fuerzas nodales CT (Carga Tipo) ............... 1 -> peso propio, 2 -> centrifuga, 3 -> cargas puntuales DIMENSIONES =========== ntype......................> Tipo de analisis (0=>Elasticidad Lineal 1D) nelem......................> No. de elementos nnode......................> No. de nodos por elemento [numeración local] npoin......................> No. de puntos nodales [numeración global] ndime......................> No. de dimensiones del espacio ndofn......................> No. de VPs por nodo nstre......................> No. de VSs por punto nmats......................> No. de materiales [distintos] nprop......................> No. de propiedades [de cada material] ngaue......................> No. de puntos de Gauss por elemento nload......................> No. de CTs [consideradas en el analisis] ncase......................> No. de hipotesis de carga [combinaciones de CTs] nvfix......................> No. de VPs prescritas nevab......................> No. de VPs por elemento [=ndofn*nnode] ntotv......................> No. total de VPs [=ndofn*npoin] ntotc......................> No. total de terminos independientes [=ncase*ntotv] nrigi......................> No. de componentes de la matriz de rigidez global [simetrica, almacenada en perfil] DEFINICION DEL PROBLEMA ======================= lnods(inode,ielem).........> numero global del nodo "inode" del elemento "ielem" matno(ielem)...............> material del elemento "ielem" props(iprop,imats).........> propiedad "iprop" del material "imats" nofix(i,ivfix).............> numero global del nodo [i=1] y grado de libertad coartado [i=2] correspondiente a la VP prescrita "ivfix" presc(ivfix)...............> valor de la VP prescrita "ivfix" ngfix(ivfix)...............> numero global de la VP prescrita "ivfix" coord(idime,ipoin).........> coordenada global "idime" del punto nodal "ipoin" NIVEL DE ELEMENTO [al procesar un elemento en concreto] ======================================================= posge(idime,igaue).........> coordenada local "idime" del punto de Gauss "igaue" weige(igaue)...............> peso del punto de Gauss "igaue" shape(inode,igaue).........> valor de la funcion de forma del nodo "inode" en el punto de Gauss "igaue" derie(idime,inode,igaue)...> derivada de la funcion de forma del nodo "inode" en el punto de Gauss "igaue" respecto a la coordenada local "idime" elcod(idime,inode).........> coordenada global "idime" del nodo "inode" gcord(idime,igaue).........> coord. glob. "idime" del punto de Gauss "igaue" xjace(idime,jdime,igaue)...> derivada de la coordenada global "idime" respecto a la coordenada local "jdime" [matriz jacobiana] en el punto de gauss "igaue" djacb(igaue)...............> det. jacobiano en el punto de gauss "igaue" xjaci(jdime,idime,igaue)...> derivada de la coordenada local "jdime" respecto a la coordenada global "idime" [inversa de la matriz jacobiana] en el punto de gauss "igaue" cartd(idime,inode,igaue)...> derivada de la funcion de forma del nodo "inode" en el punto de Gauss "igaue" respecto a la coordenada global "idime" dmatx(istre,jstre).........> componente "istre,jstre" de la matriz D bmatx(istre,ievab,igaue)...> componente "istre,ievab" de la matriz B=LN en el punto de Gauss "iguae" smatx(istre,ievab,igaue)...> componente "istre,ievab" de la matriz S=DB en el punto de Gauss "iguae" dvolu(igaue),..............> volumen de integracion en el punto de gauss "igaue" [=espesor*jacobiano*peso de integracion] estif(ievab,jevab).........> componente "ievab,jevab" de la matriz K=BtDB [matriz de rigidez del elemento] eload(ievab)...............> componente "ievab" del vector b [vector de TI del elemento] eldis(idofn,inode).........> VP "idofn" del nodo "inode" [desplazamientos del elemento] strsg(istre,igaue,icase)...> VS "istre" en el punto de Gauss "igaue" correspondiente al caso de carga "icase" [tensiones en puntos de Gauss del elemento] TRATAMIENTO DE CARGAS ===================== fload(iload,icase).........> coeficiente de seguridad de la CT "iload" correspondiente a la hipotesis de carga "icase" kload(iload)...............> tipo de CT correspondiente a la carga "iload" [1=peso propio, 2=centrifuga, 3=cargas puntuales] nilod(iload)...............> número de nodos cargados [en la CT "iload" es de tipo 3=cargas puntuales] rilod(i,iload).............> componente "i" del vector de aceleracion grav. [si la carga "iload" es de tipo 1=peso propio], o componente "i" del vector velocidad de rotacion [si la carga "iload" es de tipo 2=centrifuga] pforc(idofn)...............> componente "idofn" de la carga que actua en un nodo bforc(idofn,igaue).........> componente "idofn" de una fuerza por unidad de volumen en el punto de Gauss "igaue" REACCIONES ========== treac(ivfix,icase).........> reaccion correspondiente a la VP prescrita "ivfix" correspondiente a la hipotesis de carga "icase" totre(idofn)...............> componente "idofn" de la suma total de reacciones MATRIZ DE RIGIDEZ Y VECTOR DE TÉRMINOS INDEPENDIENTES (=VECTOR DE VPs) ====================================================================== lpntn(ipoin)...............> renumeracion global del punto nodal "ipoin" lpamc(itotv)...............> posicion en el vector "rgmat" en la que se almacena la componente "itotv,itotv" de la matriz de rigidez global simetrica [vector de punteros del almacenamiento en perfil] (¡OJO!=> Los nodos están renumerados) gload(itotv,icase).........> componente "itotv" del vector de TI correspondiente al caso de carga "icase" [antes de resolver el sistema], y componente "itotv" del vector de VPs correspondiente al caso de carga "icase" [después de resolver el sistema] (¡OJO! => Los nodos están renumerados) rgmat(irigi)...............> componente "itotv,jtotv" de la matriz de rigidez global simetrica almacenada en perfil (¡OJO! => Los nodos estan renumerados)