Example d21: Contact mechanics between a rigid block and deformable muscle

This example illustrates contact between a rigid block and deformable muscle.

Fig. 1 Before contact	Fig. 2 After contact mechanics solution

The comfile run by this example is as follows:

################################################

# Define parameters,regions,and coordinate system
fem def para;r;contact;example
fem def region;r;contact;example
fem def coord 3,1

#Read in deformable body
fem def nodes;r;rectus;example 
fem def bases;r;collapsed;example 
fem def elements;r;rectus;example

# Define fibre field (nj 4,5,6)
fem def fibre;r;rectus;example 
fem def elements;r;rectus;example fibre

# Export initial rectus muscle
fem exp nodes;rectus_ref as rectus_ref
fem exp elements;rectus_ref as rectus_ref 

# Read in target
fem def node;r;target;example region 2 
fem def element;r;target;example region 2

# Define target fibre field (nj 4,5,6)
fem def fibre;r;target;example reg 2 
fem def elements;r;target;example fibre reg 2
fem change node translate by 0,5,0 reg 2

# Export target
fem exp nodes;target as target region 2
fem exp elem;target as target region 2 

fem group nodes 29,30,31,32,40,1,2,3,4,33 as FIXALL reg 1
fem group nodes 39,27,38,24,37,19,36,15,35,11,34,7 as FIXXY reg 1

# Method
fem def equa;r;rectus;example
fem def mate;r;rectus;example 
fem def init;r;rectus;example 
  
# Method (Target)
fem def equa;r;target;example reg 2
fem def mate;r;target;example reg 2 
fem def init;r;target;example reg 2

# define contact parameter values
fem def contact;r;contact;example

# Define solve for region 1 only
fem def solve;r;rectus;example reg 1

$CONVERGED=0;
for($i=1;$i<=2;$i++){
if($i==2){ fem change nodes translate by 0,-5,0 reg 2 }
$CONVERGED=0;
while ($CONVERGED==0){ 

######################################Projection

# update XP from YP
FEM up geom from sol to 1..3 reg 1

# Set contact Xi points on specified faces as 8x8 points
fem def xi;c contact_points faces 55,58,61,66 points 8 reg 1

# Define data at xi positions 
fem def data;c from_xi

# Update slave info
fem update data field to slave

# Project onto target face (Region 2 face 91)
fem def xi;c closest_face faces 91 region 2
 
# store projection gap in data fields
fem update data field from gap region 2  

# Update master info
fem update data field to master 

# Place initial geometry YP(ny,3) into XP
FEM up geom from sol YP_index 3 to 1..3 reg 1

#####################################Mechanics problem

# Solve finite elasticity/contact problem
fem solve error 1.0D-03 iterate 22 reg 1

}
}
# Export final deformed mesh 
#fem exp nodes;rectus_def as rectus_def offset 300
#fem exp elements;rectus_def as rectus_def offset_elem 300
#fem exp nodes;target_def as target_def offset 400 reg 2
#fem exp elements;target_def as target_def offset_elem 400 reg 2

Additional testing commands:

Files used by this example are:

Name               Modified     Size

example_d21.com    23-Apr-2008  2.5k
collapsed.ipbase   21-Apr-2003  8.7k
contact.ipcont     25-Sep-2008  627 
contact.ippara     21-Feb-2008  6.0k
contact.ipregi     16-Apr-2003   93 
rectus.ipelem      16-Apr-2003   11k
rectus.ipelfb      16-Apr-2003  8.2k
rectus.ipequa      29-Sep-2004  2.2k
rectus.ipfibr      16-Apr-2003   30k
rectus.ipinit      02-May-2003   94k
rectus.ipmate      16-Apr-2003  5.6k
rectus.ipnode      16-Apr-2003   36k
rectus.ipsolv      16-Aug-2010  2.2k
rectus.ipsolv.old  13-Apr-2007  2.1k
target.ipelem      16-Apr-2003  548 
target.ipelfb      16-Apr-2003  411 
target.ipequa      29-Sep-2004  2.2k
target.ipfibr      16-Apr-2003   13k
target.ipinit      02-May-2003   37k
target.ipmate      16-Apr-2003  5.6k
target.ipnode      16-Apr-2003   15k
test_output.com    19-Apr-2003    0 

Download the entire example:

Name                      Modified     Size

examples_d_d2_d21.tar.gz  17-Aug-2010   41k

Testing status by version:

Testing status by file:

• cmiss_test.log
  

  Success	hpc_cm64_irix:	cmiss_test.log.retain.
  Success	hpc_cm_irix:	cmiss_test.log.retain.
  Success	hpc_cmo64_irix:	cmiss_test.log.retain.
  Success	hpc_cmo_irix:	cmiss_test.log.retain.
  rs6000-aix
  Success	cm:	cmiss_test.log.retain.
  Success	cm-debug:	cmiss_test.log.retain.
  Success	cm64:	cmiss_test.log.retain.
  Success	cm64-debug:	cmiss_test.log.retain.
  x86_64-linux
  Success	cm:	cmiss_test.log.retain.

  
  
• cmiss_test.out
  

  Success	hpc_cm64_irix:	ndiff test: no significant differences with generic answer.
  Success	hpc_cm_irix:	ndiff test: no significant differences with generic answer.
  Success	hpc_cmo64_irix:	ndiff test: no significant differences with generic answer.
  Success	hpc_cmo_irix:	ndiff test: no significant differences with generic answer.
  rs6000-aix
  Success	cm:	ndiff test: no significant differences with generic answer.
  Success	cm-debug:	ndiff test: no significant differences with generic answer.
  Success	cm64:	ndiff test: no significant differences with generic answer.
  Success	cm64-debug:	ndiff test: no significant differences with generic answer.
  x86_64-linux
  Success	cm:	ndiff test: no significant differences with generic answer.

  
  

Html last generated: Sun Mar 6 05:51:24 2016

Input last modified: Mon Aug 16 11:35:04 2010