Home

Descargar

1. AsIDA uE o ans 5 a Implementaci n us E 1 IMPLEMENTACION 1 1 Implementaci n El anexo del presente documento contiene todo el c digo que corresponde el programa que aqu se intenta explicar El c digo que aqu se expone est implementado de tal forma que hay un c digo fuente que es el que va llamando a las correspondientes subrutinas Esta subrutina principal se ha llamado main_Fun1000 es la que controla todo el programa y es la que hay que ejecutar desde la ventana de comandos de MATLAB para iniciar el programa Al ejecutar main_Fun1000 se crea una ventana que lee el n mero de materiales que componen la esquina y si esta es abierta o cerrada Desde aqu se llama a la subrutina intromat incluida en main_Fun90 Esta subrutina es la que lee el n mero de materiales que componen la esquina y en funci n de eso crea una nueva ventana preguntando sobre el tipo de material que es cada uno de ellos Adem s desde intronmat se llama a otras cuatro subrutinas introMaterial calculate2 calculate y getroots La subrutina introMaterial es la que lee el tipo de material que es cada uno de los que componen la esquina y en funci n de eso llama a las siguientes subrutinas que cargan los datos del material Estas son en funci n de si se trata de un material is tropo ort tropo o transversalmente 66e cc is tropo plso12 i pOrto4 i y pTransIso i r
2. entaci n u E 3 1 Esquinas abiertas y cerradas El archivo Opened crea una ventana que pregunta sobre las condiciones de contorno e interfaz de la esquina ejecuta a su vez la subrutina calcular3 que llama a las subrutinas frictionless_o i bounded_o i Free_o i Fixed_o i Simetry_o i Antisimetry_o i Ur_restricted_o i Ur_allowed_o i U3_restricted_o i U3_allowed_o i BC_Friction_o i en funci n de las condiciones que se impongan En el caso de las condiciones de contorno homog neas las subrutinas correspondientes activan una ventana por cada condici n de contorno en la que se introduce el ngulo inicial o final seg n corresponda En el caso de la condici n de interfaz pegado deslizamiento y fricci n se introduce el ngulo que posiciona la interfaz y en los casos de fricci n se introduce adicionalmente el coeficiente de fricci n Estos datos se leen y se introducen en las correspondientes matrices de condiciones de contorno e interfaz en la siguiente subrutina que se ejecuta dentro de main_Fun1000 que es la subrutina calculate La subrutina Closed crea una ventana que pregunta sobre las condiciones interfaz de la esquina y ejecuta a su vez la subrutina calculate3 que llama a las rutinas frictionless i friction i bounded_i i en funci n de las condiciones que se impongan En los casos de co
3. ente sin problemas de interferencias entre variables o en sus resultados En la siguiente figura se muestra un diagrama de flujos para intentar clarificar el proceso que se sigue en el programa Implementaci n 5o da gt nass trictionless_o frictionless BC _friction_oli friction _o i bounded_o bounded_o i Free_ofi Fixed_o i Simmetry_o i Antisimmetry_ofi Ur restricted_ofi Ur allowed_ofi U3_restricted_ofi 13 allowed_o i main_Fun1000 O Kemal4 Kima 4 0 Figu 3 2 Esquema del programa
4. es en funci n de si estamos en casos en los que exista fricci n o no En los casos en los que no exista fricci n se podr n ejecutar las subrutinas newtonroots o bisecroots en funci n de si se quiere usar el m todo de Newton o el de la bisecci n para calcular las ra ces del problema dentro del intervalo real de O a 1 Dichas subrutinas crean unas ventanas en las que se introducen los valores cerca de los cuales se encuentran las soluciones del problema y mediante las subrutinas solnewton y bisection calculan las ra ces del problema y las almacenan para mostrarlas posteriormente En los casos en los que exista fricci n se podr ejecutar la subrutina solnonsquare3 que utiliza el algoritmo de Levenberg Marquardt para calcular las ra ces Esta subrutina crea unas ventanas en las que se introducen los valores de 4 entre O y 1 y entre r y m cerca de los cuales se encuentran las soluciones del problema y mediante las subrutinas fsolnonsquare3 calcula las ra ces del problema y las almacenan para mostrarlas posteriormente Tambi n dentro de getroots se llama a la subrutina showRoots que permite visualizar las ra ces por pantalla y escribe el archivo de texto Roots txt en el que se almacenan las ra ces soluci n del problema de la esquina Adem s esta funci n borra todos los par metros anteriores usados durante el c lculo del problema lo que permite el uso del programa nuevam
5. espectivamente Estas subrutinas est n fuera del archivo principal y se guardan en archivos independientes que el c digo lee El archivo plso12 i crea una ventana que lee los datos de materiales is tropos y llama a su vez a la subrutina introE71 que calcula la matriz de transferencia del material El archivo pOrto4 i crea una ventana que lee los datos de materiales ort tropos y llama a su vez a la subrutina introE40 que calcula la matriz de transferencia del material El archivo pTransIso i crea una ventana que lee los datos de materiales transversalmente is tropos y llama a su vez a la subrutina introElso que calcula la matriz de transferencia del material Tanto las subrutinas plso12 i pOrto4 i y pTransIso i como las subrutinas introE71 introE40 y introElso se encuentran guardadas en archivos independientes de la rutina principal de forma que si se quisiera a adir alg n otro tipo de material a los c lculos se podr a hacer a adiendo una nueva subrutina y s lo habr a que llamarla desde introMaterial Una vez que se ha ejecutado introMaterial la funci n principal guarda las matrices de transferencia de cada material y es la subrutina calculate2 la que lee si la esquina es abierta o cerrada y en funci n de eso se ejecutan la subrutina Opened o Closed Dichas subrutinas tambi n se encuentran fuera del archivo principal Implem
6. ndici n de interfaz de pegado deslizamiento y fricci n se introduce el ngulo que posiciona la interfaz y adicionalmente en los casos de fricci n se introduce el coeficiente de fricci n Estos datos se leen y se introducen en las correspondientes matrices de condiciones de contorno e interfaz en la siguiente subrutina que se ejecuta dentro de main_Fun1000 que es la subrutina calculate La funci n calculate monta las matrices de condiciones de contorno e interfaz seg n corresponda lee cuales son las matrices de transferencia de cada material y monta la matriz de transferencia global del problema de la esquina multilateral en funci n de si esta es abierta o cerrada y de si hay condiciones de fricci n o no En el caso de que la matriz globla de trasferencia sea cuadrada es decir no haya fricci n dentro de esta subrutina se dibuja la parte real del determinante para valores de 4 entre 0 y 1 para de esa manera identificar los ceros de la ecuaci n dd GH En el caso de que la matriz de transferencia global sea rectangular incluya condiciones de fricci n dentro de la subrutina calculate se llamar a la subrutina procesador 7 que dibuja los valores de omin para los distintos valores de 4 entre O y 1 y entre r y m para de esa manera identificar los ceros de la ecuaci n AIDA Op gt m a Implementaci n us x o ans ES AE La subrutina getroots es la que calcula las ra c

Descargar

Contents

Download Pdf Manuals

Related Search

Related Contents