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

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

Before contact After contact mechanics solution
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:

StatusTestedReal time (s)
hpc_cm64_irixSuccessThu Apr 1 12:01:32 2004519
hpc_cm_irixSuccessFri Aug 20 02:14:09 2004465
hpc_cmo64_irixSuccessSun Jul 31 01:44:22 200592
hpc_cmo_irixSuccessSun Aug 15 01:51:19 200464
rs6000-aix
cmSuccessSun Aug 15 02:50:27 200432
cm-debugSuccessSat Aug 21 04:57:25 20042808
cm64SuccessWed Mar 4 01:21:21 200929
cm64-debugSuccessThu Apr 1 10:48:10 2004588
x86_64-linux
cmSuccessSun Mar 6 00:01:48 201610

Testing status by file:


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

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


CMISS Help / Examples / d / d2 / d21