Работа с разреженными матрицами в МКЭ

[Vassil]

Можете посмотреть сайт г-на Matthias K. Gobbert из
University of Maryland, Baltimore County (UMBC)

http://www.math.umbc.edu/~gobbert/ ,

в котором есть под-сайт

http://www.math.umbc.edu/~gobbert/femlab/

как в MATLAB ползуют МКЭ ?!

[Аспир]

Помогите составить уравнение МКЭ для уравнения типа:
a*n(x)-b*n(x)*f1(x)-b*f2(x)*dn(x)/dx=0, где a и b - постоянные, f1(x), f2(x) - известные функции (f1(x)=df2(x)/dx).

Dan_Te:

Пожалуйста, используйте тег math!

[Аспир]

Помогите составить уравнение МКЭ для задачи о дрейфе электронов (дырок) в переменном электрическом поле (уравнение непрерывности для электрического заряда).
Заранее благодарен за помощь!

[VLarin]

Для вариационной постановки? а зачем
У Вас уже есть некое уравнение, воспользуйтесь методом взвешенных невязок.
И такой вопрос - зачем для такой задачи МКЭ понадобился?
И что Вы именно хотите посчитать?

[Аспир]

Надо посчитать распределение электронов и дырок в полупроводнике под действием внешнего электрического поля, созданного пьезоэлектрическим резонатором. Система уравнений состоит из уравнений движения и Максвелла (для описания колебаний пьезоэлектрика)+уравнения непрерывности для электронов и дырок в этом поле. Уравнения МКЭ для первых двух -известны, а вот как быть с третьим???

[Аспир]

Посоветуйте, пожалуйста, как составить матрицу жесткости для расчета колебаний кострукции, которая состоит из двух слоев, между которыми скользкий контакт (нормальная компонента смещений - непрерывна, остальные две терпят разрыв на границе). Как учитывать такие условия?
Заранее благодарен.

[Zai]

Для многослойных изгибных пластин и оболочек вопрос о скольжении не имеет значения, если выполняется гипотеза Киргоффа-Лява о радиусе кривизны, так что можите рассматривать Вашу оболочку без скольжения.

[Аспир]

Для многослойных изгибных пластин и оболочек вопрос о скольжении не имеет значения, если выполняется гипотеза Киргоффа-Лява о радиусе кривизны, так что можите рассматривать Вашу оболочку без скольжения.
У меня тот случай, когда надо сравнить два типа граничных условий: первый тип- жесткая склейка двух пьезоэлектриков. При этом условия непрерывности смещений и механических напряжений на границе выполняются автоматически при построении матрицы жесткости. Для скользкого контакта непрерывной остается только нормальная компонента смещений и мех. напряжений. Тангенциальные смещения терпят разрыв. Вот с учетом этих условий у меня и возникает проблема. Непонятно, как модифицировать матрицу жесткости? И надо ли ее вообще как-то изменять?
Буду благодарен за отклики.