R/EtiquetasGLNF.R

In IzaelNava/Armaduras: Stiffness method: Truss Analysis

#' @keywords function
#' @name EtiquetasGLNF
#'
#' @title Obtención de deltas (diferencias de las coordenadas xN y xF, yN y yF y zN y zF).
#'
#'
#' @details
#' Los deltas de cada elemento se obtienen como la diferencias de las coordenas N 
#' (considerado como el inicio del elemento) y F (final del elemento).
#'
#' @param nodos Nodos del modelo como data.frame
#' @param elementos Elementos del modelos como data.frame
#' @param truss Armadura en 2D o 3D. Si truss == 2, se cosidera 2D, de otra forma es 3D.
#'
#'
#' @return
#' Devuelve los deltas de las corrdenadas x, y y/o z.
#' Por ejemplo: deltax = xF - xN ... para cada elemento.
#'
#' @example 
#' nodos <- data.frame(nodos=1:2, x=c(0, 6), y=c(0, 0), z=c(0, 0))
#' elementos <- data.frame(Elemento=1, nN=1, nF=2)
#' truss <- 3
#' deltas <- dltaxyz(nodos, elementos, truss)
#' 
#' @export
#'


EtiquetasGLNF <- function(EtiquetasGl, elementos, numElem, truss){

  EtiquetasNodosNF <- data.frame(matrix(0, numElem, (truss*2+1)))
  if(truss == 2){
    colnames(EtiquetasNodosNF) <- c("Elemento", "Nx", "Ny", "Fx", "Fy")
  } else {
    colnames(EtiquetasNodosNF) <- c("Elemento", "Nx", "Ny", "Nz", "Fx", "Fy", "Fz")
  }
  EtiquetasNodosNF[,1] <- elementos[,1]
  
  
  
  
  
  for(i in 1:numElem){
    EtiquetasNodosNF[i, 2:(truss+1)] <- EtiquetasGl[elementos[i,2], 2:(truss+1)]
    EtiquetasNodosNF[i, (truss+2):(truss*2+1)] <- EtiquetasGl[elementos[i,3], 2:(truss+1)]
  }

  
  return(EtiquetasNodosNF)
}