Metafor

ULiege - Aerospace & Mechanical Engineering

User Tools

Site Tools


commit:2018:07_17

Differences

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

Link to this comparison view

Both sides previous revisionPrevious revision
Next revision
Previous revision
Last revisionBoth sides next revision
commit:2018:07_17 [2018/07/17 17:06] – [Test case] boemercommit:2018:07_17 [2018/07/17 18:55] boemer
Line 1: Line 1:
 ====== Commit 2018-07-17 ====== ====== Commit 2018-07-17 ======
 +
 +The general feature of this commit is described in the first section, then the corresponding test case, and finally, some minor modifications.
 +
  
 ===== Description ===== ===== Description =====
Line 14: Line 17:
 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 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).
  
-{{ :commit:2018:07_17:contactoutofplanethicknessevolutiontest.png?400 |}}+{{ :commit:2018:07_17:contactoutofplanethicknessevolutiontest.png?300 |}}
  
 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: 
Line 21: Line 24:
   * ''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}) )$   * ''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.+In all three cases, the gap should be the same (besides errors due to finite precision computation). 
 + 
 + 
 +===== Minor modifications ===== 
 + 
 +  * The property ''NBR_CONTACT_INTERACTION'' was replaced by ''NUM_CONTACT_INTERACTION'' to make sure  that the specific number of a contact interaction is provided and not the total number; 
 +  * The address of clifton was replaced by the address of blueberry in ''stats.py''
 + 
 + 
 +===== Added [a] / deleted [d] / modified [m] / renamed [r] test cases ====== 
 + 
 +<code> 
 +[a] apps/gps/contThickEvol.py 
 +</code> 
 + 
 + 
 + --- // [[dominik.boemer@uliege.be|Dominik Boemer]] 2018/07/17 19:00 //
commit/2018/07_17.txt · Last modified: 2018/07/17 19:10 by boemer

Donate Powered by PHP Valid HTML5 Valid CSS Driven by DokuWiki