explosion and boundary
Brought to you by:
ampuero
Menu
▾
▴
Source time function: You are free to use any STF for your simulation. I don't have the expertise to advice which one is more apropriate for a real explosion.
Half-round internal boundary: there are two possible approaches (see section 3 of SEM2DPACK's user manual):
a. Adapted mesh: Generate a mesh that adapts to the circular interface. You can do that with a mesh generation software (CUBIT, EMC2, EZ4U, etc).
b. Non-adapted mesh: Use a heterogeneous crustal model (see input block DIST_HETE1) and a regular cartesian mesh (see input block MESH_CART). In this approach the mesh does not conform to the material boundary.
The "adapted mesh" approach is the most accurate but also the most time-consuming during mesh generation. Running SEM2DPACK takes similar time in both approaches.
Last edit: Jean-Paul Ampuero 2013-05-08
&DIST_HETE1 file col . The parameters of file is ncol ( the number of date) ,and nx ,nz ,xo,zo(bottom-left corner). I want obtain a half-round through DIST_HETE1 , how should I find the ncol ,xo,zo ? And how should I define the two parts' material properties ?
Last edit: hanjianbin 2013-05-09
What is the DIST_*** used for ? Which part does it belong ? Mesh , BC or what else ??
You should use &DIST_HETE1 to set the values of the density and elastic properties throughout the medium in &MAT_ELASTIC. There is an example in the directory EXAMPLES/InaBox. Look at the "Example 3" in the Par.inp file. It uses a file "material.tab" that contains density and velocity values of a random medium, it was created by the matlab script make_material_file.m. You can modify that script to create a medium with a semi-circular inclusion, instead of a random medium.
Thank you for your guide . I have know the &DIST_HETE1. But my model is very different from Example 3 ,it only has one material . This is my model. It is made of two kind of material . I want to ask how to define the material properties and how to mesh the element .
which variable obtain the row of date from the file "material.tab" ? such as the first row of date is distributed to the first element or other element??
Last edit: hanjianbin 2013-05-09
why does that I install Emc2 that is downloaded from the internet is always fail ?
I cannot help you install EMC2. You should ask the EMC2 developers about that.
Please tell me which of the two approaches you want to use: (1) adapted mesh or (2) non-adapted mesh? (these two are explained in my first reply)
can the material be mesh by different nx and nz ? that is to say ,in a part ,I use nx=10 , in other part nx= 5 ,so that the result is more accuracy
Last edit: hanjianbin 2013-05-16
No, nx has to be the same in all layers.
dear Jean-Paul,
why can't I build my model with QSPLINE and spline ? And the example of spline in SEM2DPACK/EXAMPLES/TestLayers is not executed. can you send me some examples about QSPLINE and spline ?