Comments (7)
akaszynski
commented on May 25, 2024
Shell281 isn’t supported yet, but I’d be happy to add it. Do you have a small sample result file that’s I can use for testing on my end?
from pymapdl.
kub77
commented on May 25, 2024
Currently the smallest rst file I have is around 800MB, but I can make a simple version and upload that in the next hour!
from pymapdl.
akaszynski
commented on May 25, 2024
Perfect. Upload it here or email me the file at [email protected]
I’ll be including this file in the unit tests, so 1 MB or smaller is better (< 1 MB). Also, please include a text output of the results that you'd like me to achieve like in #17.
from pymapdl.
kub77
commented on May 25, 2024
Here is the test file that I just made, don't know if it is too simple, but hopefully it is enough.
sample.zip
In PyAnsys I type
test_result = pyansys.ResultReader('sample.rst')
estress,elem,enode = test_result.ElementStress(0)
print(estress[23])
print(enode[23])
And get
[[nan nan nan nan nan nan]
[nan nan nan nan nan nan]
[nan nan nan nan nan nan]
[nan nan nan nan nan nan]]
[ 1 82 92 8]
And in Ansys I get
ELEMENT= 24 SHELL281
NODE SX SY SZ SXY SYZ SXZ
1 -50.863 -0.63898E-030 -215.25 -0.18465E-015 0.10251E-013 -47.847
82 13.635 -0.74815E-030 -178.71 -0.21958E-015 0.11999E-013 17.232
92 -7.1801 -0.84355E-030 -213.77 -0.56253E-015 0.13214E-013 2.0152
8 -47.523 -0.96156E-030 -204.12 -0.30574E-014 0.12646E-013 2.4081
1 107.75 0.13549E-029 454.30 0.45391E-015-0.21674E-013 100.34
82 -28.816 0.15515E-029 372.47 0.38337E-015-0.24955E-013 -35.077
92 14.454 0.16547E-029 429.43 0.80716E-015-0.26217E-013 1.2719
8 94.254 0.19148E-029 409.31 0.59899E-014-0.25281E-013 -3.5690
It would also be really useful to be able to read the Nodal Forces and Moment from the Elemental Solution using the:
ElementSolutionData(0,'ENF',sort=True)
From PyAnsys for element 24:
array([ 7.1140683e-01, 2.5775826e-06, 1.8592998e+00, 1.7531972e-03,
-5.4216904e-12, -6.6381943e-04, 7.8414015e-02, 4.7199319e-06,
-1.2074181e+00, -9.0049638e-04, -5.2645028e-12, -3.2152122e-05,
7.3660083e-02, 2.5742002e-05, -1.1951995e+00, -2.7250897e-04,
1.0039868e-12, 1.5112829e-04, -1.9362889e-01, 4.7199323e-06,
1.3849777e+00, 6.4305059e-05, 2.9884493e-12, -2.2116321e-04,
3.0604819e-02, -4.8676171e-05, -1.0389121e-01, 5.7917450e-16,
-2.7263033e-25, 4.0045388e-17, -8.5023224e-02, 2.9796447e-05,
-5.3827515e+00, -2.2202423e-03, 3.8493188e-11, 7.6806801e-04,
-8.5418850e-01, 2.4989351e-06, -3.3126956e-01, -9.2828198e-04,
6.3002242e-11, 8.8052053e-05, 2.3875487e-01, -2.1378659e-05,
4.9762526e+00, 2.5969518e-03, 5.0141464e-11, -1.8303801e-04],
dtype=float32)])
And From Ansys:
ELEM= 24 FX FY FZ
1 0.71141 0.25776E-005 1.8593
82 0.78414E-001 0.47199E-005 -1.2074
92 0.73660E-001 0.25742E-004 -1.1952
8 -0.19363 0.47199E-005 1.3850
83 0.30605E-001-0.48676E-004-0.10389
86 -0.85023E-001 0.29796E-004 -5.3828
93 -0.85419 0.24989E-005-0.33127
9 0.23875 -0.21379E-004 4.9763
ELEM= 24 MX MY MZ
1 0.17532E-002-0.54217E-011-0.66382E-003
82 -0.90050E-003-0.52645E-011-0.32152E-004
92 -0.27251E-003 0.10040E-011 0.15113E-003
8 0.64305E-004 0.29884E-011-0.22116E-003
83 0.57917E-015-0.27263E-024 0.40045E-016
86 -0.22202E-002 0.38493E-010 0.76807E-003
93 -0.92828E-003 0.63002E-010 0.88052E-004
9 0.25970E-002 0.50141E-010-0.18304E-003
from pymapdl.
akaszynski
commented on May 25, 2024
You can already get the nodal forces and moment from the element solution with:
elem, res = test_result.ElementSolutionData(0, 'ENF', sort=True)
print(res[23].reshape(8, -1))It might not have been clear that you had to reshape the array to have it match ANSYS's results. I can have it automatically reshaped in future releases, but I'll have to test it for several element types first.
I've been able to add SHELL281, but I'm getting hung up on KEYOPT(8), which is set to 2 in this solution set. According to the PRNSOL command help, ANSYS prints the solution from the top and bottom surfaces, but according to the result file header, the results are written for each layer.
I'm working out where ANSYS stores the layer information in the result file. Once I have that, it's just a matter of reading in only the top and bottom layers. Can you rerun this sample with keyopt(8) = 0 and then upload the sample file? It would be great to be able to test for both key options.
from pymapdl.
kub77
commented on May 25, 2024
Oh perfect, I didn't realize that!
I accidentally didn't save the original file, but I tried rerunning the simulation from the data just stored in the original sample rst file with KeyOpt 8 as 0 and I think it worked! If it didn't just let me know and I'll make a new file for it.
Thank you so much for the help!
from pymapdl.
akaszynski
commented on May 25, 2024
Looking over old issues and saw this is still active. Seems like it reads correctly given keyopt=8
import pyansys
# open result file and read in element results
filename = 'sample_key8=0.rst'
result = pyansys.read_binary(filename)
elem, res = result.element_solution_data(0, 'ENF', sort=True)
print(res[23].reshape(8, -1))[[ 7.1140683e-01 2.5775826e-06 1.8592998e+00 1.7531972e-03
-5.4216904e-12 -6.6381943e-04]
[ 7.8414015e-02 4.7199319e-06 -1.2074181e+00 -9.0049638e-04
-5.2645028e-12 -3.2152122e-05]
[ 7.3660083e-02 2.5742002e-05 -1.1951995e+00 -2.7250897e-04
1.0039868e-12 1.5112829e-04]
[-1.9362889e-01 4.7199323e-06 1.3849777e+00 6.4305059e-05
2.9884493e-12 -2.2116321e-04]
[ 3.0604819e-02 -4.8676171e-05 -1.0389121e-01 2.0976351e-16
2.5021406e-25 2.5852081e-17]
[-8.5023224e-02 2.9796447e-05 -5.3827515e+00 -2.2202423e-03
3.8493188e-11 7.6806801e-04]
[-8.5418850e-01 2.4989351e-06 -3.3126956e-01 -9.2828198e-04
6.3002242e-11 8.8052053e-05]
[ 2.3875487e-01 -2.1378659e-05 4.9762526e+00 2.5969518e-03
5.0141464e-11 -1.8303801e-04]]
from pymapdl.
Related Issues (20)
- Failed scheduled build HOT 1
- Remark on Documentation Notes for "PostProcessing.nodal_displacement(...)" and "PostProcessing.nodal_component_stress(...)" and maybe others? HOT 1
- Failed scheduled build HOT 3
- Failed scheduled build HOT 1
- ACCOPTION causes pymapdl run to fail HOT 7
- Failed scheduled build HOT 13
- Some examples are not plotting properly
- Better error handling
- Failed scheduled build
- Failed scheduled build
- Python versions outdated
- CLI documentation is not updated
- Array parameters retrieval fail when the size is very big HOT 1
- We need a `mapdl.screenshot` command
- Plotting nodal displacement results in one node being selected after the plot, but only when using CINT for crack calculation.
- TRNOPT command is not updated
- Eliminate redundancies regarding ``PYMAPDL_START_INSTANCE``
- Failed scheduled build HOT 1
- PyMAPDL and HPC clusters
- Failed scheduled build HOT 3
Recommend Projects
-
React
A declarative, efficient, and flexible JavaScript library for building user interfaces.
-
Vue.js
🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.
-
Typescript
TypeScript is a superset of JavaScript that compiles to clean JavaScript output.
-
TensorFlow
An Open Source Machine Learning Framework for Everyone
-
Django
The Web framework for perfectionists with deadlines.
-
Laravel
A PHP framework for web artisans
-
D3
Bring data to life with SVG, Canvas and HTML. 📊📈🎉
-
Recommend Topics
-
javascript
JavaScript (JS) is a lightweight interpreted programming language with first-class functions.
-
web
Some thing interesting about web. New door for the world.
-
server
A server is a program made to process requests and deliver data to clients.
-
Machine learning
Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.
-
Visualization
Some thing interesting about visualization, use data art
-
Game
Some thing interesting about game, make everyone happy.
Recommend Org
-
Facebook
We are working to build community through open source technology. NB: members must have two-factor auth.
-
Microsoft
Open source projects and samples from Microsoft.
-
Google
Google ❤️ Open Source for everyone.
-
Alibaba
Alibaba Open Source for everyone
-
D3
Data-Driven Documents codes.
-
Tencent
China tencent open source team.
from pymapdl.