Metafor

ULiege - Aerospace & Mechanical Engineering

User Tools

Site Tools


commit:futur:db

Differences

This shows you the differences between two versions of the page.

Link to this comparison view

Both sides previous revision Previous revision
Next revision
Previous revision
commit:futur:db [2018/07/17 11:49]
boemer [Description]
commit:futur:db [2018/07/17 12:03] (current)
boemer
Line 12: Line 12:
 ===== Test case ===== ===== Test case =====
  
-The test case in the following figure is identical to what was decribed previously.  A square is pushed against a contact tool (frictionless) by constant pressure $p$, its lateral movement is blocked and its out-of-plane thickness increases.+The test case in the following figure is identical to what was decribed previously.  A square is pushed against a contact tool (frictionless) with a pressure $p$, its lateral movement is blocked and its out-of-plane thickness increases (as a function of time). 
 + 
  
 The gap can be computed in various ways, numerically, analytically and semi-analytically/numerically.  This is what is done in three different extractors to test the good implementation of the previous feature:  The gap can be computed in various ways, numerically, analytically and semi-analytically/numerically.  This is what is done in three different extractors to test the good implementation of the previous feature: 
-  * First, the numerical computation simply consists in extracting the gaps in the Metafor computation and choosing the one at $P_{4}$. +  * ''numericalGap'': First, the numerical computation simply consists in extracting the gaps in the Metafor computation and choosing the one at $P_{4}$. 
-  * Second, the analytical computation provides the value of the gap via the following formula: $\text{gap} = p/k_{n}$, where $k_{n}$ is the normal penalty coefficient. +  * ''analyticalGap'': Second, the analytical computation provides the value of the gap via the following formula: $\text{gap} = p/k_{n}$, where $k_{n}$ is the normal penalty coefficient. 
-  * Thirdly, the semi-analytical/numerical computation consists in dividing the nodal contact force along $\mathbf{y}$, $f_{y}$, at $P_{4}$ by the nodal area $l_{x}l_{z}/n_{x}$ (where $l_{z} = \epsilon_{z}$, since the initial out-of-plane thickness is 1, and where $n_{x}$ is the number of elements along $\mathbf{x}$) and the normal penalty $k_{n}$: $\text{gap} = f_{y}/( k_{n}l_{z}l_{x}/(2n_{x}) )$+  * ''semiNumAnalyGap'': Thirdly, the semi-analytical/numerical computation consists in dividing the nodal contact force along $\mathbf{y}$, $f_{y}$, at $P_{4}$ by the nodal area $l_{x}l_{z}/n_{x}$ (where $l_{z} = \epsilon_{z}$, since the initial out-of-plane thickness is 1, and where $n_{x}$ is the number of elements along $\mathbf{x}$) and the normal penalty $k_{n}$: $\text{gap} = f_{y}/( k_{n}l_{z}l_{x}/(2n_{x}) )$ 
 + 
 +In all three cases, the gap should be the same.
commit/futur/db.1531820960.txt.gz · Last modified: 2018/07/17 11:49 by boemer