Dear @Jason.Jonkman ,
Thank you for your prompt responses.I have granted you access and updated the permissions to “Anyone with the link”: https://docs.google.com/spreadsheets/d/1Af0DmB8C5zpIFAcvqJdFG5F0RO94bSSS/edit?usp=drive_link&ouid=100534356249426457824&rtpof=true&sd=true. To further avoid access issues, I have also uploaded a compressed file My manual calculations, which contains the same Excel file with my manual calculations: https://drive.google.com/file/d/1zHvuNFr2rDinfstqcolTQzdHd1qSAS6R/view?usp=drive_link.
Best regards,
Dear @Yangyang.Li,
I spent a little time checking your math in your spreadsheet, and didn’t spot any obvious errors, but I don’t have the time to review everything and there are several complicated equations that I don’t see calculations for (only numerical values entered in), and didn’t try to calculate those myself. I’m sure there is a detail missing or small error somewhere. I’m not aware of a specific bug in the force outputs calculations of SubDyn. If you want to dig further, I would run your model through the debugger and check each term in the equations against your own to spot the difference.
Best regards,
Dear @Jason.Jonkman ,
Thank you for your prompt responses. I would like to further understand the discrepancies between my manual calculations and the SubDyn results. If the issues remain unresolved, they may continue to cause problems in the future. I hope you could take some more time to point out where my mistakes lie, and I am more than willing to provide any additional information or materials that may help.
No calculation formulas are provided (only numerical values entered in), generally sourced from the SD.sum file of SubDyn. The “several complicated equations” you referred to, I suspect, may correspond to the Shear, Jx, and Jy parameters in the “Member Parameters” section of the spreadsheet. These parameters are directly taken from the “Elements” section of the sum file, where Jx = Ixx and Jy = Iy.
Best regards,
Dear @Yangyang.Li,
I don’t really have the time to look into this further, and given that I’m not aware of an issue, it is not a high priority for me.
One thing I recalled is that because you are outputting loads in the interior of a member, SubDyn will average the load between the two elements touching that node, which I don’t see that you are computing. Again, I would suggest running it through the debugger to check each term. Alternatively, pick a simpler case (such as a single beam with a known applied load) that would be easier to identify if there is an actual issue with the calculation, which I don’t believe there is.
Best regards,