книги / САПР изделий из композиционных материалов. Моделирование процессов деформирования и разрушения в среде ANSYS
.pdfРис. П2.1. Расчетная схема
Последовательность действий через ГПИ:
– описать скалярные переменные, описывающие геометрию конст-
рукции; меню: |
|
|
|
|
Utility Menu |
Parameters |
Scalar Parameters... |
|
|
Ввести в строке ввода: pi = 3.1416 Accept |
Dr = 10 Accept |
Ar = |
||
=60 Accept |
B2r = 40 Accept → Jpol = pi*Dr**4/32 → Accept |
Close |
||
– задать тип КЭ; меню: |
|
|
|
|
Main Menu |
Preprocessor |
Element Type |
Add/Edit/Delete |
Add |
Beam → 3D elastic 4;
– задать геометрические характеристики сечения (двутавр 20а): площадь поперечного сечения (AREA = pi*dr**2/4), момент инерции (IZZ = jpol/2) и полярный момент инерции (IXX = jpol) балок в форме вещест-
венных констант (рис. П2.2); меню: |
|
|
|
|
||||
Main Menu |
Preprocessor |
Real Constants |
Add/Edit/Delete |
Add ... |
||||
– описать свойства материала; модуль упругости и модуль сдвига |
||||||||
(ÌÏà); ìåíþ: |
|
|
|
|
|
|
|
|
Main Menu |
Preprocessor |
Material Props |
Constant Isotropic |
OK |
||||
EX2e |
6 |
|
|
|
|
|
|
|
GXY8e |
5 |
|
|
|
|
|
|
|
Для более поздних версий ANSYS соответственно: |
|
|||||||
Main Menu |
Preprocessor |
Material Props |
Material Models |
|||||
Structural |
Linear |
Elastic |
Orthotropic; |
|
|
|
||
– создать узлы по координатам на плоскости; меню: |
|
|||||||
Main Menu |
Preprocessor |
|
Create |
Nodes |
In Active CS → |
|||
NODE |
Node Number |
|
1 |
|
|
|
||
X,Y,Z |
Location in active CS |
0 |
0 |
0 |
Apply |
|||
NODE |
Node Number |
|
2 |
|
|
|
||
X,Y,Z |
Location in active CS |
Ar |
0 |
0 |
Apply |
181
Рис. П2.2. Диалоговое окно описания вещественных констант
NODE |
Node Number |
3 |
|
|
|
|
X,Y,Z |
Location in active CS Ar+B2r |
0 |
0 |
Apply |
||
NODE |
Node Number |
4 |
|
|
|
|
X,Y,Z |
Location in active CS Ar+B2r*2 |
0 |
0 |
Apply |
||
NODE |
Node Number |
5 |
|
|
|
|
X,Y,Z Location in active CS |
(Ar+B2r)*2 |
0 |
0 |
OK |
||
– создать линейные элементы по парам узлов в последовательности: |
||||||
1 − 2; 2 − 3; 3 − 4; 4 − 5; ìåíþ: |
|
|
|
|
|
|
Main Menu |
Preprocessor |
Create |
Elements |
|
Auto Numbered |
|
Thru Nodes |
выбрать узлы 1 и 2 |
Apply; выбрать узлы 2 и 3 → |
||||
Apply è ò.ä. |
|
|
|
|
|
|
182
– задать условия закрепления: заделка в узлах 1 и 5 (ограничение по перемещениям и углу поворота); меню:
Main Menu Preprocessor (Solution) (Define) Loads Apply Displacements On Nodes выбрать узлы 1 и 5 All DOF OK;
–задать условия нагружения крутящие моменты в узлах 2 и 4; меню:
Main Menu Preprocessor (Solution) (Define) Loads Apply Force/Moment On Nodes выбрать узлы 2 и 4 OK в диалоговом окне в поле Lab выбрать MX, в поле VALUE внести значение: 8e4 OK;
–запустить решатель; меню:
Main Menu Solution Solve Current LS;
– просмотреть опорные реакции; варианты меню:
Main Menu General Postproc List Results Reaction Solu MX; Utility Menu List Results Reaction Solution MX;
– просмотреть угол закручивания; меню:
Main Menu General Postproc List Results Nodal Solution DOF solution Rotation ROTX;
Utility Menu List Results Nodal Solution DOF solution
Rotation ROTX. Соответствующий данному диалоговому способу командный метод приведен в листинге П2.1.
/com, |
ПРОГРАММА РЕШЕНИЯ ЗАДАЧИ СОПРОТИВЛЕНИЯ |
|
/com, |
МАТЕРИАЛОВ ПО ТЕМЕ: “КРУЧЕНИЕ” |
|
/com, |
|
|
/com, |
Сборник задач по сопротивлению материалов / |
|
/com, |
Под ред. В.К. Качурина. М.: Наука, 1970. |
|
/com, |
(по задаче 4.33. |
Ñ.94.) |
/com, |
|
|
/com, |
Итоговые результаты записываются |
|
/com, |
в стандартный файл вывода *.out |
|
/com, |
|
|
/prep7 |
! Входим в препроцессор |
|
|
! Все данные в системе СГС |
|
pi=3.1416 |
|
|
dr=10 |
! dr - диаметр стержня |
|
ar=60 |
! ar, b2r - геометрические параметры |
|
b2r=40 |
|
|
183
jpol=pi*dr**4/32 |
! константа для полярного момента |
|
! инерции |
et,1,beam4 |
! балочные КЭ ВЕАМ4 |
r,1,pi*dr**2/4,jpol/2 ! REAL CONST :
!площадь поперечного сечения AREA=pi*dr**2/4 ,
!момент инерции IZZ=jpol/2
rmore„jpol |
! полярный момент инерции |
|
! (8-ая REAL-константа для ВЕАМ4) |
Mp,ex,1,2e6 ! модуль Юнга ЕХ=2е6 |
|
mp,gxy,1,8å5 |
! модуль сдвига GXY=8е5 |
n,1 |
! Определяем узлы по координатам |
n,2,ar |
|
n,3,ar+b2r |
|
n,4,ar+b2r*2 |
|
n,5,(ar+b2r)*2 |
|
*do,i,1,4 |
! Определяем элементы ВЕАМ4 |
e,i,i+1 |
|
*enddo |
|
finish |
|
/solu |
! Входим в решатель |
antype,static |
! Статический тип анализа |
d,1,all,0„5,4 |
! Задаем условия закрепления в узлах |
f,2,mx,8å4 |
! Задаем крутящие моменты в узлах |
f,4,mx,8å4 |
|
outpr,basic,1 |
! Вывод основных результатов в стан- |
|
! дартный файл |
solve |
! Решаем СЛАУ МКЭ |
finish |
|
/prep7 |
! Вариант с КЭ РIРЕ16 |
et,1,pipe16 |
! КЭ - трубки РIРЕ16 |
r,1,dr,dr/2 |
! REAL CONST : |
|
! внешний диаметр сечения OD=dr, |
|
! толщина стенок трубы TKWALL=dr/2 |
|
! (сплошная труба) |
finish |
|
/solu |
|
solve |
|
184
finish /exit
Листинг П2.1. Текст программы расчета стержня на кручение
В результате выполнения программы в файле jobname.out (jobname имя создаваемых рабочих файлов) среди различной выводной информации найдем значения реактивных моментов МХ = 8е4 (кг·см) и угол закручивания ROTX = 0.611155е 2 (рад), что совпадает с ответом к задаче.
Повторить решение, заменив КЭ BEAM4 на КЭ PIPE16, используя добавление КЭ в список типов КЭ и переназначение атрибутов всех КЭ (см. главу 3).
Приложение 3
Расчет ферм [16]
В качестве примера расчета фермы рассмотрим задачу ¹ 4.67 из задачника [17]. В данной задаче требуется определить опорные реакции и усилия в стержнях фермы, показанной на рис. П3.1 вместе с действующими в узлах силами.
Отнесем ферму к декартовой системe координат 0XY, и по рисунку, приведенному в задачнике, рассчитаем координаты ее узлов. На рис. П3.1 узлы фермы обозначены метками n1,…,n5. Эти узлы имеют следующие координаты (в метрах): n1 − (0; 0), n2 − (2; 0), n3 − (4; 0), n4 − (2.5;3/2), n5 − (1;3).
Как известно, ферма есть конструкция, состоящая из стержней, соединенных между собой шарнирами (узлами фермы). Поскольку силы, действующие на ферму, прилагаются в шарнирах, то все стержни фермы испытывают только осевые усилия растяжения-сжатия. Следовательно, подходящими КЭ для стержней ферм являются стержневые элементы LINK1 в случае плоских ферм и LINK8 в случае пространственных ферм. Узлы конечно-элементной модели будут тогда совпадать с узлами фермы, а каждый стержень будет отдельным КЭ.
Рис. П3.1. Расчетная схема конструкции
Приведенная на рис. П3.1 ферма с количеством стержней å = 7 и количеством узлов n = 5 является статически определимой, так как для нее выполняется условие статической определимости плоских ферм: å = 2n−3. Данную задачу можно решать методами теоретической механи-
186
ки, и тогда не существенны никакие определяющие параметры стержней, кроме их длин. Однако для КЭ LINK1 (или LINK8) требуется задать, по крайней мере, одно материальное свойство (модуль Юнга ЕХ) и одну вещественную константу КЭ (площадь поперечного сечения АRЕА). Примем произвольно, что для всех стержней фермы ЕХ = 2·1011 (Í/ì2) è AREA = 1·10−4 (ì2). Отметим, что в рассматриваемой статически определимой задаче значения этих параметров не будут влиять на итоговые искомые величины, подлежащие определению.
Последовательность действий через ГПИ:
– описать скалярные переменные F_R (имя файла для вывода резуль-
татов расчета) и S3 (значение |
3); ìåíþ: |
|
|
Utility Menu Parameters |
Scalar Parameters... |
|
|
Ввести в строке ввода: F_R=’SMs1’ |
Accept |
S3=SQRT(3) |
|
Accept Close; |
|
|
|
− задать тип КЭ; меню:
Main Menu → Preprocessor → Element Type → Add/Edit/Delete → Add → Link → 2D spar 1;
– задать площадь поперечного сечения стержней в форме вещест-
венной константы; меню: |
|
|
|
|
|
||
Main Menu |
Preprocessor |
Real Constants Add/Edit/Delete |
Add |
||||
AREA |
1.e−4; |
|
|
|
|
|
|
– описать свойства материала; меню: |
|
|
|
||||
Main Menu |
Preprocessor |
Material Props ... ÅÕ |
2å11. |
|
|||
Для более поздних версий ANSYS соответственно: |
|
||||||
Main Menu |
Preprocessor |
Material Props |
Material Models |
||||
Structural |
|
Linear Elastic |
Isotropic; |
|
|
|
|
– создать узлы по координатам на плоскости; меню: |
|
||||||
Main Menu |
Preprocessor |
Create |
Nodes |
In Active CS; |
|
||
– создать линейные элементы по парам узлов в последовательности: |
|||||||
1 − 2; 2 − 3; 3 − 4; 4 − 5; 5 − 1; 5 − 2; 2 − 4; ìåíþ: |
|
|
|
||||
Main Menu |
Preprocessor |
Create |
Elements |
Auto Numbered |
Thru |
||
Nodes выбрать узлы 1 и 2 |
Apply; выбрать узлы 2 и 3 Apply è ò.ä. |
||||||
Созданные элементы будут иметь следующие |
атрибуты |
(другие |
не описаны): номер типа элемента − ТУРЕ=1 (LINK1); номер вещественной константы − REAL=1 (AREA); номер набора материальных свойств − ÌÀÒ=1 (ÅÕ);
187
– задать условия закрепления: подвижная шарнирная опора в узле 1 (ограничение по перемещениям по оси Y) и неподвижный шарнир в уз-
ле 3 (ограничение по перемещениям по осям Õ è Y); ìåíþ: |
|
|||||
Main Menu |
Preprocessor |
(Solution) |
(Define) Loads |
Apply |
||
Displacements |
On Nodes |
выбрать узлы 1 и 3 |
UY Apply âû- |
|||
брать узел 3 |
UX |
OK; |
|
|
|
|
– задать условия нагружения − описать усилия, приложенные к уз- |
||||||
ëàì; ìåíþ: |
|
|
|
|
|
|
Main Menu |
Preprocessor |
(Solution) |
(Define) Loads |
Apply |
||
Force/Moment |
On Nodes |
выбрать узлы 1 и 3 |
Apply |
в диалого- |
||
вом окне выбрать FY и в поле Value внести значение: −1å3 |
Apply |
|||||
выбрать узлы 4 и 5 |
Apply |
в диалоговом окне выбрать FY и в поле |
||||
Value внести значение: −2å3 |
OK; |
|
|
|
||
– запустить решатель; меню: |
|
|
|
|||
Main Menu |
Solution |
Solve Current LS; |
|
|
–для записи результатов в файл или окно вывода используйте меню:
Utility Menu File Switch Output to File;
Utility Menu → File → Switch Output to → Output Window;
–просмотреть опорные реакции; варианты меню:
Main MenuGeneral PostprocList ResultsReaction SoluAll items; Utility Menu → List → Results → Reaction Solution → All items;
– просмотреть усилия в стержнях; варианты меню:
Main Menu General Postproc List Results Element Solution LineElem results;
Utility Menu → List → Results → Element Solution → LineElem results. Соответствующий данному диалоговому способу командный метод
(листинг П3.1):
F_R=’SMs1’ |
! Определяем скалярные переменные |
S3=SQRT(3) |
|
/PREP7 |
! Входим в препроцессор |
ET,1,LINK1 |
! Задаем тип КЭ |
R,1,1.å-4 |
! Задаем площадь сечения |
ÌÐ,ÅÕ,1,2Å11 |
! Задаем модуль упругости |
N,1,0,0 |
! Создаем узлы |
N,2,2,0 |
|
188
N,3,4,0 |
|
N,4,2.5,S3/2 |
|
N,5,1,S3 |
|
Å,1,2 |
! Создаем элементы по узлам |
Å,2,3 |
|
Å,3,4 |
|
Å,4,5 |
|
Å,5,1 |
|
Å,5,2 |
|
Å,2,4 |
|
FINISH |
! Выходим из препроцессора |
/SOLU |
! Входим в решатель |
ÀNÒÓÐÅ,SÒÀÒIÑ |
! Тип анализа - статический |
|
! Определяем шарнирные опоры |
|
! в узлах |
D,1,UY,0„3,2 |
! UY=0 в узлах 1 и 3 |
D,3,UX |
! UX=0 â óçëå Ç |
|
! Задаем силы в узлах |
F,1,FY,-1ÅÇ„3,2 |
! FY=-1еЗ в узлах 1 и 3 |
F,4,FY,-2ÅÇ„5 |
! FY=-2еЗ в узлах 4 и 5 |
SOLVE |
! Запуск решателя |
FINISH |
! Выходим из решателя |
/ÐÎSÒ1 |
! Входим в постпроцессор |
/OUTPUT,F_R,txt |
! Направляем вывод в файл SMs1.txt |
PRRSOL |
! Печатаем опорные реакции |
PRESOL,SMISC,1 |
! Печатаем усилия в стержнях |
/OUTPUT |
! Закрываем файл для записи |
FINISH |
! Выходим из постпроцессора |
Листинг П3.1. Текст программы расчета фермы
Перечень команд удобно записать в файл текстового вида построчно (рекомендуется файл с расширением *.txt) и считать его с помощью меню:
Utility Menu File Read Input from...
Сравнение полученных результатов с ответами к задаче ¹ 4.67 [10] (см. табл. П3.1 и П3.2) позволяет оценить правильность работы программы и точность расчетов ферм с помощью ANSYS.
189
|
|
|
Ò à á ë è ö à Ï 3 . 1 |
|
Реакции в узлах |
|
|
||
|
|
|
|
|
Номер узла |
1 |
|
3 |
|
Реакции в узлах (кН) |
Ответ из задачника |
3.25 |
|
2.75 |
|
Решение ANSYS |
3.25 |
|
2.75 |
Ò à á ë è ö à Ï 3 . 2
Усилия в стержнях
Номер стержня |
1 |
2 |
3 |
4 |
5 |
6 |
7 |
|
Усилия (кН) |
Ответ из задачника |
1.3 |
3.03 |
−3.5 |
−2.5 |
−2.6 |
1.73 |
−1.73 |
|
Решение ANSYS |
1.299 |
3.0311 |
−3.5 |
−2.5 |
−2.581 |
1.732 |
−1.7321 |