[ИванМоскалев]

Добрый день. При написании программы наткнулся на одну проблему, над которой очень долго мучаюсь и не могу придумать выхода из нее.
Вот код программы. Кратко суть программы в том, что она создает стержень, закрепленный снизу,задается переменное(в тексте оно постоянное) ускорения и осуществляется динамический расчет. Но мне нужно запихать еще сюда демпфер. И вот тут начинается проблема. Когда пытаюсь присвоить переменным значения этих точек, то первая точка ReferencePoint-находится и все хорошо с ней, но вот вторую точку не могу присвоить и все тут! Подскажите, пожалуйста, что делать дальше..буду очень благодарен!
Ошибка в выделенной строке следующая:"TypeError: xVertices; found Vertex, expecting GeomSequence"


from abaqus import *
from part import *
from material import *
from section import *
from assembly import *
from step import *
from amplitude import *
from interaction import *
from abaqusConstants import *
from load import *
from mesh import *
from job import *
from sketch import *
from visualization import *
from connectorBehavior import *
from xyPlot import *
mdb.models['Model-1'].ConstrainedSketch(name='__profile__', sheetSize=20.0)
mdb.models['Model-1'].sketches['__profile__'].Line(point1=(0, -1.0),point2=(0.0,0.0))
mdb.models['Model-1'].Part(dimensionality=TWO_D_PLANAR, name='Part-1', type=DEFORMABLE_BODY)
mdb.models['Model-1'].parts['Part-1'].BaseWire(sketch=mdb.models['Model-1'].sketches['__profile__'])
mdb.models['Model-1'].Material(name='Material-1')
mdb.models['Model-1'].materials['Material-1'].Elastic(table=((100000,0.32), ))
mdb.models['Model-1'].RectangularProfile(name='Profile-1',a=1,b=1)
mdb.models['Model-1'].BeamSection(name='Section-1', integration=DURING_ANALYSIS, profile='Profile-1', material='Material-1')
et=mdb.models['Model-1'].parts['Part-1'].edges.findAt(((0.0,-0.5,0.0),))
rew = Region(edges=et)
mdb.models['Model-1'].parts['Part-1'].SectionAssignment(region=rew, sectionName='Section-1')
elemType=ElemType(elemCode=B21,elemLibrary=STANDARD)
mdb.models['Model-1'].parts['Part-1'].setElementType(elemTypes=(elemType,), regions=rew)
mdb.models['Model-1'].parts['Part-1'].seedPart(size=1)
mdb.models['Model-1'].parts['Part-1'].generateMesh()
myViewport = session.Viewport(name='Cantilever Beam Example',origin=(20, 20), width=150, height=120)
myAssembly = mdb.models['Model-1'].rootAssembly
#myViewport.assemblyDisplay.setValues(mesh=ON)
#myViewport.assemblyDisplay.meshOptions.setValues(meshTechnique=ON)
#myViewport.setValues(displayedObject=myAssembly)
mdb.models['Model-1'].ImplicitDynamicsStep(name='Step-1',previous='Initial',timePeriod=30,timeIncrementationMethod=FIXED,maxNumInc=600,initialInc=0.05)
myAssembly = mdb.models['Model-1'].rootAssembly
myInstance = myAssembly.Instance(name='beamInstance',part=mdb.models['Model-1'].parts['Part-1'], dependent=ON)
Bassse = myInstance.vertices.findAt(((0.0,-1.0,0.0),))
bassse = Region(vertices=Bassse)
mdb.models['Model-1'].EncastreBC(name='FixBase',createStepName='Step-1',region=bassse)
mass = myInstance.vertices.findAt(((0.0,0.0,0.0),))
Ma = Region(vertices=mass)
mdb.models['Model-1'].rootAssembly.engineeringFeatures.PointMassInertia(name='mass-1',region=Ma,mass2=1000)
mdb.models['Model-1'].parts['Part-1'].ReferencePoint(point=(0.0,-1.0,0.0))
fi = mdb.models['Model-1'].parts['Part-1'].referencePoints.findAt(((0.0,-1.0,0.0)))
Fi = Region(referencePoints=((fi),))
tw = mdb.models['Model-1'].parts['Part-1'].vertices.findAt(((0.0,0.0,0.0)))
Tw = Region(xVertices=tw)
myAssembly.engineeringFeatures.TwoPointSpringDashpot(axis=NODAL_LINE, name='Springs/Dashpots-1', regionPairs=(Fi,Tw), springStiffness=216.0)

[metod]

А что данная задача не может быть сделана в CAE?

[Dmitry88]

Человек сам себе усложняет решение задачи.))

[ИванМоскалев]

Сообщение от Dmitry88 Человек сам себе усложняет решение задачи.))

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


----- добавлено через ~39 мин. -----

Сообщение от metod А что данная задача не может быть сделана в CAE?

ответил выше

[metod]

ну так получается, что должно быть 500 запусков CAE?

[СергейД]

ansys apdl?

[Meknotek]

Побуду немного "кэпом", но функция(конструктор) ждет аргумент типа GeomSequence, а не объект типа Vertex (если я правильно понял из документации, Vertex - одиночная вершина, а GeomSequence - последовательность объектов (в т.ч. вершин). vertices - как раз и есть экземпляр субкласса GeomSequences, метод которого (findAt) похоже возвращает одиночный Vertex

причем парой строк выше Вы для Region задаете параметр в виде кортежа из одного элемента, может тут аналогично нужно?

[ИванМоскалев]

Так-то логично то,что Вы говорите,но я пробовал изначально vertices, но тогда он пишет что не найдено такого объекта,вот я и перепробовал кучу всяких возможных вариантов. Пробовал еще создать вторую ReferencePoint(ее можно создать),но потом когда ищу ее снова ничего не находит.

----- добавлено через 39 сек. -----

Сообщение от Meknotek Побуду немного "кэпом", но функция(конструктор) ждет аргумент типа GeomSequence, а не объект типа Vertex (если я правильно понял из документации, Vertex - одиночная вершина, а GeomSequence - последовательность объектов (в т.ч. вершин). vertices - как раз и есть экземпляр субкласса GeomSequences, метод которого (findAt) похоже возвращает одиночный Vertex причем парой строк выше Вы для Region задаете параметр в виде кортежа из одного элемента, может тут аналогично нужно?

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

----- добавлено через ~1 мин. -----

Сообщение от metod ну так получается, что должно быть 500 запусков CAE?

Именно так,а то и больше.

[СергейД]

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

[Meknotek]

Сообщение от ИванМоскалев Так-то логично то,что Вы говорите,но я пробовал изначально vertices, но тогда он пишет что не найдено такого объекта,вот я и перепробовал кучу всяких возможных вариантов. Пробовал еще создать вторую ReferencePoint(ее можно создать),но потом когда ищу ее снова ничего не находит..

Какие варианты Вы пробовали и какие ошибки Вам выдавал интерпретатор?

[ИванМоскалев]

Сообщение от СергейД все это проще сделать в ансис. была ситуация с вариантным просчетом модели в цикле около 11 тысяч раз. на выходе был итоговый массив с основными результатами.

Не спорю!Но лицензия есть только на Абакус

----- добавлено через ~6 мин. -----

Сообщение от Meknotek Какие варианты Вы пробовали и какие ошибки Вам выдавал интерпретатор?

Вот сейчас сделал так:
tw = mdb.models['Model-1'].parts['Part-1'].vertices.findAt(((0.0,0.0,0.0)))
Tw = Region(vertices=tw)
Выдает:"TypeError: vertices; found Vertex, expecting GeomSequence"

Потом через RP

myAssembly.ReferencePoint(point=(0.0,0.0,0.0))
k=myInstance.referencePoints.findAt(((0.0,0.0,0.0)))
Ki = Region(referencePoints=((k),))

Выдает:"TypeError: referencePoints[0]; found None, expecting tuple"

Так пробовал.Больше ничего не смог придумать.

[Meknotek]

Сообщение от ИванМоскалев Вот сейчас сделал так: tw = mdb.models['Model-1'].parts['Part-1'].vertices.findAt(((0.0,0.0,0.0))) Tw = Region(vertices=tw) Выдает:"TypeError: vertices; found Vertex, expecting GeomSequence"

а если попробовать:
Tw = Region(vertices=(tw,))
?

Сообщение от ИванМоскалев Потом через RP myAssembly.ReferencePoint(point=(0.0,0.0,0.0)) k=myInstance.referencePoints.findAt(((0.0,0.0,0.0))) Ki = Region(referencePoints=((k),)) Выдает:"TypeError: referencePoints[0]; found None, expecting tuple"

ну тут похоже переменная k почему-то равна None (findAt ничего не нашел?)

[ИванМоскалев]

Сообщение от Meknotek а если попробовать: Tw = Region(vertices=(tw,)) ? ну тут похоже переменная k почему-то равна None (findAt ничего не нашел?)

Получается что так..Просто создать точку вот таким образом:"mdb.models['Model-1'].parts['Part-1'].referencePoints.findAt(((0.0,-1.0,0.0)))" можно только один раз, поэтому я создаю ее через rootAssembly и т.д. и ищу ее там же,но не пойму почему он не находит ничего. А если через обычный vertices,то находит Vertex и его не взять ни xVertices ни просто vertices. Какая-то странная ситуация на самом деле

В случае
tw = mdb.models['Model-1'].parts['Part-1'].vertices.findAt(((0.0,0.0,0.0)))
Tw = Region(vertices=(tw,))

Отвечает:"TypeError: vertices; found tuple, expecting GeomSequence"

[Meknotek]

Сообщение от ИванМоскалев В случае tw = mdb.models['Model-1'].parts['Part-1'].vertices.findAt(((0.0,0.0,0.0))) Tw = Region(vertices=(tw,)) Отвечает:"TypeError: vertices; found tuple, expecting GeomSequence"

а если передать ...vertices (напрямую) что выдаст?

[ИванМоскалев]

Сообщение от Meknotek а если передать ...vertices (напрямую) что выдаст?

Не понял Вас??Как понять напрямую?

[Meknotek]

Я например попробовал так:
tw = mdb.models['Model-1'].parts['Part-1'].vertices[0:1]
(просто чтобы получить последовательность из одного Vertex)

все равно дальше не идет

[ИванМоскалев]

Сообщение от Meknotek Я например попробовал так: tw = mdb.models['Model-1'].parts['Part-1'].vertices[0:1] (просто чтобы получить последовательность из одного Vertex) все равно дальше не идет

Сделал вот таким образом еще:

myAssembly.engineeringFeatures.TwoPointSpringDashpot(axis=NODAL_LINE,dashpotBehavior=OFF, dashpotCoefficient=0.0, name='Springs/Dashpots-1', regionPairs=((Ma,Fi),),springBehavior=ON , springStiffness=216.0)

Ma-Раньше найдена.

Но отвечает:"AbaqusException: No valid region pairs have been supplied for the springs/dashpots."

[Meknotek]

У Вас fi получается "пустым" набором. Так и должно быть?
>>>print fi
({})

[ИванМоскалев]

Сообщение от Meknotek У Вас fi получается "пустым" набором. Так и должно быть? >>>print fi ({})

О!!!!Спасибо огромное!На самом деле это я с дуру начал вводить RP,хотя можно было обойтись уже найденными точками!В итоге сделал вот так:

myAssembly.engineeringFeatures.TwoPointSpringDashpot(axis=NODAL_LINE,dashpotBehavior=OFF, dashpotCoefficient=0.0, name='Springs/Dashpots-1', regionPairs=((Ma,bassse),),springBehavior=ON , springStiffness=216.0)

Спасибо Вам огромнейшее,что повозились со мной!