¿Cómo incluir cargas variables en el espacio?

En el siguiente ejemplo modelaremos una viga de acero de 5 m, que arranca en z=0, en donde
está empotrada, y se extiende hasta z=5 m. Su sección es cuadrada de canto 0.1 m. Se aplica una
carga Fy (N/m) variable lineal con la altura z, siendo Fy(z=0)=10000, y Fy(z=5)=0.

Obviaremos el procedimiento de construcción de la geometría y la malla al ser relativamente

En el módulo Geometry, debemos crear dos vértices, uno en (0,0,0), otro en (0,0,5), y luego la
línea que los une. A continuación, creamos el grupo de línea que llamamos “viga” (para asignar
las propiedades de la viga, así como la carga), y luego el grupo de vértice “empo”, para poder
aplicar la condición de empotramiento.

1
Luego, al pasar al módulo Mesh, generaremos la malla 1D, dividiendo la línea en 5 segmentos
iguales (5 elementos). Debemos también generar los grupos desde la geometría (Groups from
Geometry).

Luego, debemos pasar al módulo AsterStudy. El procedimiento para la creación del caso es
similar a lo que hemos hecho para los otros casos, excepto en la función a utilizar para aplicar la
carga variable. Recordemos:

Clicar con botón derecho en “CurrentCase”, y clicar en “Add Stage”.

Clicar con botón derecho en “Stage_1”, y clicar en Mesh > Read a mesh. Seleccionar la malla y
clicar OK.

Clicar con botón derecho en “Stage_1”, y clicar en Model definition > Assign finite element.
Seleccionar “Finite element assignment”, aplicar en “Everywhere”, con “Phenomenon”:
“Mechanic”, y “Modelisation”: “POU_D_E” (viga Euler-Bernoulli). Clicar OK, y luego OK.

Clicar con botón derecho en “Stage_1”, y clicar en Model definition > AFFE_CARA_ELEM.
Seleccionar “POUTRE”, añadir Item, y clicar en Edit. Seleccionar “Group of element”: “viga”. Tipo
de sección: “Section”: “Rectangle”. Luego en “CARA”, seleccionar “H”, y en “Value”:0.1. Clicar
OK. Luego activar “Orientation”, añadir Item, y clicar en Edit. Seleccionar “Group of element”:
“viga”. En “CARA” : VECT_Y (vector director del eje y de la sección), y en “Value” añadir hasta
tener tres casillas, e insertar: 0,1,0. Clicar OK. Clicar OK.

Clicar con botón derecho en “Stage_1”, y clicar en Material > Define a material. Activar “Linear
Isotropic”, e introducir Young’s modulus 210E9, y Poisson’s ratio=0.3. Clicar OK. Clicar OK.

Clicar con botón derecho en “Stage_1”, y clicar en Material > Assign a material. Activar “Model”.
Añadir en Material assignment un ítem, clicar Edit. Seleccionar “Everywhere”, y luego en
“Material” añadir un ítem (debe aparecer el único material por ahora definido que se llama
“mater”. Clicar OK. Clicar OK.

2
Ahora viene el primer paso diferente. Ahora vamos a crear una función. Clicar con botón
derecho en “Stage_1”, y clicar en Functions and LIsts > Define function.

3
En el menú de dicha función, en primer lugar, cambiamos el nombre del objeto en donde se
guarda, y escribimos “funcion1”. Seguidamente vamos al desplegable “Parameter name”, y
seleccionamos “Z”. Es decir, vamos a generar una función cuya variable independiente es la Z. A
continuación, seleccionamos “Coordinates” en el modo de definir la función, y por último
clicamos en “Edit” a su derecha.

En la nueva pantalla que aparece, añadimos una nueva fila clicando en “Append row”

4
Rellenamos la tabla que define la función, donde en z=0, tenemos que el valor de la función
10000 (Fz(0)=10000), y en z=5, el valor de la función es 0 (Fz(5)=0). Realmente lo que estamos
haciendo es generar una función a través de puntos de interpolación. Clicamos OK. Y luego OK
de nuevo.

Tras definir la función. Ahora vamos a definir la carga variable. Para ello clicamos botón
derecho en “Stage_1”, luego “BC and Load” > “AFFE_CHAR_MECA_F”.

5
Esta función es muy similar a la que aparece al clicar en “Assign mechanical load” (cuyo nombre
en francés es “AFFE_CHAR_MECA”, sólo que ahora en vez de asignar valores discretos, podemos
asignar valores dados por una función. En primer lugar, lo que haremos es darle un nombre
significativo, para ello ponemos en “Name” el nombre “varload”. Después, como la fuerza que
vamos a aplicar es una carga distribuida lineal (fuerza por unidad de longitud), debemos activar
“FORCE_POUTRE” (carga distribuida en vigas), luego añadir un ítem, y luego clicar en editar:

En la nueva pantalla que nos aparece, ahora activamos “Group of element” y clicamos en Edit:

6
Seleccionamos el grupo que antes llamamos “viga”, y clicamos OK.

Y luego, activamos la componente FY, y ya nos debería aparece seleccionado “funcion1”.

7
Clicamos OK, y luego OK. Si tuviésemos otros tipos de cargas distribuidas no constantes
podríamos seguir completando la lista de cargas en este menú. Si la fuerza fuese una presión,
fuerza por unidad de superficie y en dirección de la normal, pues tendríamos que activar
“PRES_REP”. O si fuerza otro tipo de carga, pues la metodología sería similar. Recordemos que
la definición de todas estar cargas se puede encontrar en la documentación de CodeAster. En
este caso, en inglés en [Link] y
en francés en [Link]

Para otras cargas que fuesen constantes, tendríamos que realizar lo que veníamos haciendo
hasta ahora en prácticas previas.

Para incluir la condición de empotramiento en la base, tendremos que clicar en “Stage_1” > BC
and Load > Assign mechanical load. Luego seleccionar “Enforce DOF”, añadir ítem, clicar en
Editar, seleccionar en “Group of node” el grupo “empo”, y poner DX=0, DY=0, DZ=0,
DRX=0,DRY=0,DRZ=0. Clicar OK, y luego OK. El nombre por defecto que nos aparece es “load2”,
pero se podría poner cualquier otro.

Una vez hecho esto, nos aparecen dos entradas en “BC and load”, una cuyo nombre es “varload”
(en donde se está aplicando la carga variable), y otra se llama “load2” (en donde se aplica la
única condición de contorno constante).

Clicar con botón derecho en “Stage_1”, y clicar en Analysis > Statics > Static mechanical analysis.
Ahora operamos como ya hemos hecho otras veces, excepto que ahora hay que añadir dos ítems
en “Loads”, en donde deberemos seleccionar por un lado la carga constante “load2”, y por otro
lado la carga variable “varload”.

8
Si quisiéramos obtener resultados en tensiones, etc., deberemos añadir cálculos posteriores al
análisis a través de “Post processing”, como ya vimos.

En este caso, para terminar, sencillamente vamos a exportar los resultados en desplazamientos
“DEPL” para visualizarlos en ParaVis. Para ello, recordemos que podemos ir directamente a clicar
botón derecho en “Stage_1”, y luego clicar en Output > Set output results, y completar como ya
hemos visto en otras prácticas.

Por último, resolvemos. Nos debería aparecer en verde, con lo cual a menos ha resuelto:

9
Para verificar la solución numérica, en este caso particular existe solución analítica para el
desplazamiento del extremo libre (z=5). Si buscamos en cualquier buscador “prontuario viga
carga triangular”, o mejor en inglés “beam deflection triangular load” podremos encontrar
referencias para dicha solución (por ejemplo, en [Link]
deflection-tables). Que resulta ser:

𝑞𝑞𝐿𝐿4
𝑤𝑤𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒 =
30𝐸𝐸𝐸𝐸
Donde q es el valor de la carga triangular (el valor máximo), L es la longitud de la viga, E el módulo
de elasticidad del material, e I es la inercia de la sección. En nuestro caso:

10000 · 54
𝑤𝑤𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒 = = 0.119 m
1
30 · 210 · 109 · · 0.13 · 0.1
12
El valor de esta deflexión la podemos observar al plotear la deformada de la viga (como la carga
está aplicada en dirección y, el único desplazamiento no nulo está en esta dirección):

Y vemos que efectivamente, el valor de esta deflexión en el extremos libre coincide con la
analítica.

10