7.6.3.1 Características del grupo

lectura en IELMNT
READ (LECG,1000) NPAR
FORMAT(20I4)

nota columna variable definición
  1-4 NPAR(1) dar el número 8
  5-8 NPAR(2) (NUME) cantidad de elementos en este grupo
NUME $>$ 0
1 9-12 NPAR(3) (INDNL) código que indica el tipo de análisis no lineal
= 0, por defecto igual a "1"
= 1, material no lineal únicamente
  13-16 NPAR(4) (IPARX) orden del elemento paraxial
= 0, únicamente
2 17-20 NPAR(5) (ITYP2D) código que indica el tipo de modelación adoptada para este grupo
= 1, sólo deformación plana
3 21-24 NPAR(6) (NDIM) cantidad de "grados de libertad locales" de los nodos de estos elementos
= 2, g.d.l. $u_y$ y $u_z$ (únicamente mecánica)
= 0, por defecto igual a 2
4 25-28 NPAR(7) (MXNODS) cantidad máxima de nodos empleados para describir un elemento del grupo:
MXNODS= 2 ó 3 unicamente
por defecto igual a 3
  29-32 NPAR(8) no utilizado
  33-36 NPAR(9) no utilizado
5 37-40 NPAR(10) (NINT) orden de la integración numérica
por defecto igual a 2
$1 \leq \mbox{\texttt{NINT}} \leq 4$
  41-44 NPAR(11) no utilizado
  45-48 NPAR(12) no utilizado
  49-52 NPAR(13) no utilizado
6 53-56 NPAR(14) (INCTYP) tipo de campo incidente
= 0, no hay campo incidente
= 1, onda plana armónica
= 2, onda plana de tipo Ricker
= 3, onda plana arbitraria
= 4, campo arbitrario
7 57-60 NPAR(15) (MODEL) número del modelo elegido para el material
= 1, únicamente
por defecto igual a 1
  61-64 NPAR(16) (NUMMAT) cantidad de "juegos de propiedades" distintas para el modelo de comportamiento seleccionado
$\texttt{NUMMAT} \geq 1$
no hay valor por defecto
  65-68 NPAR(17) no utilizado
6 69-72 NPAR(18) (NMBNO) cantidad de nodos paraxiales
(sólo si INCTYP=4)
6 73-76 NPAR(19) (NMBDT) cantidad de pasos de tiempo
(solamente si INCTYP=4)
  77-80 NPAR(20) (IEXPL) código que define si el grupo será tratado implícita o explícitamente en un cálculo dinámico con código de integración IOPE=10 (§1.7).
= 1, grupo explícito (únicamente)

notas

  1. En la versión actual emplear NPAR(3)= 0 ó 1.
    Solamente las no linealidades materiales pueden ser modeladas: los desplazamientos y las deformaciones se suponen infinitamente peque?os.
  2. No se pueden emplear elementos paraxiales en esfuerzos planos o en simetría de revolución.
  3. Los elementos paraxiales sólo pueden ser empleados en mecánica pura.
  4. NPAR(7) determina el número máximo de nodos empleados para describir cualquier elemento del grupo.
    Un mínimo de 2 nodos y un máximo de 3 nodos se emplean para describir estos elementos.
  5. Un orden de integración de 2 es suficiente.
    Si el elemento se deforma significativamente, un orden de integración más importante debe ser empleado.
  6. Para definir un campo incidente existen dos posibilidades:
    1. una onda plana:
      • armónica
      • de Ricker
      • arbitraria
    2. o bien un campo arbitrario
    1. Para las ondas planas se debe definir los materiales, las características de las ondas y las coordenadas del punto de control. El punto de control es un punto ubicado en la superficie libre (Figura 7.38).
      Los dos primeros casos corresponden a curvas analítica donde la amplitud de los desplazamientos se define en el punto de control.
      Para la ondas planas arbitrarias se debe definir las aceleraciones horizontales y verticales en el punto de control bajo la forma de curvas de carga de acuerdo a las indicaciones del anexo §[*].
      En ambos casos, el programa calcula los campos de esfuerzos y de velocidades necesarias a nivel de los paraxiales, tomando en cuenta los retardos en el tiempo de llegada de las ondas incidentes directas o reflejadas.

      Figura 7.38:
      \begin{figure}\begin{center}\epsfig{file=figures/npar8_fig1.eps,clip=,width=.4\textwidth}
\end{center} \end{figure}

      Se asume que el punto de control corresponde a la solución del semi-espacio y que no está afectada por la presencia de irregularidades geométricas o de comportamiento.
      Atención: Para INCTYP=1, 2, y 3 es la composición de las ondas $S_v$ et $P$ la que permite obtener el campo esperado en el punto de control.

    2. Para un campo arbitrario se deben calcular los campos de velocidades y de esfuerzo bajo la forma de un archivo PREFIJO.prxi con la estructura siguiente:
      NMBNO,NMBDT,T0PRX,DTPRX

      NUM(1),RHO(1),Cs(1),Cp(1)
      NUM(2),RHO(2),Cs(2),Cp(2)
      ..........................................
      ..........................................
      XNUM(n),RHO(n),Cs(n),Cp(n) n = NMBNO

      Vy(1),Vy(2),......................Vy(NMBDT) 1er nodo
      Vz(1),Vz(2),.......................Vz(NMBDT)
      SIGyy(1),SIGyy(2),........SIGyy(NMBDT)
      SIGzz(1),SIGzz(2),..........SIGzz(NMBDT)
      SIGyz(1),SIGyz(2),.........SIGyz(NMBDT)

      Vy(1),Vy(2),......................Vy(NMBDT) 2do nodo
      Vz(1),Vz(2),.......................Vz(NMBDT)
      SIGyy(1),SIGyy(2),........SIGyy(NMBDT)
      SIGzz(1),SIGzz(2),..........SIGzz(NMBDT)
      SIGyz(1),SIGyz(2),.........SIGyz(NMBDT)
      ..........................................
      ..........................................

      Vy(1),Vy(2),......................Vy(NMBDT) $n$ésimo nodo
      Vz(1),Vz(2),.......................Vz(NMBDT)
      SIGyy(1),SIGyy(2),........SIGyy(NMBDT)
      SIGzz(1),SIGzz(2),..........SIGzz(NMBDT)
      SIGyz(1),SIGyz(2),.........SIGyz(NMBDT)
      avec
      NMBNO: cantidad de nodos paraxiales
      NMBDT: cantidad de pasos de tiempo
      T0PRX: tiempo inicial
      DTPRX: paso de tiempo
      NUM($i$): número del nodo $i$
      RHO($i$) : densidad de masa en el nodo $i$
      Cs($i$): velocidad de la onda $S$ en el nodo $i$
      Cp($i$): velocidad de la onda $P$ en el nodo $i$
      Vy($j$): velocidad horizontal en el nodo $i$ en el instante $j$
      Vz($j$): velocidad vertical en el nodo $i$ en el instante $j$
      SIGyy($j$): esfuerzo $yy$ en el nodo $i$ en el instante $j$
      SIGzz($j$): esfuerzo $zz$ en el nodo $i$ en el instante $j$
      SIGyz($j$): esfuerzo $yz$ en el nodo $i$ en el instante $j$
  7. Un sólo modelo de comportamiento está disponible: medio lineal elástico (MODEL=1).

Esteban Saez 2010-12-30