MES,CAD,CAE,FEA
News Search:
StartGaleriaForumDownloadArtykułyNewsyLinkiZasoby wiedzyKontaktZasady Stosujemy pliki cookies więcej...02/28/2020 12:05:08
Reklama
Nawigacja
Galeria
Forum
Download
Artykuły
FAQ
Linki
Zasoby wiedzy
Kontakt
Szukaj
Losowa Fotka
Naprężenia w modelu 2 top
Naprężenia w modelu 2 top
Naprężenia w płytach
Użytkownicy Online
CAD, MES -metoda elementów skończonych,obliczenia inżynierskie i metody numeryczne
WITAMY:
dabeeek
jako nowego użytkownika.

Zarejestrowanch Uzytkowników: 1373

Super Administratorzy: 1
Administratorzy: 1
Użytkownicy: 1371

Użytkownicy Online:

BE-FEA01:04:08
dabeeek 1 dzień
Vinne 3 dni
mikelfr 2 tygodni
pvalue 3 tygodni

Gości Online: 14

Twoje IP to: 3.229.118.253

Kategorie Forum 53
Wątki na Forum 1248
Posty na Forum 5601
Komentarzy 83
Newsy 141
Artykuły 10
Ściągniętych plików 62
Kategorie Downloads 4
Ciekawe Strony 13
Postów w Shoutbox 404
Przykłady
Model tarczowy - wyniki

Model tarczowy - wyniki

Połączenie sworzniowe
Naprężenia redukowane model porównawczy 1

Naprężenia redukowane model porównawczy 1

Problem Oli
Skan tomograficzny zęba

Skan tomograficzny zęba

Ząb
Wartość naprężeń H-M-H At Z1

Wartość naprężeń H-M-H At Z1

Problem Oli
Paszcza

Paszcza

FEM-art, ukryte piękno obliczeń
Partnerzy

Strony wspierającej nas firmy:
Leasing,faktoring, kredyt gdy szukasz gotówki.
Opony to tanie opony letnie oraz katalog
zawierający firmy transportowe
Strony naszych przyjaciół:

Darmowe arkusz Excela
Strony sponsorów:
Zawsze się znajdzie forum transportowe na giełdzie ładunków/ na stronie giełdy ładunków
ciekawe i miłe sercu rzeczy to: fajne fotki dzieci na stronie mojebaby.pl lub też warte obejrzenia zwierzęta domowe ze strony mojzwierz.pl
-->../lll
Ostatnie Komentarze
marek_ac dnia 11.01.2017 08:14

Demo download (nowy link): (poprawiono patrz: Linki-> Przyjaciele -> Vanadis)

Zobacz Komentarze Artykułów
marek_ac dnia 21.04.2015 05:53

... mała poprawka: double wprowadzono w wersji 1.3
Zobacz Komentarze Newsów
marek_ac dnia 20.04.2015 18:59

Jeżeli do oprogramowania inżynierskiego można podpiąć własne procedury to NVIDIA udostępnia darmowy driver NVCC za pomocą którego można kompilować kod napisany w C. Można to robić z Visual Studio (od wersji 10) pod Windows lub używając gcc pod Linuxem.
Wówczas oprogramowanie optymalizuje się pod posiadaną kartę graficzną. Warto mieć kartę o CUDA capability większym niż 2.0 ze względu na obsługę typu double.
Zobacz Komentarze Newsów
marek_ac dnia 14.04.2015 00:51

Witam,

Na dysku:

TUTAJ

zamieściłem demo (wersję windowsową) modelu rozprzestrzeniania się zanieczyszczeń.

W pliku vanadis.txt zawarty jest krótki opis modelu.

Marek
Zobacz Komentarze Artykułów
steleri dnia 17.03.2015 21:54

Czy ktoś z forum korzysta może z Elmera? Bardzo ciekawi mnie jak to działa. Bo z tego tutka wynika, że jest OK, a tymczasem ja próbowałem... bezskutecznie
Zobacz Komentarze Newsów

Ostatnie aktywne Tematy rss1 rss2
Forum Wątek Odpowiedzi Ostatni Post
Sprawy studenckie Symulacja procesu ciągnienia drutu (MES) 1 BE-FEA 02/24/2020 23:18
Sprawy studenckie Sprzężenie pola mechanicznego i akustycznego - Code Aster 3 BE-FEA 02/24/2020 19:04
Teoria MES inżynieria odwrotna metody różnic skończonych 1 BE-FEA 02/05/2020 09:07
Sprawy ANSYSA Implicit a explicit 2 BE-FEA 02/03/2020 23:32
Sprawy studenckie Pomoc w zadaniu z MES 1 BE-FEA 01/29/2020 18:59
Sprawy studenckie Praca Magisterska Ansys 2 piotr 11/14/2019 17:08
Jestem kompletnie zielony Współczynnik konwekcji 1 BE-FEA 09/30/2019 19:14
Wszystko o CAE, co nie pasuje do powyższych forum Bezpłatne webinarium CAE - jak działa program do symulacj... 0 MatLit 08/28/2019 11:20
Teoria MES Łączenie elementów skończonych różnego typu 1 BE-FEA 06/10/2019 21:59
Teoria MES Macierz pochodnych funkcji kształtu [B] 14 BE-FEA 06/05/2019 14:44
Jestem kompletnie zielony Jak połączyć wolne powierzchnie w Patranie 21 BE-FEA 05/31/2019 17:52
Teoria MES Rodzaje solverów MES 14 BE-FEA 05/24/2019 19:12
Jestem kompletnie zielony Problem z importem plików do Patrana 6 Solid-ny amator 04/23/2019 23:35
Jestem kompletnie zielony [ANSYS] Wytłaczanie uszczelek 1 Rogal 04/08/2019 11:26
Sprawy studenckie Odpłatna pomoc 2 piotr 03/28/2019 21:56
Zobacz Temat
CAD, MES -metoda elementów skończonych,obliczenia inżynierskie i metody numeryczne | Oprogramowanie CAE | Sprawy Abaqusa
Autor Odczytanie napręzen
Goral
Użytkownik

Postów: 21
Data rejestracji: 19.03.13
Dodane dnia 07/02/2013 22:20
Witam,
Potrzebuje odczytac wartosci naprezen z setu w postaci linii. Prosilbym o rozwiazanie, najlepiej w postaci kodu Pythona.
Pozdrawiam
Autor RE: Odczytanie napręzen
Tomekk
Użytkownik

Postów: 242
Miejscowość: Norge
Data rejestracji: 01.10.09
Dodane dnia 07/02/2013 23:16
Witam,
Zajrzyj do "Abaqus Scripting User's Manual" (jakbyś zajrzał, to byś głupiego pytania nie zadawał) tam znajdziesz wszystkie informacje, których potrzebujesz i pewnie całą masę innych rzeczy o istnieniu których nie wiedziałeś. A jak już to przeglądniesz i pojawią Ci się konkretne pytania to najpierw zajrzyj tu: "Abaqus Scripting Reference Manual", potem prześpij się z problemem i jak dalej będzie nie jasne i TWÓJ skrypt nie będzie działać, to daj znać, wtedy pomyślimy.


pozdr
Tomek

-----------------------------
człowiekiem jestem ...
Autor RE: Odczytanie napręzen
Goral
Użytkownik

Postów: 21
Data rejestracji: 19.03.13
Dodane dnia 07/02/2013 23:50
Zajrzałem. Potrafie odczytać wartości ugięć z dowolnego seta np. w postaci linii:
import sys

from odbAccess import *
from abaqusConstants import *
from odbMaterial import *
from odbSection import *

odbpath=sys.argv[1]
outpath=sys.argv[2]

out = open(outpath,'w+')
odb = openOdb(path=odbpath)

nSet = odb.rootAssembly.instances['PART-1-1'].nodeSets['MONIT']
frame = odb.steps['Step-1'].frames[-1]
dispField = frame.fieldOutputs['U']
dispSubField = dispField.getSubset(region=nSet)
fieldValues = dispSubField.values

for v1 in fieldValues:
print >> out, '%11.10e' % (v1.data[0])

Niestety zamiana U na S nie daje skutku, bo w wezlach nie da sie odczytac naprezen i dlatego wartosci, ktore otrzymuje sa rowne 0.
Znalazlem cos innego, NFORC. Niestety podmiana rozniez nie dziala... choc w modelu zdefiniowalem zapisywanie NFORCE.
Prawdopodobnie nie dosc precyzyjnie odczytuje NFORC i otrzymuje "zmienna" zamiast skalaru jak dla U, tylko wektor.
Autor RE: Odczytanie napręzen
Tomekk
Użytkownik

Postów: 242
Miejscowość: Norge
Data rejestracji: 01.10.09
Dodane dnia 07/04/2013 19:10
Hej,
Naprężeń nie masz zapisanych w węzłach. Zmień set z węzłami na set z elementami, później zmień U na S i powinno być git.



pozdr
Tomek

-----------------------------
człowiekiem jestem ...
Autor RE: Odczytanie napręzen
Goral
Użytkownik

Postów: 21
Data rejestracji: 19.03.13
Dodane dnia 07/24/2013 12:02
Ok. Doszedłem do etapu, gdzie ptrafie odczytać naprężenia S22. Problem z tym, że tylko dla integration point 1. Chciałbym jeszcze panować nad punktem odczytu, najlepiej jakby CENTROID lub dla wszystkich punktow integration point (wtedy wyliczę ze średniej, szukaną wartość). Poniżej kod, który udało mi się wykonać.

import sys

from odbAccess import *
from abaqusConstants import *
from odbMaterial import *
from odbSection import *

odbpath=sys.argv[1]
outpath=sys.argv[2]

out = open(outpath,'w+')
odb = openOdb(path=odbpath)

nSet = odb.rootAssembly.instances['PART-1-1'].elementSets['SET-1']
frame = odb.steps['Step-1'].frames[-1]
dispField = frame.fieldOutputs['S']
comp = dispField.getScalarField(componentLabel='S11')
dispSubField = comp.getSubset(region=nSet)
fieldValues = dispSubField.values

for v1 in fieldValues:
print >> out, '%21.20e % (v1.data)

Proszę o pomoc.
Autor RE: Odczytanie napręzen
Tomekk
Użytkownik

Postów: 242
Miejscowość: Norge
Data rejestracji: 01.10.09
Dodane dnia 07/25/2013 20:15
Hej,

Popatrzyłem, trochę na Twój kod. Pobawiłem się abaqusem, po małych modyfikacjach u mnie też działa. Ale tak:
- jak miałeś output z samego S - a następnie iterowałeś po v1.data to w wyniku najprawdopodobniej ostałeś s11,
aby uzyskać misesy wystarczy zamiast v1.data wrzucić v1.mises
-on jak listuje Ci plik z wynikami, to jedzie po kolei element i wszystkie punkty całkowania dla elementu i potem kolejny element i kolejne znowu po kolei wszystkie punkty całkowania (sprawdź, sprawdź liczbę elementów pomnóż razy ilość punktów całkowania i powinno się zgadzać.

Teraz trochę wiedzy tajemnej, w wierszu poleceń w abaqusie są dwa fajne polecenia pomagające się zorientować co się dzieje i jakie pole do popisu mamy, a mianowicie polecenia print ##### oraz polecenie dir(#####) te ##### wklejasz komendę, czy tak jak masz u siebie zmienną, np. dispField.
Tylko to wszystko działa jak po prostu skopiujesz po kolei wszystkie polecenia ze skryptu łącznie z importami.


Także wygląda na to, że wszystkie cyfurki masz, tylko nie bardzo wiesz na co patrzysz :)

to polecenie:
print v1


wyswietla taki wynik:
({'baseElementType': 'C3D8', 'conjugateData': None, 'conjugateDataDouble': 'unknown', 'data': array([2.9844274520874, 17.7202224731445, 0.716619491577148, -7.66159963607788, -0.248252272605896, 0.491990804672241], 'f'), 'dataDouble': 'unknown', 'elementLabel': 1161, 'face': None, 'instance': 'OdbInstance object', 'integrationPoint': 8, 'inv3': 20.7496337890625, 'localCoordSystem': None, 'localCoordSystemDouble': 'unknown', 'magnitude': None, 'maxInPlanePrincipal': 0.0, 'maxPrincipal': 20.9967250823975, 'midPrincipal': 0.703007757663727, 'minInPlanePrincipal': 0.0, 'minPrincipal': -0.278463333845139, 'mises': 20.8018245697021, 'nodeLabel': None, 'outOfPlanePrincipal': 0.0, 'position': INTEGRATION_POINT, 'precision': SINGLE_PRECISION, 'press': -7.14042329788208, 'sectionPoint': None, 'tresca': 21.2751884460449, 'type': TENSOR_3D_FULL})

i teraz tak, jak wpiszesz:
print v1.mises
daje taką odpowiedź:
20.8018245697

podświetliłem wyżej tą cyferkę.

print v1.elementLabel
1161

print v1.integrationPoint
8
(tak, robiłem na solidach pierwszego rzędu z pełnym całkowaniem - C3D8)

także z poziomu frame.fieldOutputs['S'] i odpowiednio zapisując to:

for v1 in fieldValues:
print >> out, '%21.20e % (v1.data)

możesz śmiało wszystkie info niezbędne zapisać do pliku
np. tak:

for v1 in frame.values:
print >> out, 'El.No.: %4.0d; Int.Point: %d; VonMises stress: %3.5e' % (v1.elementLabel, v1.integrationPoint, v1.mises)



El.No.: 1; Int.Point: 1; VonMises stress: 4.92549e+02
El.No.: 1; Int.Point: 2; VonMises stress: 3.49180e+02
El.No.: 1; Int.Point: 3; VonMises stress: 5.29923e+02
El.No.: 1; Int.Point: 4; VonMises stress: 3.72549e+02
El.No.: 1; Int.Point: 5; VonMises stress: 1.02864e+03
El.No.: 1; Int.Point: 6; VonMises stress: 9.23107e+02
El.No.: 1; Int.Point: 7; VonMises stress: 1.08069e+03
El.No.: 1; Int.Point: 8; VonMises stress: 9.67219e+02
El.No.: 2; Int.Point: 1; VonMises stress: 5.99695e+02
El.No.: 2; Int.Point: 2; VonMises stress: 5.42943e+02
El.No.: 2; Int.Point: 3; VonMises stress: 6.02034e+02
El.No.: 2; Int.Point: 4; VonMises stress: 5.47491e+02
El.No.: 2; Int.Point: 5; VonMises stress: 8.79371e+02
El.No.: 2; Int.Point: 6; VonMises stress: 8.37208e+02
El.No.: 2; Int.Point: 7; VonMises stress: 8.81743e+02
El.No.: 2; Int.Point: 8; VonMises stress: 8.41092e+02
El.No.: 3; Int.Point: 1; VonMises stress: 5.19093e+02
El.No.: 3; Int.Point: 2; VonMises stress: 5.20190e+02
El.No.: 3; Int.Point: 3; VonMises stress: 5.07182e+02
El.No.: 3; Int.Point: 4; VonMises stress: 5.09875e+02
El.No.: 3; Int.Point: 5; VonMises stress: 8.06099e+02
El.No.: 3; Int.Point: 6; VonMises stress: 8.06445e+02
El.No.: 3; Int.Point: 7; VonMises stress: 7.93927e+02
El.No.: 3; Int.Point: 8; VonMises stress: 7.95243e+02





pozdr
Tomek

-----------------------------
człowiekiem jestem ...
Edytowane przez Tomekk dnia 07/25/2013 20:56
Autor RE: Odczytanie napręzen
Tomekk
Użytkownik

Postów: 242
Miejscowość: Norge
Data rejestracji: 01.10.09
Dodane dnia 07/25/2013 20:18
Duże piwo po obronie się należy

EDIT:
Zapomniałem na śmierć, na końcu, za pętlą dodaj:
out.close()


pozdr
Tomek

-----------------------------
człowiekiem jestem ...
Edytowane przez Tomekk dnia 07/25/2013 20:54
Autor RE: Odczytanie napręzen
Goral
Użytkownik

Postów: 21
Data rejestracji: 19.03.13
Dodane dnia 07/31/2013 23:26
Dzięki bardzo Tomekk za pomoc.
Rozwiązanie problemu było po innej stronie. Okazało się, że zastosowałem w modelu elementy z redukcją liczby punktów Gaussa... Wyłączenie tej funkcji rozwiązało mój problem.
Końcowa wersja wygląda następująco:

import sys

from odbAccess import *
from abaqusConstants import *
from odbMaterial import *
from odbSection import *

odbpath=sys.argv[1]
outpath=sys.argv[2]
i=sys.argv[3]

out = open(outpath,'w+')
odb = openOdb(path=odbpath)

nSet = odb.rootAssembly.instances['PART-'+str(i)+'-1'].elementSets['MONIT']
frame = odb.steps['Step-1'].frames[-1]
dispField = frame.fieldOutputs['S']
comp = dispField.getScalarField(componentLabel='S22')
dispSubField = comp.getSubset(region=nSet)
fieldValues = dispSubField.values

for v1 in fieldValues:
print >> out, '%5d, %5d, %21.20e' % (v1.elementLabel, v1.integrationPoint, v1.data)

Dodam tylko, dla osób które będą chciały skorzystać z przekazanej w temacie wiedzy, pamiętajcie o tabulatorach w skryptach Pythona.

Autor RE: Odczytanie napręzen
Tomekk
Użytkownik

Postów: 242
Miejscowość: Norge
Data rejestracji: 01.10.09
Dodane dnia 08/01/2013 07:23
[b]Goral napisał/a/b]
Dzięki bardzo Tomekk za pomoc.
Rozwiązanie problemu było po innej stronie. Okazało się, że zastosowałem w modelu elementy z redukcją liczby punktów Gaussa... Wyłączenie tej funkcji rozwiązało mój problem.

No tak, trzeba wiedzieć ile danych powinno być jak się już je sprawdzać zaczyna

[b]Goral napisał/a/b]
Dodam tylko, dla osób które będą chciały skorzystać z przekazanej w temacie wiedzy, pamiętajcie o tabulatorach w skryptach Pythona.

Dziwnie, że tu na forum jakoś strasznie ciężko wstawić taby, czy chociażby spacje aby jakoś to wyglądało.


A tak na marginesie, fajnie, że zaczęło działać.


pozdr
Tomek

-----------------------------
człowiekiem jestem ...
Skocz do Forum:
Reklama
Najtańszy sklep z dobrymi oponami

Szukaj
Google
 
Logowanie
Nazwa Użytkownika

Hasło



Nie jesteś jeszcze naszym Użytkownikiem?
Kilknij TUTAJ żeby się zarejestrować.

Zapomniane hasło?
Wyślemy nowe, kliknij TUTAJ.
Ankieta dla Użytkowników
Jakiego używasz programu MES

Nastran (wszelkiej maści)

Ansys

Abaqus

Cosmos/M

Comsol

Algor

MARC

Inny

Musisz się zalogować, żeby móc głosować w tej Ankiecie.
Reklama
Shoutbox
Tylko zalogowani mogą dodawać posty w shoutboksie.

Konrad96
02.12.2019
Witam, czy jest tu ktoś kto dobrze zna się na programie Robot struktural i pomógłby mi zamodelować belkę żelbetową podobną do belki użytej w badaniu? Proszę o kontakt jeśli ktoś chciałby pomóc

maq9633
10.11.2019
Witam, Zaczynam pisać pracę z płaskich den zbiorników ciśnieniowych i potrzebuje pomocy w programie Ansys workbench, proszę o kontakt

Vloiamwc
13.01.2019
Witam. Czy ktoś pomoże przygotować mnie do egzaminu z mes?

markos2101
30.01.2018
Dziękuję chociaż chodziło o znaczenie w MES - czyli zamiana obciążeń ciągłych na skupione

jasiu
29.01.2018
kongruentny: właściwy ; adekwatny ; zgodny ; poprawny ; godziwy ; dobry ; odpowiedni ; nadający się ; akuratny ; w sam raz ; należyty ; słuszny ; stosowny ; trafny ; przystający ;

markos2101
19.01.2018
Poszukuję definicji obciążenia kongruentnego - pomoże ktoś?

Adasiooo1
13.01.2018
witam szukam osób znających dobrze MES i abaqusa w zakresie żelbetu budownictwo

piotr
29.05.2017
Klachula nie jest do zadawania pytań od tego jest forum...

Geralt7586
28.05.2017
moźe ktoś mi pomoże ucze się femapa i mam problem jak użyję bar do rury to m=23kg a jak tube to 11,3 i wyrzuca mi blad SEKRRS przy tube a przy bar nie xd

kzurawski
08.04.2017
Jeśli chodzi o felgę to tu mamy jedną przeliczoną http://caelinux.org/wiki/
index.php/Contriblaws/C
ode_Aster/10_x_cases/plas
ticity


Archiwum
Wydarzenia
PonWtrCzwPiSoNie
1 2
3 4 5 6 7 8 9
10 11 12 13 14 15 16
17 18 19 20 21 22 23
24 25 26 27 28 29
Copyright Piotr© 2007- - 2020
Powered by Peter-Fusion © 2009 - Modified Aztec Theme