[derik]

Здраствуйте уважаемые форумчане! Просьба тех кто считает в Ансисе подсказать - где у меня ошибка в схеме! Хелп перечитал, он мне уже скоро сниться будет :-( помогите пожалуйста!
вот код:

Код:

[Выделить все]

/prep7
a=4
h=2.5
b=.4
t=.2
q=100
/view,1,1,1,1
/ang,1
/rep,fast
blc5,a/2,a/2,a,a,
k,,b/2,b/2,0
k,,a-b/2,b/2,0
k,,b/2,a-b/2,0
k,,a-b/2,a-b/2,0
k,,b/2,b/2,-h
k,,a-b/2,b/2,-h
k,,b/2,a-b/2,-h
k,,a-b/2,a-b/2,-h
lstr,5,9
lstr,6,10
lstr,7,11
lstr,8,12
nummrg,all
boptn,,keep
boptn,keep,1
boptn,nwar,0
boptn,vers,rv52
btol,1e-005,
lina,5,1
lina,6,1
lina,8,1
lina,7,1
mptemp,,,,,,,,
mptemp,1,0
mpdata,ex,1,,3e8
mpdata,prxy,1,,.3
mpdata,dens,1,,2500
et,1,shell63
r,1,t,,,,,,
et,2,beam4
keyopt,2,9,9
r,2,b*b,(b*b**3)/12,(b*b**3)/12,b,b,
lsel,s,,,5
lsel,a,,,6
lsel,a,,,7
lsel,a,,,8
latt,1,2,2,
asel,s,,,1
aatt,1,1,1
allsel,all
lesize,5,,,b,1
lesize,6,,,b,1
lesize,7,,,b,1
lesize,8,,,b,1
lmesh,5,6,
lmesh,7,8,
aesize,1,b/2,
amesh,1
LGEN,3,all, , , , ,h, ,0   
AGEN,3,all, , , , ,h, ,0
nummrg,all
allsel,all
/eshape,1
/replot
finish
/sol
nsel,s,loc,x,b/2
nsel,r,loc,y,b/2
nsel,r,loc,z,-h
D,all, , , , , ,ALL, , , , , 
nsel,s,loc,x,b/2
nsel,r,loc,y,a-b/2
nsel,r,loc,z,-h
D,all, , , , , ,ALL, , , , , 
nsel,s,loc,x,a-b/2
nsel,r,loc,y,b/2
nsel,r,loc,z,-h
D,all, , , , , ,ALL, , , , , 
nsel,s,loc,x,a-b/2
nsel,r,loc,y,a-b/2
nsel,r,loc,z,-h
D,all, , , , , ,ALL, , , , ,
allsel,all
acel,0,0,9.8,
solve
finish
/post1
set,first

aka raduntsev-nias, Kostya_PC - спасибо, подумаю над выбором.

[Cfytrr]

aka raduntsev-nias,
Насколько я понял, в этой модели бетон работает упруго даже в растянутой зоне

Потом, автор предлагает натянуть арматуру на соответствующее усилие и вполне закономерно считает, что после этого бетон действительно "трещать" не должен.
А нет ли у Вас подобных отчетов для простого железобетона, где бетон вынужден в сжатой и растянутой зоне работать по-разному?

[aka raduntsev-nias]

Cfytrr,
Вас интересует процесс развития трещин в растянутых зонах бетона? Те сравнение аналитического метода расчета жб балки с расчетом в мкэ программах и сравнимость с результатами эксперимента? те увидеть виртуальное образование трещин, разрушение балки и тд.
Такого отчета у меня нет. Поискал в сети, нашел наверное похожий
Тоже бы с удовольствием посмотрел на технику моделирования таких процессов.
Нашел. Сильно не вникал. Может не совсем то, что Вам нужно. )
Искал целенаправленно abaqus/ То же можно и в ansys найти думаю.

[derik]

Код:

[Выделить все]

fini
/clear,nostart
/outp,output,txt,C:\ansys\output
/cwd,'C:\ansys\out'
/units,mks

/prep7

*set,a,4
*set,b,.4
*set,c,.2
*set,q,5e4
*set,p,2e5

et,1,185
mptemp,,,,,,,
mptemp,1,0
mpdata,ex,1,,2.04e5
mpdata,prxy,1,,.3
tb,biso,1,1,2,
tbtemp,0
tbdata,,.84e5,3044

blc5,,,a,b,c
csys,5
vgen,2,all,,,,90,,,0
csys,0
vgen,,2,,,-b/2,,c/2,,,1
vovlap,all

mshkey,0
mshape,2,3d
esize,b/4
vsweep,all
nummrg,all

asel,s,,,5,6
asel,a,,,11,12
da,all,all,
sfa,34,1,press,q
alls
fini

/sol

/gst,1
nldiag,nrre,on
antype,0
nlgeom,1
nsubst,250,5500,75
outres,erase
outres,all,all
time,1
solve
fini

/post1

set,first
plnsol,uy,sum,0,1

Добрый вечер!
на критику первую нелинейную модель выставляю

[aka raduntsev-nias]

Сообщение от 6pa4o первую нелинейную модель выставляю

круто! у меня в вб так пока и не получилось с нелином

[st2008]

Сообщение от 6pa4o на критику первую нелинейную модель выставляю

Выдает такое сообщение...

One or more elements have become highly distorted. Excessive distortion of elements is usually a symptom indicating the need for corrective action elsewhere. Try incrementing the load more slowly (increase the number of substeps or decrease the time step size). You may need to improve your mesh to obtain elements with better aspect ratios. Also consider the behavior of materials, contact pairs, and/or constraint equations. If this message appears in the first iteration of first substep, be sure to perform element shape checking.

[derik]

вообщем нужно увеличить число подшагов или понизить временной шаг. с этой командой у меня пока туго не могу найти как подбирать эти самые шаги, их количество и мин. и макс. число шагов

[aka raduntsev-nias]

Сообщение от 6pa4o найти как подбирать эти самые шаги, их количество и мин. и макс. число шагов

или может с материалами еще что то не то

[derik]

материал брал из верификационного примера...то кол-во шагов которое там же указано, эту ошибку выдавало раз 5.

[soppr]

Маленькая контактная задачка для понимания контакта в ansys
Два кубика, один скользит по другому. Верхний кубик имеет модуль упругости 2000 и коэф. Пуассона 0.3. Нижний модуль 2e12 и Пуассона 0.3.
Верхний кубик нагружен на верхней площади распр. нагрузкой 2 и на боковой грани 0.5. Нижний кубик заделан со всех сторон (кроме контактной). Задано трение 0.3, отношение статического к динамическому трению - 1.2. Трение покоя (cohesion) - 0.4.

Код:

[Выделить все]

/COM,  Structural
!*  
/PREP7  
!*  
ET,1,SOLID186
!*  
ET,2,TARGE170
!*  
ET,3,CONTA174
!*  
!*  
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,EX,1,,2000   
MPDATA,PRXY,1,,0.3  
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,EX,2,,2e12   
MPDATA,PRXY,2,,0.3  
BLOCK,0,1,0,1,0,0.5,
BLOCK,0,1,0,1,0.5,1,
CM,_Y,VOLU  
VSEL, , , ,       1 
CM,_Y1,VOLU 
CMSEL,S,_Y  
!*  
CMSEL,S,_Y1 
VATT,       1, ,   1,       0   
CMSEL,S,_Y  
CMDELE,_Y   
CMDELE,_Y1  
!*  
CM,_Y,VOLU  
VSEL, , , ,       2 
CM,_Y1,VOLU 
CMSEL,S,_Y  
!*  
CMSEL,S,_Y1 
VATT,       2, ,   1,       0   
CMSEL,S,_Y  
CMDELE,_Y   
CMDELE,_Y1  
!*  
ESIZE,0.1,0,
FLST,5,2,6,ORDE,2
FITEM,5,1
FITEM,5,-2  
CM,_Y,VOLU  
VSEL, , , ,P51X 
CM,_Y1,VOLU 
CHKMSH,'VOLU'   
CMSEL,S,_Y  
!*  
MSHAPE,0,3d
MSHKEY,1
VMESH,_Y1   
MSHKEY,0
!*  
CMDELE,_Y   
CMDELE,_Y1  
CMDELE,_Y2  
!*  

/COM, CONTACT PAIR CREATION - START 
CM,_NODECM,NODE 
CM,_ELEMCM,ELEM 
CM,_KPCM,KP 
CM,_LINECM,LINE 
CM,_AREACM,AREA 
CM,_VOLUCM,VOLU 
/GSAV,cwz,gsav,,temp
MP,MU,1,0.3
MAT,1
MP,EMIS,1,7.88860905221e-031
R,3
REAL,3
ET,4,170
ET,5,174
R,3,,,1.0,0.1,0,
RMORE,,,1.0E20,0.0,1.0, 
RMORE,0.4,0,1.0,,1.0,0.5
RMORE,0,1.0,1.2,0.0,,1.0
KEYOPT,5,4,0
KEYOPT,5,5,0
KEYOPT,5,7,0
KEYOPT,5,8,0
KEYOPT,5,9,0
KEYOPT,5,10,2
KEYOPT,5,11,0
KEYOPT,5,12,0
KEYOPT,5,2,0
KEYOPT,4,5,0
! Generate the target surface   
ASEL,S,,,2  
CM,_TARGET,AREA 
TYPE,4
NSLA,S,1
ESLN,S,0
ESLL,U
ESEL,U,ENAME,,188,189   
NSLE,A,CT2 ! CZMESH patch (fsk qt-40109 8/2008) 
ESURF   
CMSEL,S,_ELEMCM 
! Generate the contact surface  
ASEL,S,,,7  
CM,_CONTACT,AREA
TYPE,5
NSLA,S,1
ESLN,S,0
NSLE,A,CT2 ! CZMESH patch (fsk qt-40109 8/2008) 
ESURF   
ALLSEL  
ESEL,ALL
ESEL,S,TYPE,,4  
ESEL,A,TYPE,,5  
ESEL,R,REAL,,3  
/PSYMB,ESYS,1
/PNUM,TYPE,1
/NUM,1
EPLOT   
ESEL,ALL
ESEL,S,TYPE,,4  
ESEL,A,TYPE,,5  
ESEL,R,REAL,,3  
CMSEL,A,_NODECM 
CMDEL,_NODECM   
CMSEL,A,_ELEMCM 
CMDEL,_ELEMCM   
CMSEL,S,_KPCM   
CMDEL,_KPCM 
CMSEL,S,_LINECM 
CMDEL,_LINECM   
CMSEL,S,_AREACM 
CMDEL,_AREACM   
CMSEL,S,_VOLUCM 
CMDEL,_VOLUCM   
/GRES,cwz,gsav
CMDEL,_TARGET   
CMDEL,_CONTACT  
/COM, CONTACT PAIR CREATION - END   
/MREP,EPLOT

/SOL
FLST,2,5,5,ORDE,2
FITEM,2,8
FITEM,2,-12 
!*  
/GO 
DA,P51X,ALL,
FLST,2,1,5,ORDE,1
FITEM,2,1
/GO 
!*  
SFA,P51X,1,PRES,2
FLST,2,1,5,ORDE,1
FITEM,2,5
/GO 
!*  
SFA,P51X,1,PRES,0.5

Solve
/post1

Лог файл конечно кривоват)). Но работает)
Комментарии господа)

[aka raduntsev-nias]

Сообщение от soppr Комментарии господа)

решилась. но трение 0,3 - нессиметричный решатель нужно. или с меньшим решить

[soppr]

Сообщение от aka raduntsev-nias решилась. но трение 0,3 - нессиметричный решатель нужно. или с меньшим решить

От величины трения зависит решатель? Не мог бы подсказать, в каком разделе найти эту информацию - в решателях или контактной?

Возникли еще ряд вопрсов:
1. При выборе “жесткой” и “поддатливой” поверхностей руководствовался принципами, написанными в ANSYS help, те контактной поверхностью выбрал наименее жесткую (верхний куб). При этом расчет не производился. Как только поменял наименование поверхностей местами (нижний куб с большим модулем упругости стал контактным), расчет прошел. С чем связано?
2. При вводе коэффициента трения он присваивается какому то из материалов. Влияет ли на результаты задание трения разным материалам? (как я полагаю - нет)
3. Контактная жесткость. В helpANSYSсказано, что этот параметр вычисляется автоматически (как я полагаю в зависимости от модулей упругости, геометрии и других). Тем не менее есть рекомендация по начальному заданию контактной жесткости равной 1 для объемных задач. То есть его нужно подбирать вручную? Что это за параметр и как он влияет на результаты?

[aka raduntsev-nias]

сколько вопросов. да я то откуда все это знаю.
могу только посоветовать.

Сообщение от soppr 1. При выборе “жесткой” и “поддатливой” поверхностей

решение контактных задач, CADFEM

Сообщение от soppr При вводе коэффициента трения он присваивается какому то из материалов. Влияет ли на результаты задание трения разным материалам

на САПРе это обсуждали. так понял, что задавать ДВА разных коэффициента. Но не уверен.

Сообщение от soppr Контактная жесткость. В helpANSYSсказано, что этот параметр вычисляется автоматически

да, если установить опцию. имеет и другие значения. для трения можно 0,1 поставить. тут есть про это в теме.

[soppr]

1. CADFEM читал и сделал как там написано и именно так она не решалась) Поэтому попоробовал изменить CONTACT и TARGET местами.
2. В контактном мастере задается friction только для одного ID материала. Хотя если контактный мастер не использовать то можно задать каждому материалу. Поэтому и вопрос возник.

[aka raduntsev-nias]

NROP,UNSYM нужно добавить
и субстепов поболее. в контакте сетку побейте сильней.(на контакте)
еще у тебя временного интервала нет

[soppr]

raduntsev, большое спасибо за советы.

[aka raduntsev-nias]

soppr, если помогло, хорошо. не за что.
Сравнение

[derik]

погрешности бывают и у машин это абакус, Макс?

[aka raduntsev-nias]

Сообщение от 6pa4o это абакус, Макс?

да, Ярослав. если есть интерес, вот простой пример из хелпа abaqus
твоя задачка очень помогла, выручил. всю страну научил apdl у)

Код:

[Выделить все]

"""
beamExample.py

Reproduce the cantilever beam example from the
Appendix of the Getting Started with
Abaqus: Interactive Edition Manual.
"""

from abaqus import *
import testUtils
testUtils.setBackwardCompatibility()
from abaqusConstants import *

# Create a model.

myModel = mdb.Model(name='Beam')

# Create a new viewport in which to display the model
# and the results of the analysis.

myViewport = session.Viewport(name='Cantilever Beam Example',
    origin=(20, 20), width=150, height=120)
    
#-----------------------------------------------------

import part

# Create a sketch for the base feature.

mySketch = myModel.ConstrainedSketch(name='beamProfile',
    sheetSize=250.)

# Create the rectangle.

mySketch.rectangle(point1=(-100,10), point2=(100,-10))

# Create a three-dimensional, deformable part.

myBeam = myModel.Part(name='Beam', dimensionality=THREE_D,
    type=DEFORMABLE_BODY)

# Create the part's base feature by extruding the sketch 
# through a distance of 25.0.

myBeam.BaseSolidExtrude(sketch=mySketch, depth=25.0)

#-----------------------------------------------------

import material

# Create a material.

mySteel = myModel.Material(name='Steel')

# Create the elastic properties: youngsModulus is 209.E3
# and poissonsRatio is 0.3

elasticProperties = (209.E3, 0.3)
mySteel.Elastic(table=(elasticProperties, ) )

#-------------------------------------------------------

import section

# Create the solid section.

mySection = myModel.HomogeneousSolidSection(name='beamSection',
    material='Steel', thickness=1.0)

# Assign the section to the region. The region refers 
# to the single cell in this model.

region = (myBeam.cells,)
myBeam.SectionAssignment(region=region,
    sectionName='beamSection')

#-------------------------------------------------------

import assembly

# Create a part instance.

myAssembly = myModel.rootAssembly
myInstance = myAssembly.Instance(name='beamInstance',
    part=myBeam, dependent=OFF)

#-------------------------------------------------------

import step

# Create a step. The time period of the static step is 1.0, 
# and the initial incrementation is 0.1; the step is created
# after the initial step. 

myModel.StaticStep(name='beamLoad', previous='Initial',
    timePeriod=1.0, initialInc=0.1,
    description='Load the top of the beam.')

#-------------------------------------------------------

import load

# Find the end face using coordinates.

endFaceCenter = (-100,0,12.5)
endFace = myInstance.faces.findAt((endFaceCenter,) )

# Create a boundary condition that encastres one end
# of the beam.

endRegion = (endFace,)
myModel.EncastreBC(name='Fixed',createStepName='beamLoad',
    region=endRegion)

# Find the top face using coordinates.

topFaceCenter = (0,10,12.5)
topFace = myInstance.faces.findAt((topFaceCenter,) )

# Create a pressure load on the top face of the beam.

topSurface = ((topFace, SIDE1), )
myModel.Pressure(name='Pressure', createStepName='beamLoad',
    region=topSurface, magnitude=0.5)

#-------------------------------------------------------

import mesh

# Assign an element type to the part instance.

region = (myInstance.cells,)
elemType = mesh.ElemType(elemCode=C3D8I, elemLibrary=STANDARD)
myAssembly.setElementType(regions=region, elemTypes=(elemType,))

# Seed the part instance.

myAssembly.seedPartInstance(regions=(myInstance,), size=10.0)

# Mesh the part instance.

myAssembly.generateMesh(regions=(myInstance,))

# Display the meshed beam.

myViewport.assemblyDisplay.setValues(mesh=ON)
myViewport.assemblyDisplay.meshOptions.setValues(meshTechnique=ON)
myViewport.setValues(displayedObject=myAssembly)

#-------------------------------------------------------

import job

# Create an analysis job for the model and submit it.

jobName = 'beam_tutorial'
myJob = mdb.Job(name=jobName, model='Beam',
    description='Cantilever beam tutorial')

# Wait for the job to complete.

myJob.submit()
myJob.waitForCompletion()

#-------------------------------------------------------

import visualization

# Open the output database and display a
# default contour plot.

myOdb = visualization.openOdb(path=jobName + '.odb')
myViewport.setValues(displayedObject=myOdb)
myViewport.odbDisplay.display.setValues(plotState=CONTOURS_ON_DEF)

myViewport.odbDisplay.commonOptions.setValues(renderStyle=FILLED)

[derik]

Код:

[Выделить все]

fini
/clear,nostart
/outp,output,txt,C:\ansys\output
/cwd,'C:\ansys\out'
/units,mks

/prep7

*set,w1,.18
*set,w2,.18
*set,w3,.443
*set,t1,.011
*set,t2,.011
*set,t3,.0078
*set,a,3
*set,p,5e3

et,1,188
keyopt,1,8,1
keyopt,1,9,1
mptemp,,,,,,,
mptemp,1,0
mpdata,ex,1,,2.15e8
mpdata,prxy,1,,.3
tb,biso,1,1,2,
tbtemp,0
tbdata,,1.95e8,2540

sectype,1,beam,i,,0
secoffset,cent
secdata,w1,w2,w3,t1,t2,t3,
k,,
k,,a,
lstr,1,2

lesize,all,a/40,
lmesh,all

sfbeam,all,1,pres,p,p,

ksel,s,loc,x,0
dk,all,,,,,all,
alls
fini
/eshape,1

/sol

/gst,1
nldiag,nrre,on
antype,0
nlgeom,1
nsubst,20,100,10
outres,erase
outres,all,all
time,.5
solve
fini

/post1

set,first
etable,fx,smisc,1
etable,fx,smisc,14
etable,my,smisc,3
etable,mz,smisc,24
etable,tq,smisc,4
etable,sfz,smisc,5
plls,my,my,1,0

вот еще один пример - балочка обыкновенная!