SteveDoyle2
commented
1 year ago I’d need to see a small example with the bdf/f06/op2. A small example also let’s me add it as a test. There are a few different types of CTRIA6 output formats depending on solution type/solver.
The time is for modes/time/frequency. You should not be running into that for a static case. The difference in the blocks is that the time=0 blocks sets the node/element ids and some other stuff that is common to all time steps.
On Fri, Aug 18, 2023 at 4:19 AM gornezo @.***> wrote:
Dear community,
I have encountered an isssue when trying to load the results of an OP2 model containing both CTRIA3 and CTRIA6 results. Stress output is set to CORNER, thus this forces the op2 reader to load the results of the CTRIA3 in the same directory 'ctria6_stress' that would be normally used by the CTRIA6 results. I believe this is what causes the issue.
Here is an extract of my code:
# read in model model = OP2(debug=False) model.read_op2(file_) #load results dictionaries stress_ctria3 = model.ctria3_stress stress_ctria6 = model.ctria6_stressThe obtained error message is: File "C:\Program Files\Python39\lib\site-packages\pyNastran\op2\tables\oes_stressStrain\oes.py", line 5347, in _oes_cquad4_144 obj.data[obj.itime, istart:iend, :] = results ValueError: could not broadcast input array from shape (712,8) into shape (3800,8)
In ,my test model I have used 89 CTRIA3 Elements -> 712 entries (4 nlayersnlements). Then the CTRIA6 data containing 475 CTRIA6 elements -> 3800 entries (4nlayersnlements) is attempted to be written. Strangely (or maybe because im not understanding something), the itime value when trying to write the CTRIA6 elements is set to 1 instead of 0, which causes the sizing of the data container not to be updated. If I update the code in the oes.py file line 5233 : from:
if obj.itime == 0:
to:
if obj.itime in [0, 1]:
then the program will run, but the content of the ctria6_stress dictionary will only provide the values of the CTRIA6 elements and not the CTRIA3.
My calculation is a linear calculation, thus the time increment of 1 does not seem logical to me. Is this behaviour intended? Then, it remains the question, how can I get results for the case where I have this mixed elements? it seems to work fine for CQUAD4 and CQUAD8 mixture.
I appreciate your help and of course thank you for all the hard work set into this amazing library!
Kind regards, Chris
— Reply to this email directly, view it on GitHub https://github.com/SteveDoyle2/pyNastran/issues/740, or unsubscribe https://github.com/notifications/unsubscribe-auth/AAICUWPS4PCWKSM6VYTXMTLXV5FU7ANCNFSM6AAAAAA3VMH6GY . You are receiving this because you are subscribed to this thread.Message ID: @.***>
gornezo
Dear community,
I have encountered an isssue when trying to load the results of an OP2 model containing both CTRIA3 and CTRIA6 results. Stress output is set to CORNER, thus this forces the op2 reader to load the results of the CTRIA3 in the same directory 'ctria6_stress' that would be normally used by the CTRIA6 results. I believe this is what causes the issue.
Here is an extract of my code:
The obtained error message is: File "C:\Program Files\Python39\lib\site-packages\pyNastran\op2\tables\oes_stressStrain\oes.py", line 5347, in _oes_cquad4_144 obj.data[obj.itime, istart:iend, :] = results ValueError: could not broadcast input array from shape (712,8) into shape (3800,8)
In ,my test model I have used 89 CTRIA3 Elements -> 712 entries (4nlayersnlements). Then the CTRIA6 data containing 475 CTRIA6 elements -> 3800 entries (4nlayersnlements) is attempted to be written. Strangely (or maybe because im not understanding something), the itime value when trying to write the CTRIA6 elements is set to 1 instead of 0, which causes the sizing of the data container not to be updated. If I update the code in the oes.py file line 5233 : from:
to:
then the program will run, but the content of the ctria6_stress dictionary will only provide the values of the CTRIA6 elements and not the CTRIA3.
My calculation is a linear calculation, thus the time increment of 1 does not seem logical to me. Is this behaviour intended? Then, it remains the question, how can I get results for the case where I have this mixed elements? it seems to work fine for CQUAD4 and CQUAD8 mixture.
I appreciate your help and of course thank you for all the hard work set into this amazing library!
Kind regards, Chris