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
commit:2018:07_17 [2018/07/17 17:06] – [Test case] boemercommit:2018:07_17 [2018/07/17 19:10] (current) boemer
Line 1: Line 1:
 ====== Commit 2018-07-17 ====== ====== Commit 2018-07-17 ======
  
-===== Description =====+The general feature of this commit is described in the first section, then the corresponding test case, and finally, some minor modifications. 
 + 
 + 
 +===== 1. Description =====
  
 In this commit, the out-of-plane thickness evolution in the generalized plane strain state is coherently accounted for in the computation.  Previously, the possible increase or decrease of the out-of-plane thickness via the ''depend'' command was not taken into account in the contact area computation.  This led to an incoherently changing gap in the following scenario, where ''AIC_ONCEPERSTEP'' is activated. In this commit, the out-of-plane thickness evolution in the generalized plane strain state is coherently accounted for in the computation.  Previously, the possible increase or decrease of the out-of-plane thickness via the ''depend'' command was not taken into account in the contact area computation.  This led to an incoherently changing gap in the following scenario, where ''AIC_ONCEPERSTEP'' is activated.
Line 10: Line 13:
  
  
-===== Test case =====+===== 2. 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) 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). 
 + 
 + 
 +===== 3. 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.1531839962.txt.gz · Last modified: 2018/07/17 17:06 by boemer

Donate Powered by PHP Valid HTML5 Valid CSS Driven by DokuWiki