Continue to Site

Eng-Tips is the largest engineering community on the Internet

Intelligent Work Forums for Engineering Professionals

  • Congratulations waross on being selected by the Eng-Tips community for having the most helpful posts in the forums last week. Way to Go!

ABAQUS Contact problem using SPH Method

Status
Not open for further replies.

samnaib

Materials
Jul 16, 2014
2
hello

I have been facing a Problem in obtaining a water Impact on a material box using sph method. i do not receive any error during Simulation but there is a warning message about
"Some nodes involved in general contact have penetrated their
tracked faces by more than 50.000 percent of the typical element
dimension in the general contact domain, 4.5263. Please check the
node set "InfoNodeDeepPenetFirst" in Abaqus/Viewer for more
information. Please make sure that these nodes and their
contacting faces have sufficient mass to allow the code to
calculate realistic contact penalties for enforcing the contact
constraint."

Simulation does not give anz splash instead the water ball passes through the material.

How do i solve it?

Below I give the Input file
******
******
*Heading
** Job name: watersplash_1 Model name: Model-1
** Generated by: Abaqus/CAE 6.12-1
*Preprint, echo=NO, model=NO, history=NO, contact=NO
**
** PARTS
**
*Part, name=Box
*Node
1, -10., -2.5, 20.
.....
.....
*Nset, nset=_PickedSet5, internal, generate
1, 484, 1
*Elset, elset=_PickedSet5, internal, generate
1, 300, 1
*Elset, elset=_boxsurf_S1, internal, generate
271, 300, 1
*Elset, elset=_boxsurf_S6, internal, generate
1, 298, 3
*Elset, elset=_boxsurf_S2, internal, generate
1, 30, 1
*Elset, elset=_boxsurf_S4, internal, generate
3, 300, 3
*Elset, elset=_boxsurf_S3, internal
1, 2, 3, 31, 32, 33, 61, 62, 63, 91, 92, 93, 121, 122, 123, 151
152, 153, 181, 182, 183, 211, 212, 213, 241, 242, 243, 271, 272, 273
*Elset, elset=_boxsurf_S5, internal
28, 29, 30, 58, 59, 60, 88, 89, 90, 118, 119, 120, 148, 149, 150, 178
179, 180, 208, 209, 210, 238, 239, 240, 268, 269, 270, 298, 299, 300
*Surface, type=ELEMENT, name=boxsurf
_boxsurf_S1, S1
_boxsurf_S6, S6
_boxsurf_S2, S2
_boxsurf_S4, S4
_boxsurf_S3, S3
_boxsurf_S5, S5
** Section: boxsection
*Solid Section, elset=_PickedSet5, material=Steel
1.,
*End Part
**
*Part, name=Water
*Node
1, 0., 20., 0.
2, 0., 15., 5.
........
........
*Element, type=PC3D, elset=elwater
11671,1
11672,2
11673,3
.......
*Nset, nset=Nset, generate
1, 16947, 1
*Solid section, elset=elwater, material=water, controls=sph_fac
0.00005
*End Part
**
**
** ASSEMBLY
**
*Assembly, name=Assembly
**
*Instance, name=Box-1, part=Box
0., 0., -10.
*End Instance
**
*Instance, name=Water-1, part=Water
*End Instance
**
*Nset, nset=_PickedSet7, internal, instance=Box-1, generate
1, 481, 4
*Elset, elset=_PickedSet7, internal, instance=Box-1, generate
1, 298, 3
*Nset, nset=watersurf, instance=Water-1, generate
1, 16947, 1
*Elset, elset=_boxsurf_S1, internal, instance=Box-1, generate
271, 300, 1
*Elset, elset=_boxsurf_S6, internal, instance=Box-1, generate
1, 298, 3
*Elset, elset=_boxsurf_S2, internal, instance=Box-1, generate
1, 30, 1
*Elset, elset=_boxsurf_S4, internal, instance=Box-1, generate
3, 300, 3
*Elset, elset=_boxsurf_S3, internal, instance=Box-1
1, 2, 3, 31, 32, 33, 61, 62, 63, 91, 92, 93, 121, 122, 123, 151
152, 153, 181, 182, 183, 211, 212, 213, 241, 242, 243, 271, 272, 273
*Elset, elset=_boxsurf_S5, internal, instance=Box-1
28, 29, 30, 58, 59, 60, 88, 89, 90, 118, 119, 120, 148, 149, 150, 178
179, 180, 208, 209, 210, 238, 239, 240, 268, 269, 270, 298, 299, 300
*Surface, type=ELEMENT, name=boxsurf
_boxsurf_S1, S1
_boxsurf_S6, S6
_boxsurf_S2, S2
_boxsurf_S4, S4
_boxsurf_S3, S3
_boxsurf_S5, S5
*Surface, type=NODE, name=watersurf
Water-1.Nset
*End Assembly
**
** MATERIALS
**
*Material, name=Steel
*Density
8.05,
*Elastic
210., 0.3
*Material, name=Water
*Density
1.,
*Eos, type=USUP
0.111803,0.,0.
** ----------------------------------------------------------------
**
** STEP: Step-1
**
*Step, name=Step-1
*Dynamic, Explicit
, 1.
*Bulk Viscosity
0.06, 1.2
**
** BOUNDARY CONDITIONS
**
** Name: fixedboundary Type: Symmetry/Antisymmetry/Encastre
*Boundary
_PickedSet7, ENCASTRE
**
** LOADS
**
** Name: cload Type: Concentrated force
*Cload
watersurf, 2, -5e-05
**
** CONTACT
*contact
*contact inclusions
boxsurf, watersurf
**
** OUTPUT REQUESTS
**
*Restart, write, number interval=1, time marks=NO
**
** FIELD OUTPUT: F-Output-1
**
*Output, field, variable=PRESELECT
**
** HISTORY OUTPUT: H-Output-1
**
*Output, history, variable=PRESELECT
*End Step

 
Status
Not open for further replies.

Part and Inventory Search

Sponsor