# Shell+Bolt Reaction forces

**URL:** <https://prepomax.discourse.group/t/shell-bolt-reaction-forces/2211>\
**Category:** General Questions\
**Created:** [November 24, 2024, 4:10pm UTC](https://prepomax.discourse.group/t/shell-bolt-reaction-forces/2211 "2024-11-24T16:10:37Z")\
**Posts on this page:** 13\
**Page:** 1

<div class="post-metadata">

**Author:** ![Agazzotti](https://avatars.discourse-cdn.com/v4/letter/a/3da27b/32.png) [@Agazzotti](https://prepomax.discourse.group/u/Agazzotti)\
**Post date:** [November 24, 2024, 4:10pm UTC](https://prepomax.discourse.group/t/shell-bolt-reaction-forces/2211/1 "2024-11-24T16:10:37Z")

</div>

HI  
I tried a test consisting of two plates connected by a bolt to get the forces on the bolt.  
a) In the first case I applied a surface force using a rigid body on the upper plate, but the result was not consistent as the upper plate was unloaded. I don’t understand why. In the attached model I have deactivated this behavior.  
b) To the same model I applied the forces on the upper part of the plate as surface forces (Fx= Fy = 50000 N)  
On the bolt the recation forces are in balance with the external forces but an Fz also appears on the bolt (Fz= -15780.37 N) which does not seem to me to be balanced by the contact forces on the surfaces of the two plates, I don’t understand why.  
Thank you for any clarifications you can provide me  
Agazzotti  
[Ag\_1.pmx](https://prepomax.discourse.group/uploads/short-url/mcVOojQ0jYGZLPvyCTUy4lXLL9b.pmx) (5.8 MB)

---

<div class="post-metadata">

**Author:** ![FEAnalyst](https://yyz2.discourse-cdn.com/free1/user_avatar/prepomax.discourse.group/feanalyst/32/688_2.png) [@FEAnalyst](https://prepomax.discourse.group/u/FEAnalyst)\
**Post date:** [November 24, 2024, 4:36pm UTC](https://prepomax.discourse.group/t/shell-bolt-reaction-forces/2211/2 "2024-11-24T16:36:51Z")

</div>

I can’t check your file now but if you share some screenshots of the model (to show how it’s assembled and loaded), I might be able to give you some tips. Sharing the .inp file would also be helpful.

I’m mostly wondering how you modeled the connection (shell plates and solid bolt?) and how you measure the forces.

---

<div class="post-metadata">

**Author:** ![Agazzotti](https://avatars.discourse-cdn.com/v4/letter/a/3da27b/32.png) [@Agazzotti](https://prepomax.discourse.group/u/Agazzotti)\
**Post date:** [November 24, 2024, 5:39pm UTC](https://prepomax.discourse.group/t/shell-bolt-reaction-forces/2211/3 "2024-11-24T17:39:00Z")

</div>

All elements are modeled with solid elements  
Force:N  
Dimension: mm  
Uploading: Ag\_1\_file.zip…

---

<div class="post-metadata">

**Author:** ![FEAnalyst](https://yyz2.discourse-cdn.com/free1/user_avatar/prepomax.discourse.group/feanalyst/32/688_2.png) [@FEAnalyst](https://prepomax.discourse.group/u/FEAnalyst)\
**Post date:** [November 25, 2024, 1:22am UTC](https://prepomax.discourse.group/t/shell-bolt-reaction-forces/2211/4 "2024-11-25T01:22:54Z")

</div>

It seems that the upload failed, probably due to clicking Reply too soon (before uploading was completed).

---

<div class="post-metadata">

**Author:** ![Agazzotti](https://avatars.discourse-cdn.com/v4/letter/a/3da27b/32.png) [@Agazzotti](https://prepomax.discourse.group/u/Agazzotti)\
**Post date:** [November 25, 2024, 9:10am UTC](https://prepomax.discourse.group/t/shell-bolt-reaction-forces/2211/5 "2024-11-25T09:10:55Z")

</div>

I resend the files  
[Ag\_1\_file.zip](https://prepomax.discourse.group/uploads/short-url/x8FMsH4N89kc7mm5mE4gLOmRZTi.zip) (4.5 MB)

---

<div class="post-metadata">

**Author:** ![FEAnalyst](https://yyz2.discourse-cdn.com/free1/user_avatar/prepomax.discourse.group/feanalyst/32/688_2.png) [@FEAnalyst](https://prepomax.discourse.group/u/FEAnalyst)\
**Post date:** [November 25, 2024, 10:58pm UTC](https://prepomax.discourse.group/t/shell-bolt-reaction-forces/2211/6 "2024-11-25T22:58:09Z")

</div>

Thanks. Is the bolt split into two parts to be able to apply preload to it later ? Because it might be better to try with a single/compound part for the bolt, only having an internal surface in the middle to request a surface print with SF (instead of a node print with RF). You will need keywords for that though.

---

<div class="post-metadata">

**Author:** ![ANYS](https://yyz2.discourse-cdn.com/free1/user_avatar/prepomax.discourse.group/anys/32/7737_2.png) [@ANYS](https://prepomax.discourse.group/u/ANYS)\
**Post date:** [November 25, 2024, 11:49pm UTC](https://prepomax.discourse.group/t/shell-bolt-reaction-forces/2211/7 "2024-11-25T23:49:00Z")

</div>

In addition to FEAnalyst suggestion, sum all contact pair forces component by component. I mean, look at all the contacts acting on a component and add them. They should be in balance with the external load.

NOTE: this problem is like a bottle opener.

EDITED: I have check the force balance and ccx gives the CF for the Slave surface not the Master as I assumed. I guess that’s why contact pressure is also given for the slave. Not an issue but that changes the sign criteria with respect what I said. (I was wrong sorry) .  
Negative for the master and positive for the slave.

---

<div class="post-metadata">

**Author:** ![Agazzotti](https://avatars.discourse-cdn.com/v4/letter/a/3da27b/32.png) [@Agazzotti](https://prepomax.discourse.group/u/Agazzotti)\
**Post date:** [November 26, 2024, 1:42pm UTC](https://prepomax.discourse.group/t/shell-bolt-reaction-forces/2211/8 "2024-11-26T13:42:06Z")

</div>

Thanks for your replies.  
I ran further checks and got the following results.  
If we consider the lower surface of the upper plate as the MASTER element we obtain a total contact force on the surface equal to Fz= 14082.58N which does not balance the force on the bolt Fz= 15857.60N (error = 11.2% not acceptable).  
If we consider the lower surface of the plate as a SLAVE element we obtain a total force Fz= 15857.6N which perfectly balances the Fz of the bolt. I don’t understand why. To avoid incorrect evaluations, it would be better not to provide the results of the MASTER surfaces as they would appear to be incorrect.  
Previously I had also verified that if the surface force is assigned through the rigid body the plate is not stressed and that therefore the rigid body command not only serves to distribute the loads but modifies the elastic behavior of the element. If this were the case, it would be appropriate to inform users of this aspect

---

<div class="post-metadata">

**Author:** ![ANYS](https://yyz2.discourse-cdn.com/free1/user_avatar/prepomax.discourse.group/anys/32/7737_2.png) [@ANYS](https://prepomax.discourse.group/u/ANYS)\
**Post date:** [November 26, 2024, 2:00pm UTC](https://prepomax.discourse.group/t/shell-bolt-reaction-forces/2211/9 "2024-11-26T14:00:07Z")

</div>

> [@Agazzotti](#):
>
> upper plate as the MASTER element we obtain a total contact force on the surface equal to Fz= 14082.58N which does not balance the force on the bolt Fz= 15857.60N (error = 11.2% not acceptable)

Upper plate balance involve one external action and three contacts. TIE contact must be extracted separately with node print. Contact print do not recognize it as there is no contact pair.  
I have recently update my previous post. CF extracted with the Contact print card provides the contact force at the Slave . Master CF= -Slave CF . MAster or Slave definition should not affect the overall force balance (Althought I experienced some issues with two masters defined on the same element in oposite faces).

> [@Agazzotti](#):
>
> the rigid body command not only serves to distribute the loads but modifies the elastic behavior of the element

Yep. It becomes Rigid. That should not affect the equilibrium equations ¿isn’t it?

---

<div class="post-metadata">

**Author:** ![Agazzotti](https://avatars.discourse-cdn.com/v4/letter/a/3da27b/32.png) [@Agazzotti](https://prepomax.discourse.group/u/Agazzotti)\
**Post date:** [November 26, 2024, 2:49pm UTC](https://prepomax.discourse.group/t/shell-bolt-reaction-forces/2211/10 "2024-11-26T14:49:09Z")

</div>

No. The stresses change because the contact surface between the two plates (one rigid and the other not) are different; in the specific case the action Fz increases by approximately 9% considering a plate as a rigid body. For this reason I believe that users must be informed of this aspect in order to evaluate when to insert the rigid body or not.

---

<div class="post-metadata">

**Author:** ![Agazzotti](https://avatars.discourse-cdn.com/v4/letter/a/3da27b/32.png) [@Agazzotti](https://prepomax.discourse.group/u/Agazzotti)\
**Post date:** [November 26, 2024, 2:53pm UTC](https://prepomax.discourse.group/t/shell-bolt-reaction-forces/2211/11 "2024-11-26T14:53:19Z")

</div>

I forgot… as already indicated the forces change considering the SLAVE or MASTER surface with a difference of approximately 11%. I don’t understand why.

---

<div class="post-metadata">

**Author:** ![ANYS](https://yyz2.discourse-cdn.com/free1/user_avatar/prepomax.discourse.group/anys/32/7737_2.png) [@ANYS](https://prepomax.discourse.group/u/ANYS)\
**Post date:** [November 26, 2024, 3:00pm UTC](https://prepomax.discourse.group/t/shell-bolt-reaction-forces/2211/12 "2024-11-26T15:00:53Z")

</div>

This is what I get. It matches. I have centered the model to see if I can understand the moments (They are given about the origin)

 ![imagen](https://global.discourse-cdn.com/free1/uploads/prepomax/original/2X/e/e9490dcc9b2eb19a3199d98bbace99c27c4a6779.png)

EDITED: Bolt Top and Bottom sorted in their right position.  
Section print don’t provide as accurate results. Maybe that’s the origin of your discrepancy. Marked in red colour.

---

<div class="post-metadata">

**Author:** ![FEAnalyst](https://yyz2.discourse-cdn.com/free1/user_avatar/prepomax.discourse.group/feanalyst/32/688_2.png) [@FEAnalyst](https://prepomax.discourse.group/u/FEAnalyst)\
**Post date:** [November 26, 2024, 11:22pm UTC](https://prepomax.discourse.group/t/shell-bolt-reaction-forces/2211/13 "2024-11-26T23:22:46Z")

</div>

> [@Agazzotti](#):
>
> the rigid body command not only serves to distribute the loads but modifies the elastic behavior of the element. If this were the case, it would be appropriate to inform users of this aspect

It’s explained in the CalculiX User’s Manual but IMO Abaqus documentation does it better:

> A rigid body is a collection of nodes, elements, or surfaces whose motion is governed by the motion of a single node, called the rigid body reference node. The relative positions of the nodes and elements that are part of the rigid body remain constant throughout a simulation. Therefore, the constituent elements do not deform but can undergo large rigid body motions.
