Homec4science
Diffusion LibMultiScale (lammps_2020) Browse

libmultiscale/examples/atc_xiao_Al/lm.config
95e94036595blammps_2020

Last Change
Raw File

lm.config
View Options

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 AtomsUnits
# 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 AtomsUnits
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
c4science · Help