6 AUTOMATIC REZONING AUTOMATIC REZONING Introduction
|
|
- Drusilla Henderson
- 5 years ago
- Views:
Transcription
1 AUTOMATIC REZONING AUTOMATIC REZONING 6.1 Introduction When FLAC is run in large strain mode, the grid may become distorted to such an extent that the simulation halts with a bad geometry error message. The objective of automatic rezoning logic is to map the existing stresses, velocities and displacements onto a new, more regular, grid, so that the run can continue. Because interpolation is used in the re-mapping, some smearing and accuracyloss is inevitable. In order to minimize these effects, the stresses are re-mapped for each subzone separately: thus, the overlap of a new subzone with each of the old subzones is performed in order to obtain a properly weighted estimate of the stress tensor in the new subzone. A progressively wider local search of possible overlaps is performed, in order to minimize search time and ensure that it is a linear process (time proportional to the number of zones). 6.2 Main Features In designing the automatic-rezoning logic, it is recognized that any FLAC model always uses some type of grid generator devised by the user. This may be as simple as a series of commands (e.g., GENERATE and INITIAL commands), or may be as complicated as a multi-line FISH function that produces a grid to conform to desired boundary shapes. Thus, the logic repeatedly invokes the user s grid generator within the rezone logic. When a REZONE command is issued, the following actions occur: 1. The model state is saved in a temporary area of memory. Specifically, subzone stresses, gridpoint velocities and displacements, material properties, fluid properties and pore pressure (if CONFIG gwflow is specified), temperature and thermal properties (if CONFIG thermal is specified) are saved. In addition, various conditions, such as fix flags, are saved. 2. A user-defined FISH function is called (from within the rezone logic). It is the task of this function to create a new grid. Even if the grid generator is a simple set of commands, these commands are embedded in the FISH function. (For technical reasons, a special form of embedding commands in a FISH function must be used; this is described below.) The form of the new grid is unrestricted, but normally it will be a version of the original grid that conforms to the new boundary geometry. 3. Stresses, velocities and displacements are mapped from the old grid to the new grid. In the rezoning logic, zone variables such as stress components are mapped, based on area overlapping fractions. For example, in the illustrative case shown in Figure 6.1, dash lines represent the new mesh, and solid lines define the old mesh. For the zone indicated by dash lines, its zone variables (e.g., stress components) are derived from the overlapping zones using area fractions (i.e., A1-A6) as the weighting factors. For nodal variables, such as velocity and displacement at a gridpoint in the new mesh, values are interpolated from the triangular subzone within which the gridpoint is located. The weighting factors are determined by the distance to each of the gridpoints of the host
2 6-2 Theory and Background subzone, as demonstrated in Figure 6.2. Such an interpolation scheme satisfies force and moment equilibrium. A detailed description of gridpoint weighting can be found in Section in Structural Elements. Figure 6.1 Zone variable mapping Figure 6.2 Gridpoint variable mapping 4. Advection of properties and state variables (e.g., plastic strains) is performed. Material properties are mapped as zone variables. This mapping is similar to that for stress components, except that stress components are re-mapped using subzone overlapping fractions, while material properties are re-mapped using full zone overlapping fractions.
3 AUTOMATIC REZONING 6-3 For integer-type properties, such as table number and plasticity indicator, a majority rule is adopted to determine the value for the new zone. 5. Histories used in rezoning are coordinate-based. Absolute locations (in x-y space) of gridpoint/zone histories are now stored and their indices in i-j space updated according to their absolute locations during the rezoning operation. 6. For dynamic simulations, the subzone masses are recalculated from the updated density. Simple tests show that momentum is conserved reasonably well. Normally a series of alternating STEP and REZONE commands will be issued (e.g., from a loop within a controlling FISH function). The rezoning logic can also be invoked automatically when a bad geometry error occurs by using the SET rez func command. At present, the rezone logic is not general: it does not recognize all features offered in FLAC. Support of such features will be added gradually over time. In particular, the logic excludes axisymmetry, interfaces, attach, multistepping and FISH models. Further, only one material model can be used in a rezoning run at present.
4 6-4 Theory and Background 6.3 Input Instructions for Automatic Rezoning FLAC Commands The following commands are provided to perform automatic rezoning during a FLAC run: REZONE <keyword value>... This command given without arguments invokes rezoning. The following keywords may also be given to set rezoning characteristics: set keyword history Location-based history advection is activated. method name The grid generator FISH function to be called within a rezone episode is specified. A function of the given name must exist. separate By default, each subzone receives weighted stresses from all old subzones that it overlaps. This results in some sharing (and smearing) of stress components between overlays (which are normally kept separate). This command forces subzones of a particular overlay to receive new stresses only from the same overlay in the old grid. smooth from i1,j1 to i2, j2 A surface range that will be subjected to a smoothing process is specified. The smoothing process consists of setting the coordinate of each surface gridpoint to the average of the three local gridpoints (the target gridpoint and its left and right neighbors). surface from i1,j1 to i2, j2 table n <seq ns> A surface range is specified. The table keyword specifies a table number to which the list of surface coordinates is written. This table may then be used by the user-defined FISH function to constrain the boundary of the new grid. If more than one surface is involved in the rezoning operation, the optional keyword seq is required to specify the sequence number (from 1 to 20) for each surface.
5 AUTOMATIC REZONING 6-5 SET rez func name The grid generator FISH function that is invoked automatically when a bad geometry error occurs is specified. A function of the given name must exist FISH Intrinsic Function The following FISH intrinsic function is provided to perform automatic rezoning during a FLAC run: rez exe(s) The intrinsic function rez exe (which always returns zero) causes the command contained in the given string, s, to be executed. Note that COMMAND...END COMMAND must not be used in a FISH function invoked within the rezoning logic (i.e., a function specified by the REZONE set method command). The commands that may be executed by rez exe( ) are limited to INITIAL, GENERATE, MODEL, PROPERTY and APPLY. If the string cannot fit on one line of FISH code, then it must be composed of separate concatenated sub-strings stored as a symbol (which is then given as the argument to rez exe( ) instead of being given directly as a quoted string).
6 6-6 Theory and Background 6.4 A Simple Rezoning Demonstration The automatic-rezoning capability is demonstrated for the case of a simple slope geometry. The slope is shown in Figure 6.3. The cohesive strength of the slope material is reduced to zero so that the slope fails. When run in large-strain mode, the run stops with a bad geometry error message after approximately 2000 calculation steps. JOB TITLE :. FLAC (Version 6.00) LEGEND Feb-08 13:10 step E+00 <x< 6.333E E+01 <y< 4.333E Grid plot 0 2E Itasca Consulting Group, Inc. Minneapolis, MN Figure 6.3 Simple slope initial geometry The data file for this model, including the rezoning facility, is listed in Example 6.1. The model is first stepped to an equilibrium state with a high cohesive strength assigned to the material. The cohesion is then set to zero, and the large-strain calculation mode is invoked. Automatic rezoning is implemented via the user-defined FISH function dostage slope. This function first calls the user-defined FISH function sloperezone with the command rezone set method sloperezone sloperezone is a simple grid generator composed of four commands: def sloperezone oo = rez exe( gen ) oo = rez exe( mod mohr ) oo = rez exe( gen tab 5 ) oo = rez exe( mod null reg 1 20 )
7 AUTOMATIC REZONING 6-7 The model surface is identified in dostage slope, and the table which will contain the list of surface coordinates is specified (id = 5) with the command rezone set surf from 1,11 to 61,21 tab 5 Then, a smoothing process is invoked to smooth the model surface with the command rezone set smooth from 1,11 to 61,21 The rezoning logic is implemented with the command rezone and then cycling is continued with the command step continue dostage slope is invoked automatically when a bad geometry error occurs by giving the command set rez func dostage slope Note that MOVIE commands are provided in the Example 6.1 data file in order to create a movie file of the distorting grid. The final shape of the slope after 10,000 steps is shown in Figure 6.4: JOB TITLE :. FLAC (Version 6.00) LEGEND Feb-08 13:10 step E+00 <x< 6.333E E+01 <y< 4.333E Grid plot 0 2E Itasca Consulting Group, Inc. Minneapolis, MN Figure 6.4 Simple slope geometry after 10,000 steps
8 6-8 Theory and Background Example 6.1 Automatic rezoning of a simple slope grid m mohr prop dens 2000 bulk 2e8 sh 1e8 fric 10 coh 1e20 tens 1e20 notnull table gen tab 1 mod null reg 1 20 fix x y j=1 fix x i=1 fix x i=61 set grav 10 solve save slope1.sav ; ini xv 0 yv 0 xdis 0 ydis 0 prop coh 0 set large def dostage slope command unmark rez set method=sloperezone rez set surf from 1,11 to 61,21 tab 5 rez set smooth from 1,11 to 61,21 rezone step continue Command def sloperezone oo = rez exe( gen ) oo = rez exe( mod mohr ) oo = rez exe( gen tab 5 ) oo = rez exe( mod null reg 1 20 ) set rez func dostage slope ; plot grid movie on movie step on 200 ; step save slope2.sav ret
9 AUTOMATIC REZONING Example Applications Failure of a Benched Slope A two-bench slope excavated in sensitive clay fails when the excavation level reaches 19 m depth. The initial geometry of the slope at this excavation depth is shown in Figure 6.5. A closeup view of the slope is shown in Figure 6.6. The properties of the clay are as follows: mass density (ρ) 2000 kg/m 3 bulk modulus (K) 1.7 MPa shear modulus (G) MPa initial cohesion (c) 50 kpa above 19 m depth 70 kpa below 19 m depth The clay exhibits a softening response under load. The behavior is simulated by using the strainsoftening constitutive model (MODEL ss) with the cohesion of the clay above the 19 m depth, degrading from 50 kpa to 5 kpa as a function of plastic shear strain. The data file for this model is listed in Example 6.2. An in-situ stress state is applied first, prior to excavation, with the vertical stress conforming to gravitational loading and the ratio of horizontal to vertical stress equal to Stresses are initialized using the FISH function inistress. The excavation is then made instantaneously (with the MODEL null command). The slope begins to fail. If the model is run without automatic rezoning, a bad geometry error occurs in zones at the toe of the slope after approximately 450 steps. Automatic rezoning is added to the model via the user-defined FISH function rezone, listed in Example 6.3. The rezoning operation involves four steps: 1. The gridpoints along the slope surface are stored in a table using the REZONE set surface from... to... table... command. Table 4 is designated to store the gridpoint data. 2. A new mesh is generated using the user-defined FISH function remesh. The slope surface is redefined in the new mesh by using the FISH function gentab4 to null zones above the surface defined by table 4. The slope surface is also smoothed with the REZONE set smooth command. 3. The rezoning logic is activated with the REZONE command to re-map zone and gridpoint values onto the new mesh. 4. The calculation stepping is continued. The rezone function is executed whenever the bad geometry error occurs by giving the command SET rez zone rezone before beginning the calculation.
10 6-10 Theory and Background JOB TITLE :. (*10^2) FLAC (Version 6.00) LEGEND 4-Mar-08 11:21 step E+01 <x< 5.389E E+01 <y< 1.159E+02 cohesion 5.000E E Itasca Consulting Group, Inc. Minneapolis, MN Figure 6.5 Benched slope initial slope geometry JOB TITLE :. FLAC (Version 6.00) LEGEND 4-Mar-08 12:25 step E+00 <x< 2.600E E+01 <y< 9.100E Grid plot 0 5E Itasca Consulting Group, Inc. Minneapolis, MN Figure 6.6 Benched slope closeup view
11 AUTOMATIC REZONING 6-11 Example 6.2 Benched slope grid 200,56 gen -20.0, , , ,70.0 model ss table 3-100,75 0,75 5,75 7.5, , , , ,84 100,84 gen tab 3 prop dens 2000 prop bulk 1.7e6 shear 0.625e6 prop cohesion = 45e3 ctab 99 j=21,56 table e e3 prop cohesion 70e3 ctab 100 j=1,20 table e e3 set gravity 10 def inistress loop i (1,izones) loop j (1,jzones) if model(i,j) # 1 depth = y(i,j) sy = -1. * * 10. * depth syy(i,j) = sy sxx(i,j) = sy * 0.65 szz(i,j) = sy * 0.65 if loop loop inistress fix x i 1 fix x y j 1 fix x i 201 hist 999 unbal solve save bench1.sav set large model null region 1,56 fix x y i 1 call rezone.fis save bench2.sav set rez func rezone set geom 0.21 step save bench3.sav
12 6-12 Theory and Background Example 6.3 Benched slope function def rezone command unmark rez set surfac from 1,21 to 201,57 table 4 rez set meth= remesh rez set smooth from 1,21 to 201,57 rezone step continue command ; def gentab4 loop ii (1,izones) loop jj (1,jzones) if model(ii,jj) # 1 x = (x(ii,jj)+x(ii+1,jj)+x(ii,jj+1)+x(ii+1,jj+1))/4. y = (y(ii,jj)+y(ii+1,jj)+y(ii,jj+1)+y(ii+1,jj+1))/4. yt = table(4, x) if y > yt oo = rez exe( model null i=ii j=jj ) if if loop loop ; def remesh oo=rez exe( gen -20.0, , , ,70.0 ) oo=rez exe( model ss ) gentab4 The failed slope geometry at the final equilibrium state is shown in Figure 6.7. A contour plot of cohesion at this state is given in Figure 6.8, and indicates the failure surface that has developed.
13 AUTOMATIC REZONING 6-13 JOB TITLE :. (*10^2) FLAC (Version 6.00) LEGEND 4-Mar-08 16:30 step E+01 <x< 5.389E E+01 <y< 1.159E Grid plot 0 2E Itasca Consulting Group, Inc. Minneapolis, MN Figure 6.7 Final state of failed bench slope JOB TITLE :. (*10^2) FLAC (Version 6.00) LEGEND 5-Mar-08 9:19 step E+01 <x< 5.389E E+01 <y< 1.159E+02 cohesion -1.00E E E E E E E E E+04 Contour interval= 1.00E+04 Boundary plot E Itasca Consulting Group, Inc. Minneapolis, MN Figure 6.8 Contour plot of cohesion for failed bench slope
14 6-14 Theory and Background Punch Problem A rigid punch is driven at a constant velocity into a cohesive, frictionless material. When run in large-strain calculation mode, a bad geometry error occurs after driving the punch 2.1 m into the material (approximately 10,000 steps). The automatic-rezoning logic allows the punch to be driven almost entirely through the material. The mesh is automatically re-created when a bad geometry condition is encountered, by invoking the rezoning logic. The initial grid is shown in Figure 6.9. The model is 55 m wide and 30 m deep. The punch area is 5 m wide at the top of the model. A constant velocity (denoted by the vectors in the figure) of m/step is specified to apply the loading over the punch area. The model has roller boundary conditions on the sides, and is fixed in both directions at the base. Note that the model is symmetric about a vertical line through the center of the punch, and the FLAC model could be reduced to half of the model shown in Figure 6.9. The whole model is run to illustrate the use of multiple tables to define the surface geometry when rezoning. JOB TITLE :. FLAC (Version 6.00) LEGEND 25-Feb-08 14:08 step E+00 <x< 5.806E E+01 <y< 4.556E Grid plot 0 1E Velocity vectors max vector = 2.000E E Itasca Consulting Group, Inc. Minneapolis, MN Figure 6.9 FLAC grid with constant-velocity punch
15 AUTOMATIC REZONING 6-15 The data file for this problem is listed in Example 6.4. The vertical velocity is applied and the model is first run for 4000 steps before the rezoning logic is invoked. The distorted grid at this stage is shown in Figure Example 6.4 Punch problem grid mod mohr prop dens 1000 bulk 2e8 sh 1e8 coh 1e5 tens 1e10 fix x i=1 fix x y j=1 fix x i=56 fix x y i=26,31 j=31 ini yvel -2e-4 i=26,31 j=31 save punch a.sav ; set large ncw=100 step 2000 save punch b.sav ; step 2000 save punch c.sav ; def inisetup width = 55.0 ntab1 = 1 ntab2 = 2 nz1 = jzones nz1p1 = nz1 + 1 irmax = igp jrmax = jgp jlmax = jgp inisetup ; call dostage punch.fis dostage punch save punch d.sav ; step save punch e.sav ret
16 6-16 Theory and Background JOB TITLE :. FLAC (Version 6.00) LEGEND 3-Mar-08 14:39 step E+00 <x< 5.806E E+01 <y< 4.567E+01 X-displacement contours -4.00E E E E E-01 Contour interval= 2.00E-01 Velocity vectors max vector = 2.013E E -4 Grid plot E Itasca Consulting Group, Inc. Minneapolis, MN Figure 6.10 x-displacement contours and velocity vectors after driving the punch 0.8 m before rezoning In order to continue the run when a bad geometry error occurs, a new mesh is generated by using the user-defined FISH function dostage punch, listed in Example 6.5. The rezoning process consists of four steps: 1. The gridpoints associated with the model surface are stored in tables using the REZONE set surface from... to... table... seq... command. Two tables are created: table 1 stores the gridpoints from the upper-left corner of the model to the left boundary gridpoint of the punch; and table 2 stores the gridpoints from the right boundary gridpoint of the punch to the upper-right corner. 2. A new mesh is generated using the user-defined FISH function footrezone with a new surface defined by tables 1 and 2. The mesh is also scanned to remove zones above the new surface. 3. The rezoning logic is activated with the REZONE command to re-map zone and gridpoint values onto the new mesh. 4. The punch velocity loading is fixed again and the calculation stepping is continued. The dostage punch function is executed whenever the bad geometry error occurs, by issuing the command set rez func dostage punch
17 AUTOMATIC REZONING 6-17 Example 6.5 Punch rezoning functions def dostage punch command unmark rez set surf from 1,jlmax to 26,nz1p1 tab ntab1 seq 1 rez set surf from 31,nz1p1 to irmax,jrmax tab ntab2 seq 2 rez set meth=footrezone rezone free x y i 2 55 j 2 31 fix x y i 26,31 j=nz1p1 step continue Command ; def footrezone h1 = y(26,nz1+1) tbs2 = table size(ntab2) h2 = 0.0 loop n (1,tbs2) h2 = max(h2,ytable(ntab2,n)) Loop nz1 = int(h1 * float(jzones) / h2) hinc = h1 / float(nz1) ytop = float(jzones) * hinc nz1p1 = nz1 + 1 nztop = jzones oo = rez exe( gen 0,0 0,ytop width,ytop width,0 ) oo = rez exe( mod mohr ) oo = rez exe( mod null i 26,30 j nz1p1,nztop ) oo = rez exe( gen tab ntab1 ) oo = rez exe( gen tab ntab2 ) loop i (1,26) loop j (nz1+1,jzones) if model(i,j) # 1 xcen = (x(i,j)+x(i+1,j)+x(i,j+1)+x(i+1,j+1))/4.0 ycen = (y(i,j)+y(i+1,j)+y(i,j+1)+y(i+1,j+1))/4.0 if ycen > table(ntab1,xcen) ii = i jj = j oo = rez exe( mod null i=ii j=jj ) if if Loop Loop
18 6-18 Theory and Background loop i (31,izones) loop j (nz1+1,jzones) if model(i,j) # 1 xcen = (x(i,j)+x(i+1,j)+x(i,j+1)+x(i+1,j+1))/4.0 ycen = (y(i,j)+y(i+1,j)+y(i,j+1)+y(i+1,j+1))/4.0 if ycen > table(ntab2,xcen) ii = i jj = j oo = rez exe( mod null i=ii j=jj ) if if Loop Loop jlmax = jgp section loop j (1,jzones) if model(1,j) = 1 jlmax = j exit section if Loop Section y(igp,jrmax) = table(ntab2,x(igp,jrmax)) jrmax = jgp section loop j (1,jzones) if model(izones,j) = 1 jrmax = j exit section if Loop Section y(igp,jrmax) = table(ntab2,x(igp,jrmax)) set rez func dostage punch Figures 6.11 and 6.12 display the distorted grid after driving the punch 13.8 m and 24.8 m, respectively, into the mesh.
19 AUTOMATIC REZONING 6-19 JOB TITLE :. FLAC (Version 6.00) LEGEND 25-Feb-08 15:08 step E+00 <x< 5.806E E+01 <y< 4.701E+01 X-displacement contours -2.00E E E E E E E E E+00 Contour interval= 5.00E-01 Grid plot 0 1E 1 Velocity vectors max vector = 2.063E E -4 Itasca Consulting Group, Inc. Minneapolis, MN Figure 6.11 x-displacement contours and velocity vectors after driving the punch 13.8 m JOB TITLE :. FLAC (Version 6.00) LEGEND 25-Feb-08 15:07 step E+00 <x< 5.806E E+01 <y< 4.754E+01 X-displacement contours -3.00E E E E E E E+00 Contour interval= 1.00E+00 Grid plot 0 1E 1 Velocity vectors max vector = 2.648E E Itasca Consulting Group, Inc. Minneapolis, MN Figure 6.12 x-displacement contours and velocity vectors after driving the punch 24.8 m
20 6-20 Theory and Background 6.6 Accuracy and Zoning Sensitivity The accuracy of the rezoning logic and the sensitivity of the logic to zone density is tested in the following sections Accuracy A simple 2 2 uniform square grid is used to verify the correctness of implementation of the rezoning logic. As illustrated in Figure 6.13, each zone has a size of unity, but different density (e.g., 1000 at (1,1), 2000 at (1,2), 3000 at (2,1) and 4000 at (2,2)) before rezoning. The y-displacement of the top gridpoints (j = 3) and middle gridpoints (j = 2) are also initialized to different values (1e-3 at j = 2, 2e-3 at j = 3). In the rezoning operation, the line of i = 2 is translated in the x-direction by α (= 0.2), and the line of j = 2 is translated in the y-direction by β (= 0.3) to form a new grid. The re-meshing logic is then invoked to map zone quantities and interpolate gridpoint quantities from the old grid to the new grid. The data file is listed in Example 6.6. The comparison made in FISH function check shows that the expected mapped values of density in zones (1,1) and (2,1) and vertical displacement at gridpoint (2,2) exactly match the actual values provided by FLAC. Figure 6.13 Rezoning of a 2 2 grid: solid lines are the original grid, and dashed lines are the new grid. The original and new grids share the same boundary.
21 AUTOMATIC REZONING 6-21 Example 6.6 Accuracy test new def const ; prop d11 = d12 = d21 = d22 = yd2 = 1e-3 yd3 = 2e-3 ; geom x1 = 0. x2 = 2. y1 = 0. y2 = 2. alpha = 0.2 beta = 0.3 const grid 2 2 model elastic prop dens d11 i 1 j 1 prop dens d12 i 1 j 2 prop dens d21 i 2 j 1 prop dens d22 i 2 j 2 ini ydisp yd2 j 2 ini ydisp yd3 j 3 def remesh oo=rez exe( ini x add alpha i = 2 ) oo=rez exe( ini y add beta j = 2 ) def rezone command rez set method = remesh rezone command save accu grid.sav rezone ; note that the formulations for expected values may be different ; if alpha and/or beta are changed. def check area = x2 * y2 / 4.
22 6-22 Theory and Background ; check density (1,1) frac11 = 1./ area frac21 = alpha/ area frac12 = beta/ area frac22 = alpha * beta/ area val = d11* frac11+ d21* frac21+ d12* frac12+ d22* frac22 exp d11 = val/( frac11+ frac21+ frac12+ frac22) num d11 = density(1,1) ; check density (2,1) frac21 = (2.- alpha)/ area frac22 = (2.- alpha)* beta/ area val = d21* frac21+ d22* frac22 exp d21 = val/( frac21+ frac22) num d21 = density(2,1) ; check ydisp(2,2) exp yd22 = yd3* beta+ yd2*(1.- beta) num yd22 = ydisp(2,2) check print exp d11 num d11 print exp d21 num d21 print exp yd22 num yd Zoning Sensitivity As illustrated in Figure 6.14(a), a rectangular block filled with two different materials is fixed in both directions at all boundaries. The gridpoints at the material interface are forced to displace upward following a given sine-shaped profile. Once the system evolves to the position shown in Figure 6.14(b) (which takes 5000 steps), the model is stepped for another 5000 steps with the velocities at all gridpoints reversed. Figure 6.14(c) shows the contour plot of material property at step 10,000. If a re-meshing action is made just before reversing velocities, the material property contour will be shown in Figure 6.14(d); some smearing of material property is induced in the re-meshing process. The test is performed using a grid. The smearing is found to be significantly reduced in finer grids, for example, as shown in Figure 6.14(e) for a grid, and Figure 6.14(f) for a grid.
23 AUTOMATIC REZONING 6-23 Figure 6.14 Material advection testing; dark and light areas indicate the same material property with different values. State (c) is for no rezoning; states (d), (e) and (f) are for rezoning in progressively finer grids. Example 6.7 Zone density sensitivity new def const ; prop den1 = bul1 = 1e8 she1 = 3e7 den2 = bul2 = 2e8 she2 = 6e7 ; geom (dimensions) x1 = 0. x2 = 50. y1 = 0. y2 = 30. ; geom (IJ) iz = 25;50;100 jz = 15;30;600 ig = iz+1 jg = jz+1 jm = jz/2
24 6-24 Theory and Background jm1 = jm+1 const grid iz jz gen ( x1, y1) ( x1, y2) ( x2, y2) ( x2, y1) model elastic prop dens den1 bulk bul1 shear she1 j = 1, jm prop dens den2 bulk bul2 shear she2 j = jm1, jz def setsine j = jm1 loop i(2, iz) yvel(i,j) = amp1 * sin(pi * x(i,j)/x( ig,j)) loop def genmesh oo = rez exe( gen x1, y1 x1, y2 x2, y2 x2, y1 ) def rezone command rezone set method genmesh rezone command set amp1 = 1e-3 setsine fix x y j=1 fix x y j= jg fix x y i=1 fix x y i= ig set large ncw=100 fix y j= jm1 step 5000 plot hold bulk fill fix rezone plot hold bulk fill fix ini yvel mul -1 j= jm1 step 5000 plot hold bulk fill fix
4 GRID GENERATION GRID GENERATION General Comments
GRID GENERATION 4-1 4 GRID GENERATION 4.1 General Comments Unlike many modeling programs based on the finite element method, FLAC organizes its zones (or elements ) in a row-and-column fashion, like a
More informationFISH BEGINNER S GUIDE FISH BEGINNER S GUIDE. 1.1 Introduction
FISH BEGINNER S GUIDE 1-1 1 FISH BEGINNER S GUIDE 1.1 Introduction FISH is a programming language embedded within FLAC that enables the user to define variables and functions. These functions may be used
More informationSSR Polygonal Search Area
SSR Polygonal Search Area 22-1 SSR Polygonal Search Area In this tutorial, Phase2 is used to determine the factor of safety of a slope using the shear strength reduction (SSR) method. The SSR Polygon Search
More information5 MISCELLANEOUS MISCELLANEOUS FLAC Runtime Benchmark
MISCELLANEOUS 5-1 5 MISCELLANEOUS 5.1 FLAC Runtime Benchmark FLAC has been tested on a number of different computers. The calculation rates are compared here for a 9684-zone model of Mohr-Coulomb material
More informationFLAC/Slope User s Guide Contents - 1
FLAC/Slope User s Guide Contents - 1 TABLE OF CONTENTS 1 FLAC/SLOPE 1.1 Introduction... 1-1 1.1.1 Overview... 1-1 1.1.2 Guide to the FLAC/Slope Manual... 1-2 1.1.3 Summary of Features... 1-2 1.1.4 Analysis
More information3DEC 3DEC ITASCA ITASCA VERSION 5.0 VERSION 5.0
ITASCA Consulting Group, Inc. An Itasca International Company 3DEC VERSION 5.0 Advanced, Three Dimensional Discrete Element Modeling for Geotechnical Analysis of Rock, Blocky Structures and Structural
More informationGEOSTUDIO Tutorials Results and Procedure Comparison
GEOSTUDIO Tutorials Results and Procedure Comparison Angel Francisco Martinez Application Engineer MIDAS IT Estados Unidos Integrated Solver Optimized for the next generation 64-bit platform Finite Element
More informationFISH BEGINNER S GUIDE FISH BEGINNER S GUIDE. 1.1 Introduction and Overview
FISH BEGINNER S GUIDE 1-1 1 FISH BEGINNER S GUIDE 1.1 Introduction and Overview FISH is a programming language embedded within UDEC that enables the user to define variables and functions. These functions
More information12 Blocks Bouncing down Slope
Blocks Bouncing down Slope 12-1 12 Blocks Bouncing down Slope 12.1 Problem Statement Engineers are sometimes required to design slopes that must prevent loose rocks from falling onto roadways. The slopes
More informationLab Practical - Boundary Element Stress Analysis
Lab Practical - Boundary Element Stress Analysis Part A The Basics In this example we will model a simple set of excavations using the Rocscience program EXAMINE 2D. The first step is to define the geometry
More informationExample Application 4. Undrained Cylindrical Cavity Expansion in a Cam-Clay Medium
Example Application 4 Undrained Cylindrical Cavity Expansion in a Cam-Clay Medium 1 Cylindrical Cavity Model Geometry 2 Saturated Clay Properties Cam-Clay Model 3 Modeling Procedure Step 1 Step 2 Step
More informationCHAPTER 4. Numerical Models. descriptions of the boundary conditions, element types, validation, and the force
CHAPTER 4 Numerical Models This chapter presents the development of numerical models for sandwich beams/plates subjected to four-point bending and the hydromat test system. Detailed descriptions of the
More informationA Novel Approach to High Speed Collision
A Novel Approach to High Speed Collision Avril Slone University of Greenwich Motivation High Speed Impact Currently a very active research area. Generic projectile- target collision 11 th September 2001.
More information1.992, 2.993, 3.04, 10.94, , Introduction to Modeling and Simulation Prof. F.-J. Ulm Spring FE Modeling Example Using ADINA
1.992, 2.993, 3.04, 10.94, 18.996, 22.091 Introduction to Modeling and Simulation Prof. F.-J. Ulm Spring 2002 FE Modeling Example Using ADINA H Hgρ w ργ H = B = 10 m g = 9.81 m/s 2 ρ = 2400 kg/m 3 ρ w
More informationTerrain settlement analysis
Engineering manual No. 21 Updated: 02/2018 Terrain settlement analysis Program: File: FEM Demo_manual_21.gmk This example contains the solution to terrain settlement under surcharge loading using the Finite
More informationLab Practical - Finite Element Stress & Deformation Analysis
Lab Practical - Finite Element Stress & Deformation Analysis Part A The Basics In this example, some of the basic features of a finite element analysis will be demonstrated through the modelling of a simple
More informationSimulation of Connector Assembly AA
Simulation of Connector Assembly AA Date: Tuesday, March 1, 2016 Designer: Solidworks Study name: Horizontal Stress in AA inner tab fold Analysis type: Static Table of Contents Model Information... 2 Study
More informationSETTLEMENT OF A CIRCULAR FOOTING ON SAND
1 SETTLEMENT OF A CIRCULAR FOOTING ON SAND In this chapter a first application is considered, namely the settlement of a circular foundation footing on sand. This is the first step in becoming familiar
More informationPrescribed Deformations
u Prescribed Deformations Outline 1 Description 2 Finite Element Model 2.1 Geometry Definition 2.2 Properties 2.3 Boundary Conditions 2.3.1 Constraints 2.3.2 Prescribed Deformation 2.4 Loads 2.4.1 Dead
More informationSimulation of Connector Assembly C
Simulation of Connector Assembly C Date: Sunday, March 6, 2016 Designer: Solidworks Study name: Horizontal Stress Test on C inner bend Analysis type: Static Table of Contents Model Information... 2 Study
More informationLearning Module 8 Shape Optimization
Learning Module 8 Shape Optimization What is a Learning Module? Title Page Guide A Learning Module (LM) is a structured, concise, and self-sufficient learning resource. An LM provides the learner with
More informationChapter 7 Practical Considerations in Modeling. Chapter 7 Practical Considerations in Modeling
CIVL 7/8117 1/43 Chapter 7 Learning Objectives To present concepts that should be considered when modeling for a situation by the finite element method, such as aspect ratio, symmetry, natural subdivisions,
More informationStructural static analysis - Analyzing 2D frame
Structural static analysis - Analyzing 2D frame In this tutorial we will analyze 2D frame (see Fig.1) consisting of 2D beams with respect to resistance to two different kinds of loads: (a) the downward
More informationFinite Element Method. Chapter 7. Practical considerations in FEM modeling
Finite Element Method Chapter 7 Practical considerations in FEM modeling Finite Element Modeling General Consideration The following are some of the difficult tasks (or decisions) that face the engineer
More informationVOLCANIC DEFORMATION MODELLING: NUMERICAL BENCHMARKING WITH COMSOL
VOLCANIC DEFORMATION MODELLING: NUMERICAL BENCHMARKING WITH COMSOL The following is a description of the model setups and input/output parameters for benchmarking analytical volcanic deformation models
More informationBonded Block Modeling of a Tunnel Excavation with Support
Bonded Block Modeling of a Tunnel Excavation with Support 1 Modeling Procedure Step 1 Step 2 Step 3 Step 4 Step 5 Step 6 Generate the blocks to define the tunnel geometry Assign properties to blocks and
More informationWorkshop 15. Single Pass Rolling of a Thick Plate
Introduction Workshop 15 Single Pass Rolling of a Thick Plate Rolling is a basic manufacturing technique used to transform preformed shapes into a form suitable for further processing. The rolling process
More informationTHERMAL EXPANSION OF A NAVIGABLE LOCK
THERMAL EXPANSION OF A NAVIGABLE LOCK 15 THERMAL EXPANSION OF A NAVIGABLE LOCK A navigable lock is temporarily 'empty' due to maintenance. After some time there is significant increase of the air temperature,
More informationANSYS Workbench Guide
ANSYS Workbench Guide Introduction This document serves as a step-by-step guide for conducting a Finite Element Analysis (FEA) using ANSYS Workbench. It will cover the use of the simulation package through
More informationME 475 FEA of a Composite Panel
ME 475 FEA of a Composite Panel Objectives: To determine the deflection and stress state of a composite panel subjected to asymmetric loading. Introduction: Composite laminates are composed of thin layers
More informationDeveloper s Tips. Groundwater analysis using Slide
Developer s Tips Groundwater analysis using Slide In this article we introduce the new groundwater/seepage analysis module, which was built into version 5.0 of Slide. The groundwater module uses a finite
More informationStability with FE Stresses
1 Introduction Stability with FE Stresses In the finite element stresses stability analysis method, the stresses in the ground can be computed using SIGMA/W, and then SLOPE/W uses the SIGMA/W stresses
More informationChapter 20. Finite Element Method (FEM) Introduction
Chapter 20. Finite Element Method (FEM) Introduction The objective of this chapter is to explain basic terms of the particular field of problems and practical application of GEO 5 FEM program to solve
More informationTekla Structures Analysis Guide. Product version 21.0 March Tekla Corporation
Tekla Structures Analysis Guide Product version 21.0 March 2015 2015 Tekla Corporation Contents 1 Getting started with analysis... 7 1.1 What is an analysis model... 7 Analysis model objects...9 1.2 About
More informationStructural static analysis - Analyzing 2D frame
Structural static analysis - Analyzing 2D frame In this tutorial we will analyze 2D frame (see Fig.1) consisting of 2D beams with respect to resistance to two different kinds of loads: (a) the downward
More informationModule 1.2: Moment of a 1D Cantilever Beam
Module 1.: Moment of a 1D Cantilever Beam Table of Contents Page Number Problem Description Theory Geometry Preprocessor 6 Element Type 6 Real Constants and Material Properties 7 Meshing 9 Loads 10 Solution
More informationFOUNDATION IN OVERCONSOLIDATED CLAY
1 FOUNDATION IN OVERCONSOLIDATED CLAY In this chapter a first application of PLAXIS 3D is considered, namely the settlement of a foundation in clay. This is the first step in becoming familiar with the
More informationModeling Foundations in RS
Modeling Foundations in RS 3 Piled Raft Modeling in RS 3 Deep foundation piles are commonly used to increase foundation stability and to increase the bearing capacity of structural systems. The design
More informationRevised Sheet Metal Simulation, J.E. Akin, Rice University
Revised Sheet Metal Simulation, J.E. Akin, Rice University A SolidWorks simulation tutorial is just intended to illustrate where to find various icons that you would need in a real engineering analysis.
More informationLagrangian methods and Smoothed Particle Hydrodynamics (SPH) Computation in Astrophysics Seminar (Spring 2006) L. J. Dursi
Lagrangian methods and Smoothed Particle Hydrodynamics (SPH) Eulerian Grid Methods The methods covered so far in this course use an Eulerian grid: Prescribed coordinates In `lab frame' Fluid elements flow
More informationSUBMERGED CONSTRUCTION OF AN EXCAVATION
2 SUBMERGED CONSTRUCTION OF AN EXCAVATION This tutorial illustrates the use of PLAXIS for the analysis of submerged construction of an excavation. Most of the program features that were used in Tutorial
More informationNew developments in numerical modelling of pile installation
New developments in numerical modelling of pile installation Nguyen Phuong, Frits van Tol, Alexander Rohe 18 September 2014 KIVI Geotechnical Lectures Evening TU Delft Displacement piles à installation
More informationTips about Springback and compensation with ETA/Dynaform. DYNAFORM Team June, 2015
Tips about Springback and compensation with ETA/Dynaform DYNAFORM Team June, 2015 1 Simulation Basics 2 Simulation Basics! Mesh! Implicit and Explicit! Time step! Contact! Material Model " Yielding Surfaces
More informationComputer Life (CPL) ISSN: Finite Element Analysis of Bearing Box on SolidWorks
Computer Life (CPL) ISSN: 1819-4818 Delivering Quality Science to the World Finite Element Analysis of Bearing Box on SolidWorks Chenling Zheng 1, a, Hang Li 1, b and Jianyong Li 1, c 1 Shandong University
More informationApplication of Finite Volume Method for Structural Analysis
Application of Finite Volume Method for Structural Analysis Saeed-Reza Sabbagh-Yazdi and Milad Bayatlou Associate Professor, Civil Engineering Department of KNToosi University of Technology, PostGraduate
More informationA Sliding Block Analysis
1 Introduction A Sliding Block Analysis This example simulates the sliding of a block on a rigid surface. The primary purpose of this example is to check the behavior of the slip (interface) elements in
More informationThe part to be analyzed is the bracket from the tutorial of Chapter 3.
Introduction to Solid Modeling Using SolidWorks 2007 COSMOSWorks Tutorial Page 1 In this tutorial, we will use the COSMOSWorks finite element analysis (FEA) program to analyze the response of a component
More informationCIVIL MANUFACTURING MINING OIL & GAS POWER GENERATION
CIVIL MANUFACTURING MINING OIL & GAS POWER GENERATION ICG14-BRO-UDEC600-US-3 ABOUT UDEC The Universal Distinct Element Code (UDEC) is two-dimensional numerical software that simulates the response of loading
More informationAppendix B: Creating and Analyzing a Simple Model in Abaqus/CAE
Getting Started with Abaqus: Interactive Edition Appendix B: Creating and Analyzing a Simple Model in Abaqus/CAE The following section is a basic tutorial for the experienced Abaqus user. It leads you
More informationSTATIC AND DYNAMIC 2-DIMENSIONAL FINITE ELEMENT ANALYSIS OF CONTINUA
VERSAT-S2D & VERSAT-D2D Version 2016.6.18 STATIC AND DYNAMIC 2-DIMENSIONAL FINITE ELEMENT ANALYSIS OF CONTINUA - USING WINDOWS XP, WINDOWS 7 & Windows 10 Volume 2: USER MANUAL 1998-2018.03 Dr. G. WU 1998-2018.03
More informationExample 24 Spring-back
Example 24 Spring-back Summary The spring-back simulation of sheet metal bent into a hat-shape is studied. The problem is one of the famous tests from the Numisheet 93. As spring-back is generally a quasi-static
More informationCOMPUTER AIDED ENGINEERING. Part-1
COMPUTER AIDED ENGINEERING Course no. 7962 Finite Element Modelling and Simulation Finite Element Modelling and Simulation Part-1 Modeling & Simulation System A system exists and operates in time and space.
More informationPLAXIS 2D - SUBMERGED CONSTRUCTION OF AN EXCAVATION
PLAXIS 2D - SUBMERGED CONSTRUCTION OF AN EXCAVATION 3 SUBMERGED CONSTRUCTION OF AN EXCAVATION This tutorial illustrates the use of PLAXIS for the analysis of submerged construction of an excavation. Most
More informationSet No. 1 IV B.Tech. I Semester Regular Examinations, November 2010 FINITE ELEMENT METHODS (Mechanical Engineering) Time: 3 Hours Max Marks: 80 Answer any FIVE Questions All Questions carry equal marks
More informationModule 1.7: Point Loading of a 3D Cantilever Beam
Module 1.7: Point Loading of a D Cantilever Beam Table of Contents Page Number Problem Description Theory Geometry 4 Preprocessor 6 Element Type 6 Material Properties 7 Meshing 8 Loads 9 Solution 15 General
More information3D Finite Element Software for Cracks. Version 3.2. Benchmarks and Validation
3D Finite Element Software for Cracks Version 3.2 Benchmarks and Validation October 217 1965 57 th Court North, Suite 1 Boulder, CO 831 Main: (33) 415-1475 www.questintegrity.com http://www.questintegrity.com/software-products/feacrack
More informationMSC/PATRAN LAMINATE MODELER COURSE PAT 325 Workbook
MSC/PATRAN LAMINATE MODELER COURSE PAT 325 Workbook P3*V8.0*Z*Z*Z*SM-PAT325-WBK - 1 - - 2 - Table of Contents Page 1 Composite Model of Loaded Flat Plate 2 Failure Criteria for Flat Plate 3 Making Plies
More informationThis tutorial illustrates how to set up and solve a problem involving solidification. This tutorial will demonstrate how to do the following:
Tutorial 22. Modeling Solidification Introduction This tutorial illustrates how to set up and solve a problem involving solidification. This tutorial will demonstrate how to do the following: Define a
More information3 SETTLEMENT OF A CIRCULAR FOOTING ON SAND (LESSON 1) Figure 3.1 Geometry of a circular footing on a sand layer
SETTLEMENT OF A CIRCULAR FOOTING ON SAND (LESSON 1) 3 SETTLEMENT OF A CIRCULAR FOOTING ON SAND (LESSON 1) In the previous chapter some general aspects and basic features of the PLAXIS program were presented.
More informationCreating and Analyzing a Simple Model in Abaqus/CAE
Appendix B: Creating and Analyzing a Simple Model in Abaqus/CAE The following section is a basic tutorial for the experienced Abaqus user. It leads you through the Abaqus/CAE modeling process by visiting
More informationRevision of the SolidWorks Variable Pressure Simulation Tutorial J.E. Akin, Rice University, Mechanical Engineering. Introduction
Revision of the SolidWorks Variable Pressure Simulation Tutorial J.E. Akin, Rice University, Mechanical Engineering Introduction A SolidWorks simulation tutorial is just intended to illustrate where to
More informationChapter 3 Analysis of Original Steel Post
Chapter 3. Analysis of original steel post 35 Chapter 3 Analysis of Original Steel Post This type of post is a real functioning structure. It is in service throughout the rail network of Spain as part
More informationCFD Post-Processing of Rampressor Rotor Compressor
Gas Turbine Industrial Fellowship Program 2006 CFD Post-Processing of Rampressor Rotor Compressor Curtis Memory, Brigham Young niversity Ramgen Power Systems Mentor: Rob Steele I. Introduction Recent movements
More informationTWO-DIMENSIONAL PROBLEM OF THE THEORY OF ELASTICITY. INVESTIGATION OF STRESS CONCENTRATION FACTORS.
Ex_1_2D Plate.doc 1 TWO-DIMENSIONAL PROBLEM OF THE THEORY OF ELASTICITY. INVESTIGATION OF STRESS CONCENTRATION FACTORS. 1. INTRODUCTION Two-dimensional problem of the theory of elasticity is a particular
More information1. The program has automatic generation of technical reports using customized Word templates as *.dotx-files.
SOFTWARE FOR SIMULATION AND OPTIMIZATION OF METAL FORMING PROCESSES AND PROFILE EXTRUSION QForm VX 8.2.3 new facilities and features October 2017 The new version of QForm has some amazing new features,
More informationAccuracy of the Rubicon Toolbox Finite Element Model
Accuracy of the Rubicon Toolbox Finite Element Model Introduction This document deals with the accuracy and recommended use of the Rubicon Toolbox Finite Element module. The document is intended to provide
More informationTutorial 23: Sloshing in a tank modelled using SPH as an example
Tutorial 23: Sloshing in a tank modelled using SPH as an example This tutorial gives a basic introduction to SPH modelling in Abaqus CAE. The tutorial will take you through a basic model of g forces acting
More informationALE and Fluid-Structure Interaction in LS-DYNA March 2004
ALE and Fluid-Structure Interaction in LS-DYNA March 2004 Workshop Models 1. Taylor bar impact 2. One-dimensional advection test 3. Channel 4. Underwater explosion 5. Bar impacting water surface 6. Sloshing
More informationGEO-SLOPE International Ltd, Calgary, Alberta, Canada Sheet Pile Wall
1 Introduction Sheet Pile Wall Deep excavation on a level group usually results in a very low factor of safety, unless it is properly reinforced. The purpose of this example is to illustrate how the stability
More informationTutorial 17. Using the Mixture and Eulerian Multiphase Models
Tutorial 17. Using the Mixture and Eulerian Multiphase Models Introduction: This tutorial examines the flow of water and air in a tee junction. First you will solve the problem using the less computationally-intensive
More informationModule 1.5: Moment Loading of a 2D Cantilever Beam
Module 1.5: Moment Loading of a D Cantilever Beam Table of Contents Page Number Problem Description Theory Geometry 4 Preprocessor 7 Element Type 7 Real Constants and Material Properties 8 Meshing 9 Loads
More informationUse 6DOF solver to calculate motion of the moving body. Create TIFF files for graphic visualization of the solution.
Introduction The purpose of this tutorial is to provide guidelines and recommendations for setting up and solving a moving deforming mesh (MDM) case along with the six degree of freedom (6DOF) solver and
More informationDYNAMIC ANALYSIS OF A GENERATOR ON AN ELASTIC FOUNDATION
DYNAMIC ANALYSIS OF A GENERATOR ON AN ELASTIC FOUNDATION 7 DYNAMIC ANALYSIS OF A GENERATOR ON AN ELASTIC FOUNDATION In this tutorial the influence of a vibrating source on its surrounding soil is studied.
More informationDam Construction by Stages
1 Introduction Dam Construction by Stages This simple example demonstrates the simulation of staged construction of an embankment on soft ground. The primary purposes of this example are to demonstrate
More informationOasys Frew. Copyright Oasys 2013
Oasys Frew Copyright Oasys 2013 All rights reserved. No parts of this work may be reproduced in any form or by any means - graphic, electronic, or mechanical, including photocopying, recording, taping,
More informationGuidelines for proper use of Plate elements
Guidelines for proper use of Plate elements In structural analysis using finite element method, the analysis model is created by dividing the entire structure into finite elements. This procedure is known
More informationSimulations of tunnel excavation in 2D and 3D conditions considering building loads
Simulations of tunnel excavation in D and 3D conditions considering building loads T. Nakai, E. Sung, H.M. Shahin & M. Yamamoto Nagoya Institute of Technology, Nagoya, Japan ABSTRACT: Two-dimensional model
More informationMetafor FE Software. 2. Operator split. 4. Rezoning methods 5. Contact with friction
ALE simulations ua sus using Metafor eao 1. Introduction 2. Operator split 3. Convection schemes 4. Rezoning methods 5. Contact with friction 1 Introduction EULERIAN FORMALISM Undistorted mesh Ideal for
More informationProblem description. The FCBI-C element is used in the fluid part of the model.
Problem description This tutorial illustrates the use of ADINA for analyzing the fluid-structure interaction (FSI) behavior of a flexible splitter behind a 2D cylinder and the surrounding fluid in a channel.
More informationAnalysis Steps 1. Start Abaqus and choose to create a new model database
Source: Online tutorials for ABAQUS Problem Description The two dimensional bridge structure, which consists of steel T sections (b=0.25, h=0.25, I=0.125, t f =t w =0.05), is simply supported at its lower
More information2D & 3D Semi Coupled Analysis Seepage-Stress-Slope
D & D Semi Coupled Analysis Seepage-Stress-Slope MIDASoft Inc. Angel Francisco Martinez Civil Engineer MIDASoft NY office Integrated Solver Optimized for the next generation 6-bit platform Finite Element
More informationME 442. Marc/Mentat-2011 Tutorial-1
ME 442 Overview Marc/Mentat-2011 Tutorial-1 The purpose of this tutorial is to introduce the new user to the MSC/MARC/MENTAT finite element program. It should take about one hour to complete. The MARC/MENTAT
More informationAn Introduction to SolidWorks Flow Simulation 2010
An Introduction to SolidWorks Flow Simulation 2010 John E. Matsson, Ph.D. SDC PUBLICATIONS www.sdcpublications.com Schroff Development Corporation Chapter 2 Flat Plate Boundary Layer Objectives Creating
More informationLab 9: FLUENT: Transient Natural Convection Between Concentric Cylinders
Lab 9: FLUENT: Transient Natural Convection Between Concentric Cylinders Objective: The objective of this laboratory is to introduce how to use FLUENT to solve both transient and natural convection problems.
More informationFinite Element Course ANSYS Mechanical Tutorial Tutorial 3 Cantilever Beam
Problem Specification Finite Element Course ANSYS Mechanical Tutorial Tutorial 3 Cantilever Beam Consider the beam in the figure below. It is clamped on the left side and has a point force of 8kN acting
More informationSolidWorks Flow Simulation 2014
An Introduction to SolidWorks Flow Simulation 2014 John E. Matsson, Ph.D. SDC PUBLICATIONS Better Textbooks. Lower Prices. www.sdcpublications.com Powered by TCPDF (www.tcpdf.org) Visit the following websites
More informationDevelopment of the Compliant Mooring Line Model for FLOW-3D
Flow Science Report 08-15 Development of the Compliant Mooring Line Model for FLOW-3D Gengsheng Wei Flow Science, Inc. October 2015 1. Introduction Mooring systems are common in offshore structures, ship
More informationExamine 2D 8.0 is here!
Examine 2D 8.0 is here! Rocscience is pleased to announce the release of Examine 2D version 8.0, an upgrade to our very popular FREE stress analysis program. Examine 2D uses the boundary element method
More informationTABLE OF CONTENTS 1 FISH BEGINNER S GUIDE 2 FISH REFERENCE. FISH in FLAC Contents Introduction Tutorial
FISH in FLAC Contents - 1 TABLE OF CONTENTS 1 FISH BEGINNER S GUIDE 1.1 Introduction... 1-1 1.2 Tutorial... 1-2 2 FISH REFERENCE 2.1 Introduction and Overview... 2-1 2.2 FISH Language Rules, Variables
More informationTutorial 7 Finite Element Groundwater Seepage. Steady state seepage analysis Groundwater analysis mode Slope stability analysis
Tutorial 7 Finite Element Groundwater Seepage Steady state seepage analysis Groundwater analysis mode Slope stability analysis Introduction Within the Slide program, Slide has the capability to carry out
More informationSimulation of Fuel Sloshing Comparative Study
3. LS-DYNA Anwenderforum, Bamberg 2004 Netfreie Verfahren Simulation of Fuel Sloshing Comparative Study Matej Vesenjak 1, Heiner Müllerschön 2, Alexander Hummel 3, Zoran Ren 1 1 University of Maribor,
More informationExercise 1. 3-Point Bending Using the GUI and the Bottom-up-Method
Exercise 1 3-Point Bending Using the GUI and the Bottom-up-Method Contents Learn how to... 1 Given... 2 Questions... 2 Taking advantage of symmetries... 2 A. Preprocessor (Setting up the Model)... 3 A.1
More informationPHASED EXCAVATION OF A SHIELD TUNNEL
5 PHASED EXCAVATION OF A SHIELD TUNNEL The lining of a shield tunnel is often constructed using prefabricated concrete ring segments, which are bolted together within the tunnel boring machine to form
More informationEssay 5 Tutorial for a Three-Dimensional Heat Conduction Problem Using ANSYS
Essay 5 Tutorial for a Three-Dimensional Heat Conduction Problem Using ANSYS 5.1 Introduction The problem selected to illustrate the use of ANSYS software for a three-dimensional steadystate heat conduction
More informationSlope Stability of Open Pit Mine in 2D & 3D
Slope Stability of Open Pit Mine in D & D MIDASoft Inc. Angel Francisco Martinez Civil Engineer Email : a.martinez@midasit.com Integrated Solver Optimized for the next generation64-bit platform Finite
More informationModule 1.6: Distributed Loading of a 2D Cantilever Beam
Module 1.6: Distributed Loading of a 2D Cantilever Beam Table of Contents Page Number Problem Description 2 Theory 2 Geometry 4 Preprocessor 7 Element Type 7 Real Constants and Material Properties 8 Meshing
More informationGeneral Applications
Chapter General Applications The general analysis modules can be used to calculate section properties, wind pressures on buildings and evaluate drainage systems of building roofs. General Applications
More informationEngineeringmanuals. Part3
Engineeringmanuals Part3 Engineering manuals for GEO5 programs Part 3 Chapter 1-12, refer to Engineering Manual Part 1 Chapter 13-19, refer to Engineering Manual Part 2 Chapter 20. Finite Element Method
More informationFully-Coupled Thermo-Mechanical Analysis
Fully-Coupled Thermo-Mechanical Analysis Type of solver: ABAQUS CAE/Standard Adapted from: ABAQUS Example Problems Manual Extrusion of a Cylindrical Aluminium Bar with Frictional Heat Generation Problem
More informationA Locking-free Smoothed Finite Element Formulation (Modified Selective FS/NS-FEM-T4) with Tetrahedral Mesh Rezoning for Large Deformation Problems
A Locking-free Smoothed Finite Element Formulation (Modified Selective FS/NS-FEM-T4) with Tetrahedral Mesh Rezoning for Large Deformation Problems Yuki ONISHI, Kenji AMAYA Tokyo Institute of Technology
More information