Substructure analysis of a pick-up truck model

This example illustrates the use of the substructure capability in Abaqus to efficiently simulate the vehicle dynamics of a detailed pick-up truck model traveling over road bumps.

The following topics are discussed:

ProductsAbaqus/Standard

The pick-up truck model geometry described in Inertia relief in a pick-up truck is used in this example. The model is organized as a collection of individual parts that are connected together. Twenty substructures are then created, one for each part that may undergo large motions but for which it is reasonable to assume small-strain elastic deformation (e.g., the chassis). Several parts that may deform nonlinearly (e.g., leaf springs for the rear suspension or the stabilizer bar in the front) are modeled using the usual general nonlinear modeling options. Connection points are created for each part using a coupling constraint. The parts are then attached together using appropriate connector elements. A simplified CALSPAN tire model is used (UEL) to model the radial forces in the tires. The vehicle is loaded statically by gravity, accelerated in a dynamic step on a flat road, and run over bumps. Without stress recovery in the substructures the substructure analysis runs an estimated 120 times faster than an equivalent analysis without substructures.

Geometry and materials

The pick-up truck model (1994 Chevrolet C1500) discussed here is depicted in Figure 1 riding over antisymmetric bumps. The model geometry, element connectivity, and material properties are obtained from the Public Finite Element Model Archive of the National Crash Analysis Center at George Washington University. The materials used are described in Inertia relief in a pick-up truck.

The model is organized as a collection of individual parts connected together. Most parts that undergo only small deformations in addition to a large rigid body motion are defined as substructures. Substructures are created for the following parts: the chassis, each of the four A-arms for the front suspension, each of the four wheels, the rear axle, the driveshaft, the engine/transmission, the cabin, each of the two doors, the hood, the seat, the front bumper, the truck bed, and the fuel tank.

The number of retained nodes for each substructure is determined primarily by its connection points with neighboring parts as illustrated in Figure 2 for the cabin substructure. There are twenty points associated with this substructure that are used to connect the cabin to other parts in the model. There are six retained nodes on the cabin bottom (connections to the chassis); three retained nodes for the hood connections (two hinges and the hood lock); three retained nodes for each of the two door connections (two hinges in the front and the door lock in the back); four retained nodes for the seat connections; and one node retained at the center of mass of the vehicle used for yaw, pitch, and roll measurement purposes.

Several parts deform too much to be considered substructures and are modeled using regular elements. The leaf springs in the back and the stabilizer bar in the front are both modeled with beam elements. The front brake assemblies are modeled as rigid bodies.

The connections between parts are modeled using connector elements. JOIN and REVOLUTE connectors are used to model the hinges between each of the following parts: the A-arms and chassis, the doors and cabin, the hood and cabin, the wheels and knuckles, and the leaf springs and chassis. CARTESIAN and CARDAN connections with appropriately defined connector elasticity, connector damping, and connector friction behaviors are used to define some of the bushing connections (e.g., engine mounts). Two UNIVERSAL connectors are used to model the driveshaft connections to the transmission in the front and to the differential in the back. BEAM connectors are used to model rigid connections between parts. The available connector components of relative motion applied to an AXIAL connector are used to lock (or open) the doors and the hood. The available connector components of relative motion applied to a SLOT connector are used to specify the steering by moving the steering rack. The struts are modeled using an AXIAL connection by specifying approximate nonlinear elasticity and damping. Several of the suspension-related parts are shown in Figure 3.

The radial forces in the tires are modeled approximately using a simplified CALSPAN tire model (Frik, Leister, and Schwartz, 1993) implemented via user subroutine UEL. A radial stiffness of 600 N/mm is considered.

Models

For all analyses the gravity-loaded static equilibrium configuration is found first. Since the given mesh geometry corresponds to the gravity-loaded equilibrium position and data are not available for the pre-stress in the suspension springs and tires, the pre-stress forces have to be computed. To achieve this end, a separate static stress analysis with artificial properties for the suspension springs and artificial boundary conditions is first performed, as follows. The vehicle is supported with boundary conditions in the vertical direction at the four wheel spindles (where the tire UELs will be connected) and fixed at the center of mass to prevent in-plane rigid body motion (degrees of freedom 1, 2, and 6). The stiffnesses of the suspension springs are increased artificially by a thousand times in this independent analysis to minimize deformation. The gravity load is then applied to obtain equilibrium stresses in the suspension components and reaction forces at the wheel spindles.

In the analysis of interest (with realistic properties and boundary conditions), the stresses and reaction forces obtained from the artificial static step are used as initial stresses in the suspension spring components and as pre-stress forces in the tires, respectively. A static gravity loading step is run to obtain an equilibrium configuration. This equilibrium configuration differs only slightly from the given initial geometry (the wheel spindles move laterally about two millimeters). Thus, the initial stress state in the suspension springs and tires accurately represents static equilibrium, and the vehicle is ready for dynamic loading.

The vehicle model is prescribed an initial velocity and then accelerated (0.5 g) to the desired velocity (5 m/sec or 7 m/sec) in a direct-integration implicit dynamic step. Once the “cruise” velocity is achieved, the truck model is run over symmetric or antisymmetric bumps (0.2 m high and 5.0 m long).

Results and discussion

In Figure 4 a snapshot of the truck moving forward with a velocity of 7 m/sec (25.2 km/h) and “jumping” over symmetric bumps is shown. The wheels loose contact with the ground and then land again on the road (not shown).

More results are presented for the case when the truck is riding over antisymmetric bumps (see Figure 1). Stresses are recovered for the lower left A-arm substructure and shown in Figure 5 when the front wheels have traveled 3.2 m over the bumps (for information on how to display results on the substructure, see “A script to visualize the element results for beams and shells within substructure in Abaqus” in the Dassault Systèmes Knowledge Base at http://www.3ds.com/support/knowledge-base). The radial forces on the tires are shown in Figure 6 beginning from the moment when the front tires are about to go over the bumps. A zero radial force indicates that the tire is out of contact. The yaw, pitch, and roll angles recorded using a CARDAN connector attached to the model at its center of mass are shown in Figure 7.

The advantage of using substructures instead of regular deformable elements becomes obvious when the total times needed to complete these types of analyses are compared. A full analysis using regular elements has not been performed since approximately five CPU days would be necessary to complete either of the two analyses discussed above. This total time was estimated by running a few increments, estimating the time needed per iteration, and then multiplying the time per iteration by the total number of iterations needed to complete the analysis. Using these estimates, the substructure analysis is up to 120 times faster than the regular mesh analysis, depending on the amount of recovery performed for each substructure.

The abaqus substructurecombine execution procedure can combine model and results data from two substructure output databases into a single output database. For more information, see Combining output from substructures. To combine the substructure output databases, at least one frame of field output must be requested.

Input files

tr_entire_truck_in_phase.inp

Substructure global analysis of the truck going over in-phase bumps.

tr_entire_truck_anti_phase.inp

Substructure global analysis of the truck going over antiphase bumps.

tr_road_antiphase.inp

Definition of antiphase road profile.

tr_road_inphase.inp

Definition of in-phase road profile.

tr_all_nodes.inp

All node definitions.

tr_parameters.inp

PARAMETER definitions.

tr_materials.inp

Material definitions.

tr_materials_plastic_irl.inp

Material definitions.

tr_initial_stress.inp

Initial stress definition for the leaf springs.

tr_lock_doors_and_hood.inp

CONNECTOR MOTION that keeps the doors and the hood locked.

tr_substruct_recovery.inp

Output definitions for the substructures.

tr_brake_front_left.inp

RIGID BODY definition of the front left brake assembly.

tr_brake_front_right.inp

RIGID BODY definition of the front right brake assembly.

tr_parameters_inphase.inp

PARAMETER definitions for the in-phase analysis.

tr_parameters_antiphase.inp

PARAMETER definitions for the antiphase analysis.

tr_readme.inp

Description of procedure used to run the jobs.

User subroutine

exa_tr_radial_uel.f

UEL used to define the tire model.

Substructure generation files

tr_chassis_gen.inp

Chassis.

tr_retained_chassis.inp

Retained nodes for the chassis.

tr_susp_lower_arm_left_gen.inp

Lower left A-arm.

tr_susp_lower_arm_right_gen.inp

Lower right A-arm.

tr_susp_upper_arm_left_gen.inp

Upper left A-arm.

tr_susp_upper_arm_right_gen.inp

Upper right A-arm.

tr_rear_axle_gen.inp

Rear axle.

tr_retained_rear_axle.inp

Retained nodes for the rear axle.

tr_engine_gen.inp

Engine and transmission.

tr_driveshaft_gen.inp

Driveshaft.

tr_cabin_gen.inp

Cabin and front fenders.

tr_retained_cabin.inp

Retained nodes for the cabin.

tr_hood_gen.inp

Hood.

tr_door_left_gen.inp

Left door.

tr_door_right_gen.inp

Right door.

tr_seat_gen.inp

Seat.

tr_bed_gen.inp

Truck bed.

tr_fuel_tank_gen.inp

Fuel tank.

tr_bumper_gen.inp

Front bumper.

tr_wheel_back_left_gen.inp

Rear left wheel.

tr_wheel_back_right_gen.inp

Rear right wheel.

tr_wheel_front_left_gen.inp

Front left wheel.

tr_wheel_front_right_gen.inp

Front right wheel.

Element definitions

tr_rear_susp_leaf_springs.inp

Rear leaf spring suspension.

tr_stabilizer_elts.inp

Front stabilizer bar.

tr_steering_rods_elts.inp

Steering rods and steering rack.

tr_chassis_elts.inp

Chassis.

tr_susp_lower_arm_left_elts.inp

Lower left A-arm.

tr_susp_lower_arm_right_elts.inp

Lower right A-arm.

tr_susp_upper_arm_left_elts.inp

Upper left A-arm.

tr_susp_upper_arm_right_elts.inp

Upper right A-arm.

tr_rear_axle_elts.inp

Rear axle.

tr_engine_elts.inp

Engine and transmission.

tr_driveshaft_elts.inp

Driveshaft.

tr_cabin_elts.inp

Cabin and front fenders.

tr_hood_elts.inp

Hood.

tr_door_left_elts.inp

Left door.

tr_door_right_elts.inp

Right door.

tr_seat_elts.inp

Seat.

tr_bed_elts.inp

Truck bed.

tr_fuel_tank_elts.inp

Fuel tank.

tr_bumper_elts.inp

Front bumper.

tr_wheel_back_left_elts.inp

Rear left wheel.

tr_wheel_back_right_elts.inp

Rear right wheel.

tr_wheel_front_left_elts.inp

Front left wheel.

tr_wheel_front_right_elts.inp

Front right wheel.

Multi-point constraint definitions

tr_chassis_mpc.inp

Chassis.

tr_engine_mpc.inp

Engine and transmission.

tr_cabin_mpc.inp

Cabin and front fenders.

tr_hood_mpc.inp

Hood.

tr_door_left_mpc.inp

Left door.

tr_door_right_mpc.inp

Right door.

tr_seat_mpc.inp

Seat.

tr_fuel_tank_mpc.inp

Fuel tank.

tr_bumper_mpc.inp

Front bumper.

Coupling definitions

tr_chassis_coup.inp

Chassis.

tr_susp_lower_arm_left_coup.inp

Lower left A-arm.

tr_susp_lower_arm_right_coup.inp

Lower right A-arm.

tr_susp_upper_arm_left_coup.inp

Upper left A-arm.

tr_susp_upper_arm_right_coup.inp

Upper right A-arm.

tr_rear_axle_coup.inp

Rear axle.

tr_engine_coup.inp

Engine and transmission.

tr_driveshaft_coup.inp

Driveshaft.

tr_cabin_coup.inp

Cabin and front fenders.

tr_hood_coup.inp

Hood.

tr_door_left_coup.inp

Left door.

tr_door_right_coup.inp

Right door.

tr_seat_coup.inp

Seat.

tr_bed_coup.inp

Truck bed.

tr_fuel_tank_coup.inp

Fuel tank.

tr_bumper_coup.inp

Front bumper.

tr_wheel_back_left_coup.inp

Rear left wheel.

tr_wheel_back_right_coup.inp

Rear right wheel.

tr_wheel_front_left_coup.inp

Front left wheel.

tr_wheel_front_right_coup.inp

Front right wheel.

Connector definitions

tr_conn_aarms_left.inp

Left A-arms.

tr_conn_aarms_right.inp

Right A-arms.

tr_conn_brake_left.inp

Left front brake assembly.

tr_conn_brake_right.inp

Right front brake assembly.

tr_conn_steering_rods.inp

Steering rods.

tr_conn_stabilizer.inp

Stabilizer bar.

tr_conn_leaf_springs.inp

Leaf springs.

tr_conn_engine.inp

Engine and transmission.

tr_conn_driveshaft.inp

Driveshaft.

tr_conn_cabin_to_chassis.inp

Cabin to chassis.

tr_conn_hood.inp

Hood.

tr_conn_door_left.inp

Left door.

tr_conn_door_right.inp

Right door.

tr_conn_seat.inp

Seat.

tr_conn_bed.inp

Truck bed.

tr_conn_fuel_tank.inp

Fuel tank.

tr_conn_bumper.inp

Front bumper.

tr_conn_wheels_back.inp

Rear wheels.

tr_conn_wheels_front.inp

Front wheels.

References

  1. Frik S.GLeister, and WSchwartz, “Simulation of the IAVSD Road Vehicle Benchmark Bombardier Iltis with FASIM, MEDYNA, NEWEUL, and SIMPACK,” in Multibody Computer Codes in Vehicle System Dynamics, Ed. W. Kortum and R. S. Sharp, February 1993.

Figures

Figure 1. Substructure truck model going over antisymmetric (left-right) bumps.

Figure 2. Substructure mesh for the cabin.

Figure 3. Chassis, suspension, and power-train-related parts.

Figure 4. Substructure truck model going over symmetric (left-right) bumps.

Figure 5. Recovered stresses in lower left A-arm.

Figure 6. Tire radial forces for the truck going over antisymmetric bumps.

Figure 7. Yaw, pitch, and roll for the truck going over antisymmetric bumps.