Author Archives: yhtian

About yhtian

I am an academic researcher working on Offshore Geotechnical Engineering. My blog aims to write down some work related trivial things and tricks about software, programming. It is basically a memo for me to back up some thoughts and small details. But I am more than happy if someone would visit and discuss.

Remote desktop problem with Abaqus screen resolution for python batch processing

Did not solve this problem yet. Just bookmark the below two pages for later use: https://blog.brankovucinec.com/2016/03/19/fix-remote-desktop-dpi-scaling-issues/ http://www.sqlchick.com/entries/2011/12/18/how-to-perfect-the-resolution-of-a-remote-desktop-session.html Advertisements

Posted in Uncategorized | Leave a comment

extract field data

did not get automatic extraction running just the manual method, anyway it works #manually open odb and change the belwo to extract from abaqusConstants import * myodb=’E:/CurrentProjects/LDFEUltimateState/4-CPE4R-SmoothThetaa-20X20B/UniaxialCapacity/theta=89.odb’ odb = session.openOdb(name=myodb) myfile = ‘E:/CurrentProjects/LDFEUltimateState/4-CPE4R-SmoothThetaa-20X20B/UniaxialCapacity/theta=89.txt’ session.xyDataListFromField(odb=odb, outputPosition=NODAL, variable=((‘RF’,  NODAL, ((COMPONENT, ‘RF1’), (COMPONENT, … Continue reading

Posted in Uncategorized | Leave a comment

ABAQUS display windows disappeared

Strangely ABAQUS display windows disappeared. The reason was just become drag to another screen with different resolution, where the python script windows was enlarged for coding. When drag back or restart on the main screen. All display windows disappeared. Waste … Continue reading

Posted in Uncategorized | Leave a comment

Transition mesh with equal ratio serial

lenL=xPoint[5]-C2[0] lenR=C3[0]-xPoint[5] hmax=2. hmin=0.5 ratio=hmax/hmin NL=int(np.log(ratio)/np.log((lenL-hmin)/(lenL-hmin*ratio))+1.) NR=int(np.log(ratio)/np.log((lenR-hmin)/(lenR-hmin*ratio))+1.) hminL=lenL*(ratio**(1./(NL-1))-1.)/(ratio**(NL/(NL-1.))-1.) hminR=lenR*(ratio**(1./(NR-1))-1.)/(ratio**(NR/(NR-1.))-1.) qL=ratio**(1./(NL-1.)) qR=ratio**(1./(NR-1.)) xpt=xPoint[5] for j in range(NL-1): hstep=hminL*qL**j mdb.models[‘Model-1’].sketches[‘__profile__’].Line(point1=(xpt,C2[1]), point2=(xpt-hstep,C2[1])) xpt=xpt-hstep mdb.models[‘Model-1’].sketches[‘__profile__’].Line(point1=(xpt,C2[1]), point2=C2) xpt=xPoint[5] for j in range(NR-1): hstep=hminR*qR**j mdb.models[‘Model-1’].sketches[‘__profile__’].Line(point1=(xpt,C2[1]), point2=(xpt+hstep,C2[1])) xpt=xpt+hstep mdb.models[‘Model-1’].sketches[‘__profile__’].Line(point1=(xpt,C2[1]), point2=(C3[0],C3[1]))

Posted in Uncategorized | Leave a comment

What is the “hourglass effect” in finite element analysis? How does the reduced integration, resulting in the hourglass effect, work? How can we counter the hourglass effect?

https://www.quora.com/What-is-the-hourglass-effect-in-finite-element-analysis-How-does-the-reduced-integration-resulting-in-the-hourglass-effect-work-How-can-we-counter-the-hourglass-effect … What is hourglassing? It is essentially a spurious deformation mode of a Finite Element Mesh, resulting from the excitation of zero-energy degrees of freedom. It typically manifests as a patchwork of zig-zag or hourglass like element shapes (Fig.1), … Continue reading

Posted in Uncategorized | Leave a comment

e-lnp

Posted in Uncategorized | Leave a comment

Contact problem with sharp tip

A sharp tip (CPT) penetrates into soil. ABAQUS Node-surface interaction works well both small-sliding and finite-sliding. Surface-surface interaction works in case of small-sliding but in case finite-sliding of surface-surface contact (the default option) there is no result (seem no contact … Continue reading

Posted in Uncategorized | Leave a comment