PF-feap - Preguntas Frecuentes sobre el uso de feap ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ José M.ª Goicolea, 28 marzo 2003 PREGUNTA: ~~~~~~~~ > - Como no consigo manejar bien el comando PERS (no sé que es el vector que > se introduce como segunda entrada. LA primera, el punto de vista, sí)no > consigo ver la deformada. RESPUESTA: ~~~~~~~~~ El mandato PERSpective pide tres lineas. La primera es 0 si los datos son nuevos; la segunda son las coordenadas del punto de vista; la tercera las coordenadas del vector que quieres que sea vertical en tu vista. Ten en cuenta que este último dato es necesario, una cosa es que tú pongas el ojo en un sitio, y otra cómo te orientes alrededor del ojo. He incluido un ejemplo adicional de prácticas en la página web, el bóveda de Scordelis-Lo, en el que puedes ver cómo se usa el mandato pers en un entorno batch, he incluido comentarios explicativos en el fichero de datos. En interactivo sería similar. PREGUNTA: ~~~~~~~~ > - Al definir el material, tengo que introducir los momentos de inercia de la > sección pero, respecto a que ejes?, ¿qué entiende FEAP como eje X, eje Y, > Eje Z?,¿qué son Kx y Ky? > > - Además, pide un vector de referencia o nodo de referencia. ¿Eso qué > es?Como no hay ningún ejemplo en el manual, no sé a que se refiere. > > Adjunto te envío el fichero de entrada de datos para que le eches un vistazo > y me comentes las cuestiones RESPUESTA: ~~~~~~~~~ Se trata de ejes locales que definen la orientación de tu seccion transversal, estos ejes en feap se pueden definir mediante un vector de referencia, una direccion o un nodo de referencia. Efectivamente, esto no está bien descrito en el manual de feap, por lo que a continuación lo explico. Recomiendo usar un vector de referencia, llamémosle V, cuyas coordenadas (en ejes globales) se definen mediante el mandato reference vector dentro del bloque de material. Una vez dado V, los ejes locales (x,y,z) son: e_z = eje viga; e_y = e_z ^ V; e_x = e_y ^ e_z. (e_x, e_y, e_z son los versores de cada dirección local). Lo más fácil, si la viga está en un plano, aunque sea curva, es tomar V la dirección normal a este plano. Tienes un ejemplo (puente curvo continuo) que he añadido en las páginas web del curso donde está resuelto y explicado con comentarios. Los coeficientes k_x y k_y son los de deformación a cortante, me parece que si no se dice nada toma por omisión 5/6 que corresponde a sección rectangular (ver los apuntes que os dí en el curso). PREGUNTA: ~~~~~~~~ > En alguna parte del manual he creído leer que existe un a manual con > ejemplos del FEAP.¿Tienes el fichero?¿Dónde se puede conseguir? > > ¿Existe algún manual, ejemplos o aplicaciones sencillas por Internet en > FEAP? RESPUESTA: ~~~~~~~~~ Hay un manual de ejemplos disponible en las páginas web del curso. Desgraciadamente, no hay ejemplos con vigas 3D allí, espero que te sirvan mis explicaciones. hay más referencias y otras páginas de feap a través de las páginas web de R.L. Taylor, a las que puedes acceder a través de mis páginas del curso. PREGUNTA: ~~~~~~~~ Defino el parámetro en feap Ixy=10e-5, pero feap no lo lee correctamente y me asigna mal esta propiedad de la seccion transversal de las vigas RESPUESTA: ~~~~~~~~~ Solo se pueden usar dos letras en los parámetros, por ejemplo xy=10e-5. PREGUNTA: ~~~~~~~~ Los elementos frame de feap en 3D no convergen y dan unas deformaciones de flexión ridículas, quedan bloqueados. RESPUESTA: ~~~~~~~~~ Hay un error en los elementos de feap que tiene términos de inercia cruzados (Ixy) en pequeñas deformaciones. Para solventarlo puede usarse como parche la opción de grandes deformaciones, aunque en este caso lo que hace feap es ignorar el término de inercia cruzado. En la definición del material, hay que poner una línea con «finite» justo debajo de frame. Por ejemplo: mate,1 frame finite shear,on elastic,isotropic,E,nu cross,section,A,Ix,Iy,xy,Iz refe vect 0,0,1 PREGUNTA: ~~~~~~~~ ¿Por qué no funciona la generación con boun en el siguiente ejemplo? Deseo restringir los nodos h+2 hasta el h+8. Pues bien, sólo funciona para el h+2 y el h+8, pero los generados, no. ... BOUNdary !condiciones de contorno 1,0,1,1 !nodo,inc,codigo_x,codigo_y h+2,0,1,1 h+3,1,1,1 h+8,0,1,1 ... RESPUESTA: ~~~~~~~~~ Los códigos de contorno que se desean generar deben ponerse en negativo, así: ... BOUNdary !condiciones de contorno 1,0,1,1 !nodo,inc,codigo_x,codigo_y h+2,0,1,1 h+3,1,-1,-1 ! CODIGOS -1 PARA GENERAR h+8,0,1,1 !--fin de linea ... De hecho, en tu ejemplo podrías empezar a generar desde el h+2. Esto está en el manual de usuario de feap, aunque expresado de forma críptica. PREGUNTA: ~~~~~~~~ ¿Cómo se dibujan los modos de vibración (autovalores) calculados por feap? RESPUESTA: ~~~~~~~~~ hay que emplear el mandato plot,eigv,n1,n2,n3. Por ejemplo, para el dibujo del modo 3: plot,eigv,3. Los otros parámetros son para dibujar contornos de alguna variable sobre la malla. Hay que tener en cuenta que hay que poner antes el factor de magnificación para las mallas de modos, con plot,defo,v1,n1,v2; v1 es la magnificacion para deformadas, n1 sirve para la escala (-1 para mantener y poder superponer), y v2 para el factor de magnificacion de los modos. Este último factor (v2) NO VIENE DOCUMENTADO EN EL MANUAL. PREGUNTA: ~~~~~~~~ Estoy haciendo un cálculo dinámico y parece que no funciona la integración directa en el tiempo con el metodo de Newmark, definido con el siguiente mandato: transient,newmark,0.25,0.50 RESPUESTA: ~~~~~~~~~ Parece tratarse de un error de la versión 7.3 de FEAP. En su lugar se puede obtener el mismo resultado con alpha=1 en el método HHT-alpha: transient,alpha,0.25,0.50,1.0 PREGUNTA: ~~~~~~~~ Al hacer un modelo con evolucion en el tiempo, dinamico o cuasiestático, en el cual se emplea: batch ... prop,,1 tplo ... end pdis,3,1 pdis,3,2 se come la primera historia pedida. ¿qué está pasando? RESPUESTA: ~~~~~~~~~ "prop" busca tambien la especificacion de las historias de entrada de cargas despues del "end". Aunque se utilice sin argumentos tomando el valor por omision (variacion proporcional de la carga con el tiempo) necesita una linea en blanco. Por tanto, hay que poner: batch ... prop,,1 tplo ... end pdis,3,1 pdis,3,2 PREGUNTA: ~~~~~~~~ Tengo un modelo que se deforma mucho, y deseo ir haciendo dibujos de forma que la malla inicial y las subsiguientes mallas queden todas superpuestas en la misma figura. ¿cómo puedo hacerlo sin que se salgan del area de dibujo que establecería feap en el primer dibujo? RESPUESTA: ~~~~~~~~~ Empleando los macros "plot,scale,v1" y "plot,center,v1,v2". Mientras no se haga "plot,wipe" estos valores se mantienen para los demás dibujos, que se superpondrán al primero. Por ejemplo, para la extensión de un elemento que crece mucho (de 2 a 8 en vertical), podria ponerse antes del primer dibujo podría ponerse: batch ... plot factor 0.35 plot center 0.5 0.25 plot mesh loop,time,20 time loop,tang,20 tang,,1 next,tang plot defo 0 -1 plot mesh end