8/8/2019 hpcremona

1/7

Cremona

Un programa de dimensionado de cerchas

Mariano Vazquez Esp*

Ondara, 14 de abril de 2007.

y

x

3 m

4 m 4 m

45o

100 kN

1

2 3

4

100 kN

C

A

D

B

E

F

y

x

N1 N2 N3

N4 N5

100kNR1

R2

R3

Figura 1.

Para analizar y calcular una cercha como la de laFigura 1, bastan las ecuaciones de equilibrio de la es-tatica, debido a que el numero de incognitas de fuer-za iguala exactamente el numero de ecuaciones dispo-nibles. En efecto, la articulacion del apoyo izquierdocuenta como dos bielas, cada una aportando una fuer-za ( R1, vertical; y R2, horizontal); el apoyo sobre elplano inclinado de la derecha (suponiendo que no exis-te rozamiento) cuenta como otra biela, es decir, otra

fuerza de direccion conocida ( R3, perpendicular al pla-no de apoyo); por ultimo, la solicitacion de cada barra(Ni) es otro par de fuerzas iguales y opuestas desco-nocidas; en total 8 fuerzas, igual a las ocho ecuacionesque pueden obtenerse considerando las dos ecuacionescartesianas que pueden escribirse por cada nudo.

En cerchas como esta la gran mayora las reac-ciones (Ri) pueden calcularse con ecuaciones de equili-brio global (tomando momentos en el punto 2 se obtie-ne R3 = 71 kN; por lo mismo en el punto 4, R1 = 50kN;finalmente por equilibrio horizontal, R2 = 50 kN; tam-bien podra haberse trazado un funicular con el mismo

objetivo).Una vez determinadas las reacciones, si en cada par

de ecuaciones de nudo solo apareceren dos solicitacio-

*Departamento de Estructuras de Edificacion de la Universi-dad Politecnica de Madrid

nes desconocidas, las ecuaciones pueden resolverse gra-ficamente mediante poligonos vectoriales de nudo, taly como se ha hecho en el diagrama de Cremona re-presentado en la figura (tales diagramas se deben enrealidad a Maxwell1). Por ejemplo, el polgono CBAEC representa el equilibrio de todas las fuerzasque actuan sobre el nudo 2; del mismo modo ADFEA representa el equilibrio del nudo 3. En diagramas

semejantes, todas las fuerzas aparecen representadas aescala y su valor se obtiene simplemente midiendo, sinninguna complicacion. En este caso, las solicitacionesresultan ser las indicadas en el Cuadro 1, una distri-bucion simetrica como caba esperar.

Cuadro 1:

Ni Ai,min Ai iBarra nudos kN mm2 mm

1 12 -83 463 465 4,32 13 100 556 560 2,6

3 14 -83 463 465 4,34 23 17 93 200 1,65 34 17 93 200 1,6

Ni : solicitacionAi,min: area estricta, |Ni|/(0,18 kN/mm

2)Ai : area disenadai : alargamiento

Para un acero con lmite elastico e de 260 N/mm2,tension admisible de 180 N/mm2 (coeficiente de segu-ridad = 1, 44) y modulo de Young de 210 kN/mm2,

un dimensionado teorico podra conducir a las seccio-nes Ai indicadas en el cuadro. Notese que se trata demala teora, pues los dos pares comprimidos fracasa-ran en la vida real debido a la inestabilidad inherentea la compresion (pero esta es otra historia. . . ); note-

1Los procedimientos graficos para el analisis de cerchas consti-tuan la unica herramienta practica disponible hasta la aparicionde maquinas de calculo. Las figuras recprocas, introducidas porMaxwell en 1864, constituyen el nucleo del analisis grafico desolicitaciones. Sin embargo, su exposicion resulto tan abstractaque su enorme utilidad paso desapercibida. Un lustro despues,Jenkin (1869), mostro numerosos ejemplos practicos de la aplica-cion de las figuras recprocas de Maxwell, si bien se guio para sutrazado de las reglas inventadas por W. P. Taylor, un artesano

de la epoca que traba jaba para un constructor. Paralelamente,en Suiza, la obra Die graphische Statik de K. Culmann (1866y 1875) contiene una exposicion del metodo de las seccionesy del de los nudos, aunque sin apoyarse en la formulaci on deMaxwell. La estatica grafica llego al Mediterraneo de la manodel italiano Cremona Le figure reciproche nella statica grafica,1872, razon por la cual el diagrama de Maxwell fue conocidocomo diagrama de Cremona en la Pennsula Iberica.

1

8/8/2019 hpcremona

2/7

Figura 2.

En azul la geometra inicial; en rojo, la posicion de

equilibrio, con la estructura deformada.

se tambien que las barras 23 y 34 se dimensionancon una barra cilndrica de 16 mm de diametro, aun-que podra ponerse una menos gruesa. Sea como fuere,una vez dimensionadas las barras pueden calcularse susalargamientos i, tal y como se muestra en el cuadro.

Para averiguar la deformada un metodo divertido einstructivo consiste en dibujarla mediante un programa(como AUTOCAD), cuya precision en el trazo permi-ta distinguir tan pequenos alargamientos (desde luegomediante funciones ZOOM y cosas semejantes, de otromodo no se aprecia diferencia entre la geometra ini-cial y la deformada, vease la Figura 2, arriba). Unavez trazada la deformada, pueden medirse los despla-zamientos de los nudos: por ejemplo, el apoyo derechosube por el plano inclinado 4,5 mm, lo que determinael valor del desplazamiento de la Figura 1. Se puedecomprobar el requisito de rigidez y, en su caso, dar eldiseno global por valido y pasar al diseno de detalle. . .Este es el metodo mas exacto posible.2

Calculo analtico

El calculo analtico permite escribir las dos ecuacio-nes correspondientes a cada nudo (componentes hori-zontales y verticales) y calcular las reacciones a la vezque las solicitaciones. Como siempre, hay que operarbajo un convenio de signos; tambien hay que convenirsi se manejan fuerzas sobre el nudo o las fuerzas queel nudo ejerce sobre el resto. Aqu opto por la prime-ra posibilidad. Con estas consideraciones, la ecuacionhorizontal tendra la forma:

Xl( Ql) +

i

Xl( Ni) +

j

Xl( Rj) = 0

2Aunque la aparicion de los primeros ordenadores propicio lapaulatina desaparicion de los metodos graficos, la muy posterioraparicion y divulgacion de los programas de dibujo hacen posiblevolver a ellos, pues la que siempre fue su relativa desventaja,la imprecision del trazo, ha desaparecido completamente: entreestudiantes de arquitectura e ingeniera, supuestamente proclivesal dibujo, cabr a esperar su renacimiento.

en la que el operador Xl() representa la componen-te horizontal de la fuerza implicada en el nudo l: yasea la carga Ql, las solicitaciones Ni de las barras queacometen al nudo, o las reacciones Rj de cada bielaque el nudo tenga. Notese que las bielas pueden estarorientadas segun los ejes cartesianos (como en el apoyoizquierdo del ejemplo), pero no siempre pasa eso (talocurre en el apoyo derecho). La otra ecuacion de ca-da nudo es formalmente igual, si se sustituye Xl() porYl(), la componente vertical.

Si trasladamos las componentes conocidas a la dere-cha, el aspecto de las ecuaciones cambia a:

i

nkiNi +

j

rkjRj = Bk

en la que k es el numero de ecuacion considerada. Tan-to nki como rkj representan la contribucion al equili-brio de las fuerzas en barras y bielas; si la barra o labiela no acomete al nudo que se considera en la ecua-cion k, su aportacion sera nula (nki o rkj nulo). Bkrepresenta la componente de fuerza en la ecuacion kcambiada de signo, es decir, Xl( Ql) o Yl( Ql). Lasincognitas son los modulos, Ni y Rj , de las fuerzas, esdecir, {x1 . . . xe+v} = {N1, . . . N e, R1, . . . Rv}, siendo eel numero de barras (o elementos), y v el de bielas ocondiciones de sustentacion (vnculos); en consecuen-cia podemos escribir el sistema de todas las ecuacionescomo

e+v

i=1

Akixi = Bk para k = 1 . . . 2n

siendo n el numero de nudos. O en forma matricial,mucho mas compacta:

[Aki] {xi} = {Bk} o bien Ax = B

Resuelto el sistema, se obtienen tanto el valor de lassolicitaciones como el de las reacciones. Notese que encondiciones normales 2n = e + v y si la estructura esviable el sistema tiene una unica solucion. Si por elcontrario 2n > e + v el sistema no tiene, en general,solucion: la estructura es una mecanismo capaz de mo-verse dinamicamente bajo la accion de cargas genericaspero veanse los ejemplos mas adelante. Finalmente,si 2n < e+v y la estructura es viable, existiran infini-

tas soluciones posibles. En los dos ultimos casos, y contecnicas algebraicas bien conocidas como el meto-do de los mnimos cuadrados puede extraerse algunasolucion (si existe alguna) matematicamente significa-tiva (y puede que estructuralmente tambien), pero suanalisis solo puede hacerse considerando la rotura o ladeformacion elastica, y en general es necesario fijar eldimensionado de las barras previamente : el metodo deanalisis ha de ser necesariamente hiperestatico, al con-trario que en el primer caso donde las ecuaciones de laestatica bastan (metodo isostatico).

Calculo automatico

Matrices? Sumatorios? Calculadora! Puede. . . Si-go calculando cerchas con lapiz y papel (incluso, enocasiones, sin regla, a mano alzada), pero es tal la ma-yoritaria fascinacion que ejercen hoy las maquinas de

8/8/2019 hpcremona

3/7

calculo que quiza este justificado ofrecer al publico unprograma automatico: aqu va el que hice para una delas mas, la HP48GX.

Los nudos son simplemente puntos, de manera quepueden representarse mediante una lista de numeroscomplejos. Cada numero (x;y) contiene las coordena-das en los ejes elegidos. La secuencia de puntos en lalista los numera implcitamente. Para la estructura dela Figura 1, los puntos se representan por la lista:

{ (0;0) (-4;3) (0;3) (4;3) }

objeto que puede almacenarse bajo un nombre, porejemplo, P.

Las barras conectan nudos a pares, de manera quedando el nudo inicial y final quedan determinadas: unformato conveniente es una lista de vectores de la forma[ ni nf ], siendo ni y nf la posicion de los nudos de labarra en la lista P. Las barras de la figura se representancomo:

{ [1 2][1 3][1 4][2 3][3 4] }

y quedan numeradas implcitamente por su posicionen la secuencia. Lo mas seguro es almacenarla en unavariable, por ejemplo en B.

Por razones precisas, las bielas estan emparentadastanto con barras como con cargas. Cada biela repre-senta una fuerza desconocida a medias: desconocemossu magnitud pero no su direccion. Por ejemplo, en lafigura, cualquiera que sea su valor, R3 aportara al equi-librio del nudo componentes horizontal y vertical de va-lor 0, 71R3. De hecho, el valor 0, 71 es el coeficienterk8 (k = 7 o k = 8) de las ecuaciones analticas del

nudo 4. Estos dos coeficientes pueden agruparse en unnumero complejo y antecederles con el numero del nudoen el que se situa la biela. Como la fuerza es desconoci-da, solo interesa la proporcion entre sus componentes,de manera que el apoyo sobre el plano inclinado de lafigura puede representarse con una lista { 4 (-1;-1)} (aunque { 4 (-,71;-,71) } valdra igual de bien,pero sera mas largo de escribir), y el conjunto de lasbielas por la lista de todas ellas:

{ { 2 (0;-1) } { 2 (1;0) } { 4 (-1;-

1 ) } }

Notese que cada reaccion tiene su propio convenio designos, el dibujado. Sin embargo, en la lista de bielas elsigno de sus componentes se rige por el convenio generalde los ejes xy. Si la reaccion resulta negativa, actua ensentido contrario al dibujado. La lista de vnculos sealmacena en lugar seguro, por ejemplo en V.

Por ultimo, las cargas pueden expresarse con un con-venio similar: cada carga es una lista formada por elnudo en que se aplica y por el numero complejo forma-do por sus componentes; todas ellas se agrupan en unalista. En el caso de la figura, con una unica carga:

{ { 3 (0;100) } }

que se almacena en Q. Una misma cercha puede estarsometida a distintas cargas independientes, almacena-das en variables distintas: Q1, Q2, . . .

Los tres primeros objetos, P, B y V, definen una cerchaparticular. Junto con las cargas de un caso contienen

toda la informacion relevante para el analisis de ese ca-so. Notese lo apropiado que resulta disenar programasque leen informacion, en vez de hacer estupidas pre-guntas sobre el numero de nudos de la cercha y otraszarandajas: la misma informacion puede usarse siempreque sea necesario (y con distintos propositos) sin nece-sidad de teclearla otra vez. Ademas si se cometieronerrores al escribirla por vez primera, basta corregirlospara intentarlo una segunda.

Ahora lo unico que resta es escribir el programa queforme las matrices A y B y resuelva el sistema para ob-tener x; los programas para ello (y otros utiles para lacomprobacion de la estructura) se muestran en el ane-jo final, indicando ademas las instrucciones necesariaspara almacenarlos en la memoria de la calculadora.

Uso de los programas

Una vez escritos y almacenados todos los programasy la variable CST en un directorio3, y con la ruta actu-

al apuntando a el, pulsando CST apareceran CREM ,LONGI , MULV , DIVVE , DOT y OBJ . Los ulti-

mos son utiles aunque no esenciales y se explican masadelante; es Cremona el que se encarga de escribir yresolver el sistema de ecuaciones. Con NXT , aparece-ran el resto de definiciones que se hayan almacenadoen CST.

En cualquier subdirectorio desde all, con VAR setendra acceso a las variables donde se haya ido alma-cenando la informacion de una cercha particular (en el

ejemplo: P , B , etc.).4

Cremona. Recoge de la pila una lista de puntos, una

de barras, una de vnculos y una de cargas, y deja enla pila un vector con las solicitaciones de las barras yotro con la magnitud de las reacciones. Para el ejemplode la figura, la secuencia sera:

VAR P B V Q CST CREM

que deja en la pila el vector de solicitaciones (en laposicion 2) y el de reacciones (en la 1), en todo concor-dantes con las calculadas previamente. Como las soli-citaciones son de gran utilidad, lo suyo es guardarlas,por ejemplo en N.5

Dimensionado. Para el acero considerado en el

ejemplo, pueden obtenerse las areas estrictas con la se-cuencia N 0,18 , que arroja en la pila un vectorde areas en mm2 (con los mismos valores que los indica-dos en el Cuadro 1, en la columna Ai,min). Este vectorse puede editar para hacer retoques razonables6, y depaso suprimir el irrazonable signo negativo, heredadode las solicitaciones. Como de costumbre, lo mejor esalmacenarlo, por ejemplo en A.

3Probablemente el mejor sitio es HOME si se va a usar mucho.Un lugar razonable es un directorio especfico, como CERCHAS, delque cuelguen los directorios que contengan las cerchas objetode calculo.

4Con VAR en el directorio de programas se tendra acceso

directo al resto de los programas que figuran en el Anejo: STOP ,

SOLE y OPVEC .5En lo que sigue omitire las pulsaciones de CST o VAR ,

dando por supuesto que se pulsara lo que convenga para accedera la variable que se menciona.

6Para ajustar a catalogos comerciales, por ejemplo.

8/8/2019 hpcremona

4/7

Con las areas decididas y las longitudes de las barrasse puede calcular el volumen de la cercha. Longitudessuministra el vector de longitudes de barras tras cogerde la pila la lista de puntos y la de barras. El productoescalar de esos dos vectores da el volumen de la cercha,pero en mm2 m, dado que longitudes y areas estan endistintas unidades (ciertamente puede usarse la mismaunidad en ambos casos). En definitiva, la secuencia:

P B LONGI A DOT 1E6 7850

arroja en la pila el peso de la cercha, en kilogramos eneste caso.

Requisito de rigidez

Para comprobar si ademas de ligera es rgida pue-de emplearse el principio de los trabajos virtuales, quereza as:

Fi i =

Nj j

en la que F y N son un conjunto de fuerzas exterioresen los nudos y de solicitaciones en las barras que estanen equilibrio; por su parte, y son un conjunto dedesplazamientos de los nudos y los alargamientos queproducen en las barras (alargamientos compatibles). Elsuperndice simboliza que uno y otro conjunto nonecesitan guardar mas relacion entre s que referirsea la misma cercha. (De uno de los conjuntos se sueledecir que es virtual o imaginario respecto al otro.)

Para calcular el desplazamiento de la figura em-pleamos la misma cercha, pero cargada con una fuerzaunidad en el mismo punto, sentido y direccion de . Es-to determina las fuerzas exteriores y las solicitacionesa emplear en la ecuacion anterior. Como la deforma-cion es arbitraria respecto a las fuerzas, podemos es-coger la deformacion real de la estructura inicial, de laque conocemos o podemos calcular los alargamientosde las barras y solo desconocemos los desplazamientos.En la ecuacion anterior el primer termino quedara co-mo 1kN puesto que, a parte de la fuerza unidad, soloexisten las reacciones, aplicadas en puntos y direccionessin desplazamiento, y con contribucion nula al traba-jo virtual exterior. El segundo termino quedara como

Nj j , donde Nj representa las solicitaciones en las

barras bajo la carga unitaria en la direccion de , y j

los alargamientos reales bajo la carga de 100kN, vir-tuales e imaginarios para la carga unidad pero bienreales para los 100kN. En consecuencia,

=1

1 kN

Nj j

es decir, el producto escalar del vector solicitacionesproducido por la fuerza unidad y el vector alargamien-tos bajo 100 kN (salvo el factor 1kN1 necesario parala coherencia de las unidades).

Con la maquina la secuencia de calculos es co-mo sigue. La carga unidad se representa como { { 4

(0,71;-0,71) } }, como siempre almacenada en unavariable, por ejemplo QV. La secuencia:

P B V QV CREM DROP

deja en la pila el vector {Nj }. Ahora hace falta formarel vector de los alargamientos y realizar el producto

escalar (como se hizo para calcular el peso). Pero tam-bien se pueden calcular automaticamente. La formulapara una barra es:

=N

EAL

El problema es que tanto las solicitaciones como lasareas y las longitudes de las barras est an almacena-

das en forma de vectores y necesitamos un programaque calcule un vector a partir de otros dos, de tal ma-nera que el nuevo tenga por componentes la divisiono el producto de las correspondientes componentes delos dos primeros. Esto es lo que hacen DivVec y Mul-Vec: toman dos vectores de la pila y arrojan en ella elresultado de dividir o multiplicar componente a compo-nente.7 Ambos usan OpVec, que es la version general:coge dos vectores de la pila y un programa, y arrojaen ella el resultado de evaluar el programa sobre lascomponentes correspondientes de ambos vectores, su-cesivamente.

La secuencia N A DIVVE arroja en la pila el vec-tor de tensiones; 210 calcula el vector de deforma-

ciones; por ultimo, P B LONGI MULV multiplicacada deformacion por la longitud y deja en la pila losalargamientos.

Ahora DOT ejecuta el producto escalar entre lassolicitaciones de la carga unidad (que estaban en lapila hace rato) y los alargamientos recien calculados.El resultado es cuanto sube el extremo derecho por elplano inclinado (o cuanto baja, si el valor es negativo)En que unidades? kN/mm2/(kN/mm2) m, es decir,en metros.

Resolucion de ecuaciones

Hasta ahora nada he dicho de Soleq. Se trata de unescueto programa (con una sola instruccion!), que esejecutado por Cremona casi al final, cuando las matri-ces A y B han sido ya formadas. En la HP48GX dividirmatrices tiene sentido. Si el sistema de ecuaciones seexpresa como Ax = B, formalmente x = B/A es lasolucion buscada; la division aqu simboliza la inver-sion de una matriz, lo que solo tiene sentido, en princi-pio, para matrices cuadradas. Por ello, si Soleq contu-viera >, Cremona resolvera solo aquellas cerchas

con igual numero de incognitas (barras mas bielas) queecuaciones, y solamente estas. Pero si en Soleq apare-ce LSQ en vez de / (como es el caso), Cremona sugerirauna solucion si es que existe alguna, cualquiera que seala relacion entre el numero de barras y bielas, de unaparte, y el de nudos por la otra. Lo que cambia aques simplemente el metodo de resolucion de ecuaciones.LSQ es el acronimo del ingles least squares, es decir,mnimos cuadrados; y con este metodo es posible ob-tener una solucion cuando 2n = e+v.

Por ejemplo, que ocurre si el apoyo derecho se sus-tituye por una articulacion como en la Figura 3? Losvnculos cambian, tenemos ahora cuatro bielas, es de-cir, la lista:

7Justamente lo contrario de lo que se entiende por multipli-car vectores (producto vectorial): en nuestro contexto de c alculo,objetos como las longitudes unas veces viene bien que sean vec-tores, pero en otras, es mas adecuado considerarlos como listasde numeros.

8/8/2019 hpcremona

5/7

100 kN

100 kNR1

R2

R3

R4

1

2 3 4

Figura 3.

{ { 2 (0;-1) }{ 2 (1;0) }{ 4 (0;-1)

}{ 4 (-1;0) } }

almacenada, por ejemplo, en V1. Ahora, la secuenciaP B V1 Q CREM arroja en la pila nuevas reac-

ciones y solicitaciones. Las reacciones son R1 = R3 =50kN y R2 = R4 = 33 kN. Las solicitaciones calcu-ladas son [-83 100 -83 33 33]. Puede comprobarseque globalmente o nudo a nudo hay equilibrio entretodas las fuerzas implicadas, pero tiene sentido fsico?Un teorico purista contestara no segun el siguienteargumento:

Puesto que en los nudos 2 y 4 hay articulacio-nes fijas, la distancia entre estos nudos es fijaigualmente y, en consecuencia, las barras 23y 34 no pueden experimentar alargamiento,ni deformacion ni tension ni, por tanto, solici-tacion: las solicitaciones facilitadas por Cre-mona son incorrectas, pues en esas barras esfacil ver que la solicitacion es nula.

Sin embargo, un@ disenador@ mas interesad@ en elproyecto de estructuras y en predecir como se compor-taran, podra optar por un argumento diferente:

Supongamos que el coeficiente de seguridaden el diseno de cada pieza sea = 2. Si di-seno las estructuras que soportan la cercha,a traves de las articulaciones 2 y 4, para re-sistir una fuerza horizontal de 33 = 66kNy con suficiente ductilidad, cuando la cargaen la cercha sea 100 = 200 kN, estas estruc-turas comenzaran a ceder (si no lo han he-cho antes), aumentando la distancia entre 2 y4. En consecuencia, las barras 23 y 34 co-menzaran a alargarse y entraran en traccion apartir de ese momento. Cremona me informa

que la maxima traccion en esas barras sera de33 = 66kN. Y de hecho una vez que dise-ne las barras puedo incluso calcular el alarga-miento experimentado entre 2 y 4 hasta que lacedencia de las articulaciones se detiene, gra-cias a la entrada en tension de los tirantes.

Existen las articulaciones que nunca se mueven? Si surespuesta es s es usted un purista y Cremona, conla actual version de Soleq, no le sera de utilidad: sus-tituya LSQ por /. Cuando intente resolver esta ultimacercha recibira el mensaje:8

/ Error:

INVALID DIMENSION

Si su respuesta es no, bienvenid@ al grupo! Este pro-grama, en su actual version, podra darle en ocasionespistas interesantes.

La cercha de la Figura 3 sera calificada por muchoscomo hiperestatica indicando con ello que solo la consi-deracion de las deformaciones de las barras con seccio-nes previamente definidas, permite calcular con rigory precisison las solicitaciones. Esto es cierto bajo lacarga de servicio de 100kN, pues el comportamientodebe ser elastico. Pero ya hemos visto que con solo lasecuaciones de la estatica, es posible obtener solicitacio-

nes en equilibrio que seran proporcionales a las de larotura de la cercha si y solo si ajustamos el diseno aesas solicitaciones. A fin de cuentas, Cremona permiteobtener un diseno seguro (pero no sabemos si rgido yhabra que comprobar que ductil) incluso en esta situa-cion hiperestatica. Todo depende de si nos permitimospensar isostaticamente.

Cerchas funiculares

En esencia, la suma R2 + N4 sobre el nudo 2 debeequilibrar la componente horizontal de N1, que es (ba-

jo 100kN) de 66 kN. Cremona sugiere dividir el trabajoentre la articulacion y el tirante a partes iguales. Hayinfinitas alternativas a ese reparto, pero hay algunaalternativa mejor a la de Cremona? Depende de cuantocuesten el tirante y la articulacion; no estoy pensan-do en costes monetarios muy poco interesantes. Porejemplo, si la estructura que ha de soportar la cerchaes el suelo y tiene suficiente resistencia, podemos enco-mendarle todo el trabajo: R2 = 66 y N4 = 0 sera unasolucion posiblemente mejor. Este diseno consiste, dehecho, en suprimir las dos barras del tirante, vease laFigura 4. Ahora en la lista de barras solo quedan tres:

{ [1 2][1 3][1 4] }. Con esta lista almacenada enBF, la secuencia P BF V1 Q CREM arroja en lapila nuevas reacciones y solicitaciones:

[-83 100 -83]

[50 67 50 67]

Es facil en este caso comprobar a mano que estos son losresultados correctos para esta peculiar estructura. Setrata de hecho, de una estructura funicular: si ademasde 100kN verticales hubiera algo de carga horizontal enel nudo 3, el pendolon 13 adaptara su posicion al igual

que lo hace un hilo. De hecho, un purista dira que estacercha ni siquiera es estructura, pero entonces que sonlos cables o barras de los que cuelga una lampara?

8Es lo habitual en los programas que se intercambian pordinero: pruebe a que alguno de ellos calcule (sin recurrir a sorti-legios) un cable del que cuelga una l ampara. . .

8/8/2019 hpcremona

6/7

100 kN

100 kNR1

R2

R3

R4

1

2 3 4

Figura 4.

Diseno de cerchas

En resumen, Cremona ofrece para cualquier estruc-tura de barras (de codales y tirantes), un conjunto desolicitaciones y reacciones en equilibrio estatico con lasacciones, si es que existe algun conjunto semejante. Lasolucion puede ser estable o inestable y averiguarlo estarea de quien usa el programa (piense en la cerchasimetrica respecto al eje horizontal de la ultima cer-cha considerada, pero con los 100kN dirigidos todavahacia abajo). Ademas, con areas adecuadas (iguales omayores que las que resultan de Ni/e), ese conjunto

sera proporcional a las solicitaciones en la rotura de lacercha; la razon de la proporcion entre estas y aquellassera justamente el coeficiente de seguridad . Faltaracomprobar la ductilidad y la rigidez de la estructura.Puesto que dado un esquema Cremona permite obtenersecciones bastante razonables para las barras, Cremonaes en realidad un programa de dimensionado de cer-chas.

Desde un punto de vista mas abstracto, Cremona re-suelve con su mejor criterio cualquier sistema de ecua-ciones con fuerzas que este asociado a una cercha. Tec-nicamente, su mejor criterio consiste en obtener el vec-

tor de incognitas con menor norma euleriana y que sa-tisface el sistema, si es que existe alguno. . . (Que sig-nifica norma euleriana? En esencia el modulo de unvector, pero consulte un buen libro de matematicas...)Desafortunadamente, la mnima norma euleriana delvector de solicitaciones y reacciones no puede relacio-narse positivamente con ninguna propiedad estructuralque sea deseable minimar. Por tanto, los disenos pro-puestos por el programa no seran en general optimos.

Finalmente. . .

STOPes el ultimo programa del que tengo que de-cir algo. Se trata de una version generalizada de las

rutinas de aritmetica en memoria (inspeccione

VAR ARITH ): toma de la pila un valor, una expre-sion algebraica y un programa y almacena en expresionalgebraica el resultado de evaluar el programa sobre ex-

presion algebraica y valor. Si piensa hacer uso de STOPen sus propios programas, note que la expresion alge-braica debe corresponder a una ubicacion precisa de lamemoria (es decir, ese misma expresion algebraica debe

ser admitida por STO como segundo argumento sinproducir error), y que programa debe tomar dos argu-mentos de la pila y devolver un valor (mire los ejemplosde uso en Cremona).

STOP, OpVec, MulVec, DivVec y Soleq son programasde utilidad general y su lugar es HOME o algun directo-rio especfico desde donde cuelguen los programas queoperan con vectores y matrices (y desde donde podranestos ultimos invocar a aquellos).

Para observar que hace cualquiera de los programas,y por tanto entenderlos, puede usar PRG NXT RUN

DBUG . Conviene que lo haga: toda la responsabilidadal usar un programa es suya. . .

8/8/2019 hpcremona

7/7

Anejo

Codigo fuente de los programasSitue la ruta actual en el directorio elegido para los programas y ejecute las siguientes instrucciones.

STOP

>>

NEXT q a Soleq OBJ DROP nv

+ ARRY r

>>>

>>

>>

>>

Cremona

ENTER

STO

>

STOPENTER

STO

>

Soleq

ENTER

STO

>

>>

LongitudesENTER

STO

>>

OpVec

ENTER

STO

OpVec

>>

MulVec

ENTER

STO

OpVec

>>

DivVec

ENTERSTO

{ Cremona

Longitudes MulVec

DivVec { "DOT"

> } { "OBJ "

> }

a~nada mas informacion a su gusto

}

CST

ENTERSTO

Copyleft

c

200

4,

2007,

Vazquez.

Printedwithfreesoftware:GNU/Linux

/emacs/LATEX2

/Postscript.