Uso de ansys en soporte de pendientes
Aquí hay un ejemplo de soporte de pendiente.
Ejecútelo paso a paso en ANSYS y luego compárelo y comprenda.
Ingrese al preprocesador y defina el tipo de unidad requerida. Para este problema, el modelo solo contiene un elemento, el Plano 42. Dado que este cálculo es para un problema de deformación plana bidimensional, el KEYOPT(3) del Plano 42 debe establecerse en 2. El flujo de comando correspondiente es el siguiente:
/prep7
et, 1, plane42
KEYOPT, 1, 3, 2! >
(3) Definir los parámetros del material en el modelo. Hay 7 tipos de materiales en el modelo, incluidos los parámetros del material del lecho rocoso y del cuerpo deslizante en diferentes saturaciones. El comando para definir los parámetros del material es el siguiente:
mp, ex, 1, 20e9! Modelo elástico del macizo rocoso
mp, prxy, 1, 0.22! p>
mp, dens, 1, 2600! Densidad del macizo rocoso
mp, ex, 2, 3e9!
mp,dens,2,2140
mp,ex,3,2.7e9! La saturación es 35
mp,prxy,3,0.35
mp, dens, 3, 2180
mp, ex, 4, 2.4e9! La saturación es 50
mp, prxy, 4, 0.35
mp, dens, 4, 2210
mp, ex, 5, 2.1e9! La saturación es 65
mp, prxy, 5, 0.35
mp, dens, 5, 2240
mp, ex, 6, 1.8e9! La saturación es 80
mp, prxy, 6, 0.35
mp, dens, 6 , 2270
mp, ex, 7, 1.5e9! Módulo elástico con saturación de 100
mp, prxy, 7, 0,35! p> mp, dens, 7, 2300! Densidad con saturación de 100
guardar! Base de datos de almacenamiento
(4) Establecer un modelo geométrico y dividir las unidades. Dado que la geometría del modelo es relativamente compleja, se adopta un método de modelado ascendente, es decir, primero se generan puntos, luego se generan superficies y finalmente se dividen las unidades.
k, 1, 0, 0, 0!Genera puntos clave según coordenadas
k, 2, 0, 505, 0
k, 3, 562 ,505,0
k,,621,459,0
k,,658,433,0
k,,693,393,0
k,, 802, 313, 0
k,, 850, 303, 0
k,, 892, 273, 0
k,,913,253 ,0
k,,930,243,0
k,,1034,233,0
k,,1186 ,228,0
k,,1216,223,0
k,,1232,208,0
*do,i,1,14! Crea líneas rectas de forma circular, * **Dibujó 14 líneas
l, i, i 1
*enddo
k, 16, 660, 387, 0 ¡Recrea puntos clave
k,,770,248,0
k,,827,230,0
k,,888,217 ,0
k,,930,216, 0
k,,1034,213,0
k,,1186,209,0
l, 4, 16! sobre los puntos clave
l, 16, 17
l, 17, 18
l, 18 , 19
l,19 ,20
l,20,21
l,21,22
l,22,15
nummrg, ¡todos! todos los elementos
numcmp, all! Comprime todos los números de elementos
k, 23, 1243, 198, 0! Crea puntos clave nuevamente
k,, 612, 380, 0
k,, 745, 214, 0
k,, 884, 181, 0
k,, 1242, 174, 0
l, 3, 24! Luego crea una línea recta basada en los puntos clave
l, 24, 25
l, 25, 26
l, 26, 27
l, 27, 23
l, 23, 15
nummrg, all! p> numcmp, todo! Comprime todos los números de elementos
k, 28, 1610, 198, 0 ¡Otra vez!
Crear puntos clave
k,,1610,0,0
k,,1243,0,0
k,,884,0,0 p>
k,, 745, 0, 0
k,, 612, 0, 0
l, 6, 16! los puntos clave
l, 7, 17
l, 8, 18
l, 9, 19
l, 11, 20
l, 13, 22
l, 16, 24
l, 17, 25
l, 19, 26
l, 24, 33
l, 25, 32
l, 26, 31
l, 27, 30
l, 1, 33
l, 30, 31
l, 31, 32
l, 32, 33
l, 23, 28
l, 28, 29
l, 29, 30
guardar
al, 4, 5, 15, 29! Crea polígonos basados en líneas
al, 29, 6, 16, 30
al, 7, 17, 30, 31
al, 8, 18, 31 , 32
al, 9, 10, 19, 32, 33
al, 11, 12, 20, 21, 33, 34
al, 13 , 14, 22, 34
al, 3, 15, 35, 23
al, 35, 16, 24, 36
al, 36, 17 , 18, 37, 25
al, 19, 20, 21, 22, 26, 27, 28, 37
al, 1, 2, 23, 38, 42
p>al, 38, 39, 24, 45
al, 39, 40, 25, 44
al, 40, 41, 26, 43
al, 46, 47, 48, 27, 41
allsel, all! Seleccionar todos los elementos
tipo, 1! tipo de unidad a dividir
mat, 1! Seleccione el modelo de material a dividir en unidades
mshape, 1, 2d Establezca en unidad triangular plana
ESIZE, 10, 0, ! Establecer el tamaño de la unidad
amesh, 1, 11, 1! lado 11
ESIZE,,, !
Establezca el tamaño de la unidad en el valor predeterminado
allsel, all! Seleccione todos los elementos
amesh, 12, 16, 1 ¡Divida las caras 12 a 16
guardar!
p>
terminar
!La malla de elementos finitos dividida se muestra en la Figura 9-2.
9.4 Carga y resolución
(1) Aplicar condiciones de contorno y realizar ajustes de solución.
/solu !Ingrese el solucionador
eplot !Mostrar diagrama de elementos
guardar
nsubst, 20 !Especificar subpaso de carga
p>
nlgeom, activado! Especificar deformación grande
nropt, completo! Establecer iteración newton-lapnace
habilitar predictor de solución no lineal
p>lnsrch, activado! Habilitar búsqueda lineal
outres, todo, ¡último! Muestra solo el resultado del cálculo del último paso
! , s , loc, x, 0! Selecciona el nodo
d, all, ux! Restringe el desplazamiento del nodo en la dirección X
nsel, s, loc, x, 1610! Seleccione el nodo nuevamente
d, all, ux! Restringe el desplazamiento del nodo en la dirección X
nsel, s, loc, y, 0!
d, all, all! Restringe el desplazamiento del nodo en las direcciones X e Y
allsel, all! Selecciona todos los elementos
acel, 0, 10, 0. !Aplica aceleración de gravedad
allsel, all! Selecciona todo Elemento
guardar
finalizar
(2) Resuelve el estrés de autogravedad campo
/solu! Ingrese el solucionador
tiempo, 1! Primer paso de la solución
asel, s,,, 1, 7, 1! /p>
esla, s
!Seleccione elementos de superficie
mchg, 2, todos! Cambie los parámetros de esta parte del elemento del cuerpo deslizante
allsel, todos ¡Seleccione todos los elementos
resolver! ! Resolver
guardar ! Guardar la base de datos
finalizar ! Volver al menú principal anterior o salir del solucionador
(3) Resolver la suma de tensiones bajo diferentes saturaciones Campo de desplazamiento
!******Cuando la saturación es del 35%************
/solu !Ingresa al solucionador
antype,, resto! Reiniciar
tiempo, 2! Establecer el segundo paso de solución
asel, s,,, 1, 7, 1! p>
esla, s! Seleccionar elementos de superficie
mchg, 3, all! Cambiar los parámetros de material de los elementos seleccionados
allsel, all! >
resolver !Resolver
guardar !Guardar la base de datos
finalizar !Volver al menú principal anterior o salir del solucionador
!**** ** Cuando la saturación es del 50%********
! Su anotación china es la misma que cuando la saturación es del 35%
/solu
antype, , descanso
tiempo, 3
asel, s,,, 1, 7, 1
esla, s
mchg, 4 , todos
allsel, todos
resolver
guardar
terminar
!* ***** Cuando la saturación es del 65%********
/solu
antipo, tiempo de descanso
, 4
asel, s,,, 1, 7, 1
esla, s
mchg, 5, todos
allsel, todos
resolver
guardar
terminar
!******Cuando la saturación sea del 80%******** ****
>
/solu
antipo,, descanso
tiempo, 5
asel, s,,, 1, 7, 1
esla, s
mchg, 6, all
allsel, all
resolver
guardar
terminar
!******Cuando la saturación sea 100%************
/solu
antype ,,tiempo de descanso
, 6
asel, s,,, 1, 7, 1
esla, s
mchg, 7, todos
allsel, todos
resolver
guardar
terminar
9.5 Análisis de los resultados del cálculo
9.5.1 Estabilidad de la pendiente en el estado inicial
Una vez completada la solución, lea los resultados del cálculo del primer paso y examine los diagramas de nubes de tensión, desplazamiento y tensión cortante. y deformación cortante del talud bajo la acción de su propio peso. Todos los siguientes comandos se ejecutan después de ingresar al posprocesador/POST1, por lo que los lectores deben ingresar al posprocesador general antes de ejecutar los siguientes comandos. El comando para dibujar estas imágenes de nubes es:
/POST1
SET, 1
PLNSOL, S, Y, 0, 1. en el estado inicial
PLNSOL, U, Y, 0, 1! Campo de desplazamiento en el estado inicial
PLNSOL, S, XY, 0, 1!
PLNSOL, EPTO, La estabilidad del talud se ha resuelto. Debido al espacio limitado, a continuación solo se analiza la estabilidad del talud en el estado completamente saturado. El contenido y el método del análisis se pueden aplicar al. Análisis de estabilidad de otras saturaciones.
Dado que es necesario calcular los cambios en varias cantidades físicas causadas por la saturación de agua, aún es necesario usar el comando de combinación de casos de carga en la sección 6.4.4 anterior y modificarlo en consecuencia, como se muestra a continuación.
!*******loadcase.mac*********
/POST1
LCDEF, 1, 1, , !will El primer paso de cálculo se define como caso de carga 1
LCDEF, 2, 6, , ! El sexto paso de cálculo (la saturación es 100) se define como caso de carga 2
LCASE, 2 , ! Leer en el caso de carga 2
LCOPER, SUB, 1, , , ! Reste el caso de carga 1 del caso de carga 2
Luego realice el diagrama incremental de fuerza y cortante de desplazamiento y corte. , el flujo de comando es:
PLNSOL, U, Y, 0, 1! Campo de desplazamiento en el estado inicial
PLNSOL, S, XY, 0 , 1! estado inicial
PLNSOL, EPTO, XY, 0, 1! Deformación por corte en el estado inicial
A continuación se utiliza el criterio de Mohr-Column para el estado saturado. Se calcula la distribución de la zona plástica. Los parámetros de propiedad física en el estado saturado son: C=50kpa, φ=30?, y el flujo de comando es el siguiente:
!********Mohr.mac* **. *****
/post1
asel, s,,, 1, 7, 1! Seleccionar unidad
esla, s
etable, S1_1, S, 1! Defina las tensiones principales 1.ª y 3.ª de la unidad seleccionada como tablas de unidades
etable, S3_1, S, 3
sadd ,S1,S3_1. ,,-1.00,1
sadd,S3,S1_1,,-1.00,1
sadd,smf_over,S1,S3,1.00,-1.00
cos_fai=cos(30/180*3.1415926) ! Calcula el coseno del ángulo de fricción interno
sin_fai=sin(30/180*3.1415926)
constante=2*50000 *cos_fai
sadd, S1_S3, S1, S3, 1, 1
sep, S1__S3, S1_S3,,1
sadd, smf_down, S1__S3,,sin_fai,, constante
sexp, s_m_f, smf_over, smf_down, 1, -1
esel, all
/contour
pletab, s_m_f, avg! Dibuja un diagrama de nubes del elemento de rendimiento