Section MultiScale MetalUnits

UNITCODE MetalUnits
DIMENSION 3
INCLUDE lm.config Geometries

LET WaveLength = 10
LET NBELEM = 120


###### model declaration ######
ELAST_CODE AKANTU fe1
ELAST_CODE AKANTU fe2
MD_CODE LAMMPS md 

###### couplings ######
COUPLING_CODE lbdm md fe1 XIAO GEOMETRY bridge1 BOUNDARY pad1 \
	      GRID_DIVISION 2 2 2 QUALITY 1e-3
COUPLING_CODE rbdm md fe2 XIAO GEOMETRY bridge2 BOUNDARY pad2 \
	      GRID_DIVISION 2 2 2 QUALITY 1e-3

###### stimulations ######
FILTER central_md GEOM INPUT md GEOMETRY md
STIMULATION impulse IMPULSE INPUT central_md LWAVE WaveLength*AX DIRECTION 1 INTENSITY 1.1e-3 ONESHOT 0

###### dumpers ######
DUMPER paramd PARAVIEW INPUT md DISP VEL FORCE FREQ 100 PREFIX ./ 
DUMPER parafe1 PARAVIEW INPUT fe1 DISP VEL FORCE FREQ 100 PREFIX ./ 
DUMPER parafe2 PARAVIEW INPUT fe2 DISP VEL FORCE FREQ 100 PREFIX ./ 

DUMPER unmatched-point-rbdm-boundary PARAVIEW DISP VEL FORCE ONESHOT 0 PREFIX ./ \
       				     INPUT rbdm.unmatched-point-boundary 
DUMPER unmatched-point-rbdm-bridging PARAVIEW DISP VEL FORCE ONESHOT 0 PREFIX ./ \
       				     INPUT rbdm.unmatched-point-bridging 
DUMPER unmatched-point-lbdm-boundary PARAVIEW DISP VEL FORCE ONESHOT 0 PREFIX ./ \
       				     INPUT lbdm.unmatched-point-boundary
DUMPER unmatched-point-lbdm-bridging PARAVIEW DISP VEL FORCE ONESHOT 0 PREFIX ./ \
       				     INPUT lbdm.unmatched-point-bridging

DUMPER unmatched-fe-rbdm-boundary PARAVIEW DISP VEL FORCE ONESHOT 0 PREFIX ./ \
       				  INPUT rbdm.unmatched-mesh-boundary
DUMPER unmatched-fe-rbdm-bridging PARAVIEW DISP VEL FORCE ONESHOT 0 PREFIX ./ \
       				  INPUT rbdm.unmatched-mesh-bridging 
DUMPER unmatched-fe-lbdm-boundary PARAVIEW DISP VEL FORCE ONESHOT 0 PREFIX ./ \
       				  INPUT lbdm.unmatched-mesh-boundary 
DUMPER unmatched-fe-lbdm-bridging PARAVIEW DISP VEL FORCE ONESHOT 0 PREFIX ./ \
       				  INPUT lbdm.unmatched-mesh-bridging

DUMPER matched-point-rbdm-boundary PARAVIEW DISP VEL FORCE ONESHOT 0 PREFIX ./ \
       				   INPUT rbdm.matched-point-boundary
DUMPER matched-point-lbdm-boundary PARAVIEW DISP VEL FORCE ONESHOT 0 PREFIX ./ \
       				   INPUT lbdm.matched-point-boundary

DUMPER matched-fe-rbdm-boundary PARAVIEW DISP VEL FORCE ONESHOT 0 PREFIX ./ \
       				   INPUT rbdm.matched-mesh-boundary
DUMPER matched-fe-lbdm-boundary PARAVIEW DISP VEL FORCE ONESHOT 0 PREFIX ./ \
       				   INPUT lbdm.matched-mesh-boundary


DUMPER matched-point-rbdm-bridging PARAVIEW DISP VEL FORCE ONESHOT 0 PREFIX ./ \
       INPUT rbdm.matched-point-bridging \
       ADDFIELD rbdm.weight-point \
       ADDFIELD rbdm.lambdas-point
       
DUMPER matched-point-lbdm-bridging PARAVIEW DISP VEL FORCE ONESHOT 0 PREFIX ./ \
       INPUT lbdm.matched-point-bridging \
       ADDFIELD lbdm.weight-point \
       ADDFIELD lbdm.lambdas-point 

DUMPER matched-fe-rbdm-bridging PARAVIEW DISP VEL FORCE ONESHOT 0 PREFIX ./ \
       INPUT rbdm.matched-mesh-bridging \
       ADDFIELD rbdm.weight-mesh \
       ADDFIELD rbdm.lambdas-mesh 
       
DUMPER matched-fe-lbdm-bridging PARAVIEW DISP VEL FORCE ONESHOT 0 PREFIX ./ \
       INPUT lbdm.matched-mesh-bridging \
       ADDFIELD lbdm.weight-mesh \
       ADDFIELD lbdm.lambdas-mesh

endSection

################################################################
#geometry section
################################################################
Section Geometries RealUnits

# lattice and atomic constants
LET A = 4.04525975883
LET BX = sqrt(2)
LET BY = sqrt(3)
LET BZ = sqrt(6)
LET AX = BX*A
LET AY = BY*A
LET AZ = BZ*A
LET mass = 26.9815400000000

# geometry constants
LET WIDTH = 10
LET HEIGHT = 100
LET THICKNESS = 5
LET LX = WIDTH*AX
LET LY = HEIGHT*AY
LET LZ = THICKNESS*AZ

# MD geometry constants
LET MDHEIGHT = 30
LET MDLY = MDHEIGHT*AY

# FEM geometry constants
LET FEMHEIGHT = 10
LET FEMLY = FEMHEIGHT*AY

# bridging size
LET BRIDGE = 5*AY
LET PAD = 1.5*AY

# little epsilon to match dimensions
LET eps = A*0.1

#define explicitely the geometries
GEOMETRY md CUBE BBOX 0 LX -MDLY/2 MDLY/2 0 LZ 
GEOMETRY fe1 CUBE BBOX 0 LX -LY/2 -MDLY/2+BRIDGE 0 LZ
GEOMETRY fe2 CUBE BBOX 0 LX MDLY/2-BRIDGE LY/2 0 LZ

GEOMETRY bridge1 CUBE BBOX 0 LX+eps -MDLY/2+PAD -MDLY/2+BRIDGE 0 LZ+eps
GEOMETRY pad1 CUBE BBOX 0 LX+eps -MDLY/2 -MDLY/2+PAD+eps 0 LZ+eps

GEOMETRY bridge2 CUBE BBOX 0 LX+eps MDLY/2-BRIDGE MDLY/2-PAD  0 LZ+eps
GEOMETRY pad2 CUBE BBOX 0 LX+eps MDLY/2-PAD-eps MDLY/2 0 LZ+eps

endSection

################################################################
#MD section
################################################################
Section LAMMPS:md RealUnits
LAMMPS_FILE lammps.conf
DOMAIN_GEOMETRY md
TIMESTEP 1e-3
endSection

################################################################
#FE section1
################################################################
Section AKANTU:fe1 MetalUnits
DOMAIN_GEOMETRY fe1
MESH_FILENAME mesh.msh
MATERIAL_FILENAME material.dat
TIMESTEP 1e-3
SHIFT 0 -MDLY/2-FEMLY+BRIDGE 0
PBC 1 0 1
endSection

################################################################
#FE section2
################################################################
Section AKANTU:fe2 MetalUnits
DOMAIN_GEOMETRY fe2
MESH_FILENAME mesh.msh
MATERIAL_FILENAME material.dat
TIMESTEP 1e-3
SHIFT 0 MDLY/2-BRIDGE 0
PBC 1 0 1
endSection