Products: ABAQUS/Standard ABAQUS/Explicit ABAQUS/CAE
Submodeling is the technique used in ABAQUS for analyzing a local part of a model with a refined mesh, based on interpolation of the solution from an initial global model (usually with a coarser mesh) onto the nodes on the appropriate parts of the boundary of the submodel. Shell-to-solid submodeling models a region with solid elements, when the global model is made up of shell elements.
Shell-to-solid coupling is a feature in ABAQUS by which three-dimensional shell meshes can be coupled automatically to three-dimensional solid meshes. Unlike shell-to-solid submodeling, which first performs a global analysis on a shell model followed by a submodel analysis on the continuum model, the shell and continuum meshes in a shell-to-solid coupling model are used in the same analysis.
Both shell-to-solid submodeling and shell-to-solid coupling provide cost-effective approaches to model enhancement. The purpose of this example is to demonstrate both capabilities in ABAQUS. The analysis is tested as a static process in ABAQUS/Standard and as a dynamic process in both ABAQUS/Standard and ABAQUS/Explicit. To demonstrate the shell-to-solid submodeling capability, the problem is solved quasi-statically in ABAQUS/Explicit.
In addition, an ABAQUS Scripting Interface script is included that creates a shell global model using ABAQUS/CAE. The script then uses data from the output database created by the analysis of the global model to drive a solid submodel. The script ends by displaying an overlay plot of the global model and the submodel in the Visualization module.
In this problem the joint between a pipe and a plate is analyzed. A pipe of radius 10 mm and thickness 0.75 mm is attached to a plate that is 10 mm long, 5 mm wide, and 1 mm thick. The pipe-plate intersection has a fillet radius of 1 mm. Taking advantage of the symmetry of the problem, only half the assembly is modeled. Both the pipe and the plate are assumed to be made of aluminum with 
= 69 × 103 MPa, 
= 0.3, and 
= 2740 kg/m3.
The global model for the submodeling analysis is meshed with S4R elements as shown in Figure 1.1.9–1. The fillet radius is not taken into consideration in the shell model. The static submodel is meshed using three-dimensional C3D20R continuum elements (see Figure 1.1.9–2). A coarser mesh using C3D8R elements is chosen for the dynamic tests. The shell-to-solid coupling model is meshed with S4R shell elements and C3D20R continuum elements as shown in Figure 1.1.9–3. The continuum meshes used in the static submodeling and shell-to-solid coupling analyses are identical. The continuum meshes extend 10 mm along the pipe length, have a radius of 25 mm in the plane of the plate, and use four layers through the thickness. The continuum meshes accurately model the fillet radius at the joint. Hence, it is possible to calculate the stress concentration in the fillet. The problem could be expanded by adding a ring of welded material to simulate a welded joint (for this case the submodel would have to be meshed with new element layers representing the welded material at the joint). The example could also be expanded by including plastic material behavior in the submodel while using an elastic global model solution.
A reference static solution consisting entirely of C3D20R continuum elements is also included (see Figure 1.1.9–4). The mesh of the reference solution in the vicinity of the joint is very similar to that used in the submodeling and shell-to-solid coupling analyses.
The geometry and material properties of the ABAQUS/Explicit shell-to-solid coupling model are identical to the ABAQUS/Standard models. The ABAQUS/Explicit shell model is meshed with S4R elements, and the continuum model is meshed with C3D10M elements.
The pipe is subjected to a 10 N load acting in the x-direction applied at the free end, representing a shear load on the pipe. An edge-based surface is defined at the free edge of the pipe. This surface is coupled to a reference node that is defined at the center of the pipe using the *COUPLING option in conjunction with the *DISTRIBUTING option. The 10 N load is applied to the reference point. For the dynamic cases the load is applied gradually over the entire step time by using the SMOOTH STEP parameter on the *AMPLITUDE option.
The plate is clamped along all edges. In the solid submodel, kinematic conditions are interpolated from the global model at two surfaces of the submodel: one lying within the pipe and the other within the plate. The default center zone size, equal to 10% of the maximum shell thickness, is used. Thus, only one layer of driven nodes lies within the center zone, and only these nodes have all three displacement components driven by the global solution. For the remaining driven nodes only the displacement components parallel to the global model midsurface are driven from the global model. Thus, a single row of nodes transmits the transverse shear forces from the shell solution to the solid model.
The loading and boundary conditions are such that the pipe is subjected to bending. The end of the pipe that is attached to the plate leads to deformation of the plate itself (see Figure 1.1.9–5 and Figure 1.1.9–6). From a design viewpoint the area of interest is the pipe-plate joint where the pipe is bending the plate. Hence, this area is modeled with continuum elements to gain a better understanding of the deformation and stress state.
Figure 1.1.9–7 shows the contours of the out-of-plane displacement component in the plate for both the static submodel and the shell-to-solid coupling analyses. The submodel is in good agreement with the displacement of the global shell model around the joint. The out-of-plane displacement for the shell-to-solid coupling is slightly less than that for the submodel analysis but is in good agreement with the reference solution shown in Figure 1.1.9–8.
The stress concentration in the fillet radius is obtained for the solid models. The maximum Mises stresses at the integration points and nodes for the reference solution, submodel, and shell-to-solid coupling analyses are shown in Table 1.1.9–1. As illustrated in Table 1.1.9–1 and Figure 1.1.9–9, the Mises stress computed in the shell-to-solid coupling analysis agrees very well with the reference solution. The continuity of displacements and the minimal distortion of the stress field at the shell-to-solid interface indicates that the shell-to-solid coupling has been modeled accurately. The difference in the maximum Mises stress between the submodel analysis and the shell-to-solid coupling solution (see Figure 1.1.9–10) can be partially attributed to the fact that the global shell model is more flexible than the shell-to-solid coupling model. The x-displacement at the distributing coupling reference node for the global shell model due to the 10 N load is .605 mm compared to .513 mm for the shell-to-solid coupling analysis and .512 mm for the reference solution. Thus, the static submodel mesh is subjected to slightly higher deformation. If the global shell analysis is run with an x-displacement boundary condition of .513 mm on the reference node instead of a concentrated load of 10 N, the subsequent maximum nodal Mises stress in the submodel analysis drops to 80.2 MPa, which is in better agreement with the reference solution.
The relatively large difference between the maximum Mises stresses at the integration points and the nodes in the region of the fillet (as illustrated in Table 1.1.9–1) indicates that the mesh in the fillet region is probably too coarse and should be refined. No such refinement was performed in this example.
Overall, the ABAQUS/Explicit shell-to-solid coupling analysis is in good agreement with the ABAQUS/Standard shell-to-solid coupling results. The continuity of displacements and the minimal distortion of the stress field at the shell-to-solid interface indicate that the shell-to-solid coupling has been modeled accurately. The out-of-plane displacements of the plate predicted by ABAQUS/Explicit are very close to the ABAQUS/Standard values. The maximum nodal Mises stress in the fillet region is 56 MPa, and the x-displacement at the *COUPLING reference node for the global shell model due to the 10 N load is .423 mm. The ABAQUS/Explicit analysis is solved quasi-statically by assigning a nominal density of 500 kg/m3 to the pipe-plate material and ramping up the load over 12,000 increments. Closer approximation to the static limit, achieved by reducing the density of the pipe-plate material to 50 kg/m3, results in a maximum Mises stress of 87 MPa in the fillet region, which is very close to the ABAQUS/Standard result.
The results for the submodel dynamic cases agree well with the global results. Both ABAQUS/Standard and ABAQUS/Explicit submodels read the results of the same global ABAQUS/Explicit analysis. Good agreement is also found between the ABAQUS/Explicit and ABAQUS/Standard submodel analyses.
S4R global model.
Node definitions for the S4R global model.
Element definitions for the S4R global model.
C3D20R submodel that uses the S4R global model.
Node definitions for the C3D20R submodel that uses the S4R global model.
Element definitions for the C3D20R submodel that uses the S4R global model.
S4 global model.
Node definitions for the S4 global model.
Element definitions for the S4 global model.
C3D20R submodel that uses the S4 global model.
Node definitions for the C3D20R submodel that uses the S4 global model.
Element definitions for the C3D20R submodel that uses the S4 global model.
Python script that creates an S4 global model and a C3D20R submodel using ABAQUS/CAE.
Shell-to-solid coupling model with C3D20R and S4R elements.
Node definitions for the shell-to-solid coupling model with C3D20R and S4R elements.
Node definitions for the shell-to-solid coupling model with C3D20R and S4R elements.
Node definitions for the shell-to-solid coupling model with C3D20R and S4R elements.
Element definitions for the shell-to-solid coupling model with C3D20R and S4R elements.
Element definitions for the shell-to-solid coupling model with C3D20R and S4R elements.
Element definitions for the shell-to-solid coupling model with C3D20R and S4R elements.
Shell-to-solid coupling model with C3D10 and S4R elements.
Node definitions for the shell-to-solid coupling model with C3D10 and S4R elements.
Node definitions for the shell-to-solid coupling model with C3D10 and S4R elements.
Node definitions for the shell-to-solid coupling model with C3D10 and S4R elements.
Element definitions for the shell-to-solid coupling model with C3D10 and S4R elements.
Element definitions for the shell-to-solid coupling model with C3D10 and S4R elements.
Element definitions for the shell-to-solid coupling model with C3D10 and S4R elements.
ABAQUS/Explicit shell-to-solid coupling model with C3D10M and S4R elements.
Node definitions for the ABAQUS/Explicit shell-to-solid coupling model with C3D10M and S4R elements.
Node definitions for the ABAQUS/Explicit shell-to-solid coupling model with C3D10M and S4R elements.
Node definitions for the ABAQUS/Explicit shell-to-solid coupling model with C3D10M and S4R elements.
Element definitions for the ABAQUS/Explicit shell-to-solid coupling model with C3D10M and S4R elements.
Element definitions for the ABAQUS/Explicit shell-to-solid coupling model with C3D10M and S4R elements.
Element definitions for the ABAQUS/Explicit shell-to-solid coupling model with C3D10M and S4R elements.
Reference model with C3D20R elements.
Node definitions for the reference model with C3D20R elements.
Element definitions for the reference model with C3D20R elements.
ABAQUS/Explicit S4R global model.
Node definitions for the ABAQUS/Explicit S4R global model.
Element definitions for the ABAQUS/Explicit S4R global model.
ABAQUS/Explicit C3D8R submodel.
Node definitions for the ABAQUS/Explicit C3D8R submodel.
Element definitions for the ABAQUS/Explicit C3D8R submodel.
ABAQUS/Standard C3D8R submodel.
Node definitions for the ABAQUS/Standard C3D8R submodel.
Element definitions for the ABAQUS/Standard C3D8R submodel.
Figure 1.1.9–5 Solid submodel overlaid on the shell model in the deformed state, using a magnification factor of 20.
Figure 1.1.9–6 Shell-to-solid coupling model in the deformed state, using a magnification factor of 20.
Figure 1.1.9–7 Comparison of out-of-plane displacement in the continuum mesh plate for the submodel (top) and the shell-to-solid coupling analysis (bottom).
Figure 1.1.9–8 Comparison of out-of-plane displacement in the plate for the reference solution (top) and the shell-to-solid coupling analysis (bottom).
