elements:bem2d
Differences
This shows you the differences between two versions of the page.
| Next revision | Previous revision | ||
|
elements:bem2d [2019/07/19 12:02] helene created |
elements:bem2d [2020/08/25 15:46] (current) |
||
|---|---|---|---|
| Line 14: | Line 14: | ||
| |TITLE|"BEM2D" in the first 5 columns| | |TITLE|"BEM2D" in the first 5 columns| | ||
| ^Control data (I5)^^ | ^Control data (I5)^^ | ||
| - | |NEFEQ| Number of equivalent finite elments| | + | |NEFEQ| Number of equivalent finite elements| |
| ==== Definition of an equivalent finite element ==== | ==== Definition of an equivalent finite element ==== | ||
| This section must be repeated NEFEQ times. | This section must be repeated NEFEQ times. | ||
| Line 23: | Line 23: | ||
| |NSIGP|Number of interior points where stresses should be computed| | |NSIGP|Number of interior points where stresses should be computed| | ||
| |IFTRA|Index for the calculation of nodal traction at each step: \\ = 1 Yes \\ ≠ 1 No| | |IFTRA|Index for the calculation of nodal traction at each step: \\ = 1 Yes \\ ≠ 1 No| | ||
| - | |INSOL|Index for the choice of fundamental solution (see [[elements:bem2d#Fundamental solution|explanation]]): \\ = 0 Kelvin (finite domain) \\ = 1 Kelvin (infinite domain) \\ = 2 Melan (semi-infinite)| | + | |INSOL|Index for the choice of fundamental solution (see [[elements:bem2d#Fundamental solution|explanation]] below): \\ = 0 Kelvin (finite domain) \\ = 1 Kelvin (infinite domain) \\ = 2 Melan (semi-infinite)| |
| |NRI|= 0 0 → No rotation \\ = -1 -1 → Rotation during deformation; the code computes the two most distant points \\ = I J → Rotation during deformation: I is the first node to define the angle of rotation; J is the second node to define the angle of rotation| | |NRI|= 0 0 → No rotation \\ = -1 -1 → Rotation during deformation; the code computes the two most distant points \\ = I J → Rotation during deformation: I is the first node to define the angle of rotation; J is the second node to define the angle of rotation| | ||
| |NRJ|:::| | |NRJ|:::| | ||
| + | ^Definition of the points of stress computation (2G10.0) - repeated NSIGP times^^ | ||
| + | |XP|Coordinates of the point where the stress must be computed| | ||
| + | |YP|:::| | ||
| + | ^Definition of the BEM elements (2I5/3I5)^^ | ||
| + | |NNO|Number of nodes (2 or 3 in plane state, 3 in axisymmetric)| | ||
| + | |NINTE|Number of Gauss integration points (2, 4, or 6 in plane state; 0 in axisymmetric because the program computes it automatically)| | ||
| + | |NODE(I) \\ I=1,NNO|List of nodes of the BEM element (Beware: the elements must be in order - the the first node of element i = the last node of element j)| | ||
| === Fundamental solution === | === Fundamental solution === | ||
| + | Two fundamental solutions are available: the Kelvin solution and the Melan solution. \\ | ||
| + | The Kelvin solution is used for a border defining a __closed__ contour. The considered domain is either the interior (finite domain) or the exterior (infinite domain). The order of numbering for the nodes is indicated in the figure below. \\ | ||
| + | The problem can be in plane strain state or in axisymmetric state. No volumetric force is taken into account (no specific weight, thermal dilatation, or effective stresses). \\ | ||
| + | {{:elements:bem2d_kelvin.png?600|}} | ||
| + | |||
| + | The Melan solution is used for a semi-inifinite domain with a straight border. Only the deviations from this line must be discretized. The problem is in plane strain state. No volumetric force is taken into account. \\ | ||
| + | {{ :elements:bem2d_melan.png?400 |}} | ||
| + | |||
elements/bem2d.1563530569.txt.gz · Last modified: 2020/08/25 15:33 (external edit)
Page Tools
