| Правила | Регистрация | Пользователи | Поиск | Сообщения за день | Все разделы прочитаны |  Справка по форуму | Файлообменник |

Вернуться   Форум DWG.RU > Программное обеспечение > Расчетные программы > простейшая термоупругая задача

простейшая термоупругая задача

Ответ
Поиск в этой теме
Непрочитано 03.01.2008, 17:03 #1
простейшая термоупругая задача
funny_alex
 
Регистрация: 17.12.2007
Сообщений: 37

нельзя ли решить задачу в один прием? то есть не решать сначала температурную задачу, а потом подгружать результаты в структурную задачу?
Просмотров: 13460
 
Непрочитано 06.01.2008, 23:10
#2
Разработчик

Ну типа прочнист
 
Регистрация: 12.01.2005
Москва
Сообщений: 1,649
<phrase 1=


А смысл? Есть, по большому счету, два типа задач термоупругости: связанная и несвязанная. Первая возникает, когда изменение НДС вызывает появление исочников (положительных или отрицательных) тепла или влияет на характеристики теплопроводности/теплоемкости материала. Такая задача действительно решается за один раз - т.е. в результате решения одновременно получаются и тепловые поля и НДС - но это задача большой размерности (грубо говоря, в каждом узле кроме перемещений добавляется еще одна степень свободы - температура) и, как правило, нелинейная. Поэтому, когда есть возможность (а она есть в большинстве задач) разумно последовательно решать две задачи меньшей размерности.
__________________
ZZH
Разработчик вне форума  
 
Автор темы   Непрочитано 15.01.2008, 19:19
#3
funny_alex


 
Регистрация: 17.12.2007
Сообщений: 37


Не подскажете каким образом решить связанную задачу термоупругости
funny_alex вне форума  
 
Непрочитано 16.01.2008, 11:03
#4
СергейД

расчеты МКЭ и CFD. ктн
 
Регистрация: 11.05.2005
Подмосковье
Сообщений: 2,170


изучи ансис и справишься
СергейД вне форума  
 
Непрочитано 16.01.2008, 13:43
#5
Разработчик

Ну типа прочнист
 
Регистрация: 12.01.2005
Москва
Сообщений: 1,649
<phrase 1=


Цитата:
Не подскажете каким образом решить связанную задачу термоупругости
Не подскажете каким образом спроектировать здание?
Не подскажете каким образом поставить спектакль?
Не подскажете каким образом выиграть ЧМ по футболу?
и т.п. в том же духе...
Конкретнее. Задачу в студию!
Цитата:
изучи ансис и справишься
А что, в АНСИСе есть элементы для связанной задачи? Ой сомневаюсь, уж больно экзотическая вещь для коммерческой программы.
__________________
ZZH
Разработчик вне форума  
 
Непрочитано 16.01.2008, 17:52
#6
СергейД

расчеты МКЭ и CFD. ктн
 
Регистрация: 11.05.2005
Подмосковье
Сообщений: 2,170


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

глубокоуважаемому Разработчику рекомендую для расширения кругозора
хоть часок почитать про ансис в интернете.
для начала ansys.com cadfem.ru emt.ru
СергейД вне форума  
 
Непрочитано 16.01.2008, 19:37
#7
Разработчик

Ну типа прочнист
 
Регистрация: 12.01.2005
Москва
Сообщений: 1,649
<phrase 1=


Цитата:
глубокоуважаемому Разработчику рекомендую для расширения кругозора хоть часок почитать про ансис в интернете.
Не очень понял чем это может расширить мой кругозор, с тем же успехом можно было предложить посмотреть рекламу Domestos. Я же не дилер, чтобы изучать рекламные буклеты.
Часик не потребовался, прочитал - нет там такого элемента, как я и подозревал. Связанная задача у них решается серией поочередных прогонов то тепловой, то упругой задач с перекачкой туда-сюда результатов. Сами понимаете к ответу на изначальный вопрос темы
Цитата:
нельзя ли решить задачу в один прием? то есть не решать сначала температурную задачу, а потом подгружать результаты в структурную задачу?
это никак не подходит. Так что не расширился кругозор, увы
__________________
ZZH
Разработчик вне форума  
 
Непрочитано 16.01.2008, 20:02
#8
holstenman


 
Регистрация: 04.06.2005
Сообщений: 178


Цитата:
Часик не потребовался
-читайте дальше
sequential метод не единственный способ решения связанных задач в ANSYS.
holstenman вне форума  
 
Непрочитано 16.01.2008, 20:15
#9
СергейД

расчеты МКЭ и CFD. ктн
 
Регистрация: 11.05.2005
Подмосковье
Сообщений: 2,170


ну и где же вы смотрели??????????????????????
и как???
небось,потратили две минуты а не час...
зачем же вы так самоуверенно...

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

вот хоть два элемента (их кажется, 6) в которых у элемента есть степени свободы скажем температура и перемещения и решаются они одновременно (с тзр пользователя)
вот три фрагмента скачанные только что из гугла.
наберите multifield+ansys
найдете сотни статей и диссертаций

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

The procedure for a coupled-field analysis depends on which fields are being coupled, but two distinct methods can be identified: load transfer and direct. These methods are described briefly below, and in the following chapters in detail:

"Direct Coupled-Field Analysis"

Load Transfer Methods

"The ANSYS Multi-field (TM) Solver - MFS Single-Code Coupling"

"Multi-field Analysis Using Code Coupling"

"Load Transfer Coupled Physics Analysis"

"Unidirectional Load Transfer"


ANSYS also offers the following additional coupled-field methods:

"Coupled Physics Circuit Simulation"

"Reduced Order Modeling"


1.1.1. Direct Method
The direct method usually involves just one analysis that uses a coupled-field element type containing all necessary degrees of freedom. Coupling is handled by calculating element matrices or element load vectors that contain all necessary terms. An example of a direct method coupled-field analysis is a piezoelectric analysis using the PLANE223, SOLID226, or SOLID227 elements. Another example is MEMS analysis with the TRANS126 element.

A FLOTRAN analysis using the FLOTRAN elements is another direct method. Refer to the Fluids Analysis Guide for detailed procedures for a FLOTRAN analysis.


1.1.2.1. Load Transfer Coupled Analysis - ANSYS Multi-field solver
The ANSYS Multi-field solver, available for a large class of coupled analysis problems, is an automated tool for solving load transfer coupled field problems. It supersedes the physics file-based procedure and provides a robust, accurate, and easy to use tool for solving load transfer coupled physics problems. Each physics is created as a field with an independent solid model and mesh. Surfaces or volumes are identified for coupled load transfer. A multi-field solver command set configures the problem and defines the solution sequencing. Coupled loads are automatically transferred across dissimilar meshes by the solver. The solver is applicable to static, harmonic, and transient analysis, depending on the physics requirements. Any number of fields may be solved in a sequential (or mixed sequential-simultaneous) manner.

Two versions of the ANSYS Multi-field solver, designed for different applications, offer their own benefits and different procedures:

MFS - Single code: The basic ANSYS Multi-field solver used if the simulation involves small models with all physics field contained within a single product executable (e.g., ANSYS Multiphysics). The MFS - Single code solver uses iterative coupling where each physics is solved sequentially, and each matrix equation is solved separately. The solver iterates between each physics field until loads transferred across the physics interfaces converge.

MFX - Multiple code: The enhanced ANSYS Multi-field solver used for simulations with physics fields distributed between more than one product executable (e.g., between ANSYS Multiphysics and ANSYS CFX). The MFX solver can accommodate much larger models than the MFS version. The MFX - Multiple code solver uses iterative coupling where each physics is solved either simultaneously or sequentially, and each matrix equation is solved separately. The solver iterates between each physics field until loads transferred across the physics interfaces converge.



14.226. SOLID226 - 3-D 20-Node Coupled-Field Solid
CAE Index > ANSYS Index > Release 11.0 Documentation for ANSYS







Matrix or Vector Geometry Shape Functions Integration Points
Stiffness and Stress Stiffness Matrices; and Thermal Expansion and Electrostatic Force Load Vector Brick (Equation 12–206), (Equation 12–207), and (Equation 12–208) 14
Wedge (Equation 12–183), (Equation 12–184), and (Equation 12–185) 3 x 3
Pyramid (Equation 12–168), (Equation 12–169), and (Equation 12–170) 2 x 2 x 2
Tet (Equation 12–161), (Equation 12–162), and (Equation 12–163) 4
Mass Matrix Same as stiffness matrix. 3 x 3 x 3 if brick. If other shapes, same as stiffness matrix
Pressure Load Vector Quad (Equation 12–66) and (Equation 12–67) 3 x 3
Triangle (Equation 12–43) and (Equation 12–44) 6
Thermal Conductivity Matrix and Heat Generation Load Vector Brick (Equation 12–209) 14
Wedge (Equation 12–186) 3 x 3
Pyramid (Equation 12–171) 2 x 2 x 2
Tet (Equation 12–164) 4
Specific Heat Matrix Same as thermal conductivity matrix
Convection Surface Matrix and Load Vector Quad (Equation 12–70) 3 x 3
Triangle (Equation 12–46) 6
Dielectric Permittivity and Electrical Conductivity Matrices; Charge Density, Joule Heating, and Peltier Heat Flux Load Vectors Brick (Equation 12–210) 14
Wedge (Equation 12–187) 3 x 3
Pyramid (Equation 12–172) 2 x 2 x 2
Tet (Equation 12–165) 4
Thermoelastic stiffness and damping matrices Same as combination of stiffness and thermal conductivity matrices
Piezoelectric Coupling Matrix Same as combination of stiffness matrix and dielectric matrix
Seebeck Coefficient Coupling Matrix Same as combination of electrical conductivity and thermal conductivity matrices
Surface Charge Density Load Vector Quad (Equation 12–172) 3 x 3
Triangle (Equation 12–47) 6

14.226.1. Other Applicable Sections
"Structures" describes the derivation of structural element matrices and load vectors as well as stress evaluations. General Element Formulations gives the general element formulations used by this element. "Electromagnetics" describes the derivation of dielectric and electric conduction matrices. Piezoelectrics discusses the piezoelectric capability used by the element. Piezoresistivity discusses the piezoresistive effect. Thermoelectrics discusses the thermoelectric effects. Thermoelasticity discusses the thermoelastic effects. Electroelasticity discusses the Maxwell stress electroelastic coupling.

4.13 PLANE13 2-D Coupled-Field Solid


PLANE13 has a 2-D magnetic, thermal, electrical, piezoelectric, and structural field capability with limited coupling between the fields. PLANE13 is defined by four nodes with up to four degrees of freedom per node. The element has nonlinear magnetic capability for modeling B-H curves or permanent magnet demagnetization curves. When used in a structural analysis PLANE13 has large deflection, large strain, and stress stiffening capabilities. See Section 14.13 in the ANSYS Theory Reference for more details about this element. Other coupled-field elements are SOLID5, SOLID98, and SOLID62.
Figure 4.13-1 PLANE13 2-D Coupled-Field Solid






4.13.1 Input Data
The geometry, node locations, and the coordinate system for this element are shown in Figure 4.13-1. The element input data includes four nodes and magnetic, thermal, electrical, and structural material properties. The type of units (MKS, CGS, or user defined) is specified through the EMUNIT command. EMUNIT also determines the value of MUZERO. The EMUNIT defaults are MKS units and MUZERO = 4 x 10-7 henries/meter. In addition to MUZERO, orthotropic relative permeability is specified through the MURX and MURY material property labels.
MGXX and MGYY represent vector components of the coercive force for permanent magnet materials. The magnitude of the coercive force is the square root of the sum of the squares of the components. The direction of polarization is determined by the components MGXX and MGYY. Permanent magnet polarization and orthotropic material directions correspond to the element coordinate directions. The element coordinate system orientation is as described in Section 2.3. Properties not input default as described in Section 2.4. Nonlinear magnetic B-H, piezoelectric, and anisotropic elastic properties are entered with the TB command as described in Section 2.5. Nonlinear orthotropic magnetic properties may be specified with a combination of a B-H curve and linear relative permeability. The B-H curve will be used in each element coordinate direction where a zero value of relative permeability is specified. Only one B-H curve may be specified per material.

Various combinations of nodal loading are available for this element (depending upon the KEYOPT(1) value). Nodal loads are defined with the D and the F commands. Nodal forces, if any, should be input per unit of depth for a plane analysis and on a full 360° basis for an axisymmetric analysis.

Element loads are described in Section 2.7. Pressures, convections, heat fluxes, and Maxwell force flags may be input on the element faces indicated by the circled numbers in Figure 4.13-1 using the SF and SFE commands. Positive pressures act into the element. Surfaces at which magnetic forces are to be calculated are identified by using the MXWF label on the surface load commands (no value is required). A maxwell stress tensor calculation is performed at these surfaces to obtain the magnetic forces. These forces are applied in solution as structural loads. The surface flag should be applied to "air" elements adjacent to the body for which forces are required. Deleting the MXWF specification removes the flag.

Body loads-temperature, heat generation rate, magnetic virtual displacement, and source current density-may be input at the element's nodes or as a single element value [BF, BFE]. When the temperature degree of freedom is active (KEYOPT(1)=2 or 4), applied body force temperatures [BF, BFE] are ignored. In general, unspecified nodal temperatures and heat generation rates default to the uniform value specified with the BFUNIF or TUNIF command. Heat generation from Joule heating is applied in Solution as thermal loading for static and transient analyses.

If the temperature degree of freedom is present, the calculated temperatures override any input nodal temperatures.

Air elements in which Local Jacobian forces are to be calculated may be identified by using nodal values of 1 and 0 for the MVDI label [BF]. See the ANSYS Electromagnetic Field Analysis Guide for details. These forces are not applied in solution as structural loads.

A summary of the element input is given in Table 4.13-1. A general description of element input is given in Section 2.1.

Table 4.13-1 PLANE13 Input Summary

Element Name
PLANE13

Nodes
I, J, K, L

Degrees of Freedom
AZ if KEYOPT (1) = 0
TEMP if KEYOPT (1 )= 2
UX,UY if KEYOPT (1) = 3
UX, UY, TEMP, AZ if KEYOPT (1) = 4
VOLT, AZ if KEYOPT (1) = 6
UX, UY, VOLT if KEYOPT (1) = 7

Real Constants
None

Material Properties
EX, EY, EZ, (PRXY, PRYZ, PRXZ or NUXY, NUYZ, NUXZ), ALPX, ALPY, ALPZ, DENS, GXY, DAMP, KXX, KYY, C, ENTH, MUZERO, MURX, MURY, RSVZ, MGXX, MGYY, PERX, PERY, plus BH, ANEL, and Piezoelectric data tables (see Section 2.5)

Surface Loads
Pressures:
face 1 (J-I), face 2 (K-J), face 3 (L-K), face 4 (I-L)
Convections:
face 1 (J-I), face 2 (K-J), face 3 (L-K), face 4 (I-L)
Heat Fluxes:
face 1 (J-I), face 2 (K-J), face 3 (L-K), face 4 (I-L)
Maxwell Force Flags:
face 1 (J-I), face 2 (K-J), face 3 (L-K), face 4 (I-L)

Body Loads
Temperatures:
T (I), T (J), T (K), T (L)
Heat Generations:
HG (I), HG (J), HG (K), HG (L)
Magnetic Virtual Displacements:
VD (I), VD (J), VD (K), VD (L)
Source Current Density:
spare, spare, JSZ (I), PHASE (I), spare, spare, JSZ (J), PHASE (J) spare, spare, JSZ (K), PHASE (K), spare, spare, JSZ (L), PHASE (L)

Special Features
Requires an iterative solution for field coupling (displacement, temperature, electric, magnetic, but not piezoelectric), Large deflection, Large strain, Stress stiffening, Birth and death, Adaptive descent

KEYOPT(1)
Used for degree of freedom selection
0 - AZ degree of freedom
2 - TEMP degree of freedom
3 - UX, UY degrees of freedom
4 - UX, UY, TEMP, AZ degrees of freedom
6 - VOLT, AZ degrees of freedom
7 - UX, UY, VOLT degrees of freedom

KEYOPT(2)
0 - Include extra shapes
1 - Do not include extra shapes

KEYOPT(3)
0 - Plane strain (with structural degrees of freedom)
1 - Axisymmetric
2 - Plane stress (with structural degrees of freedom)

KEYOPT(4)
0 - Element coordinate system is parallel to the global coordinate system
1 - Element coordinate system is based on the element I-J side

KEYOPT(5)
0 - Basic element printout
1 - Repeat basic solution for all integration points
2 - Nodal stress printout

KEYOPT(6)
Used only if the elastic stiffness matrix is input with the TB commands instead of with the MP commands (see Section 2.5)
0 - Input represents "stiffness" (or "C") matrix
1 - Input represents "compliance" (or "D") matrix




--------------------------------------------------------------------------------

4.13.2 Output Data
The solution output associated with the element is in two forms:

nodal degrees of freedom included in the overall nodal solution
additional element output as shown in Table 4.13-2.
Several items are illustrated in Figure 4.13-2. The element output directions are parallel to the element coordinate system. A general description of solution output is given in Section 2.2. See the ANSYS Basic Analysis Procedures Guide for ways to view results.
Figure 4.13-2 PLANE13 Element Output





The following notation is used in Table 4.13-2:

A colon ( in the Name column indicates the item can be accessed by the Component Name method [ETABLE, ESOL] (see Section 2.2.2). The O and R columns indicate the availability of the items in the file Jobname.OUT (O) or in the results file (R), a Y indicates that the item is always available, a number refers to a table footnote which describes when the item is conditionally available, and a - indicates that the item is not available.

Table 4.13-2 PLANE13 Element Output Definitions

Name
Definition
O
R

EL
Element number
Y Y
NODES
Nodes - I, J, K, L
Y Y
MAT
Material number
Y Y
VOLU:
Volume
Y Y
CENT: X, Y
Global location XC, YC
Y Y
PRES
P1 at nodes J, I; P2 at K, J; P3 at L, K; P4 at I, L
Y Y
TEMP
Input temperatures T(I), T(J), T(K), T(L)
Y Y
HGEN
Input heat generations HG(I), HG(J), HG(K), HG(L)
Y -
S:INT
Stress intensity
1 1
S:EQV
Equivalent stress
1 1
EPEL: X, Y, Z, XY
Elastic strains
1 1
EPEL: 1, 2, 3
Principal elastic strains
1 1
S: X, Y, Z, XY
Stresses (SZ=0.0 for plane stress elements)
1 1
S: 1, 2, 3
Principal stresses
1 1
TG: X, Y, SUM
Thermal gradient components and vector sum
1 1
TF: X, Y, SUM
Thermal flux (heat flow rate/cross-sectional area) components and vector sum
1 1
EF: X, Y
Electric field components (X,Y)
1 1
EF:SUM
Vector magnitude of EF
1 1
D: X, Y
Electric flux density components (X,Y)
1 1
D:SUM
Vector magnitude of D
1 1
UE, UD, UM
Elastic (UE), dielectric (UD), and electro-mechanical coupled (UM) energies
1 1
LOC
Output location (X,Y)
1 -
MUX,MUY
Magnetic permeability
1 1
H: X, Y
Magnetic field intensity components
1 1
H:SUM
Vector magnitude of H
1 1
B: X, Y
Magnetic flux density components
1 1
B:SUM
Vector magnitude of B
1 1
JSZ
Source current density, valid for static analysis only
1 1
JTZ
Total current density
1 1
JHEAT:
Joule heat generation per unit volume
1 1
FJB(X, Y)
Lorentz force components
1 1
FMX(X, Y)
Maxwell force components
1 1
FVW(X, Y)
Virtual work force components
1 1
FMAG: X, Y
Combined (FJB and FMX) force components
- 1
FACE
Face label
2 2
AREA
Face area
2 2
NODES
Face nodes
2 -
HFILM
Film coefficient at each node of face
2 -
TBULK
Bulk temperature at each node of face
2 -
TAVG
Average face temperature
2 2
HEAT RATE
Heat flow rate across face by convection
2 2
HEAT RATE/AREA
Heat flow rate per unit area across face by convection
2 -
HFLUX
Heat flux at each node of face
2 -
HFAVG
Average film coefficient of the face
2 2
TBAVG
Average face bulk temperature
- 2
HFLXAVG
Heat flow rate per unit area across face caused by input heat flux
- 2
TJB(Z)
Lorentz torque about global Cartesian +Z axis
1 1
TMX(Z)
Maxwell torque about global Cartesian +Z axis
1 1
TVW(Z)
Virtual work torque about global Cartesian +Z axis
1 1
1. Solution values are output only if calculated (based on input data).

2. Available only if a surface load is input.

Table 4.13-2a PLANE13 Miscellaneous Element Output

Description
Names of Items Output
O
R

Integration Pt. Solution
SINT, SEQV, EPEL, S, MUX, MUY, H, HSUM, B, BSUM
1 -
Nodal Solution
SINT, SEQV, S, H, HSUM, B, BSUM
2 -
1. Output at each integration point, if KEYOPT(5)=1.

2. Output at each node, if KEYOPT(5)=2.

Note-JT represents the total measurable current density in a conductor, including eddy current effects, and velocity effects if calculated.

Note-For axisymmetric solutions with KEYOPT(4)=0, the X and Y directions correspond to the radial and axial directions, respectively. The X,Y,Z, and XY stress output correspond to the radial, axial, hoop, and in-plane shear stresses, respectively.

Note-For harmonic analysis, joule losses (JHEAT), forces (FJB(X, Y), FMX(X, Y), FVW(X, Y)), and torque (TJB(Z), TMX(Z), TVW(Z)) represent time-average values. These values are stored in the "Real" data set. The macros POWERH, FMAGSUM, and TORQSUM can be used to retrieve this data.

Table 4.13-3 lists output available through the ETABLE command using the Sequence Number method. See Chapter 5 of the ANSYS Basic Analysis Procedures Guide and Section 2.2.2.2 of this manual for more information. The following notation is used in Table 4.13-3:


Name - output quantity as defined in the Table 4.13-2
Item - predetermined Item label for ETABLE command
E - sequence number for single-valued or constant element data
I,J,K,L - sequence number for data at nodes I,J,K,L
FCN - sequence number for solution items for element Face N
Table 4.13-3 PLANE13 Item and Sequence Numbers for the ETABLE and ESOL Commands
Name
Item
E
I
J
K
L

JSZ
SMISC
1 - - - -
P1
SMISC
- 4 3 - -
P2
SMISC
- - 6 5 -
P3
SMISC
- - - 8 7
P4
SMISC
- 9 - - 10
MUX
NMISC
1 - - - -
MUY
NMISC
2 - - - -
FVWX
NMISC
3 - - - -
FVWY
NMISC
4 - - - -
FVWSUM
NMISC
5 - - - -
JTZ
NMISC
7 - - - -
UE
NMISC
8 - - - -
UD
NMISC
9 - - - -
UM
NMISC
10 - - - -
TJB(Z)
NMISC
35 - - - -
TMX(Z)
NMISC
36 - - - -
TVW(Z)
NMISC
37 - - - -




FC1
FC2
FC3
FC4

AREA
NMISC
11 17 23 29
HFAVG
NMISC
12 18 24 30
TAVG
NMISC
13 19 25 31
TBAVG
NMISC
14 20 26 32
HEAT RATE
NMISC
15 21 27 33
HFLXAVG
NMISC
16 22 28 34


4.13.3 Assumptions and Restrictions
The area of the element must be positive. The element must lie in a global X-Y plane as shown in Figure 4.13-1 and the Y-axis must be the axis of symmetry for axisymmetric analyses. An axisymmetric structure should be modeled in the +X quadrants. For structural problems, the extra displacement shapes are automatically deleted for triangular elements so that a constant strain element results.
Transient magnetic analyses should be performed in a nonlinear transient dynamic analysis. A skin-effect analysis (where eddy current formation is permitted in conducting regions with impressed current loading) is performed by using KEYOPT(1)=6, specifying a resistivity, and coupling all VOLT degrees of freedom for elements in each of such regions. This is valid for both planar and axisymmetric models.

If a model has at least one element with piezoelectric degrees of freedom (displacements and VOLT) activated, then all elements where a VOLT degree of freedom is needed must be one of the piezoelectric types, and they must all have the piezoelectric degrees of freedom activated. If the piezoelectric effect is not desired in these elements, simply define very small piezoelectric material properties for them.

Note-For the axisymmetric option, the degree of freedom is VOLT*radius. See Section 14.53.3 of the ANSYS Theory Reference for details.


4.13.4 Product Restrictions
When used in the product(s) listed below, the stated product-specific restrictions apply to this element in addition to the general assumptions and restrictions given in the previous section.
ANSYS/Mechanical

Unless the Emag 2-D or Emag 3-D option is enabled, the following restrictions apply:


This element has only structural, thermal, or piezoelectric capability, and does not have magnetic capability.
The AZ degree of freedom is not active.
KEYOPT(1) defaults to 4 (UX, UY, TEMP) instead of 0, and cannot be set to 0. If set to 4 or 6, the AZ degree of freedom is not active.
The magnetic and electric material properties (MUZERO, MUR_, MG__, and the BH data table) are not allowed.
The Maxwell force flags surface loads are not applicable.
ANSYS/Emag 3-D

This element has only magnetic and electric field capability, and does not have structural, thermal, or piezoelectric capability.
The only active degrees of freedom are AZ and VOLT.
The only allowable material properties are the magnetic and electric properties (MUZERO through PERY, plus the BH data table).
The only applicable surface loads are Maxwell force flags. The heat generation body loads are not applicable. The temperature body load is only used for material property evaluation.
The element does not allow any special features.
KEYOPT(1) can only be set to 0 (default) or 6. KEYOPT(3)=2 is not applicable. KEYOPT(6) is not applicable.
ANSYS/Emag 2-D

This element has only magnetic and electric field capability, and does not have structural, thermal, or piezoelectric capability.
The only active degrees of freedom are AZ and VOLT.
The only allowable material properties are the magnetic and electric properties (MUZERO through PERY, plus the BH data table).
The only applicable surface loads are Maxwell force flags. The heat generation body loads are not applicable. The temperature body load is only used for material property evaluation.
The element does not allow any special features.
KEYOPT(1) can only be set to 0 (default) or 6. KEYOPT(3)=2 is not applicable. KEYOPT(6) is not applicable.
СергейД вне форума  
 
Непрочитано 17.01.2008, 13:12
#10
Разработчик

Ну типа прочнист
 
Регистрация: 12.01.2005
Москва
Сообщений: 1,649
<phrase 1=


Много букафф, ниасилил
Цитата:
небось,потратили две минуты а не час...
Разумеется. Смотрел, естественно, русские тексты, касающиеся термоупругости - их не так много.
Цитата:
почему вы высмеиваете мои познания в ансис? вы что его знаете не хуже меня?
Упаси Бог, я его вообще не знаю и потому Ваши познания в Ансис не оспаривал. Зато я знаю механику, в т.ч. и ее раздел, называемый термоупругостью и понимаю разницу между формулировкой связанной задачи термоупругости и возможными методами решения прикладных задач.
Обилие цитируемого материала говорит о том, что главного аргумента в пользу наличия в Ансисе элемента, позволяющего за один проход решать связанную задачу термоупругости (в чем, собственно я и сомневался) Вы там выделить не можете.
Пробежав текст по диагонали заметил один элемент для связанной пьезоэлектрической задачи
Цитата:
1.1.1. Direct Method
The direct method usually involves just one analysis that uses a coupled-field element type containing all necessary degrees of freedom. Coupling is handled by calculating element matrices or element load vectors that contain all necessary terms. An example of a direct method coupled-field analysis is a piezoelectric analysis using the PLANE223, SOLID226, or SOLID227 elements. Another example is MEMS analysis with the TRANS126 element.A FLOTRAN analysis using the FLOTRAN elements is another direct method. Refer to the Fluids Analysis Guide for detailed procedures for a FLOTRAN analysis.
Что такое MEMS и FLOTRAN не знаю, по аббревиатуре на термоупругость не похоже. Зато везде, где речь заходит о термоупругости натыкаемся на автоматизацию процесса "transfer", что говорит о раздельном решении.
Наличие у элемента дополнительных (к упругой задаче) степеней свободы, типа температуры и пр. означает лишь то, что обе задачи - тепловая и упругая - будут решаться на одной сетке, но из всего, что я углядел в текстах, на которые Вы ссылались следует, что решаться таки они будут раздельно. Если Вы, как знаток Ансис, приведете ссылку или фрагмент текста из которого следует, что это не так - буду весьма признателен. Это действительно расширит мой кругозор, будет оченно интересно узнать, как-таки ребята из Ансис разобрались с этим нехорошим функционалом связанной задачи термоупругости.
__________________
ZZH
Разработчик вне форума  
 
Непрочитано 17.01.2008, 19:42
#11
holstenman


 
Регистрация: 04.06.2005
Сообщений: 178


Цитата:
Много букафф, ниасилил
-жжоте, уважаемый
То есть если некий Разработчик читать не умеет - то и в ANSYS нет элементов для прямого решения связанных задач.
http://tigger.uic.edu/depts/accc/sof...-cou/GCOU3.htm
- хотя нафиг оно Вам надо, там букафф многа. Я уж про всякие матрицы-шматрицы молчу, вы-то о таком и не слышали по ходу.

Цитата:
Что такое MEMS и FLOTRAN не знаю
-гыг. а оно и правильно. Неча себе мозги компостировать. И ваще - "чукча не читатель, чукча писатель".
Встретимся на удаве...
holstenman вне форума  
 
Непрочитано 17.01.2008, 20:19
#12
Кулик Алексей aka kpblc
Moderator

LISP, C# (ACAD 200[9,12,13,14])
 
Регистрация: 25.08.2003
С.-Петербург
Сообщений: 39,787


Может, в более конструктивное русло беседу повернете, а?
__________________
Моя библиотека lisp-функций
---
Обращение ко мне - на "ты".
Все, что сказано - личное мнение.
Кулик Алексей aka kpblc вне форума  
 
Непрочитано 17.01.2008, 20:46
#13
holstenman


 
Регистрация: 04.06.2005
Сообщений: 178


Куда ж конструктивнее? Цитату привели, ссылку дали. Падение инженерной школы - чтение по диагонали...
holstenman вне форума  
 
Непрочитано 17.01.2008, 21:18
#14
Кулик Алексей aka kpblc
Moderator

LISP, C# (ACAD 200[9,12,13,14])
 
Регистрация: 25.08.2003
С.-Петербург
Сообщений: 39,787


Я имел в виду пригасить эмоции.
__________________
Моя библиотека lisp-функций
---
Обращение ко мне - на "ты".
Все, что сказано - личное мнение.
Кулик Алексей aka kpblc вне форума  
 
Непрочитано 18.01.2008, 01:11
#15
Разработчик

Ну типа прочнист
 
Регистрация: 12.01.2005
Москва
Сообщений: 1,649
<phrase 1=


Цитата:
имел в виду пригасить эмоции.
Нормально, обижаться - грех, эмоции от отсутствия аргументов.
Цитата:
То есть если некий Разработчик читать не умеет - то и в ANSYS нет элементов для прямого решения связанных задач
Что за нелепые наезды? Цитату - где я это сказал?!! Да и не мог, я уже признавался не раз, что в глаза не видел Ансис. Есть в математической физике куча связанных задач с удобными функционалами, которые прекрасно ложатся в русло "МКЭ системы, решающей ВСЁ". Но есть и такие связанные задачи, где сама вариационная формулировка проблематична. Уж не говоря о том, что если это удается, то говорить о возможности построения однопроходного решения в МКЭ, которое возможно только если функционал содержит квадратичную форму - не приходится.
Приведенная ссылка содержит многа букафф и опять не содержат информации о возможности Ансис за один проход решить связанную задачу термоупругости, о чем собственно и был первый пост. Я не первый раз здесь сталкиваюсь с ярыми поклонниками Ансис, которые утверждают, что "ЭТА ПРОГА МОЖЕТ ВСЕ", просто не подозревая, что принципиально нерешенные проблемы существуют не только в космофизике, но и в механике, сколь ни стара эта наука.
Оставьте эмоции, дайте конкретную ссылку и цитату, где утверждается, что Ансис или какая бы то ни было другая МКЭ программа за один проход решает связанную задачу термоупругости, и я скажу спасибо, ибо познаю новое.

ЗЫ. Ребята, остыньте. Вы - пользователи, я - разработчик. Я не владею таким обилием программ, как вы, не решаю таких сложных и оригинальных задач, используя всякие хитрые элементы, итерационные подходы и т.п. Я и своими-то программами толком пользоваться не умею: пользователи рассказывают мне о таких их возможностях, о которых я и не подозревал. Но я - механик/математик, знающий изнутри всю эту кухню и понимающий почему нельзя ни в Ансисе ни в какой-либо другой КЭ программе за один проход решить связанную задачу термоупругости или, например, расчитать балку со стесненным кручением в пластике.
__________________
ZZH
Разработчик вне форума  
 
Непрочитано 18.01.2008, 07:13
#16
Хворобьевъ


 
Сообщений: n/a


Цитата:
Сообщение от Разработчик Посмотреть сообщение
или, например, расчитать балку со стесненным кручением в пластике.
Можно - элементами beam188-189 (нарочно проверял сравнением оболочечной и балочной модели в ansys ED).
Просто результаты верны (дают хорошее совпадение) лишь на начальной стадии вступления в пластику (без ст. кручения - совпадение вплоть до развития пластических деформаций по всему сечению). Это - из-за "внутренней кухни", о чем бишь и в справке тоже написано.
Спор какой то ни о чем: вы пишите, что ANSYS в глаза не видели, но в то же время утверждаете, что там нет того, на что вам СергейД и holstenman указали по ссылке. Есть там элементы, поддерживающее одновременное решение тепловых и прочностных задач.
PS: ваши вечные "я разработчик", я "ракеты проектировал" (когда вы не находите аргументов) они не способствуют нахождению истины в споре. Я учился в университете, где каждый второй хрен (работник) тоже проектировал и проектирует до сих пор сопла ракет и прочие спутники. Не свидетельствует это ни о каких выдающихся их талантах - просто их ремесло. Как обычно - каждый специалист в своей узкой области, а есть и вообще среди них люди школьную задачу по математике не способные решить (но спутники летают прекрасно).
А в споре хороша тоько аргументация, а не аппеляция к авторитету. Даже если академик РАН со студентом спорит (я такое встречал).

Последний раз редактировалось Хворобьевъ, 18.01.2008 в 09:58.
 
 
Непрочитано 18.01.2008, 09:35
#17
kupolovich


 
Регистрация: 21.11.2007
Сообщений: 56


страсти накаляются....

P.S. Не сочтите за склочность, но всегда приятно наблюдать за "диалогом" знающих людей...
kupolovich вне форума  
 
Непрочитано 18.01.2008, 11:22
#18
СергейД

расчеты МКЭ и CFD. ктн
 
Регистрация: 11.05.2005
Подмосковье
Сообщений: 2,170


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

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

неужели, история повторяется? админ попроси разработчика все-таки сходить по указанным ссылкам...


разработчик!
если уж не умеете искать в интернете...а хотите лично от меня еще бесплатно информации...
укажите предварительно, будьте так любезны, максимальное количество одновременно Асиливаемых вами букОВ. нерусские языки (я понял) недопустимы? допустимы ли среди буков в русском тексте нерусские? каким шрифтом? какое направление диагонали вы предпочитаете?
не перевести ли вам предварительно на Олбанский?
СергейД вне форума  
 
Непрочитано 18.01.2008, 14:32
#19
Разработчик

Ну типа прочнист
 
Регистрация: 12.01.2005
Москва
Сообщений: 1,649
<phrase 1=


2Хворобьевъ
Цитата:
Можно - элементами beam188-189
Обсуждали со знатоком Ансис alle, выяснили, что нельзя и об этом наприсано в документации по Ансис:
http://dwg.ru/f/showthread.php?p=72810&#post72810
Не согласны с заключением alle, хотите продолжить тему? Мне было бы интересно получить решение простеньких тестовых задачек, начиная с той, что в посте 22 указанной темы.
Эмоциональный взрыв игнорирую - это не про меня. "Разработчик" - это не
Цитата:
не аппеляция к авторитету
а всего лишь позиционирование: не владею расчетными программами, не обладаю опытом расчета и проектирования строительных конструкций, но знаю механику, математику и зарабатываю на жизнь программированием.
Цитата:
А в споре хороша тоько аргументация
Вот и аргументируйте: дайте цитату, дайте решение задачи, и т.п., а не
Цитата:
учился в университете, где каждый второй хрен

2СергейД
А у меня были претензии?
Изначальный вопрос в теме "простейшая термоупругая задача" был:
Цитата:
нельзя ли решить задачу в один прием?
Я усомнился, что Ансис это может сделать. Вы дали ссылки, просмотрев в них русские тексты про связанную термоупругую задачу, я нашел лишь объяснения, как она решается в Ансис путем автоматизированной переброски результатов между тепловой и упругой задачами. В ответ Вы привели длиннющий английский текст в котором было упоминание про однопроходное решение связанной задачи, но не термоупругости, а пьезоэлектрической и каких-то двух еще, не похожих по аббревиатуре на термоупругость. В ответ на просьбу показать-таки где конкретно написано, что Ансис может за один проход решить связанную задачу термоупругости меня начали посылать... кто по многотомным ссылкам, а некоторые, чувствую, и дальше.
Спасибо за аргументированный диалог.
__________________
ZZH
Разработчик вне форума  
 
Непрочитано 18.01.2008, 19:59
#20
Хворобьевъ


 
Сообщений: n/a


Цитата:
Сообщение от Разработчик Посмотреть сообщение
Мне было бы интересно получить решение простеньких тестовых задачек, .
Благоволите -с...
Вложения
Тип файла: zip beam_vs_shell.zip (223.6 Кб, 286 просмотров)
 
Ответ
Вернуться   Форум DWG.RU > Программное обеспечение > Расчетные программы > простейшая термоупругая задача

Размещение рекламы
Опции темы Поиск в этой теме
Поиск в этой теме:

Расширенный поиск


Похожие темы
Тема Автор Раздел Ответов Последнее сообщение
Тестовая задача на нелинейность (МКЭ) Хворобьевъ Расчетные программы 35 28.03.2009 20:55
Практическая задача по выбору материалов для изготовления оп oo7 Прочее. Архитектура и строительство 4 21.06.2007 17:28
Задача линейная. Расчет на Лире rust-resisting Лира / Лира-САПР 19 12.03.2006 23:22
задача для LISP RSD LISP 5 22.08.2005 12:59
Задача для юзеров автокада: Vova AutoCAD 9 15.03.2004 16:34