В чем проблема?

[Alexander555 2]

Постановка задачи.

Задача на микроуровне.

Рассчитать касательное напряжение на поверхности 3D модели эндотелия при взаимодействии с 3D моделью липопротеина.

Решение:

1 шаг. Построение.

Preferences  Structural.

Построили 3D модель эндотелия, применив макросы, для создания самой поверхности и зачем из угловых точек опустив до 0 по оси Z линии, чтобы получилась объемная модель.

Затем над моделью эндотелия (первоначальная высота модели липопротеина над поверхностью эндотелия 1е-6) построили модель липопротеина (по строению липопротеин представляет собой цилиндр со средними размерами: диаметр 240Å, высота 185 Å)

После построения моделей задали свойства материалов: модуль упругости, коэффициент Пуассона и плотность (материалы - изотропные, идеально упругие).

2 шаг. Разбиение на сетку конечных элементов и другие подготовки к решению.

Во-первых, ввели тип элемента так для модели эндотелия это SHELL 281, а для липопротеина SOLID 185.

Поверхность эндотелия разбили на сетку конечных элементов (Shape: tetrahedron  free) с размером сетки 4е-8.

Поверхность липопротеина разбили тоже на тетраэдры (Shape: tetrahedron  free) с размером сетки 5е-9.

Затем задали контактные пары – узел-поверхность (node-surface), коэффициент трения задали как 0, начальное проникновение выключили.

После того модели были разбиты на сетку мы закрепили поверхность эндотелия по боковым поверхностям. И задали скорость по оси Z : -0.6м/сек.

Настройки задачи.

New analysis  Transient.

Анализ вариантов – Large Displacement Transient. Время – 1e-7. Шаги – 1. Frequency – write every substep.

3 шаг. Решение задачи.

Solve. График доходил до ¾ всего отрезка.

Замечание: после решения задачи на высоте модели липопротеина над поверхностью эндотелия 1е-6, модель липопротеина пролетела через поверхность эндотелия, затем высоту липопротеина изменяли до того чтобы, высота соответствовала скорости, т.е. чтобы липопротеин точно остановился в месте контакта с эндотелием.

Затем, чтобы получить, более точное напряжение пришлось измельчить сетку на поверхности эндотелия и поверхность липопротеина разбить более точно (Shape: hex/Wedge  sweep) с размером сетки 2е-9. И контактные пары задали как поверхность-поверхность (surface –surface).

Потом мы опять пробуем запустить задачу, задача прорешивается (Solve  Solution is done), но липопротеин находится на том же самом месте, что и раньше, т.е. он не опускается на поверхность липопротеина. Затем перешли к старой сетке и он снова не опускается.

Вопросы: Что могло послужить причиной сбоя? Как лучше решать задачи такого типа?

[Dandy 3]

А для какого момента времени вы смотрите результаты?

[Борман 2 389]

Построили 3D модель эндотелия

ввели тип элемента так для модели эндотелия это SHELL 281

Поверхность эндотелия разбили на сетку конечных элементов (Shape: tetrahedron  free)

Такого быть не может.

Поверхность эндотелия разбили на сетку конечных элементов (Shape: tetrahedron  free) с размером сетки 4е-8.

Поверхность липопротеина разбили тоже на тетраэдры (Shape: tetrahedron  free) с размером сетки 5е-9.

Разница в 10 раз. Это точно требуется ?

Затем задали контактные пары – узел-поверхность (node-surface)

Почему не поверхность-поверхность.

Анализ вариантов – Large Displacement Transient. Время – 1e-7. Шаги – 1. Frequency – write every substep.

А где размер шага по времени?

[Alexander555 2]

А для какого момента времени вы смотрите результаты?

по всему отрезку времени.

[Alexander555 2]

Такого быть не может.

Разница в 10 раз. Это точно требуется ?

Почему не поверхность-поверхность.

А где размер шага по времени?

исправил спасибо, цилиндр падает, но график до конца не просчитывается!

Поверхность эндотелия разбили на сетку конечных элементов (Shape: tetrahedron  free) с размером сетки 4е-8.

Поверхность липопротеина разбили тоже на тетраэдры (Shape: tetrahedron  free) с размером сетки 5е-9.

Разница в 10 раз. Это точно требуется ? у цилиндра очень маленькие размеры так что приходиться разбивать так!

[sergeyd 3]

у вас коэффициенты в матрице жесткости на уровне машинного нуля.

перейдите в микро, нано и тп. единицы

[Fedor 1 610]

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

[sergeyd 3]

федор! у вас абстрактные представления...обычно...

вы проверяли? я = да. влияет.

ждем ответа от автора, а не от вас.

видел расчет в фемап. размеры элементов= сотые доли мм. система СИ.

при переходе в мкм и тп. частоты изменились в семь (кажется) раз.

поправлюсь= коэфф в мж не близки к нулю, а слишком велики

как минимум пропорциональны E*L^(-2)

распечатайте суперэлемент или матрицы HBMAT= убедитесь

при их суммировании и вычитании ошибка будет слишком велика.

в нелинейности ошибка накапливается как снежный ком.

внедрение и т.п. параметры контакта на уровне машинного нуля.

[Fedor 1 610]

"вы проверяли? я = да. влияет" - это же настолько очевидный прием. Раньше были итерационные методы, а они любят, чтобы все не слишком отличалось от единицы. Никто бы не стал ухудшать свойства программы. В чем выгода то ?

Просто сам писал такие прграммульки когда-то. Чисто конкретно.

Об этом же еще Тьюарсон в Разреженных матрицах писал когда-то :)

[sergeyd 3]

сомневаюсь, что заранее будут поставлены все фильтры на все возможные ошибки (а считать эту задачу в си ошибка)

даже у вас, федор!

и от эти автоматических фильтров может быть еще хуже.

[Fedor 1 610]

Что значит считать на си ошибка? Считают же на компьютере. И что значит автоматические фильтры? Желательно раскрыть такие странные тезисы :)

Не ошибается тот, кто ничего не делает - учит народ. Что делает каждый программист? Ошибки - учили меня когда-то. Проблема не в том, чтобы ошибиться, а в том, чтобы найти ее и доказать, что это именно она.

Поскольку сейчас не программирую, то и ошибок не делаю. Это же очевидно :)

[Alexander555 2]

Спасибо большое у меня цилиндр опустился на поверхность! и напряжение проявилось )))

[sergeyd 3]

2 федор

опять общие фразы...

заведите в свою программу размеры -8 степени и убедитесь, что сделали все правильно.

заведите в ансис теперь и проверьте.

[Fedor 1 610]

Зачем? Пусть кто-нибудь проверит. Задаст размеры в мм и км и сравнит результаты. Знал людей, которые занимались микромеханикой никогда не слышал от них каких-либо проблем. Выложите APDL код который доказывал бы обратное, тогда и поговорим. Как профессионал я работаю только за хорошие деньги :)